Question: Unable to load index
0
gravatar for mohammedabeddin1
3.6 years ago by
mohammedabeddin110 wrote:

I am trying to convert .vcf to .fasta using the .fasta reference by using vcf-consensus but I am getting errors. I am unable to get the reference file as faidx location reference to the reference fasta files are generating errors. It would be nice if someone helped me by giving all the steps for converting a .vcf to .fasta file .

I am converting .vcf file to .fasta file for the location Molecular Location: base pairs 166,195,185 to 166,375,987 on chromosome 2 (Homo sapiens Annotation Release 107, GRCh38.p2) using ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/technical/reference/human_g1k_v37.fasta.gz

tool samtools bcftools • 1.5k views
ADD COMMENTlink modified 3.6 years ago by genomax74k • written 3.6 years ago by mohammedabeddin110
2
gravatar for genomax
3.6 years ago by
genomax74k
United States
genomax74k wrote:

VCF to FASTA

You should be able to index your fasta by using samtools faidx. Something along the lines of

$ samtools faidx human_g1k_v37.fasta.gz

That should generate a human_g1k_v37.fasta.gz.fai file.

ADD COMMENTlink modified 3.6 years ago • written 3.6 years ago by genomax74k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1595 users visited in the last hour