Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Exome sequencing analysis problem

    Hello All,

    I am new in next generation sequencing analysis. I am working on a vcf file which performed by GATK tools. I have a question about GT (genotype) for ChrY and ChrM. the data looks lke:
    GT : GQ : DP : PL 0/1:8.82:96:4088,0,9
    GT : GQ : DP : PL 1/1:38.92:1:390,38,0
    As my understanding, GT means genotype. Since there is only one copy for chrY or chrM, why the data showed 0/1 (heterozygous) or 1/1 (homozygous)? Or do I misunderstand the meaning of "GT" ?

    Thanks,
    Albert
    Last edited by albertyu; 11-15-2012, 04:47 PM.

  • #2
    Well for mitochondria it could be :



    As for the Y - are they in pseudoautosomal regions? Or perhaps badly genotyped indels?

    Comment


    • #3
      Originally posted by Bukowski View Post
      Well for mitochondria it could be :



      As for the Y - are they in pseudoautosomal regions? Or perhaps badly genotyped indels?

      Thanks for your reply. I checked the possibility of the pseudoautosomal genes. It seems that most of variants are not at the pseudoautosomal PAR1 and PAR2 genes location. And many of them are substitutions not indels. I would think it's because mis-mapping of BWA. ???

      Actually I originally thought this is a popular problem in Exome sequencing data since I saw several whole exome sequencing vcf data containing chrY with GT:1/0 or 1/1.

      Comment


      • #4
        Originally posted by albertyu View Post
        Thanks for your reply. I checked the possibility of the pseudoautosomal genes. It seems that most of variants are not at the pseudoautosomal PAR1 and PAR2 genes location. And many of them are substitutions not indels. I would think it's because mis-mapping of BWA. ???

        Actually I originally thought this is a popular problem in Exome sequencing data since I saw several whole exome sequencing vcf data containing chrY with GT:1/0 or 1/1.
        I just had a look at some exome data, and I think there is definitely some mismapping involved as most females in my trios a small number of chrY variations. As I suspected though in my data, most of the het calls are small indels rather than SNPs.

        I don't see what is wrong with a 1/1 call though, in Y that should just mean 'base is different to reference' surely?

        Comment


        • #5
          Originally posted by Bukowski View Post
          I just had a look at some exome data, and I think there is definitely some mismapping involved as most females in my trios a small number of chrY variations. As I suspected though in my data, most of the het calls are small indels rather than SNPs.

          I don't see what is wrong with a 1/1 call though, in Y that should just mean 'base is different to reference' surely?
          Indeed, all my female samples have chrY variations.
          From my data performed by GATK, it separated substitutions and indels. So I originally thought those are all indels and finally I realized there are a lot of SNPs in the middle of vcf files.
          If they are mostly mismapping, I am afraid to use these data. So I would filter out those data at chrY.

          Comment


          • #6
            Check out section 6 on this page:



            It explains how the genotypes on the sex chromosomes are called in GATK

            Comment


            • #7
              [QUOTE=albertyu;89712]Indeed, all my female samples have chrY variations.
              From my data performed by GATK, it separated substitutions and indels. So I originally thought those are all indels and finally I realized there are a lot of SNPs in the middle of vcf files.
              If they are mostly mismapping, I am afraid to use these data. So I would filter out those data at chrY.[/QUOTEI

              In any exome sequencing you will get some "false SNP" that have different sources: sequencing errors, aligning pseudogenes, ... That is why in an exome sequencing you can see chromosome Y in women or heterozygous SNPs in men. You can eliminate these "errors" by increasing your level of filtering in the alignement process but then the price you pay is with a tougher filtering is that you can lost true information too.

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Strategies for Sequencing Challenging Samples
                by seqadmin


                Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                03-22-2024, 06:39 AM
              • 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

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, 03-27-2024, 06:37 PM
              0 responses
              13 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-27-2024, 06:07 PM
              0 responses
              11 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-22-2024, 10:03 AM
              0 responses
              53 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-21-2024, 07:32 AM
              0 responses
              69 views
              0 likes
              Last Post seqadmin  
              Working...
              X