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!
•
written
by
acorella •
0 answers
No answers yet.
Log in to answer this question.
-
Estimating Mean Inner Distance
written by CF •
-
Automatically feeding a good mate inner distance to TopHat2
written by umn_bist •
-
tophat2 aligment ratio reduced when increase insert size
written by 52Teth •
-
Mean insert size and Mean inner distance between mate pairs
written by manekineko
-
the defination of insert size in RNA-seq, including the length of read, or not?
written by zju.whw •
-
Fragment Size: TLEN vs. isize
written by Rob
-
Manorm For Differential Histone Binding In Chip-Seq Samples
written by Diana
-
Strange Behavior In Picard Downsamplesam.Jar Algorithm?
written by jeffrey •
-
Estimated Insert Size Less Than Read Pair Sizes (180Bp For 100X100 Pe), But Seqprep Merged 0 Reads?
written by Ian Fiddes
-
Tophat Parameters Misunderstanding
written by Nicolas Rosewick
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