[BioC] missing reads in alignerBAM

Daniel.Berner at unibas.ch Daniel.Berner at unibas.ch
Fri Dec 17 09:55:44 CET 2010


hi list
When I align the N reads contained in my Illumina fastq using  
NovoAlign, convert the resulting SAM output to BAM using Samtools, and  
then upload the BAM file using readAligned (ShortRead), the resulting  
AlignedRead object contains fewer than N reads. Reconverting the BAM  
to SAM and counting the lines using wc -l, however, indicates there  
are still N reads in the BAM. Hence it appears that during upload into  
R, some reads are lost. The AlignedRead object, however, will contain  
exactly N reads when I perform the alignment using bowtie instead of  
NovoAlign.
Any idea/comment on why reads are missing?
Thanks, Daniel



More information about the Bioconductor mailing list