Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • Brittd
    Junior Member
    • Apr 2011
    • 2

    Generating vcf files with samtools

    I am new to next gen sequence analysis and I was wondering if anyone could help me with some questions that I have regarding the use of SAMtools to generate vcf files:

    I have used the following command lines to generate a vcf file:
    samtools mpileup -gf GRCh37.fa sample.bam > mpileup_sample
    bcftools view -bvcg mpileup_sample > bcf_sample
    bcftools view bcf_sample > vcf_sample
    vcfutils.pl varFilter -D 500 vcf_sample > vcf_filtered_sample

    In the resulting vcf file, the columns "ID" and "FILTER" are empty, is there any way that I can obtain values for these columns?

    Also, I would like to filter out all those samples that have a read depth less than 10 and Quality scores less than 20. Is there a way that I can do this in SAMtools?

    Lastly, I am unclear as to what the last command line (vcfutils.pl varFilter -D 500 vcf_sample > vcf_filtered_sample) does.

    Any help or suggested reading would be much appreciated!
  • mitsi
    Junior Member
    • Mar 2012
    • 1

    #2
    The ID column represents the ID of dbsnp database, i.e if the variant is reported in dbsnp then its ID will be given there, and u can reduce the read depth by setting the option -D 10 in vcfutils command.
    vcfutil command converts the binary file bcf to a more readable format of variants and also we can use various filters to get more novel variants

    Comment

    Latest Articles

    Collapse

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, 06-09-2026, 11:58 AM
    0 responses
    30 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-05-2026, 10:09 AM
    0 responses
    38 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-04-2026, 08:59 AM
    0 responses
    42 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-02-2026, 12:03 PM
    0 responses
    64 views
    0 reactions
    Last Post SEQadmin2  
    Working...