View Single Post
Old 01-14-2013, 06:01 AM   #3
Derek-C
Junior Member
 
Location: U.S.A.

Join Date: Nov 2012
Posts: 7
Default

Quote:
Originally Posted by marcowanger View Post
I remember Bam2fastq only export read pair. So, all singleton mapped will be discarded. In your case, it may be better to write a small script to parse the singleton out
Are you certain about the singletons thing? It would make sense, butl I just did some checking there with samtools flagstat and according to the numbers I'm looking at for this one accepted_hits.bam file from tophat there were more singleton reads then there were reads discarded by bam2fastq (~3 million more). So if bam2fastq doesn't export singletons, what happened to those 3 million reads.

Actually there are options in tophat for setting it so that there should be no singletons, hopefully that'll fix this.

Last edited by Derek-C; 01-14-2013 at 06:29 AM.
Derek-C is offline   Reply With Quote