Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	231,586	A→G	I155M (ATA→ATG)	yafD →	endo/exonuclease/phosphatase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	231,586	0	A	G	100.0%	237.9 / NA	63	I155M (ATA→ATG)	yafD	endo/exonuclease/phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/36);  total (27/36)

CCCTTTCGGCGGCACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAG  >  NC_000913/231456‑231714                                                                                                                                   |                                                                                                                                 cccTTTCGGCGGCACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGt                                                                                                                          <  1:1143317/139‑1 (MQ=255) cccTTTCGGCGGCACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGt                                                                                                                          <  1:445207/139‑1 (MQ=255)     tgcggcggCACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAAc                                                                                                                      <  1:86053/137‑1 (MQ=255)                    aGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTg                                                                                                       <  2:387379/139‑1 (MQ=255)                    aGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTg                                                                                                       <  2:1369931/139‑1 (MQ=255)                    aGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTg                                                                                                       <  1:717043/139‑1 (MQ=255)                     gtgtATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTgg                                                                                                      <  1:1283680/139‑1 (MQ=255)                      tgtATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGa                                                                                                     <  1:866141/139‑1 (MQ=255)                         aTTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGt                                                                                                  <  1:343256/139‑1 (MQ=255)                          ttGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTg                                                                                                       <  2:884311/133‑1 (MQ=255)                          ttGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTg                                                                                                       >  1:884311/1‑133 (MQ=255)                           tgctgcCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTg                                                                                                       >  2:999486/1‑132 (MQ=255)                           tgctgcCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTg                                                                                                       <  1:999486/132‑1 (MQ=255)                                cccGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGt                                                                                           >  2:433459/1‑139 (MQ=255)                                cccGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGt                                                                                           >  1:599842/1‑139 (MQ=255)                                  cGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTg                                                                                                       <  2:795743/125‑1 (MQ=255)                                  cGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTg                                                                                                       >  1:795743/1‑125 (MQ=255)                                   gTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAg                                                                                        <  1:1114191/139‑1 (MQ=255)                                               agaACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATgca                                                                                                                                <  2:1045343/87‑1 (MQ=255)                                                 aaCCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTg                                                                                                       >  2:1359612/1‑110 (MQ=255)                                                 aaCCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTg                                                                                                       <  1:1359612/110‑1 (MQ=255)                                                   cccATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCg                                                                                                         <  2:1175588/106‑1 (MQ=255)                                                   cccATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCg                                                                                                         >  1:1175588/1‑106 (MQ=255)                                                     cATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGc                                                                      <  2:498256/139‑1 (MQ=255)                                                      aTTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGt                                                                                                  >  2:1455350/1‑110 (MQ=255)                                                      aTTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGt                                                                                                  <  1:1455350/110‑1 (MQ=255)                                                       ttttGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGa                                                                    <  1:148761/139‑1 (MQ=255)                                                       ttttGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGa                                                                    <  1:1077479/139‑1 (MQ=255)                                                        tttGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGAt                                                                   >  1:1358944/1‑139 (MQ=255)                                                           gCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATTCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGt                                                                                                                <  2:996720/91‑1 (MQ=255)                                                           gCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGt                                                                                                                >  1:996720/1‑91 (MQ=255)                                                            cttCTGTCTAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTCGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGa                                                               <  2:844103/137‑1 (MQ=255)                                                            cGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGa                                                               <  2:1427981/139‑1 (MQ=255)                                                            cGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGa                                                               <  2:1358944/139‑1 (MQ=255)                                                                 ggCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGt                                                                                                                >  2:527161/1‑85 (MQ=255)                                                                 ggCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGt                                                                                                                <  1:527161/85‑1 (MQ=255)                                                                 ggCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCt                                                          >  1:239732/1‑139 (MQ=255)                                                                            cACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCg                                                                                                         >  2:1479229/1‑81 (MQ=255)                                                                            cACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCg                                                                                                         <  1:1479229/81‑1 (MQ=255)                                                                            cACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCgg                                               >  1:1278414/1‑139 (MQ=255)                                                                            cACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCgg                                               <  2:1272112/139‑1 (MQ=255)                                                                                           aTCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGc                                >  2:1200046/1‑139 (MQ=255)                                                                                              cATTACCTGACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCggg                             <  2:599842/139‑1 (MQ=255)                                                                                                      gACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTg                                                                                                       <  2:799490/57‑1 (MQ=255)                                                                                                      gACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTg                                                                                                       >  1:799490/1‑57 (MQ=255)                                                                                                      gACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAg                                                                                        <  2:151043/72‑1 (MQ=255)                                                                                                      gACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAg                                                                                        >  1:151043/1‑72 (MQ=255)                                                                                                      gACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCa                     >  2:1321056/1‑139 (MQ=255)                                                                                                      gACACCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCa                     >  1:173614/1‑139 (MQ=255)                                                                                                        cacCCGCCTGTTGATGGTGGTTAATATGCACGCCGTc                                                                                                                         >  2:971543/1‑37 (MQ=255)                                                                                                        cacCCGCCTGTTGATGGTGGTTAATATGCACGCCGTc                                                                                                                         <  1:971543/37‑1 (MQ=255)                                                                                                        cacCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTAc                                                                                 <  1:431293/77‑1 (MQ=255)                                                                                                        cacCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTAc                                                                                 >  2:431293/1‑77 (MQ=255)                                                                                                        cacCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGga                            >  2:229365/1‑130 (MQ=255)                                                                                                        cacCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGga                            <  1:229365/130‑1 (MQ=255)                                                                                                        cacCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAAt                   >  1:638517/1‑139 (MQ=255)                                                                                                        cacCCGCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAAt                   >  1:1326531/1‑139 (MQ=255)                                                                                                             gCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTg              <  1:1116502/139‑1 (MQ=255)                                                                                                             gCCTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTg              <  1:1162991/139‑1 (MQ=255)                                                                                                              ccTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCa                                                                 >  1:860507/1‑87 (MQ=255)                                                                                                              ccTGTTGATGGTGGTTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCa                                                                 <  2:860507/87‑1 (MQ=255)                                                                                                                         tggtTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAg  >  1:238487/1‑139 (MQ=255)                                                                                                                           gtTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAg           >  2:1017910/1‑128 (MQ=255)                                                                                                                           gtTAATATGCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAg           <  1:1017910/128‑1 (MQ=255)                                                                                                                                   |                                                                                                                                 CCCTTTCGGCGGCACATCCAGTGTATTGCTGCCCGTTACGTGAACGAGAACCCATTTTGCGTCTGGCGAAGTCGGCACTGGTGACGGTCTATCCATTACCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAG  >  NC_000913/231456‑231714

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