Filter VCF file
Hi
How do I filter the vcf file using below filters:
- number of reads per base is between 5 and 75
- base quality > or = to 30
- mapping quality >= 60
- variant quality >= 90
- distance of adjacent variant >= 5
snp
next-gen
• 4,265 views
•
link
updated
by
Ram
•
written
by
sukesh1411 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
VCFfiltering
written by Peter •Hi, A) I have a large set of samples diploid, tetraploid, and hexaploid mixed. I have used FreeBayes to build a VCF file. I want …
-
design matrix considering treatment, timepoints and batches
written by Iván •Hi, I'm analysing a microarray with 59 samples distributed across the following groups and timepoints: SampleName Group Batch Timepoint sample1 control_30 1 30 sample2 control_30 …
-
how extract our desirable data from two table?
written by starI like to find overlap regions (even with 1 bp) along with the number of counts (fragments) which have overlaps together. Input: A.bed chr1 10 …
-
Filtering VCF file
written by Prakki RamaHi all, I am trying to filter VCF file based on a publication with the following criteria. Base calls derived from these two Vcf files …
-
Why are GC% per base important in quality control of reads?
written by c.claridoHello, In quality control of reads, why do we look at the GC% per base position? I have the following result Gem. lengtes: 75 Max. …
-
VarScan.v2.3.9 filter error p-value
written by niggoK •Hi, I am getting an error message when I use the VarScan filter option. Error: Parsing Exception on line: chr1 104033 . T C . …
-
What is a good Trimmomatic parameter?
written by Arindam GhoshI am trying to clean RNAseq reads with TRIMMOMATIC using the parameters as: java -jar trimmomatic-0.36.jar PE -trimlog file_trim_log input_1.fastq.gz input_2.fastq.gz output_1P_clean.fq.gz output_1U_clean.fq.gz output_2P_clean.fq.gz output_2U_clean.fq.gz …
-
Changing samtools mpileup I16 calculations?
written by emilywynn6 •Hello, I am trying to filter some vcf files from mpileup by using the average base quality of reference bases and the average base quality …
-
rna seq data trim by trimmomatic
written by zh.khodadadi •Hi every body I am going to trim rna seq single end data(Illumina HiSeq 2000) by trimmomatics According to text: ***sets of samples were trimmed …
-
Need help understanding read depth for somaticsniper vcf's
written by JordanHi, I'm having trouble understanding read depth in vcf's called by `somatic-sniper`. Here is a sample vcf file: ``` ##fileformat=VCFv4.1 ##fileDate=20140428 ##phasing=none ##reference=human_hg19.fa ##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype"> ##FORMAT=<ID=IGT,Number=1,Type=String,Description="Genotype …
http://snpeff.sourceforge.net/SnpSift.html#filter
and it will be something like
cat variants.vcf | java -jar SnpSift.jar filter " ( QUAL >= 30 ) && ( (DP >= 5) && (DP <= 75))" > filtered.vcf
Hi sir
Thank you for your reply. I would like to know what are terms used for mapping quality, variant quality and distance of adjacent variant in vcf file?
Do you mean information about vcf specification?
http://www.1000genomes.org/wiki/Analysis/Variant%20Call%20Format/vcf-variant-call-format-version-40/
http://samtools.github.io/hts-specs/VCFv4.2.pdf
Hi sir
Can i know whether these terminologies are correct ?
DP - no of reads per base QUAL - variant quality MQ - mapping quality BQ - base quality and i cudnt find terminology for distance of adjacent variant.
Yes you are right from the header you can read
INFO=<id=dp,number=1,type=integer,description="total depth"="">
total number of sequence reads overlapping that position MQ : RMS mapping quality, similar to the version in the INFO field. (Integer) BQ : RMS base quality at this position
http://seqanswers.com/forums/showthread.php?t=9345
you mean between the reference and sample ?!
Hi sir
Thank you for your reply. I could not find adjacent variant term in .vcf file.