Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,813,482	A→G	V61A (GTG→GCG)	cedA ←	cell division modulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,813,482	0	A	G	100.0%	201.5 / NA	56	V61A (GTG→GCG)	cedA	cell division modulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/28);  total (28/28)

TTTGAATGAATGTAGCAACGAAATACAGAATTTCAGGTCATGTAACTCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCACAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCACGAAACGAATCCATCTTTATCGCATGTTCTGGCGGCGCGGGTTCCGTGCGTGGG  >  NC_000913/1813346‑1813616                                                                                                                                         |                                                                                                                                       tttGAATGAATGTAGCAACGAAATACAGAATTTCAGGTCATGTAACTCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGcgca                                                                                                                                      >  1:1555742/1‑139 (MQ=255)     aatgaatgTAGCAACGAAATACAGAATTTCAGGTCATGTAACTCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTaa                                                                                                                                  >  1:138553/1‑139 (MQ=255)     aatgaatgTAGCAACGAAATACAGAATTTCAGGTCATGTAACTCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTaa                                                                                                                                  <  1:746523/139‑1 (MQ=255)           tgTAGCAACGAAATACAGAATTTCAGGTCATGTAACTCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCg                                                                                                                            >  1:1274662/1‑139 (MQ=255)               gCAACGAAATACAGAATTTCAGGTCATGTAACTCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGa                                                                                                                        >  2:867610/1‑139 (MQ=255)                  aCGAAATACAGAATTTCAGGTCATGTAACTCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAgcg                                                                                                                     <  1:743193/139‑1 (MQ=255)                   cGAAATACAGAATTTCAGGTCATGTAACTCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAgcgc                                                                                                                    >  2:508308/1‑139 (MQ=255)                         aCAGAATTTCAGGTCATGTAACTCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAAt                                                                                                              <  2:567923/139‑1 (MQ=255)                           aGAATTTCAGGTCATGTAACTCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGa                                                                                                            <  2:1274662/139‑1 (MQ=255)                               tttCAGGTCATGTAACTCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGActct                                                                                                        <  1:824712/139‑1 (MQ=255)                                         tGTAACTCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGAct                                                                                                          >  1:386042/1‑127 (MQ=255)                                         tGTAACTCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGAct                                                                                                          <  2:386042/127‑1 (MQ=255)                                         tGTAACTCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTa                                                                                              >  1:1001153/1‑139 (MQ=255)                                          gTAACTCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTAc                                                                                             <  1:815628/139‑1 (MQ=255)                                           tAACTCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACa                                                                                            <  1:818586/139‑1 (MQ=255)                                             aCTCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACaaa                                                                                          <  2:138553/139‑1 (MQ=255)                                               tCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACaaa                                                                                          >  2:646564/1‑137 (MQ=255)                                               tCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACaaa                                                                                          <  1:646564/137‑1 (MQ=255)                                                cccGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGc                                                                                       <  2:1001153/139‑1 (MQ=255)                                                 ccGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCa                                                                                      >  2:405254/1‑139 (MQ=255)                                                   ggCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAAc                                                                                    >  2:1658011/1‑139 (MQ=255)                                                          ccGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTg                                                                             <  2:1555742/139‑1 (MQ=255)                                                                ggTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACgc                                                                       <  1:1658011/139‑1 (MQ=255)                                                                 gTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACgca                                                                      <  1:1426280/139‑1 (MQ=255)                                                                            ttACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACgc                                                                       <  2:1306255/127‑1 (MQ=255)                                                                            ttACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACgc                                                                       >  1:1306255/1‑127 (MQ=255)                                                                            ttACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACa                                                           >  2:871056/1‑139 (MQ=255)                                                                             tACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACAt                                                          >  2:157976/1‑139 (MQ=255)                                                                              aCTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATc                                                         >  1:172016/1‑139 (MQ=255)                                                                                  aGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACgc                                                                       <  2:571322/121‑1 (MQ=255)                                                                                  aGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACgc                                                                       >  1:571322/1‑121 (MQ=255)                                                                                  aGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCacga                                                     >  1:1538662/1‑139 (MQ=255)                                                                                    tCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCacgaaa                                                   >  1:852169/1‑139 (MQ=255)                                                                                    tCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCacgaaa                                                   >  1:702821/1‑139 (MQ=255)                                                                                     cACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCacgaaac                                                  >  2:716966/1‑139 (MQ=255)                                                                                          ccccTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTAc                                                                                             >  2:746549/1‑91 (MQ=255)                                                                                          ccccTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTAc                                                                                             <  1:746549/91‑1 (MQ=255)                                                                                           cccTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCACGAAACGAATcc                                            <  1:508308/139‑1 (MQ=255)                                                                                                ccTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCACGAAACGAATCCATCtt                                       >  1:671269/1‑139 (MQ=255)                                                                                                 cTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCACGAAACGAATCCATCttt                                      <  1:157976/139‑1 (MQ=255)                                                                                                    gCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCACGAAACGAATCCATCTTTATc                                   >  1:965835/1‑139 (MQ=255)                                                                                                     cGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCACGAAACGAATCCATCTTTATCg                                  >  2:618465/1‑139 (MQ=255)                                                                                                       gATCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCACGAAACGAATCCATCTTTATCGCa                                <  2:702821/139‑1 (MQ=255)                                                                                                         tCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATc                                                                                                  >  1:704948/1‑71 (MQ=255)                                                                                                         tCTGATTTGCCCAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATc                                                                                                  <  2:704948/71‑1 (MQ=255)                                                                                                                  cccAACGTTGGGCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCACGAAACGAATCCATCTTTATCGCATGTTCTggcgg                     >  2:504257/1‑139 (MQ=255)                                                                                                                            ggCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCACGAAACGAATCCATCTTTATCGCATGTTCTg                         >  1:779571/1‑125 (MQ=255)                                                                                                                            ggCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCACGAAACGAATCCATCTTTATCGCATGTTCTg                         <  2:779571/125‑1 (MQ=255)                                                                                                                            ggCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCACGAAACGAATCCATCTTTATCGCATGTTCTGGCGGCGCGGGTTcc           >  1:1467086/1‑139 (MQ=255)                                                                                                                            ggCAGATTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCACGAAACGAATCCATCTTTATCGCATGTTCTGGCGGCGCGGGTTcc           <  1:405254/139‑1 (MQ=255)                                                                                                                                 aTTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGc                                                                                       <  2:79623/58‑1 (MQ=255)                                                                                                                                 aTTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGc                                                                                       >  1:79623/1‑58 (MQ=255)                                                                                                                                 aTTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGcc                                                                           <  1:436530/70‑1 (MQ=255)                                                                                                                                 aTTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGcc                                                                           >  2:436530/1‑70 (MQ=255)                                                                                                                                 aTTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCACGAAACGAATCCATCTTTATCGCATGTTCTGGCGGCGCGGGTTCcgtgcg      <  2:671269/139‑1 (MQ=255)                                                                                                                                 aTTCAGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCACGAAACGAATCCATCTTTATCGCATGTTCTGGCGGCGCGGGTTCcgtgcg      <  1:504257/139‑1 (MQ=255)                                                                                                                                     aGGCGCAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCACGAAACGAATCCATCTTTATCGCATGTTCTGGCGGCGCGGGTTCCGTGCGTggg  <  2:1538662/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       TTTGAATGAATGTAGCAACGAAATACAGAATTTCAGGTCATGTAACTCCCGGCAAAACCGGGAGGTATGTAATCCTTACTCAGTCACTTCCCCTTCCTGGCGGATCTGATTTGCCCAACGTTGGGCAGATTCAGGCACAGTAAACGCCGGTGAGCGCAGAAATGACTCTCCCATCAGTACAAACGCAACATATTTGCCACGCAGCATCCAGACATCACGAAACGAATCCATCTTTATCGCATGTTCTGGCGGCGCGGGTTCCGTGCGTGGG  >  NC_000913/1813346‑1813616

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