Skip to content

Junction annotation question (unknown_strand) #200

@SaraAlidadiani

Description

@SaraAlidadiani

Hi,

I was using regtools junction annotation to annotate a junction bed file of junction extract output. In my results, I see a lot of significant events that are classified as unknown_strand. The junction file was extracted from STAR aligned BAM file, and I used GENCODE v42 .
As it was from Standard Illumina and unstranded, I used -s XS.

The code I used is:

for bamfile in `ls /home/leafcutter-bam/*.bam`; do
    echo Converting $bamfile to $bamfile.junc
    samtools index $bamfile
    regtools junctions extract -a 8 -m 50 -M 500000 $bamfile -o $bamfile.junc -s XS
    echo $bamfile.junc >> test_juncfiles_mayo.txt
done

Why are these events categorized as "unkown_strand"?
Could you help me with this, please?

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions