CallSNPs.py

From wiki
Revision as of 16:41, 1 April 2016 by Rf (talk | contribs) (Created page with "This script is part of "script_tools" modules. It requires an alignment in the shape of the base reference sequence and the bam files of the short read aligned to that sequen...")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

This script is part of "script_tools" modules.

It requires an alignment in the shape of the base reference sequence and the bam files of the short read aligned to that sequence (run previously, maybe with bwa or one of the bowtie programs).

Stages

1. samtools mpileup and bcftools view -cg "mpileup" collects summary information from the input BAMs and then computes and stores the likelihood of data given each possible genotype. bcftools