Question: junctions.bed file in hisat2
1
gravatar for badribio
15 months ago by
badribio230
badribio230 wrote:

Hi,

I was wondering if there is a way I can generate junctions.bed file from hisat2 just like the one in tophat.

Hisat2 documentation says use extract_splice_site.py python script but that gets splice site from gtf file, is there a way to know splice sites and no of reads associated with those sites using hisat2?

Thanks

pipeline alignment splicing hisat • 1.1k views
ADD COMMENTlink modified 15 months ago by EagleEye5.3k • written 15 months ago by badribio230
4
gravatar for EagleEye
15 months ago by
EagleEye5.3k
Sweden
EagleEye5.3k wrote:

Easiest solution will be to use featureCounts and counts reads from exon-exon junctions by providing your aligned BAM files. Check manual page 35 for detailed information.

**-J (juncCounts)**    Count the number of reads supporting each exon-exon junction.
Junctions are identified from those exon-spanning reads
(containing ‘N’ in CIGAR string) in input data. The output
result includes names of primary and secondary genes that
overlap at least one of the two splice sites of a junction. Only
one primary gene is reported, but there might be more than
one secondary gene reported. Secondary genes do not overlap
more splice sites than the primary gene. When the primary
and secondary genes overlap same number of splice sites, the
gene with the smallest leftmost base position is selected as
the primary gene. Also included in the output result are the
position information for the left splice site (‘Site1’) and the
right splice site (‘Site2’) of a junction. These include chromosome
name, coordinate and strand of the splice site. In
the last columns of the output, number of supporting reads is
provided for each junction for each library
ADD COMMENTlink modified 15 months ago • written 15 months ago by EagleEye5.3k

Thanks, does it also help with reads spanning introns?

My aim to finally use DEXSeq...

ADD REPLYlink written 15 months ago by badribio230

I have no idea about reads spanning introns. But I guess you can find explanation in the manual.

ADD REPLYlink written 15 months ago by EagleEye5.3k
4
gravatar for badribio
15 months ago by
badribio230
badribio230 wrote:

Figured out this does the job

https://regtools.readthedocs.io/en/latest/commands/junctions-extract/
ADD COMMENTlink written 15 months ago by badribio230
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: 1683 users visited in the last hour