XPRESSyourself / XPRESSpipe

An alignment and analysis pipeline for Ribosome Profiling and RNA-seq data
https://xpresspipe.readthedocs.io/en/latest/
GNU General Public License v3.0
12 stars 4 forks source link

Periodicity not working with yeast #18

Open j-berg opened 5 years ago

j-berg commented 5 years ago

riboWaltz broken for processing any yeast files

j-berg commented 4 years ago
...
Great! All reads' reference transcript IDs were found in annotation table. No reads removed.
3.09 M  (5.73 %) reads removed: mapping on negative strand.
reading W303_60mCCCP_RNA_S12_L007_R1_001_Aligned.toTranscriptome.out.bam
Good! Number of indel below indel_threshold for all reads. No reads removed.
reads: 120.72 M
Great! All reads' reference transcript IDs were found in annotation table. No reads removed.
5.98 M  (4.96 %) reads removed: mapping on negative strand.
reading W303_60mCCCP_RNA_S12_L008_R1_001_Aligned.toTranscriptome.out.bam
Good! Number of indel below indel_threshold for all reads. No reads removed.
reads: 122.59 M
Great! All reads' reference transcript IDs were found in annotation table. No reads removed.
6.06 M  (4.95 %) reads removed: mapping on negative strand.
processing W303_00m_FP_S1_L007_R1_001_Aligned.toTranscriptome.out
Error in if (extremity == "auto" & ((best_from3_tab[, perc] > best_from5_tab[,  : 
  argument is of length zero
Calls: psite
In addition: Warning messages:
1: In min(site_sub$length) :
  no non-missing arguments to min; returning Inf
2: In max(site_sub$length) :
  no non-missing arguments to max; returning -Inf
Execution halted

I think is the same issue I would get when I tried running a really small test BAM file through, even though these files should be plenty large