Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,511,751	A→G	E33G (GAA→GGA)	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,751	0	A	G	100.0%	110.2 / NA	31	E33G (GAA→GGA)	ydcS	putative ABC transporter periplasmic binding protein/polyhydroxybutyrate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/16);  total (15/16)

CAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTT  >  NC_000913/1511626‑1511879                                                                                                                              |                                                                                                                                 cAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGa                                                                                                                           <  2:129697/133‑1 (MQ=255)     taattaCGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTgg                                                                                                                       >  2:472525/1‑133 (MQ=255)                   gTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTgg                                                                                                         <  2:307461/133‑1 (MQ=255)                   gTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTgg                                                                                                         <  1:372450/133‑1 (MQ=255)                         ccTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGa                                                                                                   <  1:337218/133‑1 (MQ=255)                            tATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATAc                                                                                                >  2:374071/1‑133 (MQ=255)                            tATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATGATAGCCTGGACGGGATAc                                                                                                >  1:295818/1‑133 (MQ=255)                              tGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACAt                                                                                              >  2:440652/1‑133 (MQ=255)                                  cAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGaa                                                                                          <  2:489860/133‑1 (MQ=255)                                            gCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACaaa                                                                                <  1:440652/133‑1 (MQ=255)                                            gCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACaaa                                                                                <  1:607264/133‑1 (MQ=255)                                                  agcagcCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATa                                                                          >  2:19111/1‑133 (MQ=255)                                                  agcagcCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAAATAAACCGGGAGGGCGACTGGATATTATCGCCTGGACGGGATAAATCCACCGCGGACAAACTGATa                                                                          >  1:219888/1‑133 (MQ=255)                                                         tgtgtgCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATa                                                                   <  2:324765/133‑1 (MQ=255)                                                            gtgCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGa                                                                <  1:19111/133‑1 (MQ=255)                                                            gtgCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGa                                                                <  2:554854/133‑1 (MQ=255)                                                                gcTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGCGATACATCTAACGCGGGCAAACTGATAAACAATACGACTgg                                                            >  2:263640/1‑133 (MQ=255)                                                                 cTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAAGGCGGACAGACTGATAAACAATACGACTggg                                                           >  1:216453/1‑133 (MQ=255)                                                                      cATGACAATAGTGACTGCTCACGCGGCCTAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACg                                                      <  2:216453/133‑1 (MQ=255)                                                                                atGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCa                                                    >  1:508294/1‑125 (MQ=255)                                                                                atGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCa                                                    <  2:508294/125‑1 (MQ=255)                                                                                atGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTGGaaa                                            <  2:600085/133‑1 (MQ=255)                                                                                atGACCGCGCACGCCGCCTAACTGCCTACCTATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGaaa                                            <  1:262623/133‑1 (MQ=255)                                                                                     cGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAagag                                       >  2:387957/1‑133 (MQ=255)                                                                                     cGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCCAAAAagag                                       >  1:245451/1‑133 (MQ=255)                                                                                        tCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACa                                    >  2:143007/1‑133 (MQ=255)                                                                                               gccgAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTgcg                             <  1:374071/133‑1 (MQ=255)                                                                                               gccgAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTgcg                             <  1:143007/133‑1 (MQ=255)                                                                                                  gAACCGCCTACCAATTTAGATAAACCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCgg                          >  2:505358/1‑133 (MQ=255)                                                                                                                         aCCGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACt   >  1:502358/1‑133 (MQ=255)                                                                                                                          ccGGGAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACtt  >  1:109036/1‑133 (MQ=255)                                                                                                                              |                                                                                                                                 CAATTAATTACGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTT  >  NC_000913/1511626‑1511879

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: ‑