Hello:
I am using BWA + SAMtools/BCFtools to reference map assemble illumina reads and call SNPs in a bacterial genome respectively. I use the command:
which is supposed to filter the initial pileup and call the variant positions using bayesian inference similar to the one used in MAQ.
I'm not an expert in statistics, but I would like to know what are the priors and posteriors calculated. Thanks very much!
I am using BWA + SAMtools/BCFtools to reference map assemble illumina reads and call SNPs in a bacterial genome respectively. I use the command:
Code:
bcftools view -bvcg input.bcf > output.bcf
I'm not an expert in statistics, but I would like to know what are the priors and posteriors calculated. Thanks very much!
Comment