This is a test version of Biostars. For the public version, visit https://www.biostars.org.
rMATS can't open bam created by hisat2

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 rmats return this

gtf: 0.9844145774841309
There are 33431 distinct gene ID in the gtf file
There are 33432 distinct transcript ID in the gtf file
There are 33430 one-transcript genes in the gtf file
There are 174723 exons in the gtf file
There are 4833 one-exon transcripts in the gtf file
There are 4833 one-transcript genes with only one exon in the transcript
Average number of transcripts per gene is 1.000030
Average number of exons per transcript is 5.226220
Average number of exons per transcript excluding one-exon tx is 5.940417
Average number of gene per geneGroup is 1.978729
statistic: 0.012811422348022461
Fail to open /mnt/c/SCIProgram/trans/02hisat2/hisat2_output/LrG1.bam
/mnt/c/SCIProgram/trans/02hisat2/hisat2_output/LrG2.bam
/mnt/c/SCIProgram/trans/02hisat2/hisat2_output/LrG3.bam
Fail to open /mnt/c/SCIProgram/trans/02hisat2/hisat2_output/LrR1.bam
/mnt/c/SCIProgram/trans/02hisat2/hisat2_output/LrR2.bam
/mnt/c/SCIProgram/trans/02hisat2/hisat2_output/LrR3.bam

read outcome totals across all BAMs
USED: 0
NOT_PAIRED: 0
NOT_NH_1: 0
NOT_EXPECTED_CIGAR: 0
NOT_EXPECTED_READ_LENGTH: 0
NOT_EXPECTED_STRAND: 0
EXON_NOT_MATCHED_TO_ANNOTATION: 0
JUNCTION_NOT_MATCHED_TO_ANNOTATION: 0
CLIPPED: 0
total: 0
outcomes by BAM written to: /mnt/c/SCIProgram/trans/04rMATS/tmp_output/2022-04-10-11_18_34_323169_read_outcomes_by_bam.txt

novel: 0.005591869354248047
The splicing graph and candidate read have been saved into /mnt/c/SCIProgram/trans/04rMATS/tmp_output/2022-04-10-11_18_34_323169_*.rmats
save: 0.004040718078613281
Traceback (most recent call last):
  File "/home/ianlee/biotools/rmats_turbo_v4_1_2/rmats.py", line 595, in <module>
    main()
  File "/home/ianlee/biotools/rmats_turbo_v4_1_2/rmats.py", line 563, in main
    run_pipe(args)
  File "rmatspipeline/rmatspipeline.pyx", line 3860, in rmats.rmatspipeline.run_pipe
  File "rmatspipeline/rmatspipeline.pyx", line 3723, in rmats.rmatspipeline.split_sg_files_by_bam
  File "rmatspipeline/rmatspipeline.pyx", line 3731, in rmats.rmatspipeline.split_sg_files_by_bam
ValueError: invalid literal for int() with base 10: '/mnt/c/SCIProgram/trans/02hisat2/hisat2_output/LrG2.bam'

and this is my list file

hisat2 rmats bam
• 1,441 views • link
written by ianlee96 • ianlee96

2 answers

enter image description here

this is my list file

and what's the problem can you tell me? Can rMATS open bam created by hisat2 instead of STAR?

• 0 views • link
written by ianlee96 • ianlee96

this is my file

/mnt/c/SCIProgram/trans/02hisat2/hisat2_output$ ls
LbG1.bam      LbG3.bam.bai  LbR3.bam      LbT2.bam.bai  LrG1.bam.bai  LrR1.bam      LrR3.bam.bai  LrT3.bam
LbG1.bam.bai  LbR1.bam      LbR3.bam.bai  LbT3.bam      LrG2.bam      LrR1.bam.bai  LrT1.bam      LrT3.bam.bai
LbG2.bam      LbR1.bam.bai  LbT1.bam      LbT3.bam.bai  LrG2.bam.bai  LrR2.bam      LrT1.bam.bai  hisat2.sh
LbG2.bam.bai  LbR2.bam      LbT1.bam.bai  Lbfix.exon    LrG3.bam      LrR2.bam.bai  LrT2.bam      hisat21.sh
LbG3.bam      LbR2.bam.bai  LbT2.bam      LrG1.bam      LrG3.bam.bai  LrR3.bam      LrT2.bam.bai  stringtie.sh
• 0 views • link
written by ianlee96 • ianlee96

Log in to answer this question.

More posts like this