Calling CNV using bcftools
1
0
Entering edit mode
2.8 years ago
Raheleh ▴ 210

Hi, I used VarScan2 to call variants from whole exome sequencing data. The vcf file output looks like this:

#CHROM  POS ID  REF ALT QUAL    FILTER  INFO    FORMAT  NORMAL  TUMOR
chr1    972251  .   A   AC  .   PASS    DP=173;SOMATIC;SS=2;SSC=75;GPV=1;SPV=2.5611e-08 GT:GQ:DP:RD:AD:FREQ:DP4 0/0:.:65:63:0:0%:36,27,0,0  0/1:.:108:74:33:30.84%:46,28,18,15

I ran this command of bcftools:

bcftools cnv -c NORMAL -s TUMOR -o outdir/ -p 0 samples_concated.vcf

to get copy number variations, but it returned empty files. Can anyone help me to fix it? Is there any wrong with my vcf file or with my command? Than you!

bcftools CNV vcf • 1.8k views
ADD COMMENT
0
Entering edit mode
2.8 years ago

$ bcftools cnv

About: Copy number variation caller, requires Illumina's B-allele frequency (BAF) and Log R

Ratio intensity (LRR). The HMM considers the following copy number states: CN 2

both annotations are missing in your VCF.

ADD COMMENT
1
Entering edit mode

@Pierre Lindenbaum Dear Pierre many thanks for your reply. How can I get these two annotations for my data. Should I use another variant caller instead of VarScan2? If so, what do you recommend?

ADD REPLY
0
Entering edit mode
ADD REPLY

Login before adding your answer.

Traffic: 2479 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6