Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • How to filter out SNPs with Minor allele frequency less than 5% in VCF file

    Hello,
    I would like to filter out SNPs with allele frequencies in VCF file from multi-sample variant calling.
    I tried VCFtools with the following commands, but it doesn't work.
    vcftools --vcf input.vcf --recode --recode-INFO-all --stdout --maf 0.05 --max-maf 1 --out outpur.vcf
    Here is the vcf file format generated from samtools mpileup.
    #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT 788 712 908 544 541 581 801 381 479 487 379 570
    Chr01 11356 . C A 124 . DP=88;VDB=1.624097e-04;RPB=-2.346857e+00;AF1=0.1553;AC1=3;DP4=75,3,10,0;MQ=30;FQ=124;PV4=1,1,1,1
    GT:PLP:GQ 0/0:0,0,0:0:5 0/0:0,6,59:2:11 0/0:0,0,0:0:5 0/1:28,0,123:8:24 1/1:157,24,0:8:14 0/0:0,39,221:13:43 0/0:0,12,85:4
    :17 0/0:0,36,164:12:40 0/0:0,39,147:13:43 0/0:0,9,71:3:14 0/0:0,48,163:16:52 0/0:0,27,126:9:31
    Chr01 11384 . T C 15.9 . DP=123;VDB=5.617680e-03;RPB=1.820871e+00;AF1=0.04721;AC1=1;DP4=114,3,5,0;MQ=31;FQ=16.1;PV4=1,1,0.05,1 GT:PLP:GQ 0/0:0,12,79:4:22 0/0:0,6,60:2:16 0/0:0,0,0:0:10 0/0:0,30,163:10:40 0/0:0,27,159:9:37 0/0:0,60,244:20:70 0/0:0,12,85:4:22 0/0:0,48,163:16:58 0/1:56,0,116:15:46 0/0:0,18,122:6:28 0/0:0,69,185:23:79 0/0:0,39,135:13:49
    I wonder if anyone could help me
    Thanks

  • #2
    The DP4 part spells out how many reads aligned to reference and alternate allele. So you can parse based on that. You could do that with some finagling in Excel, if you had to.

    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, Yesterday, 06:37 PM
    0 responses
    8 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Yesterday, 06:07 PM
    0 responses
    8 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-22-2024, 10:03 AM
    0 responses
    49 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-21-2024, 07:32 AM
    0 responses
    66 views
    0 likes
    Last Post seqadmin  
    Working...
    X