Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,239,087	A→G	D485G (GAC→GGC)	yjbH →	YjbH family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,239,087	0	A	G	100.0%	162.1 / NA	46	D485G (GAC→GGC)	yjbH	YjbH family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/24);  total (22/24)

TATCAGCTGGGCGTGATGGGAACGGCAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGACTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTC  >  NC_000913/4238963‑4239220                                                                                                                             |                                                                                                                                      tATCAGCTGGGCGTGATGGGAACGGCAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACcgcg                                                                                                                        >  2:319631/1‑140 (MQ=255) tATCAGCTGGGCGTGATGGGAACGGCAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACcgcg                                                                                                                        <  1:814778/140‑1 (MQ=255) tATCAGCTGGGCGTGATGGGAACGGCAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACcgcg                                                                                                                        <  1:691618/140‑1 (MQ=255) tATCAGCTGGGCGTGATGGGAACGGCAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACcgcg                                                                                                                        >  2:13167/1‑140 (MQ=255)     aGCTGGGCGTGATGGGAACGGCAGATTTGTGGCTGATGGACCATCTGCTGACCACCGGCAGCCTGCTTGCAAATCTTGCAAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTTCCGGgcgtg                                                                                                                    <  2:229826/140‑1 (MQ=255)          ggCGTGTTGGGAACGGCAGATTTGTGGCTGACGGACCATCTGCTGCCCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTAc                                                                                                               >  1:229826/1‑140 (MQ=255)          ggCGTGATGGGAACGGCAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCgcgtgc                                                                                                                   <  2:8123/136‑1 (MQ=255)          ggCGTGATGGGAACGGCAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCgcgtgc                                                                                                                   >  1:8123/1‑136 (MQ=255)          ggCGTGATGGGAACGGCAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTAc                                                                                                               >  1:806618/1‑140 (MQ=255)          ggCGTGATGGGAACGGCAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTAc                                                                                                               >  1:851526/1‑140 (MQ=255)               gATGGGAACGGCAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGTGCGTGCGTACCCAtg                                                                                                          >  1:971545/1‑140 (MQ=255)               gATGGGAACGGCAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCAtg                                                                                                          >  2:63181/1‑140 (MQ=255)                   ggAACGGCAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTgcg                                                                                                      <  2:806618/140‑1 (MQ=255)                      aCGGCAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGa                                                                                                   <  1:758294/140‑1 (MQ=255)                        ggCAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGt                                                                                                 <  1:319631/140‑1 (MQ=255)                        ggCAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGt                                                                                                 <  2:495938/140‑1 (MQ=255)                        ggCAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGt                                                                                                 <  1:760500/140‑1 (MQ=255)                        gacAGATTTGTAGCTGACGAACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGTGCGTGCGTACCCATGTGCGCGAGt                                                                                                 <  2:971545/138‑1 (MQ=255)                          cAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTAt                                                                                               >  2:736794/1‑140 (MQ=255)                          cAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTAt                                                                                               <  2:382932/140‑1 (MQ=255)                          cAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTAt                                                                                               <  2:793404/140‑1 (MQ=255)                                  tgGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGaa                                                                                       <  2:851526/140‑1 (MQ=255)                                    gCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACg                                                                                     <  2:1006371/140‑1 (MQ=255)                                        aCGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGt                                                                                 >  1:178501/1‑140 (MQ=255)                                         cGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTc                                                                                <  2:258331/140‑1 (MQ=255)                                               aTCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTAtgtg                                                                          >  2:1018188/1‑140 (MQ=255)                                                  tgctgACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAAt                                                                       >  2:842650/1‑140 (MQ=255)                                                   gctgACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATa                                                                      <  2:178501/140‑1 (MQ=255)                                                      gACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAAcc                                                                   >  1:628978/1‑140 (MQ=255)                                                      gACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAAcc                                                                   <  1:1018188/140‑1 (MQ=255)                                                                           aaaTCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGc                                              <  1:736794/140‑1 (MQ=255)                                                                             aTCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCAt                                            <  1:13167/140‑1 (MQ=255)                                                                                    caacaaCTACGACAAGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCa                                     <  1:483615/140‑1 (MQ=255)                                                                                                  aGTTTAACTACACTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGt                       <  1:913623/140‑1 (MQ=255)                                                                                                            cacTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACg             >  2:235129/1‑140 (MQ=255)                                                                                                              cTAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACggt           >  2:659357/1‑140 (MQ=255)                                                                                                               tAATCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACggtg          <  1:659357/140‑1 (MQ=255)                                                                                                                 aTCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTTAGGTCTACggtggt        >  2:384948/1‑140 (MQ=255)                                                                                                                 aTCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGt                  <  2:215470/130‑1 (MQ=255)                                                                                                                 aTCCTCCGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGt                  >  1:215470/1‑130 (MQ=255)                                                                                                                  tcctccGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGtt       >  2:651254/1‑140 (MQ=255)                                                                                                                      ccGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAg                    <  2:544541/123‑1 (MQ=255)                                                                                                                      ccGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAg                    >  1:544541/1‑123 (MQ=255)                                                                                                                      ccGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTAtct   >  2:674795/1‑140 (MQ=255)                                                                                                                       cGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACggtggt        >  2:696348/1‑134 (MQ=255)                                                                                                                       cGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACggtggt        <  1:696348/134‑1 (MQ=255)                                                                                                                       cGCAGGGCTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTAtctc  >  2:1040481/1‑140 (MQ=255)                                                                                                                             |                                                                                                                                      TATCAGCTGGGCGTGATGGGAACGGCAGATTTGTGGCTGACGGACCATCTGCTGACCACCGGCAGCCTGTTTGCAAATCTTGCCAACAACTACGACAAGTTTAACTACACTAATCCTCCGCAGGACTCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTC  >  NC_000913/4238963‑4239220

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