Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,025,714	G→T	10.9%	Q238H (CAG→CAT) ‡	wbpZ →	glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,025,714	0	G	T	10.9%	131.1 / 7.2	46	Q238H (CAG→CAT) ‡	wbpZ	glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (21/20);  new base T (2/3);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

GCATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCT  >  NC_002947/2025572‑2025851                                                                                                                                                |                                                                                                                                           gcaTGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGgctgc‑t                                                                                                                                     <  1:142332/149‑1 (MQ=255)     gagcgagcGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGc                                                                                                                                 >  6:58441/1‑149 (MQ=255)       gcgagcgCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTC‑AGCAGCCTGCTGCTTGAGCtg                                                                                                                               <  5:18045/149‑2 (MQ=255)          agcgCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAgg                                                                                                                            <  8:214246/149‑1 (MQ=255)          agcgCGTTTTCCTGTTTGTACGCGCGTGGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAgg                                                                                                                            <  3:232293/149‑1 (MQ=255)           gcgcTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGc                                                                                                                           >  6:327049/1‑149 (MQ=255)           gcgcTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGc                                                                                                                           >  8:36194/1‑149 (MQ=255)                     cTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTC‑AGCAGCCTGCTGCTTGAGCtg                                                                                                                               >  6:319365/1‑134 (MQ=255)                     cTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTC‑AGCAGCCTGCTGCTTGAGCtg                                                                                                                               <  5:319365/135‑2 (MQ=255)                           gTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGtt                                                                                                           >  6:377375/1‑149 (MQ=255)                           gTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGt                                                                                                            >  7:210323/1‑148 (MQ=255)                                cgcgTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTgg                                                                                                      <  8:210323/149‑1 (MQ=255)                                   gTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCa                                                                                                   <  5:377375/149‑1 (MQ=255)                                          tactacAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCtt                                                                                                                         <  1:352991/120‑1 (MQ=255)                                          tactacAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCtt                                                                                                                         >  2:352991/1‑120 (MQ=255)                                             tacAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGc                                                                                            >  8:153234/1‑146 (MQ=255)                                             tacAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGc                                                                                            <  7:153234/146‑1 (MQ=255)                                               cAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGATGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGAc                                                                                        >  1:295327/1‑149 (MQ=255)                                                caGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑GGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGa                                                                                      <  4:3166/148‑1 (MQ=255)                                                   gggCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGa                                                                                   <  2:353748/149‑1 (MQ=255)                                                     gCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATg                                                                                    >  7:201976/1‑146 (MQ=255)                                                     gCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATg                                                                                    <  8:201976/146‑1 (MQ=255)                                                       tGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGACAGGAAGACGATGACa                                                                                 >  2:100005/1‑147 (MQ=255)                                                               cTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCgct                                                                       >  3:172184/1‑149 (MQ=255)                                                                     cAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTgg                                                                  <  5:58441/148‑1 (MQ=255)                                                                                tATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATg                                                      >  5:251836/1‑149 (MQ=255)                                                                                  tCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGtt                                                         >  5:189668/1‑144 (MQ=255)                                                                                  tCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGtt                                                         <  6:189668/144‑1 (MQ=255)                                                                                  tCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCa                                                    <  4:172184/149‑1 (MQ=255)                                                                                    aaCCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCAtt                                                  >  1:178711/1‑149 (MQ=255)                                                                                    aaCCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCAtt                                                  >  7:93349/1‑149 (MQ=255)                                                                                     aCCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTc                                                 <  8:331874/149‑1 (MQ=255)                                                                                       cGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCg                                                >  5:341483/1‑148 (MQ=255)                                                                                                gTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTccatc                                      >  5:89799/1‑149 (MQ=255)                                                                                                      atcatcGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCt                                <  6:89799/149‑1 (MQ=255)                                                                                                                    gCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGctgctg                                                                   <  2:361715/101‑1 (MQ=255)                                                                                                                    gCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGctgctg                                                                   >  1:361715/1‑101 (MQ=255)                                                                                                                     cTTCCTGGAAGCGCACCTGAAGCAGCAGGGTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGGAGGCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGGTTTGTTATGCATTCGTATATCCATCCCATCTGCGTTCAGAGTCGt                  <  1:210257/149‑1 (MQ=255)                                                                                                                      ttCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATg                                                      <  7:169412/112‑1 (MQ=255)                                                                                                                      ttCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATg                                                      >  8:169412/1‑112 (MQ=255)                                                                                                                         cTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCgg              >  1:233483/1‑149 (MQ=255)                                                                                                                          tCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCTGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCgg              >  3:234044/1‑148 (MQ=255)                                                                                                                            gAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGAAAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCa            <  7:36194/148‑1 (MQ=255)                                                                                                                                ctgaGCTC‑AGCAGCCTGCTGCTTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTtct       >  1:8331/4‑149 (MQ=255)                                                                                                                                     cTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGGCGTTCGGCATTTCTctgct  <  2:295327/149‑1 (MQ=255)                                                                                                                                      ctcagcagcCTGCTGCTTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTctgct  <  2:8331/146‑1 (MQ=255)                                                                                                                                                |                                                                                                                                           GCATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCT  >  NC_002947/2025572‑2025851

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