SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
TopHat --prefilter-multihits parameter id0 Bioinformatics 0 06-30-2014 12:39 PM
Merging accepted_hits and unmapped reads from tophat o/p alyamahmoud RNA Sequencing 1 11-19-2013 01:30 AM
extrat unaligned reads from unmapped.bam (tophat) vivienne_lovely Bioinformatics 1 05-26-2013 10:58 PM
de novo assembly using unmapped reads from tophat ccstaats Bioinformatics 4 04-12-2013 02:09 PM
assembly unmapped reads from tophat leifive Bioinformatics 3 03-07-2013 10:03 PM

Reply
 
Thread Tools
Old 07-01-2015, 09:29 AM   #1
imsharmanitin
Postdoc Cancer Bioinformatics
 
Location: Olso, Norway

Join Date: Dec 2014
Posts: 17
Cool TopHat -M/--prefilter-multihits; unmapped reads

hi!

i am confused when to use TopHat "-M/--prefilter-multihits and when not

"When mapping reads on the transcriptome, some repetitive or low complexity reads that would be discarded in the context of the genome may appear to align to the transcript sequences and thus may end up reported as mapped to those genes only. This option directs TopHat to first align the reads to the whole genome in order to determine and exclude such multi-mapped reads (according to the value of the -g/--max-multihits option)."

the alignment of reads first to genome and then to transcriptome seems opposite to -G option

As far i know the multi-mapped reads will end up in Unmapped.bam along with the reads which would have not mapped even if have not used -M option

so do we need to discard all unmapped.bam reads?

i am struggling to understand when to use -M option and when not and how to deal with the multi-mapped reads discarded if -M option is used
imsharmanitin is offline   Reply With Quote
Old 07-01-2015, 10:00 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,975
Default

Umapped reads file comes in handy if you find a large number of reads not mapping because you can then chase those down to see if you have a contamination or some other sort of a problem.

For some of tophat analyses you may want to run multiple versions with/without option combinations to see if one works better than the other (e.g. -G or -G with -T etc).
GenoMax is offline   Reply With Quote
Reply

Tags
multihit reads, rna seq experiment design, rna sequencing, tophat 2, unmapped reads

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:33 PM.


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