Add Strand tag to hisat2 alignments
1
0
Entering edit mode
7 months ago
RT ▴ 360

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 • 473 views
ADD COMMENT
0
Entering edit mode

Option seems to run your alignments with --dta option (see: A: Adding tag XS to HISAT2 sam files perior to StringTie )

ADD REPLY
0
Entering edit mode

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.

ADD REPLY
1
Entering edit mode
7 months ago

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/

ADD COMMENT

Login before adding your answer.

Traffic: 1621 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6