[miso-users] errors with run_events_analysis

[Maayan Kreitzman]

Hi everyone,

I'm new to MISO, and would like some help. I've successfully installed misopy-0.4.7.

I also downloaded and indexed the ensembl gff3 (no "chr" in the chromosome names).
The next step isn't working though. My library is a paired-end strand specific library.

The command I ran is this:

python misopy/run_events_analysis.py --paired-end 75 8 --compute-genes-psi /projects/mkreitzman_prj/expression_quantification_testing/testing/MISO/indexed/ /projects/mkreitzman_prj/expression_quantification_testing/testing/test_data/strand_specific/JAGuaR_rev_87647/A08586_gr3.neg.bam --output-dir /projects/mkreitzman_prj/expression_quantification_testing/testing/MISO/results --read-len 100 &> /projects/mkreitzman_prj/expression_quantification_testing/testing/MISO/results/A08586.log

the errors I get look like this (I've attached the full std out + stderr):
Traceback (most recent call last):
  File "/projects/mkreitzman_prj/expression_quantification_testing/misopy-0.4.7/misopy/run_miso.py", line 488, in <module>
    main()
  File "/projects/mkreitzman_prj/expression_quantification_testing/misopy-0.4.7/misopy/run_miso.py", line 432, in main
    run_compute_genes_from_file(options)
  File "/projects/mkreitzman_prj/expression_quantification_testing/misopy-0.4.7/misopy/run_miso.py", line 247, in run_compute_genes_from_file
    event_type=options.event_type)
  File "/projects/mkreitzman_prj/expression_quantification_testing/misopy-0.4.7/misopy/run_miso.py", line 128, in compute_gene_psi
    gene_obj)
  File "/home/mkreitzman/.local/lib/python2.7/site-packages/misopy-0.4.7-py2.7-linux-x86_64.egg/misopy/sam_utils.py", line 191, in fetch_bam_reads_in_gene
    chrom = chrom_parts[1]
IndexError: list index out of range

The bam file is position sorted and has a .bai in the same folder. The parameters for the --paired-end are just off the top of my head, but I tried without this option at all and the same thing happens. 

Thanks in advance for any ideas you may have.

Maayan
-------------- next part --------------
A non-text attachment was scrubbed...
Name: A08586.log
Type: text/x-log
Size: 13445 bytes
Desc: A08586.log
Url : http://mailman.mit.edu/pipermail/miso-users/attachments/20130320/37d9f62f/attachment.bin