Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	91,323	G→A	E98K (GAA→AAA)	ftsL →	cell division protein FtsL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	91,323	0	G	A	100.0%	173.8 / NA	48	E98K (GAA→AAA)	ftsL	cell division protein FtsL
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/25);  total (23/25)

GAGCGAGATGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGGAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATCGTAGTGCAAAAATAAGGA  >  NC_000913/91248‑91400                                                                            |                                                                              gAGCGAGATGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAgg                                                                          <  1:22013/81‑1 (MQ=255) gAGCGAGATGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAg                                                                           <  1:22012/80‑1 (MQ=255) gAGCGAGATGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAg                                                                           >  1:22071/1‑80 (MQ=255)  aGCGAGATGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAACGCGCTCGGCGACCATAGCCGGGTGAAAAGGa                                                                         >  2:22077/1‑81 (MQ=255)       gaTGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGcc                                                                    <  2:22014/81‑1 (MQ=255)        aTGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCa                                                                   <  2:22015/81‑1 (MQ=255)         tGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCAc                                                                  <  1:22020/81‑1 (MQ=255)         tGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCAc                                                                  <  2:22016/81‑1 (MQ=255)         tGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCAc                                                                  <  2:22017/81‑1 (MQ=255)         tGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCAc                                                                  <  1:22018/81‑1 (MQ=255)         tGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCAc                                                                  <  1:22019/81‑1 (MQ=255)          gCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACg                                                                 >  1:22062/1‑81 (MQ=255)            tttAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACgg                                                                >  1:22054/1‑80 (MQ=255)              tAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGaa                                                              >  2:22075/1‑80 (MQ=255)              tAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGa                                                               <  1:22021/79‑1 (MQ=255)               aGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGaaa                                                             <  2:22022/80‑1 (MQ=255)               aGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGaa                                                              >  2:22064/1‑79 (MQ=255)                   aTTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTg                                                        >  2:22056/1‑81 (MQ=255)                      gAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCa                                                      <  2:22023/80‑1 (MQ=255)                        aTGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAAt                                                   <  1:22024/81‑1 (MQ=255)                         tGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATg                                                  >  1:22066/1‑81 (MQ=255)                                aCCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATg                                           >  1:22067/1‑81 (MQ=255)                                    gATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGa                                       >  2:22070/1‑81 (MQ=255)                                     aTCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTg                                        <  2:22025/79‑1 (MQ=255)                                        cTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCg                                   >  1:22100/1‑81 (MQ=255)                                        cTTGAAGAGAATGCGATCGGCGACCATAGCCGGGTGAAAAGGATCGCCAAGGAAAAGCTGCAAATGCAGCATGTTGATcc                                    <  1:22026/80‑1 (MQ=255)                                          tGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGt                                  >  1:22081/1‑80 (MQ=255)                                             agagAATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTcac                               >  1:22083/1‑80 (MQ=255)                                              gagaATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTcac                               >  2:22063/1‑79 (MQ=255)                                              gagaATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACaa                             <  1:22027/81‑1 (MQ=255)                                               agaATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTcac                               >  2:22085/1‑78 (MQ=255)                                                gaATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAg                            >  1:22068/1‑80 (MQ=255)                                                gaATGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGa                           <  2:22028/81‑1 (MQ=255)                                                  aTGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGaaa                         <  2:22030/81‑1 (MQ=255)                                                  aTGCGCTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGaaa                         <  1:22029/81‑1 (MQ=255)                                                    gcgcTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGaaaa                        >  1:22057/1‑80 (MQ=255)                                                     cgcTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAat                       <  1:22031/80‑1 (MQ=255)                                                      gcTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAatat                     >  1:22084/1‑81 (MQ=255)                                                      gcTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAatat                     >  2:22076/1‑81 (MQ=255)                                                       cTCGGCGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATc                    <  1:22032/81‑1 (MQ=255)                                                            cGACCATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATCGTAGt               >  2:22086/1‑81 (MQ=255)                                                                cATAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATCGTAGTGCa            >  1:22069/1‑80 (MQ=255)                                                                  tAGCCGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATCGTAGTGCaaaa         <  2:22033/81‑1 (MQ=255)                                                                  tAGCCGGGTGAAAAGGATAGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATCGTAGTGCaaaa         <  2:22034/81‑1 (MQ=255)                                                                      cGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATCGTAGTGCAAAAAt       >  1:22080/1‑79 (MQ=255)                                                                      cGGGTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATCGTAGTGCAAAAATaa     <  2:22035/81‑1 (MQ=255)                                                                       gggTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATCGTAGTGCAAAAATAAg    >  2:22090/1‑81 (MQ=255)                                                                         gTGAAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATCGTAGTGCAAAAATAAGGa  <  2:22036/81‑1 (MQ=255)                                                                            |                                                                              GAGCGAGATGCTTTAGACATTGAATGGCGCAACCTGATCCTTGAAGAGAATGCGCTCGGCGACCATAGCCGGGTGGAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATCGTAGTGCAAAAATAAGGA  >  NC_000913/91248‑91400

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