Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,165,608	A→G	intergenic (+87/‑60)	yegP → / → yegQ	DUF1508 domain‑containing protein YegP/putative peptidase YegQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,165,608	0	A	G	100.0%	63.6 / NA	20	intergenic (+87/‑60)	yegP/yegQ	DUF1508 domain‑containing protein YegP/putative peptidase YegQ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/9);  total (11/9)

GCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATG  >  NC_000913/2165496‑2165743                                                                                                                 |                                                                                                                                        gCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCGCTCCCCTAACCGCTTCATCAGGTGAAg                                                                                                              <  1:555033/140‑1 (MQ=255)    aGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTGCCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCGCTCCCCTAACCGCTTCATCAGGTGAAg                                                                                                              >  1:115394/1‑137 (MQ=255)    aGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTGCCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCGCTCCCCTAACCGCTTCATCAGGTGAAg                                                                                                              <  2:115394/137‑1 (MQ=255)     gACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCGCTCCCCTAACCGCTTCATCAGGTGAAGCGGa                                                                                                          >  2:393270/1‑140 (MQ=255)           gAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCGCTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGAc                                                                                                    >  1:379553/1‑140 (MQ=255)           gAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCGCTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGAc                                                                                                    <  1:712191/140‑1 (MQ=255)                       gTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCGCTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGaa                                                                                        <  2:366791/140‑1 (MQ=255)                          aTTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCGCTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACga                                                                                     <  2:575798/140‑1 (MQ=255)                           ttACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCGCTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACgag                                                                                    >  1:30176/1‑140 (MQ=255)                                                       tCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCGCTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCttt                                                        >  2:561164/1‑140 (MQ=255)                                                                             ccACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCGCTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGaaaaa                                  >  1:708220/1‑140 (MQ=255)                                                                               acaGCGCTACAATGCCCGCCCTTAAAGTGGGGGCGCTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAata                                <  1:393270/140‑1 (MQ=255)                                                                                      tACAATGCCCGCCCTTAAAGTGGGGGCGCTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTa                         <  1:218326/140‑1 (MQ=255)                                                                                             cccgcccTTAAAGTGGGGGCGCTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTc                  >  1:492940/1‑140 (MQ=255)                                                                                             cccgcccTTAAAGTGGGGGCGCTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTc                  >  1:256461/1‑140 (MQ=255)                                                                                                            ggggCGCTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGgcg        >  2:476985/1‑135 (MQ=255)                                                                                                            ggggCGCTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGgcg        <  1:476985/135‑1 (MQ=255)                                                                                                            ggggCGCTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGAt   >  2:190760/1‑140 (MQ=255)                                                                                                            ggggCGCTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGAt   >  2:688535/1‑140 (MQ=255)                                                                                                             gggCGCTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATg  <  2:256461/140‑1 (MQ=255)                                                                                                                 |                                                                                                                                        GCCAGACGGTGAAAGACAATACGTAATTACGATGCCGGGCGTGGGTTGTTGGCGTCCGGCGCTTTACCAGAGCCCCCCACAGCGCTACAATGCCCGCCCTTAAAGTGGGGGCACTCCCCTAACCGCTTCATCAGGTGAAGCGGATCTGACCTGTCATCAGAACGAGAGAATTATGTTTAAACCGGAACTCCTTTCCCCGGCGGGAACGCTGAAAAATATGCGTTACGCTTTCGCTTATGGCGCAGATG  >  NC_000913/2165496‑2165743

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