Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,150,418	A→G	intergenic (‑153/+29)	eptC ← / ← ppc	phosphoethanolamine transferase EptC/phosphoenolpyruvate carboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,150,418	0	A	G	100.0%	60.8 / NA	21	intergenic (‑153/+29)	eptC/ppc	phosphoethanolamine transferase EptC/phosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/11);  total (10/11)

AAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGACTCTGCAATCCACGGCAGATCGGCCATCAGATGGGAATCGTTGGC  >  NC_000913/4150188‑4150662                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     aaaCCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCagaa                                                                                                                                                                                                                                  <  1:123996/251‑1 (MQ=255)                                                        gCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCa                                                                                                                                                                           <  1:193425/251‑1 (MQ=255)                                                            ggACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTa                                                                                                                                                                        >  1:80681/1‑250 (MQ=255)                                                                            cATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGCTATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGAcgcg                                                                                                                                                        >  1:159285/1‑250 (MQ=255)                                                                                                aaaaaaTATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCcttc                                                                                                                                   <  1:87242/251‑1 (MQ=255)                                                                                                     aaatCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAAAATATTGAAAACGAGGGTGTAAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGAACGAGGGTTGGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCttttt                                                                                                                              <  2:156912/249‑1 (MQ=255)                                                                                                             gAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTACCGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTgg                                                                                                                      >  2:64335/1‑251 (MQ=255)                                                                                                                       aCAATACCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGATTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCACTCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGGTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGt                                                                                                            >  1:161818/1‑251 (MQ=255)                                                                                                                              cTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTgcagca                                                                                                      >  2:199536/1‑250 (MQ=255)                                                                                                                                tAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAAcc                                                                                                   >  2:235690/1‑250 (MQ=255)                                                                                                                                 aaGGATATCGGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATAGTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCCCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACt                                                                                                   <  2:194086/250‑1 (MQ=255)                                                                                                                                              ggTATATTCAGAATTTGAAAAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACg                                                                                     <  2:87394/251‑1 (MQ=255)                                                                                                                                                        aaaTTTAAATAAAATGCAGACAAAAATATATTGAAACCGCGGGTGTAAGAACAGAAATTTTTCAGAAACCCCTCTCGCAGAAGCACGAGTTTTTGCAAAAGAGGAAGTTTAACCGTTTTTCCGCATCCCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTTCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCATTACGTTCAGCgcgt                                                                           <  2:185039/250‑4 (MQ=255)                                                                                                                                                           tttGAATAAAATGCAGACAGAAATATATTGAAAACGTGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCgg                                                                        >  1:47512/1‑251 (MQ=255)                                                                                                                                                              gAATAAAATGCAGACAGAAATATATGGAAAACGAGGGTGTGAAAACAGAAGTATTTCAAAAAACCCCCGGGCAGAAGCGCGAGGTTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGgtgt                                                                     <  2:168214/251‑1 (MQ=255)                                                                                                                                                                           aaggaataatatTAAAAACGAGGGTGTTAGAAAAGAAGTATGTCAGAAACCCCTCGCGAAGAAGCACGAGTGTTTGCAGAAGAGGAAGATTAGCCGTTATTACGCATCCCTGCCGCAATCCCGGCAATAGTGACCATTAACGCGTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAg                                                        <  2:182984/243‑1 (MQ=255)                                                                                                                                                                                   atatTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGATTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGTAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAAt                                                  >  2:18720/1‑249 (MQ=255)                                                                                                                                                                                        gAAAACGAGGGTGTGAGAACAGAAGATTTTCAGAAACCCCTCGCGCAGATCCACGAGGTTTTGCAGAAGAGGAAGATTAGCCGTTATTACGTATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCGGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGActct                                           <  2:155627/251‑1 (MQ=255)                                                                                                                                                                                             cGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGACTCTGCAAt                                      >  1:106508/1‑251 (MQ=255)                                                                                                                                                                                                       aGAACAGAAGTATTTCAAAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATCCCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGACTCTGCAATCCACGGCAGa                            <  1:211677/251‑1 (MQ=255)                                                                                                                                                                                                                                 ctcgcAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCTGCGGGTCGGTGTAAATATTCCGTCGCTGAATAGACTCTGCAATCCACGGCAGATCGGCCATCAGATGGGAATAGTTGGc  >  2:206888/3‑251 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     AAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGACTCTGCAATCCACGGCAGATCGGCCATCAGATGGGAATCGTTGGC  >  NC_000913/4150188‑4150662

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