samtools SNP calling
1
0
Entering edit mode
22 months ago
evelyn ▴ 140

I am using samtools for SNP calling from .sorted.bam file:

bcftools mpileup -f ref.fa bwa.sorted.bam | bcftools view -Ov - > bwa.bcftools.vcf

VCF file gives REF calling but does not give ALT calling: A view from vcd file:

#CHROM  POS ID  REF ALT QUAL    FILTER  INFO    FORMAT  bwa.sorted.bam
EL1.0ch00   23812   .   A   <*> 0   .   DP=1;I16=0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;QS=0,0;MQ0F=0  PL  0,0,0

Any suggestions will be helpful.

SNP • 1.5k views
ADD COMMENT
2
Entering edit mode

The <*> has a meaning - see here: What does <*> mean in a vcf file?

You should be piping mpileup into bcftools call, by the way - see here: A: Genotyping with samtools/bcftools

ADD REPLY
2
Entering edit mode
22 months ago

Hello evelyn ,

variant calling using bcftools is a 2-step process:

  1. First check each position of the reference if it contains a potential variant, using bcftools mpileup
  2. In a second step have a closer look at the positions that have potential variants and check more criteria, using bcftools call

The second step is missing in your case.

So a complete command line could look like this:

$ bcftools mpileup -Ou -f ref.fa aln.bam | bcftools call -mv -o output.vcf

fin swimmer

ADD COMMENT
0
Entering edit mode

It worked, thanks!!!

ADD REPLY

Login before adding your answer.

Traffic: 1755 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