Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,048,289	A→G	L135L (TTG→CTG)	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,289	0	A	G	100.0%	137.3 / NA	38	L135L (TTG→CTG)	gfcC	capsule biosynthesis GfcC family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/19);  total (19/19)

GCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCT  >  NC_000913/1048155‑1048425                                                                                                                                       |                                                                                                                                         gCGGGGGTGATCTTGCATATAACCGGTAACGTTACCGCCTGCCTGCCAAGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCgggg                                                                                                                                     <  2:225631/140‑1 (MQ=255) gCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCCCCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCgggg                                                                                                                                     <  1:510204/140‑1 (MQ=255) gCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCgggg                                                                                                                                     <  2:460759/140‑1 (MQ=255) gCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCgggg                                                                                                                                     <  2:136007/140‑1 (MQ=255)       gTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCGGGGGATTGc                                                                                                                               >  1:569532/1‑140 (MQ=255)          aTCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCGGGGGATTGCTGt                                                                                                                            <  1:514513/140‑1 (MQ=255)                     aaTCGGTAACGCTACGGCCTGCCTGCCAGGGGAATTGCCCGGCCCCAGCAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACt                                                                                                                 <  2:292568/140‑1 (MQ=255)                                cTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATc                                                                                                      <  2:569532/140‑1 (MQ=255)                                cTACGGCCCGCCTGCCAGGGCAATTGCCCGGCCCCCGAACCCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATc                                                                                                      <  2:342943/140‑1 (MQ=255)                                 tACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTc                                                                                                                <  1:438625/129‑1 (MQ=255)                                 tACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTc                                                                                                                >  2:438625/1‑129 (MQ=255)                                                    cAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTa                                                                                  <  2:156659/140‑1 (MQ=255)                                                          cccGGCCCCAGAAACCGCACCCTGCAGAGTAGTAGTCACTGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGAc                                                                            <  2:294765/140‑1 (MQ=255)                                                            cGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCa                                                                          <  1:81179/140‑1 (MQ=255)                                                              gCCCCAGAAACCGCAGCCAGCAGAGCAATTGTCCCCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGt                                                                        <  1:310630/140‑1 (MQ=255)                                                                         cGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATattgagat                                                             <  2:32233/140‑1 (MQ=255)                                                                                agcagAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTaa                                                      >  1:573182/1‑140 (MQ=255)                                                                                           aGTCACCGGACGCTGTAGGGTAGAGAGCGTACAGTCCCCTACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCa                                           <  2:303230/140‑1 (MQ=255)                                                                                               aCCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGAGCGGTAAGCAACCGGTGATATTGAGATTGAGTAACTGCTGGCGCACGGa                                       >  2:187517/1‑140 (MQ=255)                                                                                               aCCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGa                                       >  2:410215/1‑140 (MQ=255)                                                                                               aCCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCA‑GGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGAt                                      >  1:454129/1‑140 (MQ=255)                                                                                                ccGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGAt                                      <  1:410215/140‑1 (MQ=255)                                                                                                   gACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGCGATATTGGGATCGAGTAACTGATGGCGCACGGATTTa                                   >  1:303230/1‑140 (MQ=255)                                                                                                    aCGCTGTACGGTATAGAGCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGGACGAAATCAGGGCCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGTAGCAAGAATTTaa                                  >  2:451120/1‑140 (MQ=255)                                                                                                                  agagCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCcaca                    >  1:7783/1‑140 (MQ=255)                                                                                                                    agCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACAtc                  >  1:554679/1‑140 (MQ=255)                                                                                                                    agCGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACAtc                  >  1:174559/1‑140 (MQ=255)                                                                                                                      cGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACAtcgt                >  1:126621/1‑140 (MQ=255)                                                                                                                      cGTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACAtcgt                >  1:398843/1‑140 (MQ=255)                                                                                                                       gTATAGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACAtcgtc               >  1:402026/1‑140 (MQ=255)                                                                                                                        tataGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACAtcgtcg              >  1:517950/1‑140 (MQ=255)                                                                                                                         ataGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATcgtcgt             <  2:174559/140‑1 (MQ=255)                                                                                                                         ataGTCCCCGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATcgtcgt             >  1:129407/1‑140 (MQ=255)                                                                                                                              ccccGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCt        >  1:95408/1‑140 (MQ=255)                                                                                                                              ccccGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCt        >  1:47187/1‑140 (MQ=255)                                                                                                                                 cGACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAAGTGCTGGCGGACGAATTGAATCGTTGTCGCCGCATCGTCGTCAGCTTCt     >  2:196601/1‑140 (MQ=255)                                                                                                                                  gACCAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTg    <  2:517950/140‑1 (MQ=255)                                                                                                                                    ccAGCGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCt  <  2:668351/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                         GCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCT  >  NC_000913/1048155‑1048425

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