This is a test version of Biostars. For the public version, visit https://www.biostars.org.
[Error] getting fasta file from 1000genomes

Hi,

I'd like to get a fasta file of haplotypes from 1000 genomes.

I ran the commands below:

http://www.1000genomes.org/faq/are-there-any-fasta-files-containing-1000-genomes-variants-or-haplotypes

tabix -h ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20110521/ALL.chr17.phase1_release_v3.20101123.snps_indels_svs.genotypes.vcf.gz 17:1471000-1472000 | perl vcf-subset -c HG00098 | bgzip -c > HG00098.vcf.gz
tabix -p vcf HG00098.vcf.gz
cat ref.fa | vcf-consensus HG00098.vcf.gz > HG00098.fa

but got error massages,

[tabix] the index file is older than the vcf file. Please use '-f' to overwrite or reindex.
Can't open perl script "vcf-subset": No such file or directory

And when I tried to install vcf-tools, I got

E: Unable to locate package vcf-tools

How can I run these commands successfully?

My machine is Ubuntu 64 bit on Win7.0.

genome gene software-error
• 4,871 views • link
updated by Ram Ram • written by Mari • Mari

Can you give some more information about how you tried to install vcf-tools?

• 0 views • link
written by Matt Shirley Matt Shirley

Matt,

Thank you. I tried to run this:

marie@ubuntu:~/Downloads$ sudo apt-get install vcf-tools
Reading package lists... Done
Building dependency tree
Reading state information... Done
E: Unable to locate package vcf-tools
• 0 views • link
updated by Ram Ram • written by Mari • Mari

vcftools is not in the distrib software repo. You have to install it 'manually'. See here.

• 0 views • link
updated by Ram Ram • written by Phil S. Phil S.

Thank you Phil, maybe I move ahead a bit but got another error.

marie@ubuntu:~/Downloads/vcftools_0.1.8a/perl$ tabix -h ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20110521/ALL.chr17.phase1_release_v3.20101123.snps_indels_svs.genotypes.vcf.gz 17:1471000-1472000 | perl vcf-subset -c HG00098 | bgzip -c > HG00098.vcf.gz
[tabix] the index file is older than the vcf file. Please use '-f' to overwrite or reindex.
Broken VCF header, no column names?
 at /usr/share/perl5/Vcf.pm line 171
    Vcf::throw('Vcf4_1=HASH(0x160ee48)', 'Broken VCF header, no column names?') called at /usr/share/perl5/Vcf.pm line 845
    VcfReader::_read_column_names('Vcf4_1=HASH(0x160ee48)') called at /usr/share/perl5/Vcf.pm line 589
    VcfReader::parse_header('Vcf4_1=HASH(0x160ee48)') called at vcf-subset line 119
    main::vcf_subset('HASH(0x16051c8)') called at vcf-subset line 12
• 0 views • link
updated by Ram Ram • written by Mari • Mari

I think you can safely use the -f flag that tabix is warning you about. Currently it looks like tabix is not returning any data - you can check this by simply omitting everything after and including your first pipe character.

• 0 views • link
updated by Ram Ram • written by Matt Shirley Matt Shirley

Thank you Matt, I added -f like this and the program ran:

tabix -h ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20110521/ALL.chr17.phase1_release_v3.20101123.snps_indels_svs.genotypes.vcf.gz 17:1471000-1472000 -f

(Initially I tried to run the code below and it didn't work)

tabix -h -f ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20110521/ALL.chr17.phase1_release_v3.20101123.snps_indels_svs.genotypes.vcf.gz 17:1471000-1472000

but at the last line, I got many errors:

• 0 views • link
updated by Ram Ram • written by Mari • Mari

0 answers

No answers yet.

Log in to answer this question.

More posts like this