Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	339,376	G→A	G351E (GGG→GAG)	yahG →	DUF1116 domain‑containing protein YahG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	339,376	0	G	A	100.0%	72.3 / NA	25	G351E (GGG→GAG)	yahG	DUF1116 domain‑containing protein YahG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/12);  total (13/12)

GACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGGGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGAT  >  NC_000913/339092‑339664                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                 gACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCAACATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGa                                                                                                                                                                                                                                                                                   >  2:405508/1‑300 (MQ=255)         ggATTATTCAGGCGGGTTATTCCGTCGAGCANCAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCANGGCATCGAATACAGCACCGTGGTAACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCAGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCACGATGTTTGCCGGTTATAAGCAGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCa                                                                                                                                                                                                                                                                           >  2:459213/1‑300 (MQ=255)            ttattCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCa                                                                                                                                                                                                                                                                        >  1:103770/1‑300 (MQ=255)                    gCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAAcc                                                                                                                                                                                                                                                                <  2:657265/300‑1 (MQ=255)                                     gCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTACGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCGTGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGAGGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGAtt                                                                                                                                                                                                                                               >  2:206063/1‑300 (MQ=255)                                                  gtgtTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATTGCGAct                                                                                                                                                                                                                                  >  2:223536/1‑299 (MQ=255)                                                    gtTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACcgc                                                                                                                                                                                                                                >  1:603134/1‑300 (MQ=255)                                                    gtTTGAGTTTGTCGCCAGCAGCGACTACTTCGCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGNCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGNGACcgc                                                                                                                                                                                                                                >  2:553864/1‑300 (MQ=255)                                                              gTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTAtc                                                                                                                                                                                                                      <  1:166855/300‑1 (MQ=255)                                                                  ccAGCAGCGACTACTTCTCCGGCCCGACGTGGATGTCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGTTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGAGGGGCTATGGCGACCGCGCCCGCTAtcgtcg                                                                                                                                                                                                                  <  2:204510/300‑1 (MQ=255)                                                                                  ctcCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCAc                                                                                                                                                                                                  <  1:498943/300‑1 (MQ=255)                                                                                          cGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGgaag                                                                                                                                                                                          <  1:265645/300‑1 (MQ=255)                                                                                                                               ggcgCATGGCATCGAATACAGCACCGTGGTCACCACCATTGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGGTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCAc                                                                                                                                                     <  1:293076/300‑1 (MQ=255)                                                                                                                                  gcATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCTGGTTATAAGCCGGAAGATTCGGAGCTGGATATTGGCGACAGCGCCATCACCGAAACCTACGGTATTGGTGGATTTGCTATGGCGACTGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTACCGTCAGATGCGCGAAATCACCCt                                                                                                                                                  >  2:310788/1‑300 (MQ=255)                                                                                                                                           cGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGaaaa                                                                                                                                         <  1:262428/300‑1 (MQ=255)                                                                                                                                              aTACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCACGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTCTCGTCGAGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCCCCCTCGGTGAAAAccc                                                                                                                                      >  2:411151/1‑300 (MQ=255)                                                                                                                                                   gCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACg                                                                                                                                 <  1:342091/300‑1 (MQ=255)                                                                                                                                                          ggTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCACCGTCACCAt                                                                                                                          >  2:356146/1‑300 (MQ=255)                                                                                                                                                                              aCGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATgg                                                                                                      <  1:402492/300‑1 (MQ=255)                                                                                                                                                                                               cTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATc                                                                                     >  1:487242/1‑300 (MQ=255)                                                                                                                                                                                                                           ggTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCa                                                         >  2:558149/1‑300 (MQ=255)                                                                                                                                                                                                                                                               tGTTTGCCGGTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAg                     <  1:520591/300‑1 (MQ=255)                                                                                                                                                                                                                                                                        gTTATAAGCCGGAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCg            >  2:191328/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                  ggAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGAt  <  1:173851/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                  ggAAGATTCGGAGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTAGTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGTTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGag  <  1:489083/300‑2 (MQ=255)                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                 GACGCCGGGGATTATTCAGGCGGGTTATTCCGTCGAGCAACAGCGCGAAGTGTTTGAGTTTGTCGCCAGCAGCGACTACTTCTCCGGCCCGACGTGGATGGCGATGTGTAAAGCGGCGATGGATGCGGCGCATGGCATCGAATACAGCACCGTGGTCACCACCATGGCGCGTAACGGCGTCGAGTTCGGCCTGCGGGTCAGCGGCCTGCCGGGGCAATGGTTTACCGGCCCGGCGCAGCAGGTGATCGGCCCGATGTTTGCCGGTTATAAGCCGGAAGATTCGGGGCTGGATATCGGCGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGAT  >  NC_000913/339092‑339664

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