No announcement yet.

New Resources for 1000 Genomes

This is a sticky topic.
  • Filter
  • Time
  • Show
Clear All
new posts

  • #16
    As far as chrY

    We provide all our variation data in VCF format which serves our needs quite well, if you have a better idea for your own needs then you should be able to get all the info you need from these files to do the conversion

    Look at for streaming if you want to avoid downloading the entire data set


    • #17
      vcf file of specific sample from 1000Genome data


      Can anyone help me how to access the vcf file of a specific sample from 1000Genome data. I found the consensus file at ( but couldn't find the individual samples.

      I am trying to compare the variants found from our sequencing vs 1000Genome. if anyone has done similar analysis please let know I would to discuss wiht you offline.

      Thanks in advance


      • #18
        You should be able to get this info from our vcf files using a combination of tabix anc vcftools vcf-subset as described in our faq


        • #19

          Thanks much for your reply. I am guessing this is the example for getting the vcf of sample.

          tabix -h 17:1471000-1472000 | perl /nfs/1000g-work/G1K/work/bin/vcftools/perl/vcf-subset -c HG00098 | bgzip -c /tmp/HG00098.20100804.genotypes.vcf.gz


          • #20
            That is correct


            • #21

              how/what should I specify, if I don't have particular region to look at and want to get all genome-wide variants?

              Thanks so much for you kind help.


              • #22
                You can give tabix a whole chromosome but be aware tabix can not cope with losed network connectivity so when streaming large data volumes that can cause lossed data which means you may need to download the whole file



                • #23
                  there is another paper with 113 pages, "supplemental information"

                  with a referrence:


                  pointing back to :


                  which is 19 pages
                  Last edited by gsgs; 12-07-2012, 08:15 AM.


                  • #24

                    I tried with downloading both the vcf.gz and tbi files. but it did not work and it is difficult to interpret the error. can you see what I am doing wrong here

                    ./tabix -h /Volumes/Macintosh\ HD\ 3/1000Genome/ALL.wgs.phase1_release_v3.20101123.snps_indels_sv.sites.vcf.gz | perl ~/othertools/vcftools_0.1.10/perl/vcf-subset -c NA10851

                    [tabix] the index file exists. Please use '-f' to overwrite.
                    Broken VCF header, no column names?
                    at /Users/molpathuser1/othertools/vcftools_0.1.10/perl// line 177
                    Vcf::throw('Vcf4_1=HASH(0x7fa9d982f8d8)', 'Broken VCF header, no column names?') called at /Users/molpathuser1/othertools/vcftools_0.1.10/perl// line 869
                    VcfReader::_read_column_names('Vcf4_1=HASH(0x7fa9d982f8d8)') called at /Users/molpathuser1/othertools/vcftools_0.1.10/perl// line 604
                    VcfReader:arse_header('Vcf4_1=HASH(0x7fa9d982f8d8)') called at /Users/molpathuser1/othertools/vcftools_0.1.10/perl/vcf-subset line 122
                    main::vcf_subset('HASH(0x7fa9d98288f0)') called at /Users/molpathuser1/othertools/vcftools_0.1.10/perl/vcf-subset line 12

                    many thanks for your kind help


                    • #25
                      You need to give tabix some sort of chromosome name otherwise it doesn't know what to fetch

                      If you just want to filter the whole file you will need to use zcat

                      That being said you downloaded the sites file which contains no genotype info and no columns with individual genotypes to filter


                      • #26
                        Hi Laura,

                        I am still having trouble with extracting the variant calls of a specific sample.

                        As you pointed out earlier that I have downloaded sites file with no genotype column, so now I got this version ALL.2of4intersection.20100804.genotypes.vcf.gz vcf and tbi file from ftp site (release/20100804).

                        and I used the following command to subset the vcf file

                        tabix -fh /Volumes/Macintosh_HD_3/1000Genome/ALL.2of4intersection.20100804.genotypes.vcf.gz 1 | perl ~/othertools/vcftools_0.1.10/perl/vcf-subset -c NA10851 > NA10851/NA10851_chr1

                        but strangely the out-put file has all genotype columns. I have been following the directions given on the 1000 genome except the I don't give the range for chromosome as I want to get all variants. So I tried with giving the coordinates (see below) and result file has just the header only.

                        here is the command i used,
                        tabix -fh /Volumes/Macintosh_HD_3/1000Genome/ALL.2of4intersection.20100804.genotypes.vcf.gz 2:1-243199373 | perl ~/othertools/vcftools_0.1.10/perl/vcf-subset -c NA10851 > NA10851/NA10851_chr2

                        so now I really don't know what I am doing wrong in trying to subset the vcf file. I really appreciate for your kind help so far and would be very grateful if you could help me how to solve this.

                        thank you so much.


                        • #27
                          Hello Rama

                          Unfortunately I can not recreate your issue

                          Using your command I get a vcf file which just contains genotypes for NA10851


                          • #28
                            Thanks Laura, for trying it out.


                            • #29
                              Sorry if this is already been asked, I didn't find it..

                              I try to figure out if I can do a search by read length.
                              I am looking for reads length 101.
                              Is there a way to know this information before downloading?
                              I looked in the sequence.index but I didn't find this.

                              Thanks in advance,



                              • #30
                                I'm going to download Bam files from the Project.
                                I see two links:

                                There are some overlapping files between the two links.

                                I would like to know which one I should use?