Flexibility of inner distance parameter for Tophat
0
0
Entering edit mode
5.5 years ago
acorella ▴ 30

I downloaded ~70 SRA files for RNA-seq data and am attempting to align them. I used Picard's CollectInsertSizeMetrics to calculate insert size for two random files in this set (following this protocol: http://vinaykmittal.blogspot.com/2012/02/how-to-estimate-insert-size-for-paired.html).

However, the insert sizes from the two random files are vastly different 360, vs 180 (both with std dev around 60). The libraries were supposedly all prepared in the same way. Should I calculate insert size for each file individually to set the inner distance parameter for tophat? Or is there some flexibility in this parameter that tophat has built in?

Any help would be appreciated! Thanks!

RNA-Seq • 1.3k views
ADD COMMENT
0
Entering edit mode

I am sure there is some flexibility for this parameter, but it is probably best to work it out yourself for every data set. Can derive the inner distance with informatics.

Automatically feeding a good mate inner distance to TopHat2

ADD REPLY

Login before adding your answer.

Traffic: 1911 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6