PairwiseCallSNPs

From wiki
Revision as of 11:07, 1 December 2017 by Rf (talk | contribs) (Created page with "= Introduction = The following pipeline evolved from the callSNPs.py script developed by Miguel Pinheiro, following collaboration with Souad Alateah and operating on her myco...")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

Introduction

The following pipeline evolved from the callSNPs.py script developed by Miguel Pinheiro, following collaboration with Souad Alateah and operating on her mycobateria samples: M. Absessus and M. Tuberculosis.

The aim of the evolution were as follows:

  • High resolution, pairwise SNP calling between the samples.
  • Parallel launching of the program using Gridengine's DRMAA python library.
  • Taking advantage of recent (2017, version 1.5) improvements in the samtools/bcftools toolkit.

These aims do not however mean that earlier analyses are invalid, but rather they were less precise with lower resolution. For example one of the key problems is the high similarilty between the samples. Often they are almost identical, and the aim here was to refine and improve the method so that these samples could be more distinguishable in terms of their SNP mutations.