SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
How to get uniquely mapped reads from Tophat subeet Bioinformatics 10 11-28-2012 05:56 AM
not uniquely mapped reads unidodo RNA Sequencing 2 04-22-2011 02:07 PM
Uniquely mapped reads with bowtie mapper Bioinformatics 2 11-22-2010 10:44 PM
BWA Uniquely Mapped Reads NF_seq Bioinformatics 0 09-06-2010 03:32 AM
cufflinks and non-uniquely mapped reads clariet Bioinformatics 1 05-08-2010 11:13 AM

Reply
 
Thread Tools
Old 11-10-2010, 08:31 PM   #1
mrfox
Senior Member
 
Location: USA

Join Date: Aug 2010
Posts: 103
Default Tophat uniquely mapped reads

I have just posted two questions but this time I would make it simpler and clearer. If I do not use "-g1" when running Tophat(i.e., allow multi-hits), then how can I identify uniquely mapped reads? My understanding is they are mapped reads with "NH:i:1", where NH is "Number of reported alignments that contains the query in the current record". Am I right?
Thank you.
mrfox is offline   Reply With Quote
Old 08-28-2012, 02:03 AM   #2
wisense
Member
 
Location: Shanghai,China

Join Date: Sep 2011
Posts: 30
Default

Yes,you are right.When the 5 column of the bam file is 255,it also means it is an uniquely mapped read.
wisense is offline   Reply With Quote
Old 05-23-2013, 05:58 AM   #3
rahilsethi
Member
 
Location: Pittsburgh, PA

Join Date: May 2010
Posts: 22
Default NH:i:1 does not give uniquely mapped reads when TopHat run under default settings

No, when running TopHat "NH:i:1" will not report uniquely mapped reads when run with default settings. This is because under default settings max no of reported primary alignments are 20 even when it is actually greater than 20. When I looked for max value of NH:i: in the output of whole transcriptome sequence it is 20. So some reads, in reality may contain more than 1 alignment NH:i: will give only the no. of reported alignment and in this case only primary alignments are reported. NH:i:1 for some reads may still give you 1 if only primary alignments are 1, even if the total alignments is greater than 1
rahilsethi 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:15 AM.


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