Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,428,122	A→G	L199P (CTG→CCG)	ytfB ←	OapA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,428,122	0	A	G	100.0%	141.8 / NA	40	L199P (CTG→CCG)	ytfB	OapA family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/20);  total (20/20)

CCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGTTTTGCAAATTACTCAGCGGCTTACCCGCGCCTTCTACCTGCGCCATTGCATAGACATCGGT  >  NC_000913/4427978‑4428267                                                                                                                                                 |                                                                                                                                                  ccTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCt                                                                                                                                             <  2:62037/151‑1 (MQ=255)    ggTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACTGCGCACGAATAAAACTGCCATCCGGCTGGCGGGGGCACGGGACCTGtt                                                                                                                                          >  2:246105/1‑151 (MQ=255)            acaTCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCGTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTc                                                                                                                                  >  1:122798/1‑151 (MQ=255)            acaTCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTc                                                                                                                                  <  1:351749/151‑1 (MQ=255)                 aGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGt                                                                                                                             >  1:399317/1‑151 (MQ=255)                    aTGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTAt                                                                                                                           <  1:69217/150‑1 (MQ=255)                       gTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAAt                                                                                                                       <  2:435526/151‑1 (MQ=255)                        tttCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGCGAACGGGACCTGTTGATTATTCCCCGTATCAAt                                                                                                                       >  1:421999/1‑150 (MQ=255)                           cACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGt                                                                                                                    <  1:130419/150‑1 (MQ=255)                                   gCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGt                                                                                                           <  1:206854/151‑1 (MQ=255)                                     aCTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTca                                                                                                         <  2:381937/151‑1 (MQ=255)                                        ggAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTcacc                                                                                                       >  1:6347/1‑150 (MQ=255)                                                     aaGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGtt                                                                                         >  1:255773/1‑151 (MQ=255)                                                     aaGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGt                                                                                          <  1:612677/150‑1 (MQ=255)                                                      aGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGtt                                                                                         >  2:539928/1‑150 (MQ=255)                                                         tcatcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGGGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTgg                                                                                     >  2:321628/1‑151 (MQ=255)                                                            tcaAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCgg                                                                                  >  2:284881/1‑151 (MQ=255)                                                               aaGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGAtt                                                                               >  1:74983/1‑151 (MQ=255)                                                                 gagaCAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGAtttt                                                                             <  1:507674/151‑1 (MQ=255)                                                                 gagaCAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGAtttt                                                                             <  2:387299/151‑1 (MQ=255)                                                                  agaCAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGAtttt                                                                             <  1:45018/150‑1 (MQ=255)                                                                   gaCAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGAtttt                                                                             >  2:29287/1‑149 (MQ=255)                                                                   gaCAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCa                                                                           <  1:284888/151‑1 (MQ=255)                                                                    aCAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCAc                                                                          <  1:271470/151‑1 (MQ=255)                                                                         ccTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCAACACGCCGCTGGCGTTCTGGAGGATTTTCACCATCt                                                                     >  2:629151/1‑151 (MQ=255)                                                                             aGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGccc                                                                 <  2:126669/151‑1 (MQ=255)                                                                                 ggCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGtt                                                              <  1:544225/150‑1 (MQ=255)                                                                                   cTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGtttt                                                            <  2:100122/150‑1 (MQ=255)                                                                                                               aCGAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGTTTTGCAAATTACTCAGCGGCTTACCCGCGcc                                >  1:165857/1‑150 (MQ=255)                                                                                                                 gAATAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGTTTTGCAAATTACTCAGCGGCTTACCCGCGCCTTc                             >  2:145754/1‑151 (MQ=255)                                                                                                                  aaTAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGTTTTGCAAATTACTCAGCGGCTTACCCGCGCCTTCt                            <  2:501371/151‑1 (MQ=255)                                                                                                                   aTAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGTTTTGCAAATTACTCAGCGGCTTACCCGCGCCTTCt                            >  1:51318/1‑150 (MQ=255)                                                                                                                    tAAAACTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCAACACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGTTTTGCAAATTACTCAGCGGCTTACCCGCGCCTTCTa                           >  2:508655/1‑150 (MQ=255)                                                                                                                     aaaaCTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGTTTTGCAAATTACTCAGCGGCTTACCCGCGCCTTCTAcc                         >  2:598351/1‑151 (MQ=255)                                                                                                                        aCTGCCATCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGTTTTGCAAATTACTCAGCGGCTTACCCGCGCCTTCTACCTgc                      <  1:145757/151‑1 (MQ=255)                                                                                                                              aTCCGGCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGTTTTGCCAATTACTCAGCGGCTTACCCGCGCCTTCTACCTGCGCCAtt                >  2:386240/1‑151 (MQ=255)                                                                                                                                  ggCTGGCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGTTTTGCAAATTACTCAGCGGCTTACCCGCGCCTTCTACCTGCGCCATTGCAt            >  2:390262/1‑151 (MQ=255)                                                                                                                                      ggCGGGTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGTTTTGCAAATTACTCAGCGGCTTACCCGCGCCTTCTACCTGCGCCATTGCATAGa         <  1:545555/150‑1 (MQ=255)                                                                                                                                          ggTGAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGTTTTGCAAATTACTCAGCGGCTTACCCGCGCCTTCTACCTGCGCCATTGCATAGACATCg    <  2:441601/151‑1 (MQ=255)                                                                                                                                             gAACGGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGTTTTGCAAATTACTCAGCGGCTTACCCGCGCCTTCTACCTGCGCCATTGCATAGACATCGGt  >  2:167589/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                  CCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGCACGAATAAAACTGCCATCCGGCTGGCGGGTGAACAGGACCTGTTGATTATTCCCCGTATCAATGGTTAAGCCCGTCACCACGCCGCTGGCGTTCTGGCGGATTTTCACCATCTGCCCGTTTTGCAAATTACTCAGCGGCTTACCCGCGCCTTCTACCTGCGCCATTGCATAGACATCGGT  >  NC_000913/4427978‑4428267

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