Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,279,276	C→G	F121L (TTC→TTG)	kbaZ →	putative tagatose‑1,6‑bisphosphate aldolase 1 chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,279,276	0	C	G	93.1%	82.4 / ‑3.8	29	F121L (TTC→TTG)	kbaZ	putative tagatose‑1,6‑bisphosphate aldolase 1 chaperone
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (13/14);  minor base T (2/0);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCT‑GTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGAAGCC  >  NC_000913/3279043‑3279522                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          ttccTGATGGAAGCACCCGCCAATCAGGTAGCCCGGTTCGGCGGTTATACCGAAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCtt                                                                                                                                                                                                                                           <  2:187629/246‑1 (MQ=255)            agccaCCTCCAATCAGGTAGACCAGTTCGGCGGCTATACCGGAATGCCGCCCGCCGATTTTCGCGGCGTTGTTTGTCAGCTCGCCGACTCGTTGGATTTCCCGCAGGATGCGTTGTTTCTGGGTGGTGACCCTCTAGGGCCAAACCGCTTGCCAAACCTGCCGGCCGCTCAGGCAATGGCCACTGCCGCTGTTTTGTTTAAAAGCCACTGTGCTGCAGGGTTTAAAACAAACCCACCttattggcccctgg                                                                                                                                                                                                                               >  2:19527/5‑239 (MQ=255)                  cTCCAATCAGGTAGACCAGTTCGGCGGGTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTTGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGCTGATTTGATTAAAGGCTACGTTGCGGCAGGAGTGAAAA‑AAATCCACCTTGTTTTCAACATGttccgtc                                                                                                                                                                                                                        >  2:189687/1‑245 (MQ=255)                            gTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGGTTTTCGCGGCTTTTTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGCTTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCCAAACCCGCCGGCCGCCCCGGCAATGGCCAATGCCCATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGCTcc                                                                                                                                                                                                             >  2:46840/1‑251 (MQ=255)                              aGACCAGTTCGGCGGTTATAACGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCCGCTCGCCGACCCGTTGCAATTCCCGCAGGATTCGGTGATTCTGGGTGGGGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCCCCTGCCATGGCCAATGCCGATGATTTGATTTAAAGCTACGCTGCGGCAGGACTTAAAA‑AAATACACCTTTCTTTCAGCACGTCCTTGTCAGGagctccc                                                                                                                                                                                                             >  2:204936/1‑244 (MQ=255)                               gACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGAt                                                                                                                                                                                                          <  1:178861/251‑1 (MQ=255)                                      gcggcgtTTATCCCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTa                                                                                                                                                                                                   <  1:172160/250‑1 (MQ=255)                                                     aaaggACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCGTAACTGATGACATCGTg                                                                                                                                                                                    <  2:191382/247‑1 (MQ=255)                                                      aaTGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTgg                                                                                                                                                                                   >  2:83843/1‑251 (MQ=255)                                                                                      gtCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGAATCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCg                                                                                                                                                   >  2:89784/1‑251 (MQ=255)                                                                                                     cGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCT‑GTc                                                                                                                                     <  2:200299/249‑1 (MQ=255)                                                                                                         gAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCT‑GTCTTGAac                                                                                                                               <  1:40457/251‑1 (MQ=255)                                                                                                                                 tCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCT‑GTCTTGAACACTTTGGCGAAGCCGATCTGGAg                                                                                                        <  1:149249/250‑1 (MQ=255)                                                                                                                                     ggtggtGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCT‑GTCTTGAACACTTTGGCGCAGCCGATCTggtggatgt                                                                                                   >  2:94082/1‑244 (MQ=255)                                                                                                                                                tCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCT‑GTCTTGAACACCTTTGCGAAGCCGATCTGGAggatttcctttgtaccc                                                                                        >  2:179646/1‑235 (MQ=255)                                                                                                                                                    gggCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCT‑GTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGt                                                                                    <  1:59751/251‑1 (MQ=255)                                                                                                                                                          aaCCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCT‑GTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGt                                                                              <  1:112095/251‑1 (MQ=255)                                                                                                                                                                        ccTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCCGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTTGT‑TTGAACACTTTGGCGCAGCCGGCTTGGGGGATGTCAATGGGACCGAAGTGCCGGGACCTGTCGgcccgc                                                                >  2:202733/1‑251 (MQ=255)                                                                                                                                                                                ccgCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCCTGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCT‑GTCTTGAAAACTTTGGCGAAGCCGATCTGGAGTATGGAATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAAcc                                                        >  2:119787/1‑251 (MQ=255)                                                                                                                                                                                       ggCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCT‑GTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCg                                                 <  1:38305/251‑1 (MQ=255)                                                                                                                                                                                                cAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCT‑GTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCgg                                        <  1:88562/251‑1 (MQ=255)                                                                                                                                                                                                  aTGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCT‑GTCTTGAACACTTTTGCGAAGCCCATCTGGAGTATGTCATTGGGACCGCAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCgggc                                      >  2:51015/1‑249 (MQ=255)                                                                                                                                                                                                  aTGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCT‑GTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGGTc                                      >  1:166469/1‑251 (MQ=255)                                                                                                                                                                                                        atgatTTGATTAAAAGCTACGTTGCGGAAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGACCGTTTGGCGAAAGTGGCGGAAGAAACCT‑GTCTTGAACACTTTGGCGAAGCCGATCTGGGGTCTTTCATTGGGACCGAATGGCCGGTACCTGACGGCGCGCAAGAAGCCTTAAGAGAGCTGGCCGTccaca                                  >  1:57703/1‑245 (MQ=255)                                                                                                                                                                                                              tgattAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATAGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCAGTCTGGCGAAAGTGGCGGAAGAAAACT‑GTGTTGAACAATTTGGCGAAGCCGATATGGAGTATGTAATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGGGCACCACGcc                              >  1:38462/1‑247 (MQ=255)                                                                                                                                                                                                                attAAAAGCTACGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCT‑GTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATgc                        <  1:25799/251‑1 (MQ=255)                                                                                                                                                                                                                           cGTTGCGGCAGGATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAATCCT‑GTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAATGCCGGGACCCGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGGGCACCACGCCGGATGccgcccgccccc             >  2:167702/1‑247 (MQ=255)                                                                                                                                                                                                                                     ggATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCT‑GTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGAAGc   <  1:5989/251‑1 (MQ=255)                                                                                                                                                                                                                                      gATTGAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCT‑GTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGAAGcc  <  2:82643/251‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          TTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAA‑AAATCCACCTTGATTGCAGCATGTCCT‑GTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCT‑GTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGAAGCC  >  NC_000913/3279043‑3279522

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