This is a test version of Biostars. For the public version, visit https://www.biostars.org.
rMATS reads used is zero (0)

I've tried running rMATS with either fastq or bam files for my experiment (paired reads) and I keep getting zero used reads. I checked my fastq files and the read length is 150 but a lot of my reads are under "NOT_EXPECTED_READ_LENGTH". Can someone help me troubleshoot this? Here is my code for running with bam files:

python rmats.py --b1 /home/data/Megan/AS_analysis/rMATS/ScrVeh.txt --b2 /home/data/Megan/AS_analysis/rMATS/ScrIns.txt --gtf /home/data/Megan/AS_analysis/rMATS/gencode.v41.annotation.gtf -t paired --readLength 150 --nthread 4 --od /home/data/Megan/AS_analysis/rMATS/ScrVeh_ScrIns --tmp /home/data/Megan/AS_analysis/rMATS

output:

read outcome totals across all BAMs
USED: 0
NOT_PAIRED: 10463908
NOT_NH_1: 756430308
NOT_EXPECTED_CIGAR: 6469470
NOT_EXPECTED_READ_LENGTH: 540038736
NOT_EXPECTED_STRAND: 0
EXON_NOT_MATCHED_TO_ANNOTATION: 0
JUNCTION_NOT_MATCHED_TO_ANNOTATION: 0
CLIPPED: 137075821
total: 1450478243
splicing python rmats alternative
• 3,379 views • link
updated by benformatics benformatics • written by Megan Noonan • Megan Noonan

1 answer

Use the --variable-read-length option?

https://github.com/Xinglab/rmats-turbo/blob/v4.1.2/README.md

EDIT: You might have to download the newest version of rmats for this param

• 0 views • link
written by benformatics benformatics

Okay, I did that and now here is my output. Definitely have used reads now but am not sure if the other parameters are within a decent range (specifically it seems like exons and junctions not matched are high).

read outcome totals across all BAMs
USED: 513516620
NOT_PAIRED: 10463908
NOT_NH_1: 756430308
NOT_EXPECTED_CIGAR: 6469470
NOT_EXPECTED_READ_LENGTH: 0
NOT_EXPECTED_STRAND: 0
EXON_NOT_MATCHED_TO_ANNOTATION: 25505241
JUNCTION_NOT_MATCHED_TO_ANNOTATION: 1016875
CLIPPED: 137075821
total: 1450478243
• 0 views • link
written by Megan Noonan • Megan Noonan

are you sure your reads were aligned to the same reference + GTF as the one used for rmats? if not it could be problematic... either re-align with the same annotation as rmats or trying to build or own rmats reference. or maybe you aligned allowing for novel splice junctions??? anyway you can at least check out your results from the reads that did work and see if they make sense

• 0 views • link
written by benformatics benformatics

Good question, I'll have to check. If I start with fastq files this wouldn't matter, would it?

• 0 views • link
written by Megan Noonan • Megan Noonan

If rmats does the alignment then yes you should be fine

• 0 views • link
written by benformatics benformatics

Log in to answer this question.

More posts like this
  • rMATS reports zero evenets?
    written by PK • PK
    I'm running rmats/4.1.2 to detect the splicing events. python /software/rmats/4.1.2/rmats.py --gtf /reference_annotation/Homo_sapiens.GRCh38.86.gtf -t paired \ --b1 ko_1.txt --b2 wt_1.txt --od /Rmats_test \ --tmp /Rmats_test --libType …
  • rMATS can't open bam created by hisat2
    written by ianlee96 • ianlee96
    this is my code : python /home/ianlee/biotools/rmats_turbo_v4_1_2/rmats.py --b1 /mnt/c/SCIProgram/trans/04rMATS/list/04LrG.txt --b2 /mnt/c/SCIProgram/trans/04rMATS/list/04LrR.txt --gtf /mnt/c/SCIProgram/genome/Lb/Lbfix.gtf -t paired --readLength 50 --nthread 11 --od /mnt/c/SCIProgram/trans/04rMATS/output --tmp /mnt/c/SCIProgram/trans/04rMATS/tmp_output --novelSS and …