A weird variant call from VarScan2 mpileup2snp

Hi there,

I used VarScan.v2.4.3.jar mpileup2snp to call mutations with -min-coverage=100 and --min-var-freq= 0.01 and I got some weird variants. For example, one mutation for the output is :

21 42848562 . G C . PASS ADP=866;WT=0;HET=1;HOM=0;NC=0 GT:GQ:SDP:DP:RD:AD:FREQ:PVAL:RBQ:ABQ:RDF:RDR:ADF:ADR 0/1:0:866:866:6:5:45.45%:9.8E-1:36:34:4:2:3:2

Usually, DP=RD+AD=RDF+RDR+ADF+ADR and FREQ=AD/DP=(ADF+ADR)/(RDF+RDF+ADF+ADR), but in this case RD+AD=11, not equal to DP=866. Frequency should be 5/866=0.6%, also not equal to 45.45%.

How to explain this variant? Which value should I choose to filter? Does anyone have any idea?

Thanks so much! Rui

varscan2 mpileup2snp variant calling

• 1,796 views • link

written 8.2 years ago by rayliu0115 •

My first option would be to visually inspect the BAM with IGV or similar and see if I can make sense of the numbers.

• 1 views • link

written 8.2 years ago by bernatgel

0 answers

No answers yet.

More posts like this

htsjdk.tribble.TribbleException: The provided VCF file is malformed at approximately line number 92…
written 5.4 years ago by DareDevil

I run variantsToTable module of gatk as follows: gatk VariantsToTable -V 001.vcf -F CHROM -F POS -F ID -F REF -F ALT -F QUAL -F …
VarScan: How is the DP score in a .vcf calculated?
written 7.5 years ago by James Reeve

I'm trying to understand the formats of different SNP callers. I have made test bam files with a subset of two stickleback WGS samples. I …
VarScan.v2.3.9 filter error p-value
written 7.9 years ago 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 . …
The different results for the same tumor sample between VarScan2 somatic and mpileup2snp
written 8.3 years ago by rayliu0115 •

Dear all, I have a question about the different results for one same tumor sample between VarScan2 somatic and mpileup2snp. I used VarScan2 somatic to …
vcf-annotate not adding annotations
written 9.5 years ago by wes3985 •

I would like to annotate a vcf file with the output of a filter that detects false positives. I know a very similar question to …
VCFtools filter - genotype quality score changes
written 9.9 years ago by AW

Hi, I ran the following filtering command vcftools --vcf test.vcf --out test.test.vcf --remove-indels --remove-filtered-all --max-alleles 2 --max-alleles 2 --recode --recode-INFO-all However, when I compare the …
Varscan variant calll error on pacbio data
written 10.2 years ago by mark.rose •

Hello I'm trying to call variants on an alignment of pacbio data to a reference thusly: blasr sample.bas.h5 ref.fa -sam -clipping soft -nproc 4 -minSubreadLength …
vcf format genotypes called
written 10.4 years ago by AW

Hi, I ran VARSCAN2 mpileup2snp to call SNPs and generated a vcf file. However, I noticed a strange thing about the NC flag. This is …
Fill in missing VCF information with bcftools or vcftools or Vcf.pm
written 11.8 years ago by Lee Katz

Hi Biostars, I was wondering how I might backfill some information into the VCF format. Here are the first three lines of a VCF pooled …
Questions In Output Of Varscan 2.3.2
written 13.6 years ago by rahilsethi •

<p>Hi,</p> <p>I have some questions in the output of mpileup2snp and mpileup2indel. For substitutions obtained from mpileup2snp when viewed in VCF format, whether they have …