Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,905,434	G→A	intergenic (‑16/‑123)	queE ← / → yqcG	putative 7‑carboxy‑7‑deazaguanine synthase QueE/cell envelope stress response protein YqcG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,905,434	0	G	A	100.0%	23.1 / NA	13	intergenic (‑16/‑123)	queE/yqcG	putative 7‑carboxy‑7‑deazaguanine synthase QueE/cell envelope stress response protein YqcG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/4);  total (9/4)

GCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAGAATGCTGAAAAGGGAGACTTCCCGATCCTCAAGCTTTTCCCAGGTGTGTTTGGTGTCGCACCAGGCACAGCCAACCGGGCATCCCTGTAAACGAATAAAAATGGCGGGAACGCCGGTAAAGTAACCCTCACCTTGCAGGGTCTGGAACATCTCGTTAATCGGGTACTGCATAGCATTCTCTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTTAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCAAAATTATATTTATTTTTTCTGTATTAGTTTTTAATGATA  >  NC_000913/2905213‑2905647                                                                                                                                                                                                                              |                                                                                                                                                                                                                      nccccccACTTATCACTCTCTTTGGTCTTCGCCAGAATGCTGAAAAGGGAGACTTCCCGATCCTCAAGCTTTTCCCAGGTGTGTTTGGTGTCGCACCAGGCACAGCCAACCGGGCATCCCTGTAAACGAATAAAAATGGCGGGAACGCCGGTAAAGTAACCCTCACCTTGCAGGGTCTGGAACATCTCGTTAATCGGGTACTGCATAGCATTCTCTGTGAAATGGGTAATTGTTAATTATTGCAGATCCTg                                                                                                                                                                                          >  1:129492/2‑251 (MQ=255)                                                                                          ntCGCACCAGGCACAGCCAACCGGGCATCCCTGTAAACGAATAAAAATGGCGGGAACGCCGGGAAAGTAACCCTCACCTTGCAGGGTCTGGAACATCTCGTTAATCGGGTACTGCATAGCATTCTCTGTGAAATGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTTAATAATTAAAGGGAGTAACg                                                                                                 >  1:62701/2‑251 (MQ=255)                                                                                             gCACCAGGCACAGCCAACCGGGCATCCCTGTAAACGAATAAAAATGGCGGCAACGCCGGTAAAGTAAACCTCACCTTGCAGGGTCTGGAACATCTCGTTAATCGGGTAATGCATCGCATTCTCTGTGAAATGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGaaaaa                                                                                                                                     >  1:151263/1‑212 (MQ=255)                                                                                                   ggCACAGCCAACCGGGCATCCCTGTACACGAATAAAAAAGACGGGAACGCCGCCAAAGTAACCCTCACCTTGCCGGGGCTGGAACATCTCGTTAATCGGGTACTACCTAGCATTCTCTGTGCAATGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTTAATAATTAAAGGGAGTAACGTATTATGTc                                                                                        >  1:107328/1‑251 (MQ=255)                                                                                                   ggCACAGCCAACCGGGCATCCCTGTAAACGAATAAAAATGGCGGGAACGCCGGTAAAGTAACCCTCACCTTGCCGGGTCTGGGACATCTCGTTAATCGGGTCCTGCGTAGCATTCTCTGTGAAATGGATCATTGTTACTTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTTAATAATTAAAGGGAGTAACGTCTTATGTc                                                                                        >  1:81883/1‑251 (MQ=255)                                                                                                     cacaGCCAACCGGGCATCCCTGTAAACGAATAAAAATGGCGGGAACGCCGGTAAAGTAACCCTCACCTTGCAGGGTCTGGAACATCTCGTTAATCGGGTACTGCATAGCATTCTCTGTGAAATGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGaaaaa                                                                                                                                     >  1:2030/1‑204 (MQ=255)                                                                                                                               cGAATAAAAATGGCGGGACCGCCGGCCAAGTAACCCTCACCTTGCAGGGACTTGAACCTATCGTTAATCTGGTACCGCCTAGCATTCACTGTGACATGGATAAATGTTAATTATTGCAGATCCTGCCACAACAATCCTGTCTTATTAACATTCTGTTAAAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTTAATAATTAAAGGGAGTACCGTATTATGTCAGAAGAAAATACAGAAAATGGATTTAAt                                                            >  1:162319/1‑251 (MQ=255)                                                                                                                                                       nnnnnntANNNNNCACNTTGCAGGGTCTGGAACATCTCGTTAATCGGGTACTGCATAGCATTCTCTGTGAAATGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGa                                                                                                                                         <  2:129492/144‑1 (MQ=255)                                                                                                                                                                                   ggAACATCTCGTTAATCGGGTACTGCATAGCATTCTCTGTGACATGGATGATTGATAATTAATGCAGATCATGCCACAACACTCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACCTCTCGTTCATAATTAAAGGGAGTAACGTATTATGTCAGAAGCAAATAAAGAAACTGGATTTAATCATGTCAAAACATTCACCAAAATTATATTTATTTTTTCTGTATTAGttttt         >  1:199434/1‑250 (MQ=255)                                                                                                                                                                                   ggAACATCTCGTTAATCGGGTAATGCATAGCATTCTCTGCGAAATGGATAATTGTTAATTATTGCAGATCCTGACACACCAATCATGTCTTATTAACATACTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTTCATAATTAAAGGGAGTAACGTATTATGTAAGAAGAAAATAAAGAAAATGGATTTAATAATGTCAAAAAATTCACCAAAATTATATTTATTTTTTCTGTAtt                >  1:161499/1‑243 (MQ=255)                                                                                                                                                                                            ctttAATCGGGTACTGCATAGCATTCTCTGTTAAATTGATAAGTGTTAATTAGTGCATAGCCTGCCACAACAATCATGTCTTAGTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTTAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGTATTTAATCATGGCAAACCATTCACCAAAATTGTATGGATATTTTCGGGAGTAGGTTTTAATGATa  <  1:68674/246‑1 (MQ=255)                                                                                                                                                                                                             anacnnttCTNTGTGAAATGGATAATTGTTAATTANTGCAGATCCTGCCACAACAATCATGTCTTATNANCATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTTAANAATTAAAGGGNGTANNGTATNANGNCAGAAGAAAATAAAGAAAATGGATTTAATNATGTCAAAACATTCACCAAAATTNNATNNNTTNNTTCTg                    <  2:156133/207‑1 (MQ=255)                                                                                                                                                                                                                 cATTCTNNGTGAAATGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGaaaaa                                                                                                                                     <  2:2030/96‑1 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                                      GCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAGAATGCTGAAAAGGGAGACTTCCCGATCCTCAAGCTTTTCCCAGGTGTGTTTGGTGTCGCACCAGGCACAGCCAACCGGGCATCCCTGTAAACGAATAAAAATGGCGGGAACGCCGGTAAAGTAACCCTCACCTTGCAGGGTCTGGAACATCTCGTTAATCGGGTACTGCATAGCATTCTCTGTGAAGTGGATAATTGTTAATTATTGCAGATCCTGCCACAACAATCATGTCTTATTAACATTCTGTTACAGGCAGGTTTAAAGAAAAATGCGAAACATATCGTTAATAATTAAAGGGAGTAACGTATTATGTCAGAAGAAAATAAAGAAAATGGATTTAATCATGTCAAAACATTCACCAAAATTATATTTATTTTTTCTGTATTAGTTTTTAATGATA  >  NC_000913/2905213‑2905647

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