Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,863,111	C→A	C259* (TGC→TGA)	ygbK →	putative 3‑oxo‑tetronate kinase YgbK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,863,111	0	C	A	100.0%	113.9 / NA	33	C259* (TGC→TGA)	ygbK	putative 3‑oxo‑tetronate kinase YgbK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (15/18);  total (15/18)

CGCGATGCCCCACTGGTAACGGGCGGTTCTGGTCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAA  >  NC_000913/2862974‑2863249                                                                                                                                          |                                                                                                                                           cgcgATGCCCCACTGGTAACGGGCGGTTCTGGTCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGatc                                                                                                                                          >  2:442120/1‑140 (MQ=255)   cgATGCCCCACTGGTAACGGGCGGTTCTGGTCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGAtctc                                                                                                                                        >  2:64018/1‑140 (MQ=255)                gTAACGGGCGGTTCTGGTCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCg                                                                                                                           <  1:847241/140‑1 (MQ=255)                gTAACGGGCGGTTCTGGTCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCg                                                                                                                           <  1:491693/140‑1 (MQ=255)                   aCGGGCGGTTCTGGTCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCa                                                                                                                        >  1:588926/1‑140 (MQ=255)                   aCGGGCGGTTCTGGTCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCATCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGATTCATGATCTCAAATGACCAACCGCCa                                                                                                                        >  1:744245/1‑140 (MQ=255)                     gggCGGTTCTGGTCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAgg                                                                                                                      >  2:50211/1‑140 (MQ=255)                     gggCGGTTCTGGTCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAgg                                                                                                                      <  1:164517/140‑1 (MQ=255)                                     gCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCa                                                                                                      >  1:519745/1‑140 (MQ=255)                                     gCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCa                                                                                                      >  1:612938/1‑140 (MQ=255)                                       gATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCaaa                                                                                                    <  2:393702/140‑1 (MQ=255)                                          tggcctggccCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTg                                                                                                 >  1:462682/1‑140 (MQ=255)                                             cctggccCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCAc                                                                                              <  1:24369/140‑1 (MQ=255)                                                 gccCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGc                                                                                          <  1:442120/140‑1 (MQ=255)                                                    cGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCg                                                                                       <  2:621929/140‑1 (MQ=255)                                                    cGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCg                                                                                       <  1:599174/140‑1 (MQ=255)                                                                    aaaaCGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCa                                                                       <  1:682362/140‑1 (MQ=255)                                                                       aCGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGc                                                                    <  1:448346/140‑1 (MQ=255)                                                                         ggTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCATGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTg                                                                  <  2:99327/140‑1 (MQ=255)                                                                             aCCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCctc                                                              >  2:128835/1‑140 (MQ=255)                                                                             aCCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCctc                                                              <  2:588926/140‑1 (MQ=255)                                                                                    tCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTg                                                       >  1:902431/1‑140 (MQ=255)                                                                                         aaGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAAct                                                  <  1:64018/140‑1 (MQ=255)                                                                                              gggCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTATCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGc                                             <  2:834084/140‑1 (MQ=255)                                                                                                gCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCg                                           >  1:149850/1‑140 (MQ=255)                                                                                                    ccGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTa                                       >  1:47096/1‑140 (MQ=255)                                                                                                        tcgctGGGCGCGGCGTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGca                                   <  2:612938/140‑1 (MQ=255)                                                                                                                  cggcgTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTgg                         <  1:515082/140‑1 (MQ=255)                                                                                                                    gcgTAGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCC‑GGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCag                      >  1:329670/1‑140 (MQ=255)                                                                                                                        aGTGCTCTCCGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGt                   >  1:408916/1‑140 (MQ=255)                                                                                                                              ctcCGATTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTc             <  2:744245/140‑1 (MQ=255)                                                                                                                                 cGGTTCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTgg          >  1:507245/1‑140 (MQ=255)                                                                                                                                     tCATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCa      <  2:462682/140‑1 (MQ=255)                                                                                                                                      cATGATCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAg     >  2:617477/1‑140 (MQ=255)                                                                                                                                        tactctAAAATGCCCAACCGCC‑GGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGaa  <  2:329670/138‑1 (MQ=255)                                                                                                                                          |                                                                                                                                           CGCGATGCCCCACTGGTAACGGGCGGTTCTGGTCTGGCGATTGGCCTGGCCCGGCAGTGGGCGCAAGAAAACGGTAACCAGGCTCGCAAAGCAGGGCGTCCGCTCGCTGGGCGCGGCGTAGTGCTCTCCGGTTCATGCTCTCAAATGACCAACCGCCAGGTAGCACATTACCGTCAAATTGCACCAGCCCGTGAAGTTGATGTGGCACGCTGCCTCTCAATTGAAACTCTGGCCGCTTATGCACACGAACTGGCAGAGTGGGTTCTGGGCCAGGAA  >  NC_000913/2862974‑2863249

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