Comments (4)
Hi @JPFinnigan,
This could be happening because your mm10.fa
and GTF file do not match in the chromosome names. Can you check that they match?
from hisat2.
Hello,
The FASTA file and the snp_file have the same naming convention for chromosomes (chr_). I think this problem may stem from the fact that my snp_file is a .vcf. I was looking through the some other posted issues and it seems as though its not yet possible to pass a .vcf to extract_snps.py
. But perhaps this is not the case?
Either way I wanted to explain myself. I'm sure this has come up in your internal discussions, but we're excited about the potential of using HISAT2 to support alignment against custom reference-graphs created from SNV/INDEL call-sets derived from joint mutation calling from tumor-normal pairs. (e.g a graph reference containing the reference allele, plus germline SNP/INDEL, plus somatic SNP/INDEL). We (cc: @iskandr) would be eager to hear your thoughts.
from hisat2.
Hi @JPFinnigan,
There isn't a script to do it from a VCF file, but I think it wouldn't be too hard to write one. I think you'd have to decompose multiallelic variants and normalize them so they are all left shifted first. You can do both of those operations with vt decompose
and vt normalize
(http://genome.sph.umich.edu/wiki/Vt) @infphilo, if we did that and called SNP/indel/deletion via just looking at length differences between ref and alt, would that be enough?
from hisat2.
HISAT2 now supports VCF file - see the HISAT2 webpage. I will update the usage information to make clear which files and formats are used.
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.