frattalab / rna_seq_snakemake

Snakemake Pipeline For RNA-Seq in UCL Cluster
9 stars 4 forks source link

rseq.smk - BAM index file paths are incorrect (causing occasional crashes with junction_saturation.py) #59

Open SamBryce-Smith opened 6 months ago

SamBryce-Smith commented 6 months ago

rseqc.smk

rule read_distribution:
    input:
        bam = os.path.join(STAR_OUTDIR, "{sample}.Aligned.sorted.out.bam"),
        idx = os.path.join(STAR_OUTDIR, "{sample}.Aligned.sorted.out.bam")

This is the case for all RSeQC rules. Might be cause of slow run times for other tools that do not error out when index is not present. Very quick fix

SamBryce-Smith commented 6 months ago
image

Example errors