Question: paired-end alignment with tophat2 expected (mean) inner distance between mate pair
0
gravatar for kpr
15 months ago by
kpr60
kpr60 wrote:

How do I find the expected (mean) inner distance between mate pair?

I am having an issue with tophat and I think it is because I have 150bp paired end data and I didn't change the -r argument for tophat.

tophat -r ???? ref_genome file-1.fastq file-2.fastq

I've been looking at this post:

paired-end alignment with tophat2

I have my QC report, but I am not sure what I am looking at.

rna-seq • 480 views
ADD COMMENTlink modified 15 months ago by Carlo Yague4.4k • written 15 months ago by kpr60

You should know that the old 'Tuxedo' pipeline of Tophat(2) and Cufflinks is no longer the "advisable" tool for RNA-seq analysis. The software is deprecated/ in low maintenance and should be replaced by HISAT2, StringTie and ballgown. See this paper: Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. (If you can't get access to that publication, let me know and I'll -cough- help you.) There are also other alternatives, including alignment with STAR and bbmap, or pseudo-alignment using salmon.

ADD REPLYlink written 15 months ago by WouterDeCoster38k
0
gravatar for Carlo Yague
15 months ago by
Carlo Yague4.4k
Belgium
Carlo Yague4.4k wrote:

You can infer the mean inner distance between mates in different ways:

1. Experimental: The inner distance depends directly on the fragments sizes of the cDNA library. This can be measured by running a bioanalyzer chip on it. You can then estimate the inner distance as Inner_Dist = Fragment_Size - (Adapters + Index) - (2x Read_Length)

2. In silico: After mapping the reads (using default -r value for instance),you can infer the mean inner distance as the average genomic distance between two concordantly mapped reads. This value can be extracted from the 9th field of the SAM files (TLEN) or using Picard CollectInsertSizeMetrics.

ADD COMMENTlink written 15 months ago by Carlo Yague4.4k
Please log in to add an answer.

Help
Access

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