This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Add Strand tag to hisat2 alignments

Hi All,

Hisat2 does not report XS tag for each alignment like tophat2. How do I extract Strand info and add XS tag for each alignment?

RT

rna hisat2
• 2,163 views • link
updated by Istvan Albert Istvan Albert • written by GR GR

Thanks for your response. I think --dta adds XS tag only for reads spliced reads alignments. I need this XS tag for all the reads.

• 0 views • link
written by GR GR

1 answer

You will need to use the --rna-strandness parameter:

--rna-strandness <string>

Specify strand-specific information: the default is unstranded.

For single-end reads, use F or R.

‘F’ means a read corresponds to a transcript. ‘R’ means a read corresponds to the reverse complemented counterpart of a transcript.

For paired-end reads, use either FR or RF. With this option being used, every read alignment will have an XS attribute tag:

’+’ means a read belongs to a transcript on ‘+’ strand of genome. ‘-‘ means a read belongs to a transcript on ‘-‘ strand of genome.

(TopHat has a similar option, –library-type option, where fr-firststrand corresponds to R and RF; fr-secondstrand corresponds to F and FR.)

this info can be found in the manual:

http://daehwankimlab.github.io/hisat2/manual/

• 0 views • link
written by Istvan Albert Istvan Albert

Log in to answer this question.

More posts like this