NEW VERSION: Commands to use samtools with a bam file, input.bam,Latest version on Fourierseq: 0.1.18 installed 9/26/11 by SPHS
Useful samtools utilies:
1. samtools idxstats : This tool will provide statistics about how many reads have aligned to each sequence/chromosome in the reference genome. The input bam file must be sorted and indexed.
Code Block |
---|
samtools idxstats <in.bam>
|
2. samtools flagstat : Simple stats about how many reads mapped to the reference, how many reads were paired properly etc. The input bam file must be sorted and indexed.
Code Block |
---|
samtools flagstat <in.bam>
|
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
...