Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,103,891	G→A	P207S (CCA→TCA)	wbbK ←	putative glycosyltransferase WbbK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,103,891	0	G	A	100.0%	73.8 / NA	23	P207S (CCA→TCA)	wbbK	putative glycosyltransferase WbbK
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/14);  total (9/14)

GCAATCCCCATGTTTCTAACCTAGAGGGAAAACAAACTATATCTGAAATATTATAACAATGATCGATTTTTTCTTTATCCAAGTACCCGAGGAAATGAACATTATCCAGTCCTTCTGCAAGACTGATAATATATTTTGCATACGCATTTTCTGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTGATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCATT  >  NC_000913/2103645‑2104111                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                             gCAATCCCCATGTTTCTAACCTAGAGGGAAAACAAACTATATCTGAAATATTATAACAATGATCGATTTTTTCTTTATCCAAGTACCCGAGGAAATGAACATTATCCAGTCCTTCTGCAAGACTGATAATATATTTTGCATACGCATTTTCTGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGAAAca                                                                                                                                                                                                                          <  1:421704/251‑1 (MQ=255)  cAATCCCCATGTTTCTAACCTAGAGGGAAAACAAACTATATCCGAAATATTATAACAATGATCGATTTTTTCTTTATCCAAGTACCCGAGGAAATGAACATTATCCAGTCCTTCTGCAAGACTGATAATATATTTTGCATACGCATTTTCTGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGAAAc                                                                                                                                                                                                                           >  1:61583/1‑249 (MQ=255)        ccATGTTTCTAACCTAGAGGGACAACACACTATATCTGAAATATTATAACCCTGATCGATTTTTTCTTTATCCACGTACCCGAGGAAATGACCATTATCCAGTCCTTCTGCAAGACTGATAATATATTTTGCATACGCATTTTCTGTCCCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGAAAcagcagtgt                                                                                                                                                                                                                   >  2:237407/1‑248 (MQ=255)                    ccTAGAGGGAAAACAAACTATATCTGAAATATTATAACAATGATCGATTTTTTCTTTATCCAAGTACCCGAGGAAATGAACATTATCCAGTCCTTCTGCAAGACTGATAATATATTTTGCATACGCATTTTCTGAACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGAAACAGCAGGGTAAAATATTGTCa                                                                                                                                                                                                       <  2:100603/251‑1 (MQ=255)                                 caaacTATATCTGAAATATTATAACAATGTTCGATTTTTTCTTTATCCAAGTACCCGAGGAAATGAACATTATCAAGTCCTTCTGCAAGACTGATAATATATTTTGCATACGCATTTTCTGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGAAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGtt                                                                                                                                                                                          >  2:514631/1‑251 (MQ=255)                                             gAAATATTATAACAATGATCGATTTTTTCTTTATCCAAGTACCCGAGGAAATGAACATTATCCAGTCCTTCTGCAAGACTGATAATATATTTTGCATACGCATTTTCTGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGAAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAAtt                                                                                                                                                                               <  2:70216/250‑1 (MQ=255)                                                              aTCGATTTTTTCTTTATCCAAGTACCCGAGGAAATGAACATTATCCAGTCCTTCTGCAAGACTGATAATATATTTTGCATACGCATTTTCTGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGAAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAAGTGAGAATCATCATCAGTaa                                                                                                                                                             <  1:203731/251‑1 (MQ=255)                                                                             aTCCAAGTACCCGAGGAAATGAACATTATCCAGTCCGTCTGCAAGACTGATAATATATTTTGCATACGCATTTTCTGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGAAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAg                                                                                                                                              <  1:125866/251‑1 (MQ=255)                                                                                cAAGTACCCGAGGAAATGAACATTATCCAGTCCTTCTGCAAGACTGATAATATATTTTGCATACGCATTTTCTGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGAAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGAt                                                                                                                                            >  1:257391/1‑250 (MQ=255)                                                                                                   aCATTATCCAGTCCTTCTGCAAGACTGATAATATATTTTGCATACGCATTTTCTGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGAAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACt                                                                                                                        >  1:569173/1‑251 (MQ=255)                                                                                                        aTCCAGTCCTTCTGCAAGACTGATAATATATTTTGCATACGCATTTTCTGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGAAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACa                                                                                                                    <  1:402502/250‑1 (MQ=255)                                                                                                             gTCCTTCTGCAAGACTGATAATATATTTTGCATACGCATTTTCTGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAAGACTCGTGAAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAatgatg                                                                                                              <  2:294750/251‑1 (MQ=255)                                                                                                               ccTTCTGCAAGACTGATAATATATTTTGCATACGCATTTTCTGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGAAACATCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTGGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGtt                                                                                                            <  2:535957/251‑1 (MQ=255)                                                                                                                             gATAATATATTTTGCATACGCATTTTCTGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGAAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATAtt                                                                                              <  2:61583/251‑1 (MQ=255)                                                                                                                              ataataTATTTTGCATACGCATTTTCTGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGAAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATAttt                                                                                             >  1:368829/1‑251 (MQ=255)                                                                                                                                                        tGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGAAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTGATAAATTTTTCTTTCATcc                                                                      >  1:34648/1‑248 (MQ=255)                                                                                                                                                                        aGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGAAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTGATAAATTTTTCTTTCATCCAGAATTGTTGAACAAacac                                                   <  1:502817/251‑1 (MQ=255)                                                                                                                                                                             aaaTTTAATATTGGATTGTTCTTTCAATTTCCTGGCTGCACTAATAATAAGCTCGTAAGTTTTAAATACTCGTGAAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGGAAGTTGGCTTTTATCAGATAATTTAATTTCTGTCCGACTGACAATGATGTTATTTATAGAATATTTCCTGATAAATTTTTCGTTCGTCCAGAATTGTTGATCGAACGCTGCAg                                              <  2:395130/251‑1 (MQ=255)                                                                                                                                                                                                           ccTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGAAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTGATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCa                <  1:391509/251‑1 (MQ=255)                                                                                                                                                                                                              tgctgcACTAATAATAAGCTCGTAATTTTTAAATACTCGTGAAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTAGTTATAGAATATTTCTTGATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAg               <  2:26314/249‑1 (MQ=255)                                                                                                                                                                                                                  gcACTAATAATAAGCTCGTAATTTTTAAATACTCGTGAAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTGATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGt          >  1:186526/1‑250 (MQ=255)                                                                                                                                                                                                                        ataataaGCTCGTAATTTTTAAATACTCGTGAAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATAAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTGATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCAt   >  1:566375/1‑251 (MQ=255)                                                                                                                                                                                                                          aataaGCTCGTAATTTTTAAATACTCGTGAAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTGATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCAtt  <  2:307614/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                             GCAATCCCCATGTTTCTAACCTAGAGGGAAAACAAACTATATCTGAAATATTATAACAATGATCGATTTTTTCTTTATCCAAGTACCCGAGGAAATGAACATTATCCAGTCCTTCTGCAAGACTGATAATATATTTTGCATACGCATTTTCTGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTGATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCATT  >  NC_000913/2103645‑2104111

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