Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,716,526	G→A	intergenic (‑77/‑228)	grcA ← / → ung	stress‑induced alternate pyruvate formate‑lyase subunit/uracil‑DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,716,526	0	G	A	100.0%	125.5 / NA	36	intergenic (‑77/‑228)	grcA/ung	stress‑induced alternate pyruvate formate‑lyase subunit/uracil‑DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/19);  total (17/19)

AAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACATCAACTTATGCGGGTGTGAAATTTTACCAATTTACATTTT  >  NC_000913/2716394‑2716659                                                                                                                                     |                                                                                                                                      aaaGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTatatatat                                                                                                                                <  2:559247/140‑1 (MQ=255)               aTCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACAcc                                                                                                                  >  1:545795/1‑140 (MQ=255)                         gCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCC‑GTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCttt                                                                                                       >  1:151139/1‑140 (MQ=255)                          cGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACGATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCtt                                                                                                        >  2:607538/1‑139 (MQ=255)                          cGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACGATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCtt                                                                                                        <  1:607538/139‑1 (MQ=255)                                  gTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCC‑GTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAAt                                                                                              <  2:151139/140‑1 (MQ=255)                                       cTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAAt                                                                                          >  1:17009/1‑140 (MQ=255)                                           aTACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATaa                                                                                        >  1:307852/1‑138 (MQ=255)                                           aTACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATaa                                                                                        <  2:307852/138‑1 (MQ=255)                                             aCCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAAc                                                                                    <  2:1359872/140‑1 (MQ=255)                                             aCCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAAc                                                                                    <  1:775116/140‑1 (MQ=255)                                             aCCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAAc                                                                                    <  1:720786/140‑1 (MQ=255)                                               cTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGTCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACca                                                                                  <  2:972848/140‑1 (MQ=255)                                               cTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACca                                                                                  <  2:17009/140‑1 (MQ=255)                                                  tAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACcacac                                                                               >  2:340761/1‑140 (MQ=255)                                                   aaTCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACcacaca                                                                              <  1:340761/140‑1 (MQ=255)                                                    aTCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACAt                                                                             <  2:1199546/140‑1 (MQ=255)                                                         gtgATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACca                                                                                  >  2:1249258/1‑130 (MQ=255)                                                         gtgATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACca                                                                                  <  1:1249258/130‑1 (MQ=255)                                                         gtgATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGt                                                                        >  2:1055243/1‑140 (MQ=255)                                                         gtgATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGt                                                                        >  2:1325595/1‑140 (MQ=255)                                                         gtgATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAAAAAAAACCACACATCAAGt                                                                        >  1:1211765/1‑140 (MQ=255)                                                          tgATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACca                                                                                  >  1:382460/1‑129 (MQ=255)                                                          tgATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACca                                                                                  <  2:382460/129‑1 (MQ=255)                                                              gCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGt                                                                   <  2:545795/140‑1 (MQ=255)                                                               ccTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTc                                                                  >  1:28655/1‑140 (MQ=255)                                                                 tCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGc                                                                >  2:1251346/1‑140 (MQ=255)                                                                     cGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAAt                                                            <  2:1032095/140‑1 (MQ=255)                                                                                   gatttgTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGAAAATGGATTTTATTGttt                                              <  2:28655/140‑1 (MQ=255)                                                                                        gTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACAt                                         <  1:888762/140‑1 (MQ=255)                                                                                                  tCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACATCAACTTATGc                               <  2:191571/140‑1 (MQ=255)                                                                                                   cgcgGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACAt                                         <  1:376149/129‑1 (MQ=255)                                                                                                   cgcgGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACAt                                         >  2:376149/1‑129 (MQ=255)                                                                                                   cgcgGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACATCAACTTATGCg                              >  1:1322050/1‑140 (MQ=255)                                                                                                     cgGCCAATTGGTAAAACCATTGTTGCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACATCAACTTATGCggg                            <  1:1251346/140‑1 (MQ=255)                                                                                                                             gCTTGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACATCAACTTATGCGGGTGTGAAATTTTACCAATTTACAtt    >  2:390041/1‑140 (MQ=255)                                                                                                                               ttGAGTATATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACATCAACTTATGCGGGTGTGAAATTTTACCAATTTACAtttt  >  2:1442594/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                      AAAGAGTTCAGCAGATCGTCGTTAGCGGCTTTAGTAATCTGGATACCTGTAATCATGTGATGCCTCCCCGGCAAAATTATTTGATTTGTTCAGCCTGTCGCGGCCAATTGGTAAAACCATTGTTGCTTGAGTGTATATATACTCCTCAAACACCCTTGAATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACATCAACTTATGCGGGTGTGAAATTTTACCAATTTACATTTT  >  NC_000913/2716394‑2716659

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