Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,136,954	A→G	V97A (GTC→GCC)	wza ←	colanic acid export protein; outer membrane auxillary lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,136,954	0	A	G	100.0%	114.5 / NA	33	V97A (GTC→GCC)	wza	colanic acid export protein; outer membrane auxillary lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/15);  total (18/15)

CGAGCGTTTTCCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGACGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACTTTCCAGCGTCATATTGGGGCGCGCTACATTCGGGCGTGGGCGTAATTGGTCAATCAGGCGCGGGGTCAGCGGATAAACAT  >  NC_000913/2136820‑2137092                                                                                                                                       |                                                                                                                                           cgagcgTTTTCCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAAc                                                                                                                            >  1:659637/1‑151 (MQ=255)      gTTTTCCCGACTACGTGGACCTTGCCGTTATACGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTccc                                                                                                                       >  2:585187/1‑151 (MQ=255)        tttCCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTccccc                                                                                                                     >  2:237434/1‑151 (MQ=255)          tCCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCgg                                                                                                                   >  2:506940/1‑151 (MQ=255)                   cGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGAcgcg                                                                                                          <  2:82937/151‑1 (MQ=255)                    gTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGGGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGAcg                                                                                                            >  1:345762/1‑148 (MQ=255)                      ggACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATa                                                                                                       <  1:190511/151‑1 (MQ=255)                            tGCCGATATACGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGt                                                                                                  <  2:474385/150‑1 (MQ=255)                               cGATATACGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAAtt                                                                                              <  2:447986/151‑1 (MQ=255)                               cGATATACGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAAtt                                                                                              >  2:252125/1‑151 (MQ=255)                                gATATACGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTc                                                                                             <  2:87000/151‑1 (MQ=255)                                     aCGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGa                                                                                         <  2:312110/150‑1 (MQ=255)                                     aCGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGAt                                                                                        >  2:600620/1‑151 (MQ=255)                                      cGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGAtc                                                                                       >  2:148884/1‑151 (MQ=255)                                           aaaaCATAGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACt                                                                                  >  1:549032/1‑151 (MQ=255)                                               cATAGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACTTTcc                                                                              <  1:655730/151‑1 (MQ=255)                                                  aGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACTTTCCAGc                                                                           >  2:189347/1‑151 (MQ=255)                                                  aGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGGCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACTTTCCAGc                                                                           >  1:248690/1‑151 (MQ=255)                                                      ccGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACTTTCCAGCGTCa                                                                       <  1:626883/151‑1 (MQ=255)                                                       cGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACTTTCCAGCGTCat                                                                      >  2:649200/1‑151 (MQ=255)                                                           aGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACTTTCCAGCGTCatat                                                                    >  1:320277/1‑149 (MQ=255)                                                                 tACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACTTTCCAGCGTCATATTGGGgcgc                                                            >  1:263903/1‑151 (MQ=255)                                                                            cGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACTTTCCAGCGTCATATTGGGGCGCGCTACATTCgg                                                 >  2:406819/1‑151 (MQ=255)                                                                             ggTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACTTTCCAGCGTCATATTGGGGCGCGCTACATTCggg                                                <  1:628166/151‑1 (MQ=255)                                                                                     tGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACTTTCCAGCGTCATATTGGGGCGCGCTACATTCgggcgtgggcg                                        >  2:577080/1‑151 (MQ=255)                                                                                      ggAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACTTTCCAGCGTCATATTGGGGCGCGCTACATTCgggcgtgggcgt                                       <  2:248691/151‑1 (MQ=255)                                                                                         gCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGAGCTCACTTTCCAGCGTCATATTGGGGCGCGCTACATTCGGGCGTGGGCGTAAt                                    <  2:359536/151‑1 (MQ=255)                                                                                                     aCCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACTTTCCAGCGTCATATTGGGGCGCGCTACATTCGGGCGTGGGCGTAATTGGTCAATCAgg                        <  1:111575/151‑1 (MQ=255)                                                                                                        ggCTGGCGTGGTGAGTTCCGGGGGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACTTTCCAGCGTCATATTGGGGCGCCCTACATTCGGGCGTGGGCGTAATTGGGCAATCAGgcg                      >  1:143166/1‑150 (MQ=255)                                                                                                           tGGCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACTTTCCAGCGTCATATTGGGGCGCGCTACATTCGGGCGTGGGCGTAATTGGTCAATCAGGCGCggg                  >  1:310556/1‑151 (MQ=255)                                                                                                            ggCGTGGTGAGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACTTTCCAGCGTCATATTGGGGCGCGCTACATTCGGGCGTGGGCGTAATTGGTCAATCAGGCGCgg                   <  1:520315/149‑1 (MQ=255)                                                                                                                     aGTTCCGGGTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACTTTCCAGCATCATATTGGGGCGCGCTACATTCGGGCGTGGGCGTAATTGGTCAATCAGGCGCGGGGTCAGCGGAt        <  1:531524/151‑1 (MQ=255)                                                                                                                           gggTGATCCCAGGCGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACTTTCCAGCGTCATATTGGGGCGCGCTACATTCGGGCGTGGGCGTAATTGGTCAATCAGGCGCGGGGTCAGCGGATAAACAt  <  2:421811/151‑1 (MQ=255)                                                                                                                                       |                                                                                                                                           CGAGCGTTTTCCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGGTGTCGCTGGAGCTGCGGTACTGACCGGCTGGCGTGGTGAGTTCCGGGTGATCCCAGACGGTGACATTAAGAACGTCCCCCGGCCCGACGCGATACTGGTAATTCGCGATCTCACTTTCCAGCGTCATATTGGGGCGCGCTACATTCGGGCGTGGGCGTAATTGGTCAATCAGGCGCGGGGTCAGCGGATAAACAT  >  NC_000913/2136820‑2137092

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