[miso-users] Error when running miso

Rosannah C Cameron rosannah.cameron at einstein.yu.edu
Mon Mar 30 10:15:07 EDT 2015 

Hello,

I’ve been having some trouble getting MISO to work with my RNA-seq data. I’m relatively new to RNA-seq analysis and I’m probably just missing something relatively simple. I installed MISO using the latest version from the GitHub repository and Python’s pip install misopy. I checked that all the packages were installed correctly and was able to Test MISO and it seemed to be working correctly. I was also able to compute the insert length distribution for my BAM file.

When I try to run MISO with my BAM file I get the following errors:


miso --run /Volumes/Sherlock/RNA-Seq_analysis_RCC_Jan_2015/Reference_genome/indexed/ /Volumes/Sherlock/RNA-Seq_analysis_RCC_Jan_2015/Mut_1/Mut_1_sorted.bam --output-dir Mut_1/ --read-len 125 --paired-end 167 33

MISO (Mixture of Isoforms model)

Probabilistic analysis of RNA-Seq data for detecting differential isoforms

Use --help argument to view options.


Using MISO settings file: /Library/Python/2.7/site-packages/misopy/settings/miso_settings.txt

Computing Psi values...

  - GFF index: /Volumes/Sherlock/RNA-Seq_analysis_RCC_Jan_2015/Reference_genome/indexed

  - BAM: /Volumes/Sherlock/RNA-Seq_analysis_RCC_Jan_2015/Mut_1/Mut_1_sorted.bam

  - Read length: 125

  - Output directory: /Users/bowmanlab/Mut_1

Checking your GFF annotation and BAM for mismatches...

Checking if BAM has mixed read lengths...

Found reads of length 125 in BAM.

Mapping genes to their indexed GFF representation, using /Volumes/Sherlock/RNA-Seq_analysis_RCC_Jan_2015/Reference_genome/indexed

Searching for /Volumes/Sherlock/RNA-Seq_analysis_RCC_Jan_2015/Reference_genome/indexed/genes_to_filenames.shelve..

  - Found shelved file.

Preparing to run 4 batches of jobs...

Running batch of 8434 genes..

  - Executing: python /Library/Python/2.7/site-packages/misopy/run_miso.py --compute-genes-from-file "/Users/bowmanlab/Mut_1/batch-genes/batch-0_genes.txt" /Volumes/Sherlock/RNA-Seq_analysis_RCC_Jan_2015/Mut_1/Mut_1_sorted.bam /Users/bowmanlab/Mut_1 --read-len 125  --paired-end 167.0 33.0 --settings-filename /Library/Python/2.7/site-packages/misopy/settings/miso_settings.txt

  - Submitted thread batch-0

Running batch of 8434 genes..

  - Executing: python /Library/Python/2.7/site-packages/misopy/run_miso.py --compute-genes-from-file "/Users/bowmanlab/Mut_1/batch-genes/batch-1_genes.txt" /Volumes/Sherlock/RNA-Seq_analysis_RCC_Jan_2015/Mut_1/Mut_1_sorted.bam /Users/bowmanlab/Mut_1 --read-len 125  --paired-end 167.0 33.0 --settings-filename /Library/Python/2.7/site-packages/misopy/settings/miso_settings.txt

  - Submitted thread batch-1

/Library/Python/2.7/site-packages/misopy/miso_sampler.py:306: RuntimeWarning: Long inserted alignment at pysplicing/src/solve.c:291

  stop_cond)

Running batch of 8434 genes..

  - Executing: python /Library/Python/2.7/site-packages/misopy/run_miso.py --compute-genes-from-file "/Users/bowmanlab/Mut_1/batch-genes/batch-2_genes.txt" /Volumes/Sherlock/RNA-Seq_analysis_RCC_Jan_2015/Mut_1/Mut_1_sorted.bam /Users/bowmanlab/Mut_1 --read-len 125  --paired-end 167.0 33.0 --settings-filename /Library/Python/2.7/site-packages/misopy/settings/miso_settings.txt

  - Submitted thread batch-2

/Library/Python/2.7/site-packages/misopy/miso_sampler.py:306: RuntimeWarning: Long deleted alignment at pysplicing/src/solve.c:281

  stop_cond)

Running batch of 8435 genes..

  - Executing: python /Library/Python/2.7/site-packages/misopy/run_miso.py --compute-genes-from-file "/Users/bowmanlab/Mut_1/batch-genes/batch-3_genes.txt" /Volumes/Sherlock/RNA-Seq_analysis_RCC_Jan_2015/Mut_1/Mut_1_sorted.bam /Users/bowmanlab/Mut_1 --read-len 125  --paired-end 167.0 33.0 --settings-filename /Library/Python/2.7/site-packages/misopy/settings/miso_settings.txt

  - Submitted thread batch-3

/Library/Python/2.7/site-packages/misopy/miso_sampler.py:306: RuntimeWarning: Long deleted alignment at pysplicing/src/solve.c:281

  stop_cond)

Waiting on 4 threads...

/Library/Python/2.7/site-packages/misopy/miso_sampler.py:306: RuntimeWarning: Long inserted alignment at pysplicing/src/solve.c:291

  stop_cond)

