Question: How many should multiple alignment of tophat-out in Human RNA-seqs?
0
gravatar for Leila Kian
11 months ago by
Leila Kian10
Leila Kian10 wrote:

Hi

I'm checking my tophat output data and I would like to know the rate of this multiple alignments are acceptable or not

whether is dependent this rate of multiple alignments to type of genome or not? this RNA-seqs have belonged to Human.

Reads: Input : 36238973 Mapped : 30511468 (84.2% of input) of these: 14351857 (47.0%) have multiple alignments (15253 have >20) 84.2% overall read mapping rate.

thanks in advance

rna-seq tophat alignment • 664 views
ADD COMMENTlink modified 11 months ago by Jennifer Hillman Jackson25k • written 11 months ago by Leila Kian10
0
gravatar for Jennifer Hillman Jackson
11 months ago by
United States
Jennifer Hillman Jackson25k wrote:

Hello,

The multiple alignment rate seems high. But, this could be how you set up the job parameters with Tophat/Tophat2 (I'm not sure which are you using, and that may not matter going forward). Tophat/Tophat2 have been deprecated. Please use HISAT2 instead.

To retain more or fewer alignments per query sequence (primary alignments, secondary alignments), the parameters can be tuned. For HISAT2, please see: https://ccb.jhu.edu/software/hisat2/manual.shtml#reporting.

More help:

  • Galaxy NGS QA/QC plus RNA-seq Tutorials: https://galaxyproject.org/learn/

  • Search prior Q&A at https://www.biostars.org where there is much discussion about mapping tools, reported alignments, fastq QA (trimming, adaptor removal), and how different parameters can impact mapping rates/reported alignments.

Thanks! Jen, Galaxy team

ADD COMMENTlink written 11 months ago by Jennifer Hillman Jackson25k
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 16.09
Traffic: 169 users visited in the last hour