Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	333,763	A→G	L224P (CTG→CCG)	yahB ←	putative DNA‑bindng transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	333,763	0	A	G	100.0%	234.2 / NA	65	L224P (CTG→CCG)	yahB	putative DNA‑bindng transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/36);  total (29/36)

GGCGCGGATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCAGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCATAATATTAGGGTAGAGACGCAGTTGGCTTTCGGCGATGGGTTCCGGCACA  >  NC_000913/333632‑333882                                                                                                                                    |                                                                                                                        ggCGCGGATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGa                                                                                                                  >  1:1369555/1‑139 (MQ=255)          tATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGccc                                                                                                         <  1:137209/139‑1 (MQ=255)                  tGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTgcagca                                                                                                 >  2:1457751/1‑139 (MQ=255)                    tCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGccc                                                                                                         >  2:1141339/1‑129 (MQ=255)                    tCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGccc                                                                                                         <  1:1141339/129‑1 (MQ=255)                    tCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAAc                                                                                               >  2:1240085/1‑139 (MQ=255)                           aCCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGcc                                                                                          >  2:1141950/1‑137 (MQ=255)                           aCCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTCAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGcc                                                                                          <  1:1141950/137‑1 (MQ=255)                             cagcagGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACt                                                                                      <  1:836618/139‑1 (MQ=255)                              agcagGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACt                                                                                                                >  1:1000600/1‑112 (MQ=255)                              agcagGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACt                                                                                                                <  2:1000600/112‑1 (MQ=255)                                cagGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACtttt                                                                                   >  2:542174/1‑139 (MQ=255)                                     aTTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTttatt                                                                              >  1:997550/1‑139 (MQ=255)                                      ttGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTttatta                                                                             >  2:774015/1‑139 (MQ=255)                                       tGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTaa                                                                            >  2:33375/1‑139 (MQ=255)                                               cGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTgc                                                                                                     >  2:1306962/1‑106 (MQ=255)                                               cGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTgc                                                                                                     <  1:1306962/106‑1 (MQ=255)                                                gCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATgct                                                                   <  2:1158995/139‑2 (MQ=255)                                                gCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATgcg                                                                   <  1:1240085/139‑1 (MQ=255)                                                gCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATgcg                                                                   <  1:1215720/139‑1 (MQ=255)                                                gCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATgcg                                                                   <  1:1039509/139‑1 (MQ=255)                                                  cTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGcgcg                                                                 <  1:356979/139‑1 (MQ=255)                                                  cTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGcgcg                                                                 <  1:33375/139‑1 (MQ=255)                                                  cTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGcgcg                                                                 <  1:774015/139‑1 (MQ=255)                                                   tCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGTACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCgg                                                                <  1:1077852/139‑1 (MQ=255)                                                         aCCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGtcct                                                          <  1:1457751/139‑1 (MQ=255)                                                              gTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGtcct                                                          <  2:168824/134‑1 (MQ=255)                                                              gTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGtcct                                                          >  1:168824/1‑134 (MQ=255)                                                                 aTCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGgtg                                                              <  1:705503/127‑1 (MQ=255)                                                                 aTCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGgtg                                                              >  2:705503/1‑127 (MQ=255)                                                                    cGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCat                                                   <  1:616181/135‑1 (MQ=255)                                                                    cGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCat                                                   >  2:616181/1‑135 (MQ=255)                                                                     ggCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATgcg                                                                   >  1:675857/1‑118 (MQ=255)                                                                     ggCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATgcg                                                                   <  2:675857/118‑1 (MQ=255)                                                                     ggCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCATAatat                                              <  2:461234/139‑1 (MQ=255)                                                                       cAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTAt                                                                      >  2:415493/1‑113 (MQ=255)                                                                       cAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTAt                                                                      <  1:415493/113‑1 (MQ=255)                                                                       cAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGgtg                                                              <  2:470629/121‑1 (MQ=255)                                                                       cAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGgtg                                                              >  1:470629/1‑121 (MQ=255)                                                                                       tCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGcc                                                                                          >  2:516418/1‑77 (MQ=255)                                                                                       tCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGcc                                                                                          <  1:516418/77‑1 (MQ=255)                                                                                        cACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAAt                                                                                                                     >  2:445859/1‑49 (MQ=255)                                                                                        cACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAAt                                                                                                                     <  1:445859/49‑1 (MQ=255)                                                                                        cACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAAc                                                                                               <  2:228160/71‑1 (MQ=255)                                                                                        cACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAAc                                                                                               >  1:228160/1‑71 (MQ=255)                                                                                                    ggCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACt                                                                                                                <  1:1026360/42‑1 (MQ=255)                                                                                                    ggCACTGACATTTGGTGTTGAAATCTGGCACCGGAAATGACt                                                                                                                >  2:1026360/1‑42 (MQ=38)                                                                                                     gCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGtcctc                                                         <  1:330684/96‑1 (MQ=255)                                                                                                     gCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGtcctc                                                         >  2:330684/1‑96 (MQ=255)                                                                                                     gCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCat                                                   >  2:1528464/1‑102 (MQ=255)                                                                                                     gCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCat                                                   <  1:1528464/102‑1 (MQ=255)                                                                                                     gCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCat                                                   >  1:336352/1‑102 (MQ=255)                                                                                                     gCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCat                                                   <  2:336352/102‑1 (MQ=255)                                                                                                     gCACTGACATTTGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCATAATATTAGGGTAGAGACGCAGTTGGCTTTCGGCTATg              <  1:542174/139‑1 (MQ=255)                                                                                                               ttGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAAc                                                                                               >  2:1495807/1‑48 (MQ=255)                                                                                                               ttGGTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAAc                                                                                               <  1:1495807/48‑1 (MQ=255)                                                                                                                 ggTGTTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCATAATATTAGGGTAGAGACGCAGTTGGCTTTCGGCGATGGGTTCCGGcaca  <  2:997550/139‑1 (MQ=255)                                                                                                                   tgtTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTccac                                                      <  2:1108900/85‑1 (MQ=255)                                                                                                                   tgtTGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTccac                                                      >  1:1108900/1‑85 (MQ=255)                                                                                                                     ttGAAATGTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGcc                                                                                          <  1:1065296/47‑1 (MQ=255)                                                                                                                     ttGAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGcc                                                                                          >  2:1065296/1‑47 (MQ=255)                                                                                                                       gAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGcc                                                                                          >  2:412295/1‑45 (MQ=255)                                                                                                                       gAAATCTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGcc                                                                                          <  1:412295/45‑1 (MQ=255)                                                                                                                            cTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGgtg                                                              <  2:1232728/68‑1 (MQ=255)                                                                                                                            cTGGCACCGGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGgtg                                                              >  1:1232728/1‑68 (MQ=255)                                                                                                                                    |                                                                                                                        GGCGCGGATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCAGAATTGACTCCTGCCCGTGCAGCAACCAGCCCACTTTTTTATTAATCGTATGCGCGGTGTCCTCCACCATAATATTAGGGTAGAGACGCAGTTGGCTTTCGGCGATGGGTTCCGGCACA  >  NC_000913/333632‑333882

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