Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,295,693	C→T	8.5%	A18T (GCG→ACG)	ccmH ←	heme lyase, CcmH subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,295,693	0	C	T	8.5%	146.0 / 3.7	47	A18T (GCG→ACG)	ccmH	heme lyase, CcmH subunit
Reads supporting (aligned to +/- strand):  ref base C (26/17);  new base T (2/2);  total (28/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

TCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGC  >  W3110S.gb/2295626‑2295752                                                                    |                                                            tcttcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAg                                                          >  1:182545/1‑71 (MQ=255)  cttcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgc                                                         >  1:409491/1‑71 (MQ=255)  cttcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgc                                                         >  1:3596346/1‑71 (MQ=255)  cttcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgc                                                         >  1:1748375/1‑71 (MQ=255)  cttcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgc                                                         >  1:1845508/1‑71 (MQ=255)  cttcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgc                                                         >  1:3282898/1‑71 (MQ=255)  cttcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgc                                                         >  1:2923199/1‑71 (MQ=255)  cttcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAg                                                          >  1:1241724/1‑70 (MQ=255)  cttcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAg                                                          >  1:1724975/1‑70 (MQ=255)  cttcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAg                                                          >  1:282543/1‑70 (MQ=255)  ctccAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAgc                                                         >  1:1925166/4‑71 (MQ=255)     cAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgcg                                                        <  1:1331797/69‑1 (MQ=255)     cAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgcg                                                        <  1:2005012/69‑1 (MQ=255)     cAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCt                                                      >  1:994533/1‑71 (MQ=255)     cAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCt                                                      >  1:3395416/1‑71 (MQ=255)              gACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCg                                                >  1:452404/1‑68 (MQ=255)              gACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCg                                                >  1:3460926/1‑68 (MQ=255)               aCGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGGAACACATCGATGGTCGCCAGCGCTGAGCCGgaga                                            <  1:2560681/71‑1 (MQ=255)               aCGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGgaga                                            <  1:3592446/71‑1 (MQ=255)               aCGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGgaga                                            <  1:1887326/71‑1 (MQ=255)               aCGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGgag                                             <  1:3503518/70‑1 (MQ=255)               aCGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGgag                                             <  1:1895298/70‑1 (MQ=255)                cGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCTACACATCGATGGTCGTCAGCGCTGAGCCGCAGat                                           <  1:3670997/71‑1 (MQ=255)                cGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGat                                           <  1:3032507/71‑1 (MQ=255)                      tgttgtTCCTGTGCTTCATCTTTAAACTTCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATcatcag                                     >  1:1313046/1‑71 (MQ=255)                       gttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATcatcag                                     >  1:537544/1‑70 (MQ=255)                       gttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATcatcag                                     >  1:254655/1‑70 (MQ=255)                       gttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATcatcag                                     >  1:2580090/1‑70 (MQ=255)                           ttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATcatcagcatc                                 <  1:2134204/70‑1 (MQ=255)                             ccTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATcatcagcatcag                               <  1:2177002/70‑1 (MQ=255)                             ccTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATcatcagcatcag                               <  1:1632927/70‑1 (MQ=255)                               tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCAc                            >  1:2971511/1‑71 (MQ=255)                               tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCAc                            >  1:1908924/1‑71 (MQ=255)                                      cATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCaa                      >  1:103507/1‑70 (MQ=255)                                      cATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCa                       >  1:2261908/1‑69 (MQ=255)                                       aTCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGatcat                                        >  1:426166/1‑51 (MQ=255)                                       aTCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGatcat                                        >  1:271266/1‑51 (MQ=255)                                          tttAAACTGCAACACATCGATGGGCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                  >  1:694436/1‑70 (MQ=255)                                            tAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                  <  1:3586605/68‑1 (MQ=255)                                            tAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                  <  1:2751961/68‑1 (MQ=255)                                            tAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                  <  1:646259/68‑1 (MQ=255)                                                      cacaTCGATGGTCGCCAGCGCTGAGCCGGAGATCAtcagcatcagca                             <  1:2971696/47‑1 (MQ=255)                                                      cacaTCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTc            <  1:3487539/64‑1 (MQ=255)                                                      cacaTCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTc            <  1:2730348/64‑1 (MQ=255)                                                      cacaTCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTc            <  1:573326/64‑1 (MQ=255)                                                        caTCGATGGTCGCCAGCGCTGAGCCGGAGATCATTAGCATCAGCACGCCCAATAAAAACCTCATtgtg      >  1:2966605/1‑68 (MQ=255)                                                         aTCGATGGTCGTCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGc  >  1:2080044/1‑71 (MQ=255)                                                                    |                                                            TCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGC  >  W3110S.gb/2295626‑2295752

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