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
_______________________________________________ miso-users mailing list [email protected] http://mailman.mit.edu/mailman/listinfo/miso-users
