Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,491,137	T→C	E109E (GAA→GAG)	fic ←	stationary phase‑induced protein, putative toxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,491,137	0	T	C	100.0%	157.5 / NA	44	E109E (GAA→GAG)	fic	stationary phase‑induced protein, putative toxin
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/24);  total (20/24)

GAAGATCCGCTGTGCCAGACCACTTCCCACCCGGAAGGGATGCAGCACGTTGATTTCACAATAGTAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCTTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCAGAACGGCGTATCACCCTGATAAATATCAACTTCACGCAGTTGCCCTGCCCAGTCGAAAATATCCTGATA  >  NC_000913/3491020‑3491265                                                                                                                      |                                                                                                                                 gaagaTCCGCTGTGCCAGACCACTTCCCACCCGGAAGGGATGCAGCACGTTGATTTCACAATAGTAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATc                                                                                                             >  1:923876/1‑139 (MQ=255)            gtgCCAGACCACTTCCCACCCGGAAGGGATGCAGCACGTTGATTTCACAATAGTAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGcc                                                                                                  <  2:103197/139‑1 (MQ=255)               ccAGACCACTTCCCACCCGGAAGGGATGCAGCACGTTGATTTCACAATAGTAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATc                                                                                                                    <  2:1013459/118‑1 (MQ=255)               ccAGACCACTTCCCACCCGGAAGGGATGCAGCACGTTGATTTCACAATAGTAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATc                                                                                                                    >  1:1013459/1‑118 (MQ=255)                    ccACTTCCCACCCGGAAGGGATGCAGCACGTTGATTTCACAATAGTAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACcctct                                                                                                                               <  1:292999/102‑1 (MQ=255)                    ccACTTCCCACCCGGAAGGGATGCAGCACGTTGATTTCACAATAGTAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACcctct                                                                                                                               >  2:292999/1‑102 (MQ=255)                       cTTCCCACCCGGAAGGGATGCAGCACGTTGATTTCACAATAGTAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCa                                                                                                            >  1:459455/1‑118 (MQ=255)                       cTTCCCACCCGGAAGGGATGCAGCACGTTGATTTCACAATAGTAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCa                                                                                                            <  2:459455/118‑1 (MQ=255)                          cccACCCGGAAGGGATGCAGCACGTTGATTTCACAATAGTAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATc                                                                                                                    >  2:168409/1‑107 (MQ=255)                          cccACCCGGAAGGGATGCAGCACGTTGATTTCACAATAGTAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATc                                                                                                                    <  1:168409/107‑1 (MQ=255)                             aCCCGGAAGGGATGCAGCACGTTGATTTCACAATAGTAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGATATAAg                                                                                 >  2:991320/1‑139 (MQ=255)                               ccGGAAGGGATGCAGCACGTTGATTTCACAATAGTAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCa                                                                                                            >  1:74379/1‑110 (MQ=255)                               ccGGAAGGGATGCAGCACGTTGATTTCACAATAGTAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCa                                                                                                            <  2:74379/110‑1 (MQ=255)                               ccGGAAGGGATGCAGCACGTTGATTTCACAATAGTAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGATATAAGCa                                                                               <  2:312903/139‑1 (MQ=255)                                                                 tAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCAt                                                                                                              <  1:217142/74‑1 (MQ=255)                                                                 tAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCAt                                                                                                              >  2:217142/1‑74 (MQ=255)                                                                 tAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCAGAACGGCGTATCACCCTGATAAatat                                             <  1:814998/139‑1 (MQ=255)                                                                 tAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCAGAACGGCGTATCACCCTGATAAatat                                             >  2:995864/1‑139 (MQ=255)                                                                   aTGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCa                                                                       >  2:942362/1‑111 (MQ=255)                                                                   aTGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCa                                                                       <  1:942362/111‑1 (MQ=255)                                                                    tGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACcctctt                                                                                                                              <  2:954838/55‑1 (MQ=255)                                                                    tGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACcctctt                                                                                                                              >  1:954838/1‑55 (MQ=255)                                                                      cgcCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCAt                                                                                                              <  2:827287/69‑1 (MQ=255)                                                                      cgcCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCAt                                                                                                              >  1:827287/1‑69 (MQ=255)                                                                        ccAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCct                                                                                                >  2:283462/1‑81 (MQ=255)                                                                        ccAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCct                                                                                                <  1:283462/81‑1 (MQ=255)                                                                          aGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGcc                                                                                                  >  1:760594/1‑77 (MQ=255)                                                                          aGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGcc                                                                                                  <  2:760594/77‑1 (MQ=255)                                                                           gCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACcctct                                                                                                                               >  2:418964/1‑47 (MQ=255)                                                                           gCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACcctct                                                                                                                               <  1:418964/47‑1 (MQ=255)                                                                           gCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCAGAACGGCGTATCACCCTGATAAATATCAACTTCACg                                   >  2:627600/1‑139 (MQ=255)                                                                           gCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGAGCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCAGAACGGCGTATCACCCTGATAAATATCAACTTCACg                                   >  2:21591/1‑139 (MQ=255)                                                                            ccGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGGTCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCAGAACGGCGTATCACCCTCATAAATATCAACTTCACGc                                  <  2:923876/139‑1 (MQ=255)                                                                              gCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCAGAACGGCGTATCACCCTGATAAATATCAACTTCACGCAg                                <  1:995864/139‑1 (MQ=255)                                                                              gCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCAGAACGGCGTATCACCCTGATAAATATCAACTTCACGCAg                                <  1:643236/139‑1 (MQ=255)                                                                              gCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCAGAACGGCGTATCACCCTGATAAATATCAACTTCACGCAg                                <  1:564406/139‑1 (MQ=255)                                                                                        cTTCGCTTTCTCCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGAGCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCAGAACGGCGTATCACCCTGATAAATATCAACTTCACGCAGTtgccctgcc                      <  1:21591/139‑1 (MQ=255)                                                                                                 ctcCAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCAGAACGGCGTATCACCCTGATAAATATCAACTTCACGCAGTTGCCCTGCCCAGTCGaaa             <  1:627600/139‑1 (MQ=255)                                                                                                    cAGGCCAACCAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCAGAACGGCGTATCACCCTGATAAATATCAACTTCACGCAGTTGCCCTGCCCAGTCGAAAata          <  1:492601/139‑1 (MQ=255)                                                                                                            ccAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGc                                                                                                   >  1:65913/1‑42 (MQ=255)                                                                                                            ccAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGc                                                                                                   <  2:65913/42‑1 (MQ=255)                                                                                                            ccAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCct                                                                                                >  2:427086/1‑45 (MQ=255)                                                                                                            ccAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCct                                                                                                <  1:427086/45‑1 (MQ=255)                                                                                                            ccAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGat                                                                                      >  2:191189/1‑55 (MQ=255)                                                                                                            ccAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGat                                                                                      <  1:191189/55‑1 (MQ=255)                                                                                                            ccAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCAGAACGGCGTATCACCCTGa                                                    <  2:371336/89‑1 (MQ=255)                                                                                                            ccAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCAGAACGGCGTATCACCCTGa                                                    >  1:371336/1‑89 (MQ=255)                                                                                                            ccAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCAGAACGGCGTATCACCCTGATAAATATCAACTTCACg                                   >  2:1108165/1‑106 (MQ=255)                                                                                                            ccAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCAGAACGGCGTATCACCCTGATAAATATCAACTTCACg                                   <  1:1108165/106‑1 (MQ=255)                                                                                                            ccAGATAACCCTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCAGAACGGCGTATCACCCTGATAAATATCAACTTCACGCAGTTGCCCTGCCCAGTCGAAAATATCCTGATa  >  1:430824/1‑139 (MQ=255)                                                                                                                      |                                                                                                                                 GAAGATCCGCTGTGCCAGACCACTTCCCACCCGGAAGGGATGCAGCACGTTGATTTCACAATAGTAATGCGCCAGCCGCTCGACGAACTTCGCTTTCTCCAGGCCAACCAGATAACCTTCTTCCTCCAGATCCTGCATCAGGGCATTGCCCTCTTTTTCGATATAAGCAAAGTGGCAGAACGGCGTATCACCCTGATAAATATCAACTTCACGCAGTTGCCCTGCCCAGTCGAAAATATCCTGATA  >  NC_000913/3491020‑3491265

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