/Library/Python/2.7/site-packages/misopy/miso_sampler.py:306: RuntimeWarning: Long deleted alignment at pysplicing/src/solve.c:281

  stop_cond)

/Library/Python/2.7/site-packages/misopy/miso_sampler.py:306: RuntimeWarning: Long inserted alignment at pysplicing/src/solve.c:291

  stop_cond)

/Library/Python/2.7/site-packages/misopy/miso_sampler.py:306: RuntimeWarning: Long deleted alignment at pysplicing/src/solve.c:281

  stop_cond)

/Library/Python/2.7/site-packages/misopy/miso_sampler.py:306: RuntimeWarning: Long inserted alignment at pysplicing/src/solve.c:291

  stop_cond)

python(23371,0x7fff736a4310) malloc: *** error for object 0x7faf43321c30: incorrect checksum for freed object - object was probably modified after being freed.

*** set a breakpoint in malloc_error_break to debug

/bin/sh: line 1: 23371 Abort trap: 6           python /Library/Python/2.7/site-packages/misopy/run_miso.py --compute-genes-from-file "/Users/bowmanlab/Mut_1/batch-genes/batch-3_genes.txt" /Volumes/Sherlock/RNA-Seq_analysis_RCC_Jan_2015/Mut_1/Mut_1_sorted.bam /Users/bowmanlab/Mut_1 --read-len 125 --paired-end 167.0 33.0 --settings-filename /Library/Python/2.7/site-packages/misopy/settings/miso_settings.txt >> "/Users/bowmanlab/Mut_1/batch-logs/batch-3-03-30-15_09:23:29.log"

/bin/sh: line 1: 23369 Segmentation fault: 11  python /Library/Python/2.7/site-packages/misopy/run_miso.py --compute-genes-from-file "/Users/bowmanlab/Mut_1/batch-genes/batch-2_genes.txt" /Volumes/Sherlock/RNA-Seq_analysis_RCC_Jan_2015/Mut_1/Mut_1_sorted.bam /Users/bowmanlab/Mut_1 --read-len 125 --paired-end 167.0 33.0 --settings-filename /Library/Python/2.7/site-packages/misopy/settings/miso_settings.txt >> "/Users/bowmanlab/Mut_1/batch-logs/batch-2-03-30-15_09:23:28.log"

python(23366,0x7fff736a4310) malloc: *** error for object 0x7fd52157b420: incorrect checksum for freed object - object was probably modified after being freed.

*** set a breakpoint in malloc_error_break to debug

/bin/sh: line 1: 23366 Abort trap: 6           python /Library/Python/2.7/site-packages/misopy/run_miso.py --compute-genes-from-file "/Users/bowmanlab/Mut_1/batch-genes/batch-1_genes.txt" /Volumes/Sherlock/RNA-Seq_analysis_RCC_Jan_2015/Mut_1/Mut_1_sorted.bam /Users/bowmanlab/Mut_1 --read-len 125 --paired-end 167.0 33.0 --settings-filename /Library/Python/2.7/site-packages/misopy/settings/miso_settings.txt >> "/Users/bowmanlab/Mut_1/batch-logs/batch-1-03-30-15_09:23:28.log"

python(23364,0x7fff736a4310) malloc: *** error for object 0x7ff814359150: incorrect checksum for freed object - object was probably modified after being freed.

*** set a breakpoint in malloc_error_break to debug

/bin/sh: line 1: 23364 Abort trap: 6           python /Library/Python/2.7/site-packages/misopy/run_miso.py --compute-genes-from-file "/Users/bowmanlab/Mut_1/batch-genes/batch-0_genes.txt" /Volumes/Sherlock/RNA-Seq_analysis_RCC_Jan_2015/Mut_1/Mut_1_sorted.bam /Users/bowmanlab/Mut_1 --read-len 125 --paired-end 167.0 33.0 --settings-filename /Library/Python/2.7/site-packages/misopy/settings/miso_settings.txt >> "/Users/bowmanlab/Mut_1/batch-logs/batch-0-03-30-15_09:23:27.log"

03/30/2015 09:28:15 AM - miso_main - WARNING - Thread batch-0 might have failed...

03/30/2015 09:28:15 AM - miso_main - WARNING - Thread batch-1 might have failed...

03/30/2015 09:28:15 AM - miso_main - WARNING - Thread batch-2 might have failed...

03/30/2015 09:28:15 AM - miso_main - WARNING - Thread batch-3 might have failed…



I used the STAR aligner to create my SAM files and converted them to sorted indexed BAM files with samtools (outside of Miso). I’m working with zebrafish and downloaded my gtf file from Ensembl (Danio_rerio.Zv9.78.gtf) and then used the Galaxy (GTF-to-GFF3) to convert from gtf to a gff3 file (I used Galaxy because I had some trouble with the gtf2gff3.pl script). I wanted to include a small version of my BAM file but I wasn’t sure what the best way to do this would be.


Any suggestions would be greatly appreciated!


Thank you,

Rosannah





-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.mit.edu/pipermail/miso-users/attachments/20150330/e10296f6/attachment-0001.htm

More information about the miso-users mailing list