Question: Hisat2 alignment results can be ued for TEtoolkit analysis
0
gravatar for lkianmehr
9 months ago by
lkianmehr30
France
lkianmehr30 wrote:

Hello,

I have aligned RNA-seq reads by Hisat2, then ran TEtoolkit on sam files of Hisat2, and it works. can I trust these results? or is better to align reads again with STAR; then use TEtoolkit?

tetoolkit hisat2 • 254 views
ADD COMMENTlink modified 9 months ago by h.mon26k • written 9 months ago by lkianmehr30

Not sure what the TEtoolkit does but I don't think there will be an issue with swapping STAR for Hisat2 (or vice versa). Different flavors of mapping will give different results though but I assume the majority will be in common.

ADD REPLYlink written 9 months ago by lieven.sterck5.5k
1
gravatar for h.mon
9 months ago by
h.mon26k
Brazil
h.mon26k wrote:

It should be fine, but you probably have to set -k 200, similarly to what the TEToolkit manual suggests for STAR.

edit: I didn't see you already mapped the reads. Check the settings, the issue is, as you are trying to quantify features with potentially hundreds of copies in the genome, you want to be sure the aligned is considering them. STAR, by default, considers as unmapped reads mapping to more than 20 locations, hence you have to provide an additional parameter when you want to ten TEToolkit downstream. You should check how HISAT2 works in this regard, I am not familiar with it.

ADD COMMENTlink modified 9 months ago • written 9 months ago by h.mon26k

I checked Hisat2 in this case and faced this sentence: HISAT2 is not designed with large values for -k in mind, and when aligning reads to long, repetitive genomes large -k can be very, very slow.

ADD REPLYlink written 9 months ago by lkianmehr30
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: 616 users visited in the last hour