Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,110,033	G→A	Q233Q (CAG→CAA)	opgG →	osmoregulated periplasmic glucans biosynthesis protein G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,110,033	0	G	A	91.3%	50.3 / ‑1.4	23	Q233Q (CAG→CAA)	opgG	osmoregulated periplasmic glucans biosynthesis protein G
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (12/9);  total (13/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

TGCTCGGGGCCAGCTATTTCCGCGTGATTGGTGCAGGTCAGGTTTATGGCCTTTCTGCCCGCGGCCTGGCAATTGATACCGCCTTGCCATCGGGTGAAGAATTTCCGCGCTTCAAAGAGTTCTGGATCGAGCGTCCAAAACCGACTGATAAACGTTTAACCATCTATGCATTGCTTGACTCGCCGCGTGCGACAGGTGCTTACAAATTCGTGGTTATGCCAGGGCGTGACACGGTTGTGGATGTGCAGTCGAAAATCTATCTGCGCGATAAAGTCGGCAAACTGGGGGTTGCACCGTTAACCAGTATGTTCCTGTTTGGGCCGAACCAACCGTCGCCTGCAAATAACTATCGTCCGGAGTTGCACGACTCTAACGGTCTCTCTATCCATGCCGGTAATGGCGAATGGATCTGGCGTCCGTTGAATAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGAAAACCCGCAAGGCTTTGGTCTG  >  NC_000913/1109786‑1110270                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              tGCTCGGGGCCAGCTATTTCCGCGTGATTGGTGCAGGTCAGGTTTATGGCCTTTCTGCCCGCGGCCTGGCAATTGATACCGCCTTGCCATCGGGTGAAGAATTTCCGCGCTTCAAAGAGTTCTGGATCGAGCGTCCAAAACCGACTGATAAACGTTTAACCCTCTATGCATTGCTTGACTCGCCGCGGGCGACAGGTGATTACAAATTCGTGGTTATGCCAGGGCGTGACCCGGTTGTGGATGTGCAATCg                                                                                                                                                                                                                                            >  1:276945/1‑251 (MQ=255)                                  gagtgagggtTTATGGCCTTTATGCCCGCGGCCTGGCACTTGATACCGCCTTGCCATCGGGTGAAGAATTTCCGCGCTTCAAAGAGTTCTGGATCGAGCGTCCAAAACCGACTGATAAACGTTTAACCATCTATGCATTGCTTGACTCGCCGCGTGCGACAGGTGCTTACAAATTCGTGGTTATGCCAGGGCGTGACACGGTTGTGGATGTGCAATCGAAAATCTATCTGCGCGATAAAGTCGGCAAACTg                                                                                                                                                                                                           <  2:88540/244‑1 (MQ=255)                                       caggtTTATGGCCTTTCTGCCCGCGTCCTGGCAATTGATACCGCCTTGCCATCGGGTGAAGAATTTCCGCGCTTCAAAGAGTTCTGTATCGAGCGTCCAAAACCGACTGATAAACGTTTAACCATCTATGCATTGCTTGACTCGCCGCGTGCGACAGTTGCTTACAAATTCGTGGTTATGCCAGGGCGTGACACGGTTGTGGATGTGCAATCGAAAATCTATCTGCGCGATAAAGTCGGCAATCTggggg                                                                                                                                                                                                       <  1:241626/250‑1 (MQ=255)                                             tATGGCCTTTCTGCCCGCGGCCTGGCAATTGATACCGCCTTGCCATCGGGTGAAGAATTTCCGCGCTTCAAAGAGTTCTGGATCGAGCGTCCAAAACCGACTGATAAACGTTTAACCATCTATGCATTGCTTGACGCGCCGCGTGCGACAGGTGCTTACAAATTCGTGGTTATGCCAGGGCGTGACACGGTTGTGGATGTGCAATCGAAAATCTATCTGCGCGATAAAGTCGGCAAACTGGGGGTTGCAc                                                                                                                                                                                                 <  1:352956/250‑1 (MQ=255)                                               tGGCCTTTCTGCCCGCGGCCTGGCAATTGATACCGCCTTGCCATCGGGTGAAGAATTTCCGCGCTTCAAAGAGTTCTGGATCGAGCGTCCAAAACCGACTGATAAACGTTTAAACATCTATGCATTGCTTGACTCGCCGCGTGCGACAGGTGCTTACAAATTCGTGGTTATGCCAGGGCGTGACACGGTTGTGGATGTGCAATCGAAAATCTATCTGCGCGATAAAGACGGCAAACTGGGGGTTGCAACGt                                                                                                                                                                                              >  2:226489/1‑251 (MQ=255)                                                         gCCCGCGGCCTGGCAATTGATACCGCCTTGCCATCGGGTGAAGAATTTCCGCGCTTCAAAGAGTTCTGGATCGAGCGTCCAAAACCGACTGATAAACGTTTAACCATCTATGCATTGCTTGACTCGCCGCGTGCGACAGGTGCTTACAAATTCGTGGTTATGCCAGGGCGTGACACGGTTGTGGATGTGCAATCGAAAATCTATCTGCGCGATAAAGTCGGCAAACTGGGGGTTGCACCGTTAACCAGTa                                                                                                                                                                                     <  1:226489/250‑1 (MQ=255)                                                            cgcgGCCTGGCAATTGATACCGCCTTGCCATCGGGTGAAGAATTTCCGCGCTTCAAAGAGTTCTGGATCGAGCGTCCAAAACCGACTGATAAACGTTTAACCATCTATGCATTGCTTGACTCGCCGCGTGCGACAGGTGCTTACAAATTCGTGGTTATGCCAGGGCGTGACACGGTTGTGGATGTGCAATCGAAAATCTATCTGCGCGATAAAGTCGGCAAACTGGGGGTTGCACCGTTAACCAGTATGtt                                                                                                                                                                                 >  1:141949/1‑251 (MQ=255)                                                                 ccTGGCAATTGATACCGCCTTGCCATCGGGTGAAGAATTTCCGCGCTTCAAAGCGTTCTGGCTCGAGCGTCCAAAACCGACTGATAAACGTTTAACCATCTATGCATTGCTTGACTCGCCGCGTGTGACAGGTGCTTACAAATTCGTGGTTATGCCAGGGCGTGACACGGTTGTGGCTGTGCCATCGCAAATCTATCTGCGCGATAAAGTCGGCAAACTGGGGGTTGCACCGTTACCCAGTATGTTCCTGt                                                                                                                                                                            >  2:65387/1‑251 (MQ=255)                                                                   tGGCAATTGATACCGCCTTGCCATCGGGTGAAGAATTTCCGAGCTTCAAAGAGATCTGGATCGAGCGTCCAAAACCGACTGATAAACGTTTAACCATCTATGCATTGCTTGACTCGCCGCGTGCGACAGGTGCTTACAAATTCGTGGTTATGCCAGGGCGTGACCCGGTTGTGGCTGTGCAATCGAAAATCTATCTGCGCGATAAAGTCGGCAAACTGGGGGTGGCACCGTTAACCAGTATGTTCCTGttt                                                                                                                                                                          >  1:361769/1‑251 (MQ=255)                                                                               ccgactttccgtCGGGTGAAGAATTTCCGCGCTTCAACGCGATCTGGATCGCGCGTCCAAAACCGACTGATAAACGTTTCACAATCTATGCATTGCTTGACTCGCCGCGTGCGACAGGTGCTTACAAATTCGTGGTTATGCCAGGGCGTGACACGGTTGTGGGTGTGCAATCGGAAATCTATCTGCGCGATAAAGTCGGCAAACTTGGGGTTGCACCGTTAAACAGTTTATTCATTTTTGGGCCGcaccc                                                                                                                                                               >  2:246755/12‑249 (MQ=255)                                                                                                                                cGAGCGTCCAAAACCGCCTGATAAACGTTTAACCATCTATGCATTGCTTGACTCGCCGCGTGCGACAGGTGCTTACAAATTCGTGGTTATGCCAGGGCGTGACACGGTTGTGGATGTGCAATCGAAAATCTATCTGCGCGATAAAGTCGGCAAACTGGGGGTTGCACCGTTAACCAGTATGTTCCTGTTTGGGCCGAACCAACCGTCGCCTGCAAATAACTATCGTCCGGAGTTGCACGACTCTAACGGtc                                                                                                             >  1:285578/1‑251 (MQ=255)                                                                                                                                   gCGTCCAAAACCGACTGATAAACGTTTAACCATCTATGCATTGCTTGACTCGCCGCGTGCGACAGGTGCTTACAAATTCGTGGTTATGCCAGGGCGTGACACGGTTGTGGATGTGCAATCGAAAATCTATCTGCGCGATAAAGTCGGCAAACTGGGGGTTGCACCGTTAACCAGTATGTTCCTGTTTGGGCCGAACCAACCGTCGCCTGCAAATAACTATCGTCCGGAGTTGCACGACTCTAACGGtctct                                                                                                          <  1:101353/251‑1 (MQ=255)                                                                                                                                    cGTCCAAAACCGACTGATAAACGTTTAACCATCTATGCATTGCTTGACTCGCCGCGTGCGACAGGTGCTTACAAATTCGTGGTTATGCCAGGGCGTGACACGGTTGTGGATGTGCAATCGAAAATCTATCTGCGCGATAAAGTCGGCAAACTGGGGGTTGCACCGTTAACCAGTATGTTCCTGTTTGGGCCGAACCAACCGTCGCCTGCAAATAACTATCGTCCGGAGTTGCACGACTCTAACGGtctctc                                                                                                         <  1:65387/251‑1 (MQ=255)                                                                                                                                         aaaaCCGACTGATAAACGTTTAACCATCTATGCATTGCTTGACTCGCCGCGTGCGACAGGTGCTTACAAATTCGTGGTTATGCCAGGGCGTGACACGGTTGTGGATGTGCAGTCGAAAATCTATCTGCGCGATAAAGTCGGCAAACTGGGGGTTGCACCGTTAACCAGTATGTTCCTGTTTGGGCCGAACCAACCGTCGCCTGCAAATAACTATCGTCCGGAGTTGCACGACTCTAACGGTCTCTCTATcc                                                                                                    >  1:7790/1‑251 (MQ=255)                                                                                                                                            aCCGACTGATAAACGTTTAACCATCTATGCATTGCTTGACTCGCCGCGTGCGACAGGTGCTTACAAATTCGTGGTTATGCCAGGGCGTGACACGGTTGTGGATGTGCAATCGAAAATCTATCTGCGCGATAAAGTCGGCAAACTGGGGGTTGCACCGTTAACCAGTATGTTCCTGTTTGGGCCGAACCAACCGTCGCCTGCAAATAACTATCGTCCGGAGTTGCACGACTCTAACGGTCTCTCTATCCATg                                                                                                 <  1:372759/251‑1 (MQ=255)                                                                                                                                                     tAAACGTTTAACCATCTATGCATTGCTTGACTCGCCGCGTGCGACAGGTGCTTACAAATTCGTGGTTATGCCAGGGCGTGACACGGTTGTGGATGTGCAATCGAAAATCTATCTGCGCGATAAAGTCGGCAAACTGGGGGTTGCACCGTTAACCAGTATGTTCCTGTTTGGGCCGaacca                                                                                                                                                               <  2:102218/180‑1 (MQ=255)                                                                                                                                                     tAAACGTTTAACCATCTATGCATTGCTTGACTCGCCGCGTGCGACAGGTGCTTACAAATTCGTGGTTATGCCAGGGCGTGACACGGTTGTGGATGTGCAATCGAAAATCTATCTGCGCGATAAAGTCGGCAAACTGGGGGTTGCACCGTTAACCAGTATGTTCCTGTTTGGGCCGaacca                                                                                                                                                               >  1:102218/1‑180 (MQ=255)                                                                                                                                                        aCGTATAACCATCTCTGCATTGCTTGACTCGCCGCGTGCGACAGGGGCTTACAAATTCGTGGTTATGCCAGGGCGTGACACGGTTGTGGATGGGCAATCGAAAATCTATCTGCGCAATAAAGCCGGCAAACTGGGGGTTGCACCGTTAACCAGTATGTTCCTGTTTGGGCCGAACCAACCGTCGCCTGCAAATAACTATCGTCCGGAGTTGCACGACTCTAACGGTCTCTCTATCCATGCCGGTAATGGCg                                                                                     <  2:132731/251‑1 (MQ=255)                                                                                                                                                               accaccTATGCATTGCTTGACTCGCCGCGTGCGACAGGTGCTTACAAAATCGTGGTTATGCCAGGGCGTGACACGGTTGTGGATGTGCAGTCGAAAATCTATCTGCGCGATAAAGTCGGCAAACTGGGGGTTGCACCGTTAACCAGTATGTTCCTGTTTGGGCCGAACCAACCGTCGCCTGCAAATAACTATCGTCCGGAGTTGCACGACTCTAACGGTCTCTCTATCCATGCCGGTAATGGCGAATGGAt                                                                              <  2:7790/251‑1 (MQ=255)                                                                                                                                                                    cTATGCATTGCTTGACTCGCCGCGTGCGACAGGTGCTTACAAATTCGTGGTTATGCCAGGGCGTGACACGGTTGTGGATGTGCAATCGAAAATCTATCTGCGCGATAAAGTCGGCAAACTGGGGGTTGCACCGTTAACCAGTATGTTCCTGTTTGGGCCGAACCAACCGTCGCCTGCAAATAACTATCGTCCGGAGTTGCACGACTCTAACGGTCTCTCTATCCATGCCGGTAATGGCGAATGGATCTGGc                                                                         >  1:91307/1‑251 (MQ=255)                                                                                                                                                                                 gACTCGCAGCGTGAGACAGGAGCTTACCAATTCGTGGATATGCAATGGCGTGACACGGTTGTGGATGTACAATCGAAAAGCTATCTGCGCGATAAAGTAGGAAAACTAGGGGTTGCACCGTTAACAAGTATGTTCCTGTTTGGGCCGAACCAACCGTCGCCTGCAACTAACTATCGTCCGGAGTTGCACGACACTAACGGACTCTCTATCCATCCCGGTAATGGCGacttgctagggcctccgttgtcta                                                             >  2:330332/1‑227 (MQ=255)                                                                                                                                                                                                                                    gACACGGTTGTGGATGTGCAATCGAAAATCTATCTGCGCGATAAAGTCGGCAAACTGGGGGTTGCACCGTTAACCAGTATGTTCCTGTTTGGGCCGAACCAACCGTCGCCTGCAAATAACTATCGTCCGGAGTTGCACGACTCTAACGGTCTCTCTATCCATGCCGGTAATGGCGAATGGATCTGGCGTCCGTTGAATAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGAAAACCCGCAAGGCtt         >  1:190753/1‑251 (MQ=255)                                                                                                                                                                                                                                           ttGTGGATGTGCAATCGAAAATCTATCTGCGCGATAAAGTCGGCAAACTGGGGGTTGCACCGTTAACCAGTATGTTCCTGTTTGGGCCGAACCAACCGTCGCCTGCAAATAACTATCGTCCGGAGTTGCACGACTCTAACGGTCTCTCTATCCATGCCGGTAATGGCGAATGGATCTGGCGTCCGTTGAATAACCCGAACCATTTAGCGGTCAGCAGCTTCTCCATGGAAAACCCGCAAGGCTTTGGTCTg  >  1:89119/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              TGCTCGGGGCCAGCTATTTCCGCGTGATTGGTGCAGGTCAGGTTTATGGCCTTTCTGCCCGCGGCCTGGCAATTGATACCGCCTTGCCATCGGGTGAAGAATTTCCGCGCTTCAAAGAGTTCTGGATCGAGCGTCCAAAACCGACTGATAAACGTTTAACCATCTATGCATTGCTTGACTCGCCGCGTGCGACAGGTGCTTACAAATTCGTGGTTATGCCAGGGCGTGACACGGTTGTGGATGTGCAGTCGAAAATCTATCTGCGCGATAAAGTCGGCAAACTGGGGGTTGCACCGTTAACCAGTATGTTCCTGTTTGGGCCGAACCAACCGTCGCCTGCAAATAACTATCGTCCGGAGTTGCACGACTCTAACGGTCTCTCTATCCATGCCGGTAATGGCGAATGGATCTGGCGTCCGTTGAATAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGAAAACCCGCAAGGCTTTGGTCTG  >  NC_000913/1109786‑1110270

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