Question about Tophat2-parameter"-r <int> 50 as default"and "--mate-std-dev <int> 20 as default"
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.
• 1,886 views
•
link
written
by
Yunlong JIA •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Weird results in TopHat2 paired-end alignment
written by Aspire
-
Are there any tools that can create a very basic GTF file from contig sequences (no annotations rea…
written by O.rka
-
Getting unmapped paired-end reads as fastq from unmapped.bam obtained with tophat2
written by cris.kgs •
-
RNA-seq High Variation Between and Within Groups
written by jacobsen.jeremy •
-
optimizing the -r/--mate-inner-dist in Tophat2
written by candida.vaz •
-
Optimization of RNA-Seq data mapping with tophat2
written by Denis
-
Tophat2 and Bowtie2 : determining appropriate values for mate-inner-dist and mate-std-dev
written by irritable_phd_syndrome •
-
low percentage mapping... tophat2.1.1
written by mra8187 •
-
Estimated Fragment Mean In Cuffdiff
written by apt.university
-
Tophat Parameters Misunderstanding
written by Nicolas Rosewick
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.Thank you, maybe it's some bug for tools, thanks for your advise, I will try this BB.~