Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,417,189	A→G	A66A (GCT→GCC)	recE ←	Rac prophage; exonuclease VIII, 5' to 3' specific dsDNA exonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,417,189	0	A	G	100.0%	190.0 / NA	53	A66A (GCT→GCC)	recE	Rac prophage; exonuclease VIII, 5' to 3' specific dsDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/26);  total (27/26)

CCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCA  >  NC_000913/1417041‑1417330                                                                                                                                                     |                                                                                                                                              cctcAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGc                                                                                                                                             <  1:58797/151‑1 (MQ=255)       tctcAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGc                                                                                                                                         >  1:579622/1‑149 (MQ=255)            gTCATGTCTTCGCCGTTGACGTTGTTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCAt                                                                                                                                  >  2:246849/1‑151 (MQ=255)            gTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCAt                                                                                                                                  <  2:145293/151‑1 (MQ=255)            gTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCAt                                                                                                                                  >  2:32844/1‑151 (MQ=255)                 gTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACg                                                                                                                             <  2:420059/151‑1 (MQ=255)                         cGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGaa                                                                                                                      <  2:474121/150‑1 (MQ=255)                           ttGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGaaaa                                                                                                                    <  2:541893/150‑1 (MQ=255)                                     tATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAAGCGATCTCACCTTCGGCGGGCAGGTCATTAACGCCAGGAAAATTCGTGGcg                                                                                                           <  1:149021/149‑2 (MQ=255)                                      aTTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAAc                                                                                                         <  2:386294/150‑1 (MQ=255)                                      aTTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAAc                                                                                                         <  2:481039/150‑1 (MQ=255)                                         gCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGc                                                                                                      <  2:366336/150‑1 (MQ=255)                                         gCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCt                                                                                                     <  1:271568/151‑1 (MQ=255)                                          ccTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCtt                                                                                                    <  1:454267/151‑1 (MQ=255)                                            ttGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCttt                                                                                                   >  2:714/1‑150 (MQ=255)                                              gATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTaaa                                                                                                <  1:182597/151‑1 (MQ=255)                                                 aGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAAt                                                                                              >  1:132939/1‑150 (MQ=255)                                                       cGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAgctgct                                                                                       <  1:427681/151‑1 (MQ=255)                                                                  gctgctCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTAc                                                                            >  1:48471/1‑151 (MQ=255)                                                                        cccGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACt                                                                      <  1:539148/151‑1 (MQ=255)                                                                         ccGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACtt                                                                     <  2:753780/151‑1 (MQ=255)                                                                          cGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTa                                                                    <  2:307739/151‑1 (MQ=255)                                                                               ttAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGa                                                                <  1:240309/150‑1 (MQ=255)                                                                               ttAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGa                                                                <  2:445838/150‑1 (MQ=255)                                                                               ttAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGAt                                                               >  1:421209/1‑151 (MQ=255)                                                                                   ttCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAgg                                                           >  1:184421/1‑151 (MQ=255)                                                                                       cATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGt                                                       >  1:396439/1‑151 (MQ=255)                                                                                           tCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTcc                                                     >  1:344660/1‑149 (MQ=255)                                                                                           tCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTcc                                                     >  2:521519/1‑149 (MQ=255)                                                                                            cATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCaga                                                  <  1:400424/151‑1 (MQ=255)                                                                                            cATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGACCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTcc                                                     >  2:566648/1‑148 (MQ=255)                                                                                              tGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCagag                                                 >  2:668473/1‑150 (MQ=255)                                                                                                     tttGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGcac                                           >  1:215493/1‑149 (MQ=255)                                                                                                       tGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGt                                       >  2:703794/1‑151 (MQ=255)                                                                                                        gCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGcacacg                                        <  1:644264/149‑3 (MQ=255)                                                                                                          tGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCg                                     >  2:742033/1‑150 (MQ=255)                                                                                                              ttGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTc                                 >  1:491075/1‑150 (MQ=255)                                                                                                                gATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCaa                               >  1:72004/1‑150 (MQ=255)                                                                                                                 aTAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCaaa                              >  2:234679/1‑150 (MQ=255)                                                                                                                     cGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGt                         >  2:226718/1‑151 (MQ=255)                                                                                                                         cACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGttgct                     >  2:683048/1‑151 (MQ=255)                                                                                                                          aCTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGttgctt                    <  2:104835/151‑1 (MQ=255)                                                                                                                           atCCAGGTAAAATCGATCTCACCTTCGGAGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTAACTGACTTAACGAAAAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGttgct                     >  1:207130/2‑149 (MQ=255)                                                                                                                             ccAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGttgcttgc                  <  1:662898/150‑1 (MQ=255)                                                                                                                             ccAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGttgcttg                   <  1:297562/149‑1 (MQ=255)                                                                                                                               aGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGccc                >  1:569690/1‑150 (MQ=255)                                                                                                                               aGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGccc                >  2:576774/1‑150 (MQ=255)                                                                                                                               aGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCca               <  1:63479/151‑1 (MQ=255)                                                                                                                                   aaaaTCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAgg           >  1:665572/1‑151 (MQ=255)                                                                                                                                     aaTCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGgac         <  2:343505/151‑1 (MQ=255)                                                                                                                                         gATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGacgacg     >  1:19281/1‑151 (MQ=255)                                                                                                                                          aTCTCACCTTCGGCGGGCCGGTCATTAACGACAGGAAACTTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGt    >  1:396931/1‑151 (MQ=255)                                                                                                                                              ctccTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCa  <  2:461967/147‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              CCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCA  >  NC_000913/1417041‑1417330

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