Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Help figuring out variant call format (vcf)

    Hi,

    I know that there are quite a few threads already on this topic but I can't seem to find any specific to the questions I have. Anyway, I have a vcf file and can't figure out what is going on with certain lines. The file only contains biallelic SNPs and does not contain indels. Here's an example of three lines from my file:

    1089 5357 S1089_5357 G A 20 PASS . GT:ADP:GQ:PL 0/0:127,0:127:99:0,255,255 0/0:1
    27,10:137:99:0,52,255 0/0:43,0:43:99:0,129,255 0/1:127,19:146:99:245,0,255 0/0:127,0:127:99:0,255,255 0/
    0:126,0:126:99:0,255,255 0/0:127,0:127:99:0,255,255 0/0:127,0:127:99:0,255,255 0/1:127,12:139:96:14,0,255
    1089 5359 S1089_5359 - C 20 PASS . GT:ADP:GQ:PL 0/0:127,0:127:99:0,255,255 0/0:6
    6,0:66:99:0,198,255 0/0:13,0:13:99:0,39,255 0/0:126,0:126:99:0,255,255 0/0:106,0:106:99:0,255,255 0/0:30,0:3
    0:99:0,90,255 0/0:124,0:124:99:0,255,255 0/0:37,0:37:99:0,111,255 0/0:57,0:57:99:0,171,255
    1089 5360 S1089_5360 G A 20 PASS . GT:ADP:GQ:PL 0/1:127,18:145:99:212,0,255 0/0:1
    27,0:127:99:0,255,255 0/0:43,0:43:99:0,129,255 0/0:127,0:127:99:0,255,255 0/0:127,0:127:99:0,255,255 0/
    0:126,0:126:99:0,255,255 0/0:127,0:127:99:0,255,255 0/0:127,0:127:99:0,255,255 0/0:127,0:127:99:0,255,255

    My questions are as follows:
    1) For the second line, what does the "-" mean in the column for the reference allele?
    2) Should I be suspicious of SNPs that are found in close proximity to one-and-other? There are many instances in my file, such as in the above example, where numerous snps are positioned only single nts apart.
    3) Lastly, many of my snps have extremely large read depths (ca. 127-254 reads per locus per sample). Should I be concerned about these snps, since they could be caused from alignment of paralogous sequences?

    Any help on these questions would be great.

    Thanks,
    Bryan

  • #2
    1) ? Reference has nothing - so it's a C insertion (?) I have never seen this in the documentation or in output... Can't help - sorry.

    2) YES! That looks extremely suspicious. Check the raw alignment of the reads in IGV or similar

    3) Probably mapping errors due to repeat seq or similar - remove them - You should make a histogram of the mapping depth and see if anything looks weird - but mapping errors (with mismatches) causes "SNPs/indels" to arise - and you don't want them.

    My main suggestion: look at the raw data for these positions.

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Essential Discoveries and Tools in Epitranscriptomics
      by seqadmin


      The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist on Modified Bases...
      Yesterday, 07:01 AM
    • seqadmin
      Current Approaches to Protein Sequencing
      by seqadmin


      Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
      04-04-2024, 04:25 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 04-11-2024, 12:08 PM
    0 responses
    39 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 10:19 PM
    0 responses
    41 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 09:21 AM
    0 responses
    35 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-04-2024, 09:00 AM
    0 responses
    55 views
    0 likes
    Last Post seqadmin  
    Working...
    X