Samtools vcf file filtering
Entering edit mode
7 days ago

Hi, I am using samtools inbuilt to filter my vcf file. So before setting up my parameters, I wanted to know what is the default value set in the script for these options:


-Q INT minimum RMS mapping quality for SNPs [$opts{Q}]

     -d INT    minimum read depth [$opts{d}]

   -D INT    maximum read depth [$opts{D}]

     -a INT    minimum number of alternate bases [$opts{a}]

   -w INT    SNP within INT bp around a gap to be filtered [$opts{w}]

   -W INT    window size for filtering adjacent gaps [$opts{W}]

   -1 FLOAT  min P-value for strand bias (given PV4) [$opts{1}]

   -2 FLOAT  min P-value for baseQ bias [$opts{2}]

   -3 FLOAT  min P-value for mapQ bias [$opts{3}]

   -4 FLOAT  min P-value for end distance bias [$opts{4}]

     -e FLOAT  min P-value for HWE (plus F<0) [$opts{e}]

   -p        print filtered variants

Thank you!

samtools RNAseq sequencing • 54 views

Login before adding your answer.

Traffic: 2394 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6