Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

Example:

1. samtools mpileup -uf Euf reference.fna aln1.bam aln2.bam | bcftools view -bvcg - > var.raw.bcf

...

aln1.bam, aln2.bam : BAM files containing alignment results. You can use 1 or more alignment flies at a time.  Note that as of late 2011, the new BAQ filter seems to aggressively remove SNPs unless you "extend" it with the "-E" option.

2. bcftools view var.raw.bcf | vcfutils.pl varFilter -D10 > var.filtered.vcf

...