Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,047,955	A→G	V246A (GTT→GCT)	gfcC ←	capsule biosynthesis GfcC family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,047,955	0	A	G	100.0%	137.1 / NA	41	V246A (GTT→GCT)	gfcC	capsule biosynthesis GfcC family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/16);  total (25/16)

GTGCCGCCAAAGTCTGATTGCGACGAACCCAGCGGATCCGGGTAGGTTAATACTTCAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAACACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAG  >  NC_000913/1047827‑1048082                                                                                                                                 |                                                                                                                                gtgCCGCCAAAGTCTGATTGCGACGAACCCAGCGGATCCGGGTAGGTTAATACTTCAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAGGATAAGAATTCTTCTTCATCAGTCAGGAGCACg                                                                                                                             >  1:31085/1‑133 (MQ=255) gtgCCGCCAAAGTCTGATTGCGACGAACCCAGCGGATCCGGGTAGGTTAATACTTCAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACg                                                                                                                             >  2:512204/1‑133 (MQ=255) gtgCCGCCAAAGTCTGATTGCGACGAACCCAGCGGATCCGGGTAGGTTAATACTTCAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACg                                                                                                                             >  2:480735/1‑133 (MQ=255)  tgCCGCCAAAGTCTGATTGCGACGAACCCAGCGGATCCGGGTAGGTTAATACTTCAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGt                                                                                                                            <  2:47057/133‑1 (MQ=255)  tgCCGCCAAAGTCTGATTGCGACGAACCCAGCGGATCCGGGTAGGTTAATACTTCAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGt                                                                                                                            >  1:47057/1‑133 (MQ=255)    ccgccAAAGTCTGATTGCGACGAACCCAGCGGATCCGGGTAGGTTAATACTTCAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTg                                                                                                                          >  1:127349/1‑133 (MQ=255)    ccgccAAAGTCTGATTGCGACGAACCCAGCGGATCCGGGTAGGTTAATACTTCAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTg                                                                                                                          >  2:297692/1‑133 (MQ=255)       ccAAAGTCTGATTGCGACGAACCCAGCGGATCCGGGTAGGTTAATACTTCAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAGAAGATAAGGATTGTACGGCATCAGTAAGGAGAACGTTGCGt                                                                                                                       >  2:250548/1‑133 (MQ=255)            gTCTGATTGCGACGAACCCAGCGGATCCGGGTAGGTTAATACTTCAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTaaaa                                                                                                                  <  1:411058/133‑1 (MQ=255)             tCTGATTGCGACGAACCCAGCGGATCCGGGTAGGTTAATACTTCAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAAc                                                                                                                 <  1:512204/133‑1 (MQ=255)                       acgaACCCAGCGGATCCGGGTAGGTTAATACTTCAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATc                                                                                                       >  1:37176/1‑133 (MQ=255)                          aaCCCAGCGGATCCGGGTAGGTTAATACTTCAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTGCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTgg                                                                                                    <  2:387113/133‑1 (MQ=255)                                     tCCGGGTAGGTTAATACTTCAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATc                                                                                         >  1:441702/1‑133 (MQ=255)                                               ttAATACTTCAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATAtttt                                                                               <  1:297633/133‑1 (MQ=255)                                                   tACTTCAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCgg                                                                           >  2:104360/1‑133 (MQ=255)                                                        cAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACATACACAATCTGGTGATTGAAATCAGCATATTCTTCGGGCAGt                                                                      >  2:133954/1‑133 (MQ=255)                                                        cAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGc                                                                         <  1:206121/130‑1 (MQ=255)                                                        cAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGc                                                                         <  2:30800/130‑1 (MQ=255)                                                        cAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGc                                                                         >  1:30800/1‑130 (MQ=255)                                                        cAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGc                                                                         >  2:206121/1‑130 (MQ=255)                                                        cAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGt                                                                      >  2:468994/1‑133 (MQ=255)                                                          gCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTAc                                                                    <  2:127349/133‑1 (MQ=255)                                                             tGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATg                                                                 <  1:297692/133‑1 (MQ=255)                                                                caagCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATATGGTCATTAAAAGCAGAATATTTTTGGGGCAGTACATgcgc                                                              >  2:145237/4‑133 (MQ=255)                                                                     aCTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTGATTTAAATGAGCGTATGTTTCGGGCGGTACTTGCGCAGaaa                                                         >  2:141751/1‑133 (MQ=255)                                                                              ggCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCTGTACATGCTCAGAAAAGCccaacc                                                >  1:377739/1‑133 (MQ=255)                                                                              ggCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCccaacc                                                >  1:324179/1‑133 (MQ=255)                                                                              ggCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGACATTTAAATCAGCATATTTTTCGGGCAGTACAGGCTAAGAAAAGCccaacc                                                >  1:253643/1‑133 (MQ=255)                                                                                          gCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGc                                    >  2:496082/1‑133 (MQ=255)                                                                                                tAAAAGATCAGAATTGTTTTTCGTCAGTCAAGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCg                              <  1:250548/133‑1 (MQ=255)                                                                                                 aaaaGATATGAATTCCTCCTCACCAGTCAGCAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCgg                             <  1:141751/133‑1 (MQ=255)                                                                                                    aGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCggtgg                          >  1:367448/1‑133 (MQ=255)                                                                                                    aGATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCggtgg                          <  2:37176/133‑1 (MQ=255)                                                                                                     gATAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCggtggt                         <  2:31085/133‑1 (MQ=255)                                                                                                       tAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTc                       <  1:480735/133‑1 (MQ=255)                                                                                                       tAAGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTc                       <  1:104360/133‑1 (MQ=255)                                                                                                        aaGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTcc                      >  2:35763/1‑133 (MQ=255)                                                                                                        aaGAATTCTTCTTCATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTcc                      >  2:107777/1‑133 (MQ=255)                                                                                                                ttcttcATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGGCAGGCGCTGGTTCAACATGTCt              >  1:169561/1‑133 (MQ=255)                                                                                                                 tcttcATCAGTCAGGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCtt             <  2:441702/133‑1 (MQ=255)                                                                                                                            cagGAGCACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTGCCAGATGTCTTTTATTCCACAg  >  1:178285/1‑133 (MQ=255)                                                                                                                                 |                                                                                                                                GTGCCGCCAAAGTCTGATTGCGACGAACCCAGCGGATCCGGGTAGGTTAATACTTCAGCATGACAGGCACTGGAGACGGCAATGGCCAGGCAGCTTAAAAGATAAGAATTCTTCTTCATCAGTCAGGAACACGTTGCGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAG  >  NC_000913/1047827‑1048082

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