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 06:56 AM
not uniquely mapped reads unidodo RNA Sequencing 2 04-22-2011 03:07 PM
Uniquely mapped reads with bowtie mapper Bioinformatics 2 11-22-2010 11:44 PM
BWA Uniquely Mapped Reads NF_seq Bioinformatics 0 09-06-2010 04:32 AM
cufflinks and non-uniquely mapped reads clariet Bioinformatics 1 05-08-2010 12:13 PM

Thread Tools
Old 11-10-2010, 09:31 PM   #1
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, 03:03 AM   #2
Location: Shanghai,China

Join Date: Sep 2011
Posts: 30

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, 06:58 AM   #3
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

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 05:26 PM.

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