Question: Short Variants calling error when using bcftools mpileup
0
gravatar for mike229lin
2.6 years ago by
mike229lin20
mike229lin20 wrote:

Hello,

I want to call gvcf with bcftools mpileup containing variant site and reference homologous site blocks The command line shows below

nohup sh -c "bcftools mpileup -C50 -q25 -g 4 -t DP -f reference.fa sample.bam | bcftools call -c | vcfutils.pl varFilter -d4 -D28 > output.g.vcf" &

I also want to filtrate the variant site out which is 2 times of and one third of the average coverage of each sample. However, it only printed the header out and didn't show any variant site or non variant blocks in my output. How should I edit my command line to manipulate vcf files I want?

Thanks.

snp alignment genome • 691 views
ADD COMMENTlink modified 2.6 years ago by finswimmer14k • written 2.6 years ago by mike229lin20

Hello mike229lin,

you should first check the output after bcftools call -c, if it contains what you expected before starting filtering. Having the header but no variants in the final output sounds like there are nor variants that fulfill the criteria.

Furthermore I wonder if this is correct:

-t DP

The manual says:

-t, --targets REG[,...] similar to -r but streams rather than index-jumps
-r, --regions REG[,...] comma separated list of regions in which pileup is generated

fin swimmer

ADD REPLYlink written 2.6 years ago by finswimmer14k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2214 users visited in the last hour
_