Question: How many should multiple alignment of tophat-out in Human RNA-seqs?
0
gravatar for Leila Kian
5 weeks ago by
Leila Kian0
Leila Kian0 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 • 84 views
ADD COMMENTlink modified 5 weeks ago by Jennifer Hillman Jackson23k • written 5 weeks ago by Leila Kian0
0
gravatar for Jennifer Hillman Jackson
5 weeks ago by
United States
Jennifer Hillman Jackson23k 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 5 weeks ago by Jennifer Hillman Jackson23k
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: 111 users visited in the last hour