Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,279,087	T→C	G58G (GGT→GGC)	kbaZ →	tagatose‑1,6‑bisphosphate aldolase 1 subunit KbaZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,279,087	0	T	C	100.0%	83.0 / NA	26	G58G (GGT→GGC)	kbaZ	tagatose‑1,6‑bisphosphate aldolase 1 subunit KbaZ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/10);  total (16/10)

GCGGGCAAAACAAATGGAATTTATGCCGTTTGTTCCGCACATCCGCTGGTGCTGGAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCC  >  NC_000913/3278950‑3279219                                                                                                                                          |                                                                                                                                        gCGGGCAAAACAAATGGAATTTATGCCGTTTGTTCCGCACATCCGCTGGTGCTGGAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGCta                                                                                                                                       <  2:768456/140‑1 (MQ=255)                      tATGCCGTTTGTTCCGCACATCCGCTGGTGCTGGAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACACCCGCCGa                                                                                                                 >  2:715135/1‑140 (MQ=255)                               tGTTCCGCACATCCGCTGGTGCTGGAAGCTGCAATCCGCTACGCCAGTGCAACCCAAACCCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCATTTCGGCGGCTATACC‑GGAAGGACGCCCGCCGATTTTCGCgg                                                                                                        <  2:487692/140‑1 (MQ=255)                               tGTTCCGCACATCCGCTGGTGCTGGAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCgg                                                                                                        >  2:455435/1‑140 (MQ=255)                                                   gCTGGAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCg                                                                                   >  2:295480/1‑140 (MQ=255)                                                   gCTGGAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCg                                                                                   >  2:326205/1‑140 (MQ=255)                                                      ggAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACt                                                                                >  1:713064/1‑140 (MQ=255)                                                           cTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTg                                                                          >  2:366319/1‑140 (MQ=255)                                                                  ccGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTATACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGCTTGTCAGCTCGCCGACT‑CGTTGAATTTcc                                                                   >  1:613970/1‑140 (MQ=255)                                                                        cGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAgg                                                             >  1:592066/1‑140 (MQ=255)                                                                             gTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCg                                                        <  1:243819/140‑1 (MQ=255)                                                                             gTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCg                                                        >  2:593853/1‑140 (MQ=255)                                                                              tGCAACCCAAACGCCTTTACTGATTGAAGCACCCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGt                                                       <  1:188101/140‑1 (MQ=255)                                                                              tGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGt                                                       >  1:350982/1‑140 (MQ=255)                                                                               gCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCg                                                        >  2:597023/1‑138 (MQ=255)                                                                               gCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCg                                                        <  1:597023/138‑1 (MQ=255)                                                                               gCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGtt                                                      >  2:168807/1‑140 (MQ=255)                                                                                   accaaacGCCGTTACTGATTGAAGCAACCT‑CAATCAGGGAGGACAGGTCGGGGGCTATACCGGGAAGAACCCCCGCCG‑TTTTTGCGGCCTTGGTTTGCCAG‑TCCCCGGCTCCGTTG‑ATTTCCCGGAGGGTGCGGTGGttt                                                 >  2:491706/1‑139 (MQ=255)                                                                                     caaacGCCGTTACTGATTGAAGCACCCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGTTGATTCt                                                <  2:337442/140‑1 (MQ=255)                                                                                        ccGCCGTTACTGATTGAAGCACCCTCCAATCGGGTATACCAGTTCGGCGGCTATCCC‑GGAATGCCGCCCGCCGATTTTCGCGG‑CTTTGCTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGTTGATTCTggg                                             <  2:613970/139‑1 (MQ=255)                                                                                        ccGCCGTTACTGATTGAAGCACCCCCCAATCAGGTAGACCAGTCCGGCGGCTATCCC‑GGAAGGCCGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGTTGATTCTggg                                             <  2:592066/139‑1 (MQ=255)                                                                                        acGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGTTGATTCTggg                                             <  1:168807/140‑1 (MQ=255)                                                                                              ttACTGATTGAAGCACCCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGa                                       <  2:350982/140‑1 (MQ=255)                                                                                               tACTGATTGAAGCACCCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGAc                                      <  2:511414/140‑1 (MQ=255)                                                                                                      tGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATc                                  <  1:709050/137‑1 (MQ=255)                                                                                                      tGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATc                                  >  2:709050/1‑137 (MQ=255)                                                                                                                        ggTAGACCAGTTTGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCaaa             >  1:571556/1‑140 (MQ=255)                                                                                                                        ggTAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCaaa             <  1:455435/140‑1 (MQ=255)                                                                                                                          tAGACCAGTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAAc           >  2:141750/1‑140 (MQ=255)                                                                                                                                 gtgcgggggTTATCCG‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTgccgg    <  2:399431/137‑1 (MQ=255)                                                                                                                                 gTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTgccgg    <  1:284314/140‑1 (MQ=255)                                                                                                                                 gTTCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTgccgg    <  1:491706/140‑1 (MQ=255)                                                                                                                                   tCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTgccggcc  >  1:284519/1‑140 (MQ=255)                                                                                                                                   tCGGCGGCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTgccggcc  >  2:379509/1‑140 (MQ=255)                                                                                                                                       cggCTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTgccgg    <  1:525896/134‑1 (MQ=255)                                                                                                                                          |                                                                                                                                        GCGGGCAAAACAAATGGAATTTATGCCGTTTGTTCCGCACATCCGCTGGTGCTGGAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACC‑GGAATGACGCCCGCCGATTTTCGCGG‑CTTTGTTTGTCAGCTCGCCGACT‑CGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCC  >  NC_000913/3278950‑3279219

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