Question: Question about Tophat2-parameter"-r <int> 50 as default"and "--mate-std-dev <int> 20 as default"
gravatar for Yunlong JIA
4.3 years ago by
Yunlong JIA20
Yunlong JIA20 wrote:

I ran my data using Tophat2, for the "-r" and "--mate-std-dev", I left them as default value.

Next, I watch the paired-end reads in Savant Genome Browser, I check the insetsize for the paired reads, they are so variant.

As default, insersize should in the range of 50(+- 20), from 30~70. But I found many insersize bigger than 200nt, or overlapping.

Why the default value didn't work? The default is randomized insertsize when pairing the reads?

Hope your help, thanks.

rna-seq tophat2 parameter • 947 views
ADD COMMENTlink written 4.3 years ago by Yunlong JIA20

I just asked a different poster to do this and will ask you as well.

You can use BBMap (if your reads don't overlap) or BBMerge (if you know that R1/R2 overlap) to get an idea of insert size distribution. Both programs are part of BBMap suite. You will want to use ihist= flag to create a file with insert size distribution.

ADD REPLYlink modified 4.3 years ago • written 4.3 years ago by genomax89k

Thank you, maybe it's some bug for tools, thanks for your advise, I will try this BB.~

ADD REPLYlink written 4.3 years ago by Yunlong JIA20
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1563 users visited in the last hour