Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,025,706	Δ1 bp	8.5%	coding (706/1110 nt)	wbpZ →	glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,025,706	0	A	.	8.5%	142.7 / 7.8	47	coding (706/1110 nt)	wbpZ	glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (22/21);  new base . (2/2);  total (24/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 = 8.67e-01

TACTGGCGCGCACGCATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCT  >  NC_002947/2025559‑2025851                                                                                                                                                     |                                                                                                                                                   tACTGGCGCGCACGCATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCaa                                                                                                                                                   <  4:63020/149‑1 (MQ=255)        gcgcACGCATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAagc                                                                                                                                                 >  8:137556/1‑144 (MQ=255)        gcgcACGCATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAagc                                                                                                                                                 <  7:137556/144‑1 (MQ=255)              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                                                                                                                           >  8:36194/1‑149 (MQ=255)                        gcgcTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGc                                                                                                                           >  6:327049/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                                                                                            <  7:153234/146‑1 (MQ=255)                                                          tacAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGc                                                                                            >  8:153234/1‑146 (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                                                                                    <  8:201976/146‑1 (MQ=255)                                                                  gCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATg                                                                                    >  7:201976/1‑146 (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                                                  >  7:93349/1‑149 (MQ=255)                                                                                                 aaCCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCAtt                                                  >  1:178711/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                                                                   >  1:361715/1‑101 (MQ=255)                                                                                                                                 gCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGctgctg                                                                   <  2:361715/101‑1 (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)                                                                                                                                                     |                                                                                                                                                   TACTGGCGCGCACGCATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGG‑TGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCT  >  NC_002947/2025559‑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: ‑