Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	260,826	G→A	G147S (GGC→AGC)	proB →	glutamate 5‑kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	260,826	0	G	A	100.0%	78.6 / NA	25	G147S (GGC→AGC)	proB	glutamate 5‑kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/16);  total (9/16)

CCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTATGGCATTCACGTCGGGCAAATGCTGCTGACCCGTGCTGATATGGAAGACCGTGAACGCTTCCTGAACGCCCGCGACACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAAGGTTTGTATACCGCTGACCCGCGCAGCAATCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGT  >  NC_000913/260590‑261071                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      ccATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTATGGCATTCACGTCGGGCAAATGCTGCTGACCCGTGCTGATATGGAAGACCGTGAACGCTTCCTGAACGCCCGCGACACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCgataacgataac                                                                                                                                                                                                                                         >  2:147471/1‑251 (MQ=255)                       ggccgcggTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTATGGCATTCACGTCGGGCAAATGCTGCTGACCCGTGCTGATATGGAAGACCGTGAACGCTTCCTGAACGCCCGCGACACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGAt                                                                                                                                                                                                                     <  2:250154/245‑1 (MQ=255)                               tAGGGCAGAGTCGACTGATTCGACTGTGGGAACAGCTGTTTTCGATTTATGGCATTCACGTCGGGCAAATGCTGCTGACCCGTGCTGATATGGAAGACCGTGAACGCTTCCTGAACGCCCGCGACACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCGATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCgggg                                                                                                                                                                                                           >  1:225637/1‑250 (MQ=255)                                           gACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTATGGCATTCACGTCGGGCAAATGCTGCTGACCCGTGCTGATATGGAAGACCGTGAACGCTTCCTGAACGCCCGCGACACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTg                                                                                                                                                                                               <  2:236748/251‑1 (MQ=255)                                           gACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTATGGCATTCACGTCGGGCAAATGCTGCTGACCCGTGCTGATATGGAAGACCGTGAACGCTTCCTGAACGCCCGCGACACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTg                                                                                                                                                                                               >  2:104477/1‑251 (MQ=255)                                               gATTCAACTGTGGGAACAGCTGTTTTCGATTTATGGCATTCACGTCGGGCAAATGCTGCTGACCCGTGCTGATATGGAAGACCGTGAACGCTTCCTGAACGCCCGCGACACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTtgc                                                                                                                                                                                           <  2:242757/251‑1 (MQ=255)                                                                     ttttCGATTTATGGCATTCACGTCGGGCAAATGCTGCTGACCCGTGCTGATATGGAAGACCGTGAACGCTTCCTGAACGCCCGCGACACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAAGGttt                                                                                                                                                                     >  1:149190/1‑251 (MQ=255)                                                                           aTTTATGGCATTCACGTCGGGCAAATGCTGCTGACCCGTGCTGATATGGAAGACCGTGAACGCTTCCTGAACGCCCGCGACACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAAGGTTTGTATAc                                                                                                                                                               <  1:282551/251‑1 (MQ=255)                                                                               aTGGCATTCACGTCGGGCAAATGCTGCTGACCCGTGCTGATATGGAAGACCGTGAACGCTTCCTGAACGCCCGCGACACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAAGGTTTGTATACCGCt                                                                                                                                                           <  1:132471/251‑1 (MQ=255)                                                                                       cACGTCGGGCAAATGCTGCTGACCCGTGCTGATATGGAAGACCGTGAACGCTTCCTGAACGCCCGCGACACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAAGGTTTGTATACCGCTGACCcgcg                                                                                                                                                   <  1:105858/251‑1 (MQ=255)                                                                                                      ctgctgACCCGTGCTGATATGGAAGACCGTGAACGCTTCCTGAACGCCCGCGACACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAAGGTTTGTATACCGCTGACCCGCGCAGCAATCCgcaggc                                                                                                                                    >  1:149647/1‑251 (MQ=255)                                                                                                        gctgACCCGTGCTGATATGGAAGACCGTGAACGCTTCCTGAACGCCCGCGACACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAAGGTTTGTATACCGCTGACCCGCGCAGCAATCCgcaggc                                                                                                                                    <  1:115705/249‑1 (MQ=255)                                                                                                                  gCTGATATGGAAGACCGTGAACGCTTCCTGAACGCCCGCGACACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAAGGTTTGTATACCGCTGACCCGCGCAGCAATCCGCAGGCAGAACTGATTaa                                                                                                                        <  1:33539/251‑1 (MQ=255)                                                                                                                     gATATGGAAGACCGTGAACGCTTCCTGAACGCCCGCGACACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAAGGTTTGTATACCGCTGACCCGCGCAGCAATCCGCAGGCAGAACTGATTaa                                                                                                                        >  2:170802/1‑248 (MQ=255)                                                                                                                      atatGGAAGACCGTGAACGCTTCCTGAACGCCCGCGACACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAAGGTTTGTATACCGCTGACACGCGCAGCAATCCGCAGGCAGAACTGATTaa                                                                                                                        >  1:144340/1‑247 (MQ=255)                                                                                                                                         cTTCCTGAACGCCCGCGACACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAAGGTTTGTATACCGCTGACCCGCGCAGCAATCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACg                                                                                                 >  1:52190/1‑251 (MQ=255)                                                                                                                                               gAACGCCCGCGACACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAAGGTTTGTATACCGCTGACCCGCGCAGCAATCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACt                                                                                             <  1:104477/249‑1 (MQ=255)                                                                                                                                                                          gCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAAGGTTTGTATACCGCTGACCCGCGCAGCAATCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTc                                                                  >  2:266435/1‑249 (MQ=255)                                                                                                                                                                          gCTCGATAACAATATCGTTCCGGGAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCAATCAAAAAGGTTTGTATACCGCTGACCCGCGCAGCAATCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGAGGACGCACTGCGCGCGATTGCCTGTGACAGTGTTTCa                                                                 <  2:168429/250‑1 (MQ=255)                                                                                                                                                                                                     cAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAAGGTTTGTATACCGCTGACCCGCGCAGCAATCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCa                                     <  1:163815/251‑1 (MQ=255)                                                                                                                                                                                                             aCGATGCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAAGGTTTGTATACCGCTGACCCGCGCAGCAATCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAg                             <  1:103831/251‑1 (MQ=255)                                                                                                                                                                                                                  gCTGTCGCTACGGCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAAGGTTTGTATACCGCTGACCCGCGCAGCAATCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGTCCg                         <  2:42906/250‑1 (MQ=255)                                                                                                                                                                                                                             ggCAGAGATTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAAGGTTTGTATACCGCTGACCCGCGCAGCAATCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCtt             <  2:225637/251‑1 (MQ=255)                                                                                                                                                                                                                                    aTTAAGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAAGGTTTGTATACCGCTGACCCGCGCAGCAATCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGc      <  1:89632/251‑1 (MQ=255)                                                                                                                                                                                                                                        aGGTCAGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAATGTTTGTATACCGCTGACCCGCGCAGCAATCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACGGGCGGCATGAGTACCAAATTTCAGGCCGCTGACGTGGCTTGCCGTGCGGGt  <  1:206638/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      CCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTATGGCATTCACGTCGGGCAAATGCTGCTGACCCGTGCTGATATGGAAGACCGTGAACGCTTCCTGAACGCCCGCGACACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGATAAACTGTTGCTGCTGACCGATCAAAAAGGTTTGTATACCGCTGACCCGCGCAGCAATCCGCAGGCAGAACTGATTAAAGATGTTTACGGCATTGATGACGCACTGCGCGCGATTGCCGGTGACAGCGTTTCAGGCCTCGGAACTGGCGGCATGAGTACCAAATTGCAGGCCGCTGACGTGGCTTGCCGTGCGGGT  >  NC_000913/260590‑261071

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