Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,511,968	A→G	K105K (AAA→AAG)	ydcS →	putative ABC transporter periplasmic binding protein/polyhydroxybutyrate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,511,968	0	A	G	100.0%	124.6 / NA	35	K105K (AAA→AAG)	ydcS	putative ABC transporter periplasmic binding protein/polyhydroxybutyrate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/22);  total (13/22)

TTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACGGCACACCTTACCAATGGGGGCCGAACCTGCTG  >  NC_000913/1511831‑1512100                                                                                                                                          |                                                                                                                                     ttCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGcg                                                                                                                                    >  2:899289/1‑140 (MQ=255)  tCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAgcgc                                                                                                                                   >  2:75495/1‑140 (MQ=255)      aaaaGAGACGGGCGGCGAGGTGTATGTAAAAACCGGCGCTACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGc                                                                                                                               <  2:357353/140‑1 (MQ=255)      aaaaGAGACAGGCTGTGCGGTGCATCTGATAACCGCCCCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGc                                                                                                                               <  1:286148/140‑1 (MQ=255)        aaGAGAGAGGTTGCGCGGTGAATGTAAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAg                                                                                                                             <  2:237514/140‑1 (MQ=255)             agagggggcgcGGTAAATGGAAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGAt                                                                                                                        <  2:606046/133‑1 (MQ=255)             acgggCTGCGGGGTAAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGAt                                                                                                                        <  1:747599/137‑1 (MQ=255)             aCAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATACGGCGATGGCTCCCTGGGTTTGATTTTGGGAAAGCGCGTGCAGCCGAt                                                                                                                        >  2:318614/1‑140 (MQ=255)                ggCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGATTaa                                                                                                                     >  1:606046/1‑140 (MQ=255)                 gCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGc                                                                                                                               >  1:766456/1‑129 (MQ=255)                 gCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGc                                                                                                                               <  2:766456/129‑1 (MQ=255)                  cTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGATTAATa                                                                                                                   >  1:626211/1‑140 (MQ=255)                    gcgcGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGATTAATAcc                                                                                                                 <  1:685545/140‑1 (MQ=255)                     cgcgGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGc                                                                                                                               >  1:850518/1‑125 (MQ=255)                     cgcgGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGc                                                                                                                               <  2:850518/125‑1 (MQ=255)                               gtgAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGATTAATACCGCATTGATTcc                                                                                                      <  1:75495/140‑1 (MQ=255)                                        gccgcGACTTCCGATTAAATGGCCAATCTCATGACCACAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGaa                                                                                             <  1:318614/140‑1 (MQ=255)                                                        aaaTGGTCAGTCTGATGGCCAAGGGGGGTTACGATCTGGGTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCcgcgc                                                                             <  2:385352/140‑1 (MQ=255)                                                        aaaTGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCcgcgc                                                                             <  1:782941/140‑1 (MQ=255)                                                        aaaTGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCcgcgc                                                                             <  2:983092/140‑1 (MQ=255)                                                                      atgaCCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGa                                                               <  2:626211/140‑1 (MQ=255)                                                                       tgaCCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAACCGCACGATCCGCGCGTGGTTAAAGGTGAc                                                              >  1:496921/1‑140 (MQ=255)                                                                            aaaGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGtt                                                         >  2:717293/1‑140 (MQ=255)                                                                            aaaGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGATTAATACCGCATTGAGTCCCAACTGGAAGACGATCGATCCGGGCGTGGTTAAAGGCGACTGGtt                                                         >  2:274636/1‑140 (MQ=255)                                                                                   gTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGAc                                                              >  2:18582/1‑128 (MQ=255)                                                                                   gTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGAc                                                              <  1:18582/128‑1 (MQ=255)                                                                                               ttACGGCATCCGGCGCTGCCTCTCTCCGTATGATTATTGGTAAGCGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGtt                                      <  1:274636/140‑1 (MQ=255)                                                                                               ttACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGtt                                      <  1:881066/140‑1 (MQ=255)                                                                                                  cGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTAc                                   >  1:901916/1‑140 (MQ=255)                                                                                                    gCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACgg                                 <  1:659500/140‑1 (MQ=255)                                                                                                      aTCCGGCGGTGCCTCGTCGCGTTTGATTATGAGCAAGCGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACGGca                               <  2:493005/140‑1 (MQ=255)                                                                                                        ccGGCGATGCCTCGCTGCGTTTGATTATGGGTAAGCGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACGGcaca                             <  1:899289/140‑1 (MQ=255)                                                                                                                                 tATGGGTAAGCGCGTGCAGCCGGTTAATACTGCATGGATTCCCAACTGAAAAACGCTCGATCCGCGCGTGGTTAAAGGTGACTGGTTTAATGTTGGCGGCAAAGTTTACGGCACACCTTACCAATGGGGGCCGAACctgc    <  2:496921/140‑1 (MQ=255)                                                                                                                                   tGGGTAAGCGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACGGCACACCTTACCAATGGGGGCCGAACctgctg  >  2:178468/1‑140 (MQ=255)                                                                                                                                   tGGGTAAGCGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACGGCACACCTTACCAATGGGGGCCGAACctgctg  >  2:822894/1‑140 (MQ=255)                                                                                                                                   tGGGTAAGCGCGTGCAGCCGATTAACACCTCATTGATCCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACGGCACACCTTACCAATGGGGGCCGAACctgctg  <  1:284087/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                     TTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACGGCACACCTTACCAATGGGGGCCGAACCTGCTG  >  NC_000913/1511831‑1512100

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