Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • How to generate VCF from HISAT2 pre-built SNP index?

    My ultimate goal is to get a (n= samples, m= SNPs) data matrix. My plan was to use HISAT2 for the mapping, VCF tools for the vcf file, and then parse it to generate the data matrix I can actually mine.

    I'm using the pre-built SNP index file for H. sapiens, Ensembl GRCh38 ftp://ftp.ccb.jhu.edu/pub/infphilo/h...h38_snp.tar.gz . I have HISAT2 running smoothly for all of my samples and started reading the downstream pipeline for generating VCF files (https://ccb.jhu.edu/software/hisat2/manual.shtml).

    Code:
    samtools mpileup -uf $HISAT2_HOME/example/reference/22_20-21M.fa eg2.sorted.bam | bcftools view -bvcg - > eg2.raw.bcf
    How do I get the original fasta file or build a VCF file using this index and my sam/bam files? I was going to just download the hg38 Ensemble annotated genome but I don't think that's what I need. . . I went into the `make_grch38_snp.sh`file from the tar ball when downloading the SNPs db. I think it's building the SNP index from `Homo_sapiens.GRCh38.dna.primary_assembly.fa. Is this the file that needs to be used? (ftp://ftp.ensembl.org/pub/release-84...assembly.fa.gz)

    Also, if anyone has any insight on how to generate a data matrix from the vcf files, it would be greatly appreciated (but first I need to generate the vcf files)

    Thanks in advance
    Attached Files
    Last edited by jol.espinoz; 04-30-2017, 11:26 AM.

Latest Articles

Collapse

  • seqadmin
    Techniques and Challenges in Conservation Genomics
    by seqadmin



    The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

    Avian Conservation
    Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
    03-08-2024, 10:41 AM
  • seqadmin
    The Impact of AI in Genomic Medicine
    by seqadmin



    Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
    02-26-2024, 02:07 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 03-14-2024, 06:13 AM
0 responses
34 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-08-2024, 08:03 AM
0 responses
72 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-07-2024, 08:13 AM
0 responses
82 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-06-2024, 09:51 AM
0 responses
68 views
0 likes
Last Post seqadmin  
Working...
X