Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,130,464	T→C	T139A (ACC→GCC)	wcaD ←	putative colanic acid polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,130,464	0	T	C	100.0%	206.5 / NA	57	T139A (ACC→GCC)	wcaD	putative colanic acid polymerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/23);  total (34/23)

AGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATAT  >  NC_000913/2130323‑2130610                                                                                                                                              |                                                                                                                                                   agaTTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAAt                                                                                                                                           >  1:608625/1‑151 (MQ=255)        tGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGcc                                                                                                                                    >  1:870878/1‑151 (MQ=255)          ccAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCAt                                                                                                                                  >  1:609246/1‑151 (MQ=255)           cAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCAt                                                                                                                                  >  2:990157/1‑150 (MQ=255)                 gAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCAcc                                                                                                                           >  2:251813/1‑151 (MQ=255)                         ccGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAg                                                                                                                   <  2:8360/151‑1 (MQ=255)                           ggTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAgcg                                                                                                                 >  1:349784/1‑151 (MQ=255)                            gTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAgcgc                                                                                                                >  1:416709/1‑151 (MQ=255)                            gTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAgcgc                                                                                                                <  2:348710/151‑1 (MQ=255)                                gAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGccacc                                                                                                            >  1:688877/1‑151 (MQ=255)                                 aaaTAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGccacc                                                                                                            >  2:870572/1‑150 (MQ=255)                                   aTAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGccaccacc                                                                                                         <  1:859461/151‑1 (MQ=255)                                    tAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGccaccacc                                                                                                         >  2:221147/1‑150 (MQ=255)                                     aaaGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGccaccacca                                                                                                        <  1:827679/150‑1 (MQ=255)                                             gTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGacc                                                                                                <  1:94772/150‑1 (MQ=255)                                             gTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGacc                                                                                                >  2:871073/1‑150 (MQ=255)                                              tCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGaccac                                                                                              >  2:265169/1‑151 (MQ=255)                                               cGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGaccacc                                                                                             >  1:661780/1‑151 (MQ=255)                                               cGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGaccac                                                                                              <  2:748908/150‑1 (MQ=255)                                                gTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGaccacc                                                                                             >  1:833062/1‑150 (MQ=255)                                                   cgcTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGAACCATCTCCACCGCCTCCAGCGCCACCTCCAGCCCGACCACCAGAt                                                                                         <  2:797928/151‑1 (MQ=255)                                                    gcTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGAt                                                                                         >  1:727815/1‑150 (MQ=255)                                                    gcTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGAt                                                                                         <  2:840613/150‑1 (MQ=255)                                                    gcTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGAt                                                                                         <  2:749816/150‑1 (MQ=255)                                                       tGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATaa                                                                                       <  2:336072/149‑1 (MQ=255)                                                            ccgAATTTAATGAAATTCAG‑‑CATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGt                                                                              <  1:974754/151‑1 (MQ=255)                                                             cgAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACg                                                                               >  2:627348/1‑151 (MQ=255)                                                               aaTTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTa                                                                             <  1:268518/151‑1 (MQ=255)                                                                aTTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTaa                                                                            <  2:803621/151‑1 (MQ=255)                                                                        aaaTTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGc                                                                    >  1:126705/1‑151 (MQ=255)                                                                        aaaTTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGc                                                                    <  2:688801/151‑1 (MQ=255)                                                                          aTTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCg                                                                   >  2:192361/1‑150 (MQ=255)                                                                          aTTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGa                                                                  <  2:334618/151‑1 (MQ=255)                                                                           ttCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGAt                                                                 >  1:780107/1‑151 (MQ=255)                                                                           ttCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGAt                                                                 >  2:816773/1‑151 (MQ=255)                                                                            tCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATg                                                                >  1:560748/1‑151 (MQ=255)                                                                                cacaTAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCgatggatg                                                            >  2:140497/1‑151 (MQ=255)                                                                                            gCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAAtt                                                >  1:4411/1‑151 (MQ=255)                                                                                               gTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGtt                                             >  2:768097/1‑151 (MQ=255)                                                                                                tAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTg                                               <  1:574610/148‑1 (MQ=255)                                                                                                    tCAGATATTTACAAATCGACTCCAAAATACTGCTGCCACCGGCGAGGAAAATTTGGGCAATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATACCCTGCCATGGCTGAGGGGCGAAATTGTTTTa                                          >  2:831683/1‑149 (MQ=255)                                                                                                     cAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTaaaa                                       <  1:251841/151‑1 (MQ=255)                                                                                                     cAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTaaa                                        >  2:210128/1‑150 (MQ=255)                                                                                                        atatTTCGAAATCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGa                                    <  2:168318/151‑1 (MQ=255)                                                                                                                aaaTCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCca                            <  1:342649/151‑1 (MQ=255)                                                                                                                 aaTCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCcac                           >  1:478049/1‑151 (MQ=255)                                                                                                                 aaTCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCcac                           <  2:793723/151‑1 (MQ=255)                                                                                                                 aaTCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCcac                           <  2:932422/151‑1 (MQ=255)                                                                                                                  aTCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCcac                           >  1:215834/1‑150 (MQ=255)                                                                                                                  aTCGACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCca                            >  2:178091/1‑149 (MQ=255)                                                                                                                     gACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACAcc                        >  1:725492/1‑150 (MQ=255)                                                                                                                     gACTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACAcc                        >  1:666771/1‑150 (MQ=255)                                                                                                                      aCTCCATAATACTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCAt                      >  1:279116/1‑151 (MQ=255)                                                                                                                             aataCTGCTGCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCaat                <  1:286970/150‑1 (MQ=255)                                                                                                                                 ctgctgCCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATg            >  1:49222/1‑150 (MQ=255)                                                                                                                                   gctgcCACCGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGa           <  2:233068/149‑1 (MQ=255)                                                                                                                                           cGGCGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAatat  >  2:109186/1‑150 (MQ=255)                                                                                                                                              |                                                                                                                                                   AGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATAT  >  NC_000913/2130323‑2130610

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