Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,316,959	C→A	S87I (AGC→ATC)	phnM ←	RPnTP hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,316,959	0	C	A	100.0%	199.0 / NA	55	S87I (AGC→ATC)	phnM	RPnTP hydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (26/29);  total (26/29)

CGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCC  >  NC_000913/4316824‑4317095                                                                                                                                        |                                                                                                                                         cGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGc                                                                                                                                      <  2:908615/140‑1 (MQ=255) cGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGc                                                                                                                                      <  2:1146373/140‑1 (MQ=255)  ggTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGcc                                                                                                                                     <  1:1068374/140‑1 (MQ=255)  ggTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGcc                                                                                                                                     <  1:693000/140‑1 (MQ=255)    tGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGccac                                                                                                                                   >  2:650214/1‑140 (MQ=255)    tGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGccac                                                                                                                                   >  2:1430322/1‑140 (MQ=255)    tGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGccac                                                                                                                                   >  1:932322/1‑140 (MQ=255)          gcgcGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGc                                                                                                                                      >  2:686646/1‑131 (MQ=255)          gcgcGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGc                                                                                                                                      <  1:686646/131‑1 (MQ=255)          gcgcGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAg                                                                                                                             <  1:1094840/140‑1 (MQ=255)              ggTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGcgcg                                                                                                                         >  2:418184/1‑140 (MQ=255)               gTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGcgt                                                                                                                        <  1:418184/140‑1 (MQ=255)                      ccgcGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGgct                                                                                                                 <  2:1284183/140‑1 (MQ=255)                       cgcgTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGgctg                                                                                                                >  2:1119544/1‑140 (MQ=255)                          gTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGgctgct                                                                                                              <  2:1031255/139‑1 (MQ=255)                          gTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGgctgct                                                                                                              >  1:1031255/1‑139 (MQ=255)                                  gTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGa                                                                                                     <  1:1430322/140‑1 (MQ=255)                                  gTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGa                                                                                                     <  1:272427/140‑1 (MQ=255)                                  gTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGa                                                                                                     <  1:47918/140‑1 (MQ=255)                                       ttCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTggg                                                                                                <  2:932322/140‑1 (MQ=255)                                         cGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCa                                                                                              <  1:31111/140‑1 (MQ=255)                                           aTGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAgg                                                                                            <  1:1462194/140‑1 (MQ=255)                                            tGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGc                                                                                           <  1:59787/140‑1 (MQ=255)                                            tGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGc                                                                                           <  2:1587793/140‑1 (MQ=255)                                                        tcttctCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTc                                                                               >  2:466292/1‑140 (MQ=255)                                                         cttctCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCg                                                                              >  1:95963/1‑140 (MQ=255)                                                            ctcCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGgc                                                                           >  2:543594/1‑140 (MQ=255)                                                              ccAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGgcgc                                                                         >  2:817426/1‑140 (MQ=255)                                                               cAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGgcg                                                                          <  1:1482990/138‑1 (MQ=255)                                                               cAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGgcg                                                                          >  2:1482990/1‑138 (MQ=255)                                                                aGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGGACGGTGGTGATGCCGATCGACACCATCAGCGCGTCGTAGCTCCTCATCCCCGAGAGGGCAGGCCAGACCACTTTCGCGCGCgg                                                                       >  2:84901/1‑140 (MQ=255)                                                                   ttCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGt                                                                    >  1:939099/1‑140 (MQ=255)                                                                     ctcCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTga                                                                  <  2:667041/140‑1 (MQ=255)                                                                     ctcCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTga                                                                  >  2:1341364/1‑140 (MQ=255)                                                                          gCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAAt                                                             <  2:331678/140‑1 (MQ=255)                                                                            cGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAAttt                                                           >  2:831805/1‑140 (MQ=255)                                                                               tCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATc                                                        <  1:394436/140‑1 (MQ=255)                                                                                         cgcgCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCg                                              <  2:939099/140‑1 (MQ=255)                                                                                             cACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTAt                                          <  1:466292/140‑1 (MQ=255)                                                                                              aCGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATg                                         <  1:650214/140‑1 (MQ=255)                                                                                              aCGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATg                                         <  2:963835/140‑1 (MQ=255)                                                                                                  cGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGAt                                                   <  2:505697/126‑1 (MQ=255)                                                                                                  cGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGAt                                                   >  1:505697/1‑126 (MQ=255)                                                                                                  cGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGc                                     >  2:523414/1‑140 (MQ=255)                                                                                                        ttGCCACGGCATCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATc                               <  2:920519/140‑1 (MQ=255)                                                                                                             aCGGCTTCCAGTACGGTGGTGATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGccc                          >  1:655158/1‑140 (MQ=255)                                                                                                                             tggtgATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCTGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAgccg          >  1:748181/1‑140 (MQ=255)                                                                                                                             tggtgATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAgccg          >  2:737192/1‑140 (MQ=255)                                                                                                                             tggtgATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAgccg          >  2:813603/1‑140 (MQ=255)                                                                                                                             tggtgATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAgccg          >  2:1476639/1‑140 (MQ=255)                                                                                                                             tggtgATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAgccg          >  1:666443/1‑140 (MQ=255)                                                                                                                             tggtgATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAgccg          >  1:273414/1‑140 (MQ=255)                                                                                                                             tggtgATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAgccg          >  2:870669/1‑140 (MQ=255)                                                                                                                             tggtgATGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAgccg          >  1:50140/1‑140 (MQ=255)                                                                                                                                  aTGCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTc     >  1:66529/1‑140 (MQ=255)                                                                                                                                    gCCGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGc   <  1:1119544/140‑1 (MQ=255)                                                                                                                                     ccGATCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGcc  <  1:523414/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                         CGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCC  >  NC_000913/4316824‑4317095

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