SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
HELP!! XS tag to filter unique reads Do NOT match bowtie2 report momo Bioinformatics 24 10-11-2019 09:11 AM
how to extract unique hits from a sam file gfmgfm Bioinformatics 10 07-07-2019 08:33 AM
Exonerate: filter low quality hits masterpiece RNA Sequencing 0 01-27-2013 08:17 PM
unique hits arrchi Bioinformatics 7 06-01-2011 08:53 AM
Unique and repeat Hits pfranchini Bioinformatics 1 07-15-2009 09:51 AM

Reply
 
Thread Tools
Old 05-05-2015, 05:23 AM   #1
MattB
Member
 
Location: Norway

Join Date: Aug 2008
Posts: 35
Default Filter unique hits in Bowtie2 for stranded RNA-seq

Hi all,

We have PE stranded RNA-seq data, and we want to filter the alignment files for uniquely mapping reads (ie. no multi-hit reads).

We have done this before for 'unstranded' data based on the XS tag in the alignment file, however all I seem to get in the file now is:

XS:A:+ and XS:A:-

Any tips on how to filter multi-hit reads with stranded data? I saw that this thread (http://seqanswers.com/forums/showthread.php?t=39206) suggested calculating the MAPQ score...

Matt
MattB is offline   Reply With Quote
Old 05-05-2015, 05:34 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

You appear to be using tophat2 rather than bowtie2, which is good given that you have RNAseq reads. In that case, just filter by MAPQ, namely those having MAPQ > 3 are "uniquely aligned"*.

*For one of the many definitions of that term.
dpryan is offline   Reply With Quote
Old 05-05-2015, 06:37 AM   #3
MattB
Member
 
Location: Norway

Join Date: Aug 2008
Posts: 35
Default

Thanks!

That is correct, we are using TopHat2.

I have also been debating what 'uniquely' means, I guess it comes down to where you set a threshold
MattB 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 03:20 PM.


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