Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,048,596	G→A	G32G (GGC→GGT)	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,048,596	0	G	A	100.0%	73.2 / NA	26	G32G (GGC→GGT)	gfcC	capsule biosynthesis GfcC family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/15);  total (11/15)

TGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGTTTTTCTGGCCGCG  >  NC_000913/1048355‑1048818                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                               tGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATAGTTGCCGCCACATCGTCGTCAGCTTCTGCTCCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCACAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGCCAACATTTTCCACGGGGCCAACGGACAGTGTCTGTTGTTCACCAGGc                                                                                                                                                                                                                          >  1:183548/1‑248 (MQ=255)                     tgctgcCGCACGGATTTAATCGTTGCCGCCACATCGTCGTTCGCTTCTGCTTCACAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTCAGCGCTTTAGCTTTCGCCGCAGAGTAGGTGAGCAAGGCCCCAGGCCACCAGAGGCGCTCGCGCAGTTGCGGTTGTGTCACCAAATGGACAACATTTTCCACGGGGCCAACGGAAAGCATCTGTTGTTAAACAGGCAGCTCAATAGTCCCCGTTTCTTg                                                                                                                                                                                                   >  2:144928/1‑251 (MQ=255)                      gctgGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACCTGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGCGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACCACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTgc                                                                                                                                                                                                  >  1:134213/1‑251 (MQ=255)                           cgcACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCaa                                                                                                                                                                                              >  2:150556/1‑250 (MQ=255)                            gcACGGATTTAAGCGTTGCCGCCACATCGTCGTCAGCTTCTGCTGCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGGTCACCAGGCAGGTAAATAGTCACCGTTCCTTGCGCAAAg                                                                                                                                                                                            <  1:264417/251‑1 (MQ=255)                                   tttAATCGTTGCCCCGACATCGTCGACAGCTTCTGCTACCCAGGAAGCACGCCGTGCCATGACATGTTGGTAGCCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGAGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGCTCACCAGGCAGATAAATAGTCACCGTTCCTTTCGCAAAGGCATggg                                                                                                                                                                                     >  2:176478/1‑251 (MQ=255)                                     tAATCGTTGCCGCCACAGCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATgggg                                                                                                                                                                                    <  1:30123/250‑1 (MQ=255)                                          gtggccgccACATCGTCGTCAGCTGCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTGGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCTACGGAAAGTGTCTGTGGTTCACCAGGCAGATAAATAGTCACCTTTCCGTGCGCAAAGGCATGCGGTGTCAt                                                                                                                                                                              <  2:35419/248‑1 (MQ=255)                                           ttGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAACGTGTCTGTTGTTCACCAGGCAGATAAATAGACACCGTTCCTTGCGCAAAGGCATGGGGTGTCAtt                                                                                                                                                                             >  2:172093/1‑251 (MQ=255)                                             gccgccGCATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTGAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCGGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCGTGCGCGAAGGCATGGTGTGTCATTAc                                                                                                                                                                           <  2:132643/251‑1 (MQ=255)                                                                                        gtgCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATg                                                                                                                                <  1:144928/251‑1 (MQ=255)                                                                                             aTGACATGGTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAGGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCTTGGGGTGTCATTAAGTAAAGTACGCTGGCAATGAAATACGCATGTAATTTATTCATGGTGCg                                                                                                                           >  2:474/1‑251 (MQ=255)                                                                                                    gTTGGTAGGCTTTAAGCTCTTTAGCTTTCGCCTCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTACACCAGGCAGATAAATAGTCTCCGTTCCGTGCGCAAAGGCATGGGGTGTCAGTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTAGTCAGGGTGCGGGTTTa                                                                                                                     <  2:183548/250‑1 (MQ=255)                                                                                                               ttAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCa                                                                                                           <  1:84913/249‑1 (MQ=255)                                                                                                               ttAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCa                                                                                                           <  1:70871/249‑1 (MQ=255)                                                                                                                                cgccgcAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAAccc                                                                                         <  1:113972/250‑1 (MQ=255)                                                                                                                                   cgcAGAGTCGGTCAGCAAGGCACCAGGCAACCAGAGGCGATCGCGCAGATGCGGTTGTGTCACCACCTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGCTAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTCCGGGAATAACCCCTgc                                                                                     >  2:292561/1‑251 (MQ=255)                                                                                                                                      agagTCGGGCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGcc                                                                                  <  1:187761/251‑1 (MQ=255)                                                                                                                                          tCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGATGCGGTTGTGTCACCAACTCGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGCTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAAc                                                                              >  1:290614/1‑251 (MQ=255)                                                                                                                                           cGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACCTTTTCCCAGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTACGAAACGTCATTTCTACGGGAATAACCCCTGCGCCtcac                                                                              >  1:194664/1‑247 (MQ=255)                                                                                                                                                                            cgcgCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTGGCCGACTTTGACGAACCTGGCCAGTAGcc                                             <  2:302537/250‑1 (MQ=255)                                                                                                                                                                                            ggtCACCAACTGGACACCATTTTCCACGGGGCCAACGGAAAGTGGCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCGTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAAtgt                            <  2:171607/250‑1 (MQ=255)                                                                                                                                                                                                    aCTGGACAACATTTTCCACGGGGCCAACGGAAAGTGGCTGTTGTTCACCAGGCAGATAAAGATGCACAGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGGACGCTGGCAATTAAATACGACTGTAATTTATTCATGGTGCGGGTGTAAGAAACTTCAGTTGTACGGGAGTAACCCGTGCTCCTATCAGTTGCCGACGTTGATGAACCTGGCGAGTAGCCGAATCGACCCAGGATGTGTTTTGcc                    <  2:197999/251‑1 (MQ=255)                                                                                                                                                                                                                tttCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttg         <  1:292561‑M1/250‑8 (MQ=255)                                                                                                                                                                                                                  tCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgac       <  1:150556‑M1/250‑10 (MQ=255)                                                                                                                                                                                                                       ggggCCAACGGACAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAACGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagt  >  2:109961‑M1/1‑236 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                               TGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGTTTTTCTGGCCGCG  >  NC_000913/1048355‑1048818

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