SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem generating consensus files EmmaB Bioinformatics 2 11-09-2012 01:30 PM
Questions about generating consensus sequences from the VCF file sulicon Bioinformatics 1 07-24-2012 06:47 AM
merging VCF files kjaja Bioinformatics 2 12-20-2011 06:50 AM
Generating SFF files Xterra 454 Pyrosequencing 8 10-31-2011 01:07 PM
Generating a BAM file from Illumina export files in CASAVa 1.7 nirav99 Bioinformatics 1 09-10-2010 01:20 AM

Reply
 
Thread Tools
Old 08-24-2011, 06:07 AM   #1
Brittd
Junior Member
 
Location: South Africa

Join Date: Apr 2011
Posts: 2
Default 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!
Brittd is offline   Reply With Quote
Old 06-05-2012, 11:41 PM   #2
mitsi
Junior Member
 
Location: india

Join Date: Mar 2012
Posts: 1
Default

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
mitsi is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 02:10 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO