Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,279,461	A→G	H183R (CAT→CGT)	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,461	0	A	G	100.0%	59.6 / NA	24	H183R (CAT→CGT)	kbaZ	putative tagatose‑1,6‑bisphosphate aldolase 1 chaperone
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/12);  total (12/12)

ACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGAAGCCCATCGTCACGCCTTTGAAAAGCAAGGTTTGAATGCCATCTGGCCACGCATCATTGCCCTGGTGGTTCAACCCGGCGTCGAATTCGATCACACCAACGTTATTGATTATCAGCCCGCCAAAGCGAGCGCCTTAAGCCAGATGGT  >  NC_000913/3279260‑3279665                                                                                                                                                                                                          |                                                                                                                                                                                                             aCGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCGTGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCgcgc                                                                                                                                                              <  1:23156/250‑1 (MQ=255)                                                 gTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCGTGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGAAGCCCATCGTCACGCCTTTGAAAAGCAAGGTTTGAATGc                                                                                                              <  1:17057/250‑1 (MQ=255)                                                                   ccTTAACTGATGACATCGAGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACATGTCTTGAAAACTTTGGCGAAGCCGATCTTGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCGTGCAACCTTAAGCGAGCTAGCGGACACCACGCCGGATGCCGCCCGCGCCCCGACGGAAGCCCATCGGCACGCCTTTGAAAAGCCAGGGTTGAATGCCATCTGGTCACGCATCAtt                                                                                           >  1:403911/1‑251 (MQ=255)                                                                         cTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACCCTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCGTGANACCNNNNNCgnnnnnn                                                                                                                                                                                          >  2:298602/1‑144 (MQ=255)                                                                            atgaCATCGTTGCTGAACGCGCCGCCCTTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTTGCGAAGCCGAGCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCGTGAAAcc                                                                                                                                                                                                       <  1:368574/134‑1 (MQ=255)                                                                                    gTGGCTGAACGCGCCGCCCGTCGGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTTGCGAAGCCGATCTGGAGTATGTCATTGGTACCGGAGTGCCGGTACCTGGCGGCGCGCGTTAAACCGTAAGCGAGCTGTCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGAAGCCCATCGTCa                                                                                                                                         <  1:205752/188‑1 (MQ=255)                                                                                      ggCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTGGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCGTGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATgn                                                                                                                                                                       <  1:100251/156‑2 (MQ=255)                                                                                       ncTGAACGCGCCGCCCGTCTGGCGAAAGCGGCGGAAGAAACCTGTCTTGAACACTTCGGGGCAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTCCCTCGCGGCGCGCGTGAAACCTTAAGCGAGCTTGCGGTCCCCACGCCGGATGCCCCCAGCGCCACGCCGGAAGCCCATCGTCACGCCTTTGAAAAGCAAGGTTTGAATGCCATCTGGCCACGCATCATTACACTGGTGGTTCAACCCgg                                                                       >  1:50454/2‑251 (MQ=255)                                                                                           aaCGCGCCGCCCGTCTGGCGCAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCGTGAAAGCTTCAGCGAGCTGACGGTCACCACGCCGGCCGCCGCCCGCGCCACGCTGGAAGCCAATCGTAACGCCTTTGAAAAGCAAGGTTTGCATGCCATCTGGCCACGCATCCTTGCCCTGGTGGTTCAACCCGGCGTc                                                                   >  1:490670/1‑251 (MQ=255)                                                                                             cgcgCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCGGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCGTGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGccgcccgc                                                                                                                                                                <  1:82962/156‑1 (MQ=255)                                                                                              gcgcCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTTCCGGTACCTGGCGGCGCGCGTGAAACCTTAAGCGAGCTGGCGGGNACCNNNNNggnnnnnnnnnn                                                                                                                                                                 >  2:303043/1‑144 (MQ=255)                                                                                              gcgcCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCGTGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGAAGCCCATCGTCACGCCTTTGAAAAGCAAGGTTTGAATGCCATCTggn                                                                                                      <  1:303043/213‑2 (MQ=255)                                                                                                 ccgccCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTCAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGCAGTGCCGGTACCTGGCGGCGCGCGTGACACCTTAAGCGCGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGAAGCCCCTCGTCACGCCTTTGaaaa                                                                                                                             >  1:252406/1‑187 (MQ=255)                                                                                                   gccCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCGTGAAACCTTAAGCGAGCTGCCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGAAGcc                                                                                                                                                 >  1:271937/1‑165 (MQ=255)                                                                                                           ggCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCGTGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGAAGCCCATCGTCACGCCTTTGAAAAGCAAGGTTTGAATGCCATCTGGCCACGCATCATTGCCCTGGTGGTTCAACCCGGCGTCGAATTCGATCacacn                                                    <  1:298602/250‑2 (MQ=255)                                                                                                                     gcggacgcaacatGTCTTGAACACTTTGGCGAAGCAGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCGTGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGACGCCCATCGTCACGCCTTTGAAAAGCAAGGTTTGAATGc                                                                                                              >  1:389301/13‑182 (MQ=255)                                                                                                                                                   naaGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCGTGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGAAGCCCATCGTCACGCCTTTGAAAAGCAAGGTTTGAATGCCATCTGGCCACGCATCATTGcc                                                                                        >  1:146617/2‑174 (MQ=255)                                                                                                                                                              ggAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCGTGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGAAAGCCCAGCGTCACGCCTTTGAAAAGCAAGGTTTGAATGCCATCTGGCCACGCATCATTGCCCTGGTGGTTCAACCCGGCGTCGAATTCGATCACACCAACGTTATTGATTATCAGCCCGCCAAAGCGAGCGCCTTAAGCCAGATGGt  <  1:503934/249‑1 (MQ=255)                                                                                                                                                                                      gCCGGTACNTNNNGNNGCGCGTGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCNGCCAGCGCCACGCTGNAANCCCANCGTNAAGCCTTTGAAAAGCAAGGTTTGAATg                                                                                                               >  2:359870/1‑116 (MQ=255)                                                                                                                                                                                      gCCGGTACCTGTCTGCGGGCGTGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTTGAAGCCCATCGTCACGCCTTTGAAAAGCAAGGTTTGAATg                                                                                                               <  1:359870/116‑1 (MQ=255)                                                                                                                                                                                       ccGGTCCCTGGCGGCGCGCGTGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGAAGCCCATCGTCACGCCTTTGAAAAGc                                                                                                                           <  1:7177/103‑1 (MQ=255)                                                                                                                                                                                       ccGGTACCTGGCGGCGCGCGTGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGAAGCCCATCGTCNNGCCTTTNna                                                                                                                               >  2:7177/1‑99 (MQ=255)                                                                                                                                                                                       ccGGTACCTGGCGGCGCGCGTGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGAAGCCCATCGTCACGCCTTTGAAAAGCAAGGTTTGAATGCCATCTGGCCACGCATCATTGCCCTGGTGGTTCAACCCGGCGTCGAATTCGATca                                                        >  1:388395/1‑170 (MQ=255)                                                                                                                                                                                             ccTGGCGGCGCGCGTGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGAAGCCCATCGTCACGCCTTTGAAAAGCAAGGTTTGAATGCCATCTGGCCACGCATCATTGCCCTGGTGGTTCAACCCGGCGTCGAATTCGATCACACCAACGTGATTGATTATCAGCCCGCCAAagcgagc                   <  1:389963/201‑1 (MQ=255)                                                                                                                                                                                                          |                                                                                                                                                                                                             ACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGAAGCCCATCGTCACGCCTTTGAAAAGCAAGGTTTGAATGCCATCTGGCCACGCATCATTGCCCTGGTGGTTCAACCCGGCGTCGAATTCGATCACACCAACGTTATTGATTATCAGCCCGCCAAAGCGAGCGCCTTAAGCCAGATGGT  >  NC_000913/3279260‑3279665

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