How to make Smatools to output both SNP and INDEL variants?
I used bwa-sw to align my Sanger reads with the reference genome. Then I used Samtools mpileup to do variant calling (see commands below). However, it did NOT show INDELs variants.
<command>
$samtools mpileup -uf Sorbil.fasta 1A.bam > 1A.pileup
$bcftools view -bvcg 1A.pileup > 1A.bcf
$bcftools view 1A.bcf > 1A.vcf
<result>
##fileformat=VCFv4.0
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT outdir/02A.bam
chr_8 50759327 . T C 4.77 . DP=1;AF1=0.9999;CI95=0.5,1;DP4=0,0,0,1;MQ=60 PL:GT:GQ 33,3,0:1/1:41
chr_8 50759520 . T C 10.4 . DP=1;AF1=0.9999;CI95=0.5,1;DP4=0,0,0,1;MQ=60 PL:GT:GQ 40,3,0:1/1:41
When I tried "pileup", instead of "mpileup", it only showed INDELs, no SNPs.
<command>
$samtools pileup -vcf Sorbil.fasta 1A.bam > 1A.pileup
<result>
chr_8 50759483 a A 30 0 60 1 ,+1t M
chr_8 50759483 * +T/+T 40 0 60 1 +T * 1 0 0 0 0
I used bwa-sw to align my Sanger reads with the reference genome. Then I used Samtools mpileup to do variant calling (see commands below). However, it did NOT show INDELs variants.
<command>
$samtools mpileup -uf Sorbil.fasta 1A.bam > 1A.pileup
$bcftools view -bvcg 1A.pileup > 1A.bcf
$bcftools view 1A.bcf > 1A.vcf
<result>
##fileformat=VCFv4.0
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT outdir/02A.bam
chr_8 50759327 . T C 4.77 . DP=1;AF1=0.9999;CI95=0.5,1;DP4=0,0,0,1;MQ=60 PL:GT:GQ 33,3,0:1/1:41
chr_8 50759520 . T C 10.4 . DP=1;AF1=0.9999;CI95=0.5,1;DP4=0,0,0,1;MQ=60 PL:GT:GQ 40,3,0:1/1:41
When I tried "pileup", instead of "mpileup", it only showed INDELs, no SNPs.
<command>
$samtools pileup -vcf Sorbil.fasta 1A.bam > 1A.pileup
<result>
chr_8 50759483 a A 30 0 60 1 ,+1t M
chr_8 50759483 * +T/+T 40 0 60 1 +T * 1 0 0 0 0
Comment