![]() |
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
what is DP4 in VCF (Variant Call Format) and AF1? | gcrdb | Bioinformatics | 8 | 01-20-2015 09:37 AM |
strand bias variant call | yinshe | Bioinformatics | 2 | 10-10-2011 06:01 AM |
Variant call using BAM from Tophat | emilyjia2000 | Bioinformatics | 5 | 08-01-2011 04:37 PM |
ERANGE to call peaks from ChIP-Seq data | sonja | Bioinformatics | 3 | 08-31-2010 08:47 AM |
Using Galaxy to Call Peaks? | Giles | Bioinformatics | 4 | 02-05-2010 10:43 AM |
![]() |
|
Thread Tools |
![]() |
#1 |
Member
Location: Toronto, ON Join Date: Dec 2009
Posts: 16
|
![]()
I have reference mapped paired end illumina reads and called variants using BWA and Samtools respectively. The resulting vcf was treated to remove high coverage SNPs with
Code:
vcfutils.pl varFilter -D30 Code:
'($3=="*"&&$6>=50)||($3!="*"&&$6>=20)' Last edited by elfuser; 07-31-2012 at 12:02 PM. |
![]() |
![]() |
![]() |
Tags |
samtools vcf bwa |
Thread Tools | |
|
|