Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,502,361	A→G	C709R (TGC→CGC)	fryA ←	putative PTS multiphosphoryl transfer protein FryA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,502,361	0	A	G	100.0%	116.3 / NA	33	C709R (TGC→CGC)	fryA	putative PTS multiphosphoryl transfer protein FryA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/22);  total (11/22)

ACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCA  >  NC_000913/2502224‑2502496                                                                                                                                          |                                                                                                                                        aCGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGcgca                                                                                                                                       >  1:512145/1‑140 (MQ=255)    gATCCACTGAGATTTTGTGTGCTGGCTCGTTACGCCACAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGa                                                                                                                                    <  1:175317/140‑1 (MQ=255)    gATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGa                                                                                                                                    <  1:584440/140‑1 (MQ=255)       ccACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAGACGGATGTGCAGAGCGACCGTTAACGCCGAAGTTGCCGCGCAGGAACt                                                                                                                                 >  2:175317/1‑140 (MQ=255)         aCTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTgg                                                                                                                               <  2:114344/140‑1 (MQ=255)            gagaTTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATc                                                                                                                            <  2:470438/140‑1 (MQ=255)                    gtgtGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGctc                                                                                                                    >  2:513973/1‑140 (MQ=255)                         cGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCtt                                                                                                                  <  2:664962/137‑1 (MQ=255)                         cGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCtt                                                                                                                  >  1:664962/1‑137 (MQ=255)                         cGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTAt                                                                                                               >  2:435270/1‑140 (MQ=255)                              tCGCTACGCCAAAACCCACGCCAGTGGTAACAACCGCTTCCCGCTGCCTGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTa                                                                                                          <  2:312336/140‑1 (MQ=255)                                   aCGCCAAAACCAACGCCTGTGGCAACAATCTCTTCCCGCCGCCAGACACCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATc                                                                                                     <  2:171202/140‑1 (MQ=255)                                      ccAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTg                                                                                                  <  2:606336/140‑1 (MQ=255)                                      ccAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTg                                                                                                  <  1:470686/140‑1 (MQ=255)                                        aaaaCCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGAt                                                                                                <  1:374461/140‑1 (MQ=255)                                          aaCCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATca                                                                                              <  1:3444/140‑1 (MQ=255)                                              aaCGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATcaacaa                                                                                          <  1:123773/140‑1 (MQ=255)                                              aaCGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATcaacaa                                                                                          <  1:577121/140‑1 (MQ=255)                                                      tggtAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAAGTGGATCGCGTGCACTTGATTGCTAATATCCTGATCAACAAAGATATtc                                                                                  >  2:271411/1‑140 (MQ=255)                                                         tAACAATCTCTTCCCGCTGCCACACATCTCCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTcc                                                                               <  1:247219/140‑1 (MQ=255)                                                               tctcTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCATGCTCTTTGTTGATAAGATCCTGATTAACAAAGATATTCTCCAGCGcc                                                                         >  2:205432/1‑140 (MQ=255)                                                               tctcTTCCCGCTGCCAGACATCCTCCTCCGGCTCAAAAGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGcc                                                                         <  1:205432/140‑1 (MQ=255)                                                                            ccAGACATCTTCTTCCATCTCACACGGATGTCCAGTGCGGCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGaa                                                            <  2:270667/140‑1 (MQ=255)                                                                                            aGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTC‑ACAGCGCCCGCAGTGCGCGAACTTCTTCTTCCGGCGTa                                           >  1:226756/1‑140 (MQ=255)                                                                                             gCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTa                                           <  1:513973/140‑1 (MQ=255)                                                                                                           tCAGTGCGCCCGTTATCGCCGACGTTGCCCCGCAGGACCTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTaa                             <  1:195481/140‑1 (MQ=255)                                                                                                               tGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCt                         >  2:13006/1‑140 (MQ=255)                                                                                                               tGCGCCCGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCt                         >  2:345593/1‑140 (MQ=255)                                                                                                                    ccGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAAt                    <  2:147707/140‑1 (MQ=255)                                                                                                                    ccGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAAt                    >  1:47443/1‑140 (MQ=255)                                                                                                                     cGTTAACGCCGAGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAAtt                   <  1:435270/140‑1 (MQ=255)                                                                                                                                aGGTTGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACt        >  2:435358/1‑140 (MQ=255)                                                                                                                                    tGCCGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCgg    <  2:512145/140‑1 (MQ=255)                                                                                                                                      ccGCGCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCa  <  1:63129/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                        ACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCA  >  NC_000913/2502224‑2502496

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