Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,048,196	C→T	A166T (GCA→ACA)	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,196	0	C	T	100.0%	132.5 / NA	39	A166T (GCA→ACA)	gfcC	capsule biosynthesis GfcC family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/21);  total (18/21)

CTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATC  >  NC_000913/1048069‑1048325                                                                                                                                |                                                                                                                                  cTTTTATTCCACAGGGCAACAGAGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCaa                                                                                                                       >  1:568520/1‑140 (MQ=255)      aTTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGcc                                                                                                                  <  2:260280/140‑1 (MQ=255)      aTTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGcc                                                                                                                  <  1:813983/140‑1 (MQ=255)      aTTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGcc                                                                                                                  <  1:517292/140‑1 (MQ=255)            caGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGccc                                                                                                            >  1:380203/1‑140 (MQ=255)            caGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGccc                                                                                                            >  2:71469/1‑140 (MQ=255)            caGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGccc                                                                                                            >  2:468529/1‑140 (MQ=255)               ggCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAg                                                                                                         >  2:855907/1‑140 (MQ=255)               ggCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAg                                                                                                         <  2:656742/140‑1 (MQ=255)                gCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGa                                                                                                        <  1:71469/140‑1 (MQ=255)                  aaCAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGaaa                                                                                                      >  2:679254/1‑140 (MQ=255)                  aaCAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGaaa                                                                                                      <  1:461907/140‑1 (MQ=255)                    cAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAAcc                                                                                                    >  1:497757/1‑140 (MQ=255)                       gggCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCa                                                                                                 <  2:555396/140‑1 (MQ=255)                       gggCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCa                                                                                                 <  2:380203/140‑1 (MQ=255)                         gCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCAcc                                                                                               >  2:60389/1‑140 (MQ=255)                          cgacgaCAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCAccc                                                                                              <  1:381261/140‑1 (MQ=255)                            acgacAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCcag                                                                                            <  2:385338/140‑1 (MQ=255)                                             aGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCg                                                                           >  2:461967/1‑140 (MQ=255)                                             aGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCg                                                                           >  2:528075/1‑140 (MQ=255)                                             aGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCg                                                                           >  2:366950/1‑140 (MQ=255)                                             aGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCg                                                                           >  1:603383/1‑140 (MQ=255)                                             aGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATATTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCg                                                                           >  2:836837/1‑140 (MQ=255)                                               gCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGa                                                                         <  2:603383/140‑1 (MQ=255)                                                 gTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACg                                                                       <  2:327480/140‑1 (MQ=255)                                                     tCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGt                                                                   <  1:321255/140‑1 (MQ=255)                                                      cACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTa                                                                  <  2:497757/140‑1 (MQ=255)                                                                ttATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGtat                                                            >  1:67687/1‑136 (MQ=255)                                                                ttATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGtat                                                            <  2:67687/136‑1 (MQ=255)                                                                        gTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGt                                                <  1:463074/140‑1 (MQ=255)                                                                               ccTGCAAGGCGGGGGTGATATTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGAc                                         <  1:836837/140‑1 (MQ=255)                                                                                        cGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACggggg                                >  1:385383/1‑140 (MQ=255)                                                                                         gggggTGATCTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGa                               <  1:679254/140‑1 (MQ=255)                                                                                                  cTTGCAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTtc                      <  1:461967/140‑1 (MQ=255)                                                                                                      cAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACgggggg                               <  1:692842/127‑2 (MQ=255)                                                                                                      cAGATAATCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACgggggg                               >  2:692842/1‑126 (MQ=255)                                                                                                             tCGGTAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCg           >  2:510434/1‑140 (MQ=255)                                                                                                                 tAACGCTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACa       <  1:510434/140‑1 (MQ=255)                                                                                                                      cTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATc  >  2:628458/1‑140 (MQ=255)                                                                                                                      cTACGGCCTGTCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATc  <  1:60389/140‑1 (MQ=255)                                                                                                                                |                                                                                                                                  CTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATC  >  NC_000913/1048069‑1048325

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