SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Tophat2: bam file has no quality sequences for a lot of reads duhaimj Bioinformatics 4 04-28-2012 07:47 AM
BWA - unmapped reads Adamo Bioinformatics 23 02-22-2012 10:12 PM
why retain unmapped reads? MeganS Bioinformatics 5 01-13-2012 09:42 PM
What are the unmapped reads beelu Illumina/Solexa 1 09-09-2010 06:18 AM
% unmapped reads bioinfosm Illumina/Solexa 8 07-05-2010 01:36 AM

Reply
 
Thread Tools
Old 06-21-2012, 02:37 PM   #1
manianslab
Junior Member
 
Location: USA

Join Date: Jun 2012
Posts: 9
Default tophat2/bowtie2 inconsistency in number of unmapped reads

We are using the tophat-cufflinks pipeline for RNA seq. With tophat1-bowtie1 mapping, the reads are mapped initially to the predicted transcriptome and then unmapped reads are mapped to the genome to identify splice junctions. So the number of unmapped reads from the initial mapping (see bowtie.left_kept_reads.fixmap.log below) are carried to the next mapping (see bowtie.left_kept_reads.m2g_um.fixmap.log below; I believe m2g_um stands for map to gtf unmapped)

bowtie.left_kept_reads.fixmap.log
# reads processed: 1761894
# reads with at least one reported alignment: 1334963 (75.77%)
# reads that failed to align: 413755 (23.48%)
# reads with alignments suppressed due to -m: 13176 (0.75%)
Reported 2240278 alignments to 1 output stream(s)

bowtie.left_kept_reads.m2g_um.fixmap.log
# reads processed: 413755
# reads with at least one reported alignment: 169629 (41.00%)
# reads that failed to align: 241180 (58.29%)
# reads with alignments suppressed due to -m: 2946 (0.71%)
Reported 184162 alignments to 1 output stream(s)

However, when we try the same dataset with tophat2-bowtie2, the read counts don't seem to match. Mapping parameters are identical (zero mismatches, zero indels). However, the unmapped reads in the first mapping (bowtie.left_kept_reads.fixmap.log) is different from what the second mapping starts with (bowtie.left_kept_reads.m2g_um.fixmap.log).

bowtie.left_kept_reads.fixmap.log
1761894 reads; of these: 1761894 (100.00%) were unpaired; of these:
223293 (12.67%) aligned 0 times
468278 (26.58%) aligned exactly 1 time
1070323 (60.75%) aligned >1 times
87.33% overall alignment rate

bowtie.left_kept_reads.m2g_um.fixmap.log
413755 reads; of these: 413755 (100.00%) were unpaired; of these:
72334 (17.48%) aligned 0 times
190748 (46.10%) aligned exactly 1 time
150673 (36.42%) aligned >1 times
82.52% overall alignment rate

Any of you have any explanations? Thanks for your time!!
manianslab is offline   Reply With Quote
Old 07-13-2012, 01:47 PM   #2
liueagle
Junior Member
 
Location: NYC

Join Date: Mar 2012
Posts: 2
Default

I am also looking forward to the explanation. Recently, I also face this problem: http://seqanswers.com/forums/showthread.php?t=21663
liueagle is offline   Reply With Quote
Old 07-13-2012, 01:56 PM   #3
manianslab
Junior Member
 
Location: USA

Join Date: Jun 2012
Posts: 9
Default

@liueagle, it appears that since bowtie2 by default allows less stringent matching (e.g. gaps), the number of mapped reads is higher, but tophat perhaps checks them and makes sure they meet the criteria. Due to the repeated release of betas for cufflinks2, we have put on hold the tophat2-cufflinks2 pipeline and instead are using the stable release tophat1-cufflinks1.
manianslab is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 09:40 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO