Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
This topic is closed.
X
This is a sticky topic.
X
X
 
  • 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 http://www.1000genomes.org/faq/how-d...your-vcf-files for streaming if you want to avoid downloading the entire data set

    Comment


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

      Hi,

      Can anyone help me how to access the vcf file of a specific sample from 1000Genome data. I found the consensus file at (ftp://ftp-trace.ncbi.nih.gov/1000genomes/ftp/release) 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
      Rama

      Comment


      • #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

        1000genomes.org is your first and best source for all of the information you’re looking for. From general topics to more of what you would expect to find here, 1000genomes.org has it all. We hope you find what you are searching for!

        Comment


        • #19
          Laura,

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

          tabix -h ftp://ftp-trace.ncbi.nih.gov/1000gen...804/ALL.2of4in... 17:1471000-1472000 | perl /nfs/1000g-work/G1K/work/bin/vcftools/perl/vcf-subset -c HG00098 | bgzip -c /tmp/HG00098.20100804.genotypes.vcf.gz

          Comment


          • #20
            That is correct

            Comment


            • #21
              Laura,

              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.

              Comment


              • #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

                Comment


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


                  with a referrence:

                  ...
                  38!
                  Garrison,!E.!K.!vcflib$K$a$simple$C++$library$for$parsing$and$manipulating$
                  VCF$files,!<https://github.com/ekg/vcflib>!(2012).!

                  pointing back to :

                  1000genomes.org is your first and best source for all of the information you’re looking for. From general topics to more of what you would expect to find here, 1000genomes.org has it all. We hope you find what you are searching for!


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

                  Comment


                  • #24
                    Laura,

                    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//Vcf.pm line 177
                    Vcf::throw('Vcf4_1=HASH(0x7fa9d982f8d8)', 'Broken VCF header, no column names?') called at /Users/molpathuser1/othertools/vcftools_0.1.10/perl//Vcf.pm line 869
                    VcfReader::_read_column_names('Vcf4_1=HASH(0x7fa9d982f8d8)') called at /Users/molpathuser1/othertools/vcftools_0.1.10/perl//Vcf.pm 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

                    Comment


                    • #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

                      Comment


                      • #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.
                        Rama

                        Comment


                        • #27
                          Hello Rama

                          Unfortunately I can not recreate your issue

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

                          Comment


                          • #28
                            Thanks Laura, for trying it out.

                            Comment


                            • #29
                              Hi,
                              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,
                              Pap

                              Thanks,

                              Comment


                              • #30
                                I'm going to download Bam files from the Project.
                                I see two links:
                                ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/data/
                                and
                                ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/phase1/data/

                                There are some overlapping files between the two links.

                                I would like to know which one I should use?

                                Thanks,

                                Comment

                                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
                                32 views
                                0 likes
                                Last Post seqadmin  
                                Started by seqadmin, 03-08-2024, 08:03 AM
                                0 responses
                                71 views
                                0 likes
                                Last Post seqadmin  
                                Started by seqadmin, 03-07-2024, 08:13 AM
                                0 responses
                                80 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