Question: samtools mpileup to bcftools mpileup
2
gravatar for Ric
16 months ago by
Ric290
Australia
Ric290 wrote:

I got

samtools mpileup option `u` is functional, but deprecated. Please switch to using bcftools mpileup in future.

How could I convert the following commands to bcftools mpileup ?

samtools mpileup -uf QMg-NbQ3P-RN.fasta aln_P1O1.sorted.bam | bcftools view -bvcg - > P1O1_var.raw.bcf
bcftools view P1O1_var.raw.bcf | vcfutils.pl varFilter -D100 > P1O1_var.flt.vcf

Thank you in advance.

snp samtools bcftools • 4.5k views
ADD COMMENTlink modified 16 months ago by zx87549.1k • written 16 months ago by Ric290
1
gravatar for finswimmer
16 months ago by
finswimmer13k
Germany
finswimmer13k wrote:

Hello Ric ,

first make sure you have the latest version of samtools and bcftools.

$ bcftools mpileup -Ou QMg-NbQ3P-RN.fasta aln_P1O1.sorted.bam | bcftools call -Ou -mv  | bcftools norm -Ou -f QMg-NbQ3P-RN.fasta | bcftools view -e 'FORMAT/DP > 100' > P1O1_var.flt.vcf
  • The bcftools view command you showed was only used for bcftools version < 1 and is replaced by bcftools call
  • I recommend to normalize your indel variants, that's why I added a bcftools norm step
  • For removing variants with read depth greater 100 I used bcftools view

fin swimmer

ADD COMMENTlink written 16 months ago by finswimmer13k
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: 1576 users visited in the last hour