Comments (7)
Thanks, especially given the detailed information!, which should help us a lot to shoot the problem.
Did you use the official version of HISAT2 1.0.1-beta, available on the HISAT2 website? If so, that version didn't work with --mm option. I fixed this issue later, with the fix included in the master branch, here.
Thanks,
Daehwan
from hisat2.
I just confirmed this was fixed in the GitHub master version - I'll try to release the next release sometime this week, thanks!
from hisat2.
I also replicate this problem in the latest hisat2 build on Mac OS X Mavericks (v10.9.5). Additionally, a different error occurs on linux systems (hisat2 terminates, rather than running indefinitely like on Mac OS X):
(ERR): hisat2-align died with signal 11 (SEGV)
when administered with the following command:
~/tools/hisat2-master/hisat2 -x /mnt/ramdisk/Mus_musculus/grcm38_snp_tran/genome_snp_tran -1 input_1.fastq -2 input_2.fastq -S output_mm.sam --mm
Perhaps my command notation is incorrect?
from hisat2.
I think the problem may be with the current "grcm38_snp_tran" index on the hisat2 website. I just tried building my own index with transcripts (no SNP data) using the following command, but hisat2 now appears to execute much slower compared to the previous index on the hisat2 site:
../hisat2-2.0.2-beta/hisat2-build --ss ../splicesites.txt --exon ../exons.txt -p 28 ../Mus_musculus.GRCm38.dna.primary_assembly.fa mus_base
I used the following files from ensembl:
Primary genome - ftp://ftp.ensembl.org/pub/release-84/fasta/mus_musculus/dna/Mus_musculus.GRCm38.dna.primary_assembly.fa.gz
GTF file - ftp://ftp.ensembl.org/pub/release-84/gtf/mus_musculus/Mus_musculus.GRCm38.84.gtf.gz
The machine I used for the index building had 208GB of RAM - I suspect that should be enough to facilitate optimal '--dcv' values? Is there any other secret sauce to the reference indices?
from hisat2.
I just found out this is an issue in index building. There is nothing wrong with --mm option. I'll try to rebuild grcm38_snp_tran index. The other mouse indexes are fine.
from hisat2.
For what it's worth, I just tried running hisat2 with indices that I built using the parameters above and the "--mm" command, and the following error resulted on a Linux box:
Error while flushing and closing output
terminate called after throwing an instance of 'int'
(ERR): hisat2-align died with signal 6 (ABRT)
Parallel execution seems to work fine though.
UPDATE: Nevermind, sorry - I was just being stupid and had ran out of storage space. Everything works perfectly with the home-rolled index.
from hisat2.
I fixed this issue in the HISAT2 master branch and rebuilt genome_snp_tran indexes for the human and the mouse genomes.
from hisat2.
Related Issues (20)
- hisat2 hangs aligning axolotl reads HOT 1
- Output files(.snp, .haplotype) of hisat2_extract_snps_haplotypes_*.py are empty
- Please add the pbat option of hisat-3n
- A question about methylation information extraction
- Any plans to support Apple Silicon architecture? HOT 1
- Installation Issue Error 1 - make HOT 1
- -np argument seemingly not working
- ERR): "fastq file.fastq" does not exist. Exiting now ...
- [Bug Report] hisat2-align exited with value 137, space complexity of hisat2
- hisat2 location does not exist
- Hisat-3N mapping quality
- hisat2-build index for circRNA-seq
- hisat2-build failed for Segmentation fault
- [Future request] hisat-3n table option to report conversions summarized to genomic feature or reads counts
- Issue with hisatgenotype HOT 1
- Mapping using different parameters --very-sensitive and default
- (ERR): "ref.genome" does not exist Exiting now ...
- --directional-mapping-reverse vs. --rna-strandness on HISAT3N
- Question about calculation of base counts in hisat3Ntable
- mkfifo failed error and change $temp_dir HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from hisat2.