Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,025,719:1 +T 11.1% coding (719/1110 nt) wbpZ → glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,025,7191.T11.1% 156.3 / 15.9 45coding (719/1110 nt)wbpZglycosyl transferase
Reads supporting (aligned to +/- strand):  ref base . (21/19);  new base T (2/3);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

GCATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTGGAAGCGGCAATGTA  >  NC_002947/2025572‑2025866
                                                                                                                                                     |                                                                                                                                                   
gcaTGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGgctgct                                                                                                                                                    <  1:142332/149‑1 (MQ=255)
    gagcgagcGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGc                                                                                                                                                >  6:58441/1‑149 (MQ=255)
      gcgagcgCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAGCAGCCTGCTGCTTGAGCtg                                                                                                                                              <  5:18045/149‑2 (MQ=255)
         agcgCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAgg                                                                                                                                           <  8:214246/149‑1 (MQ=255)
         agcgCGTTTTCCTGTTTGTACGCGCGTGGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAgg                                                                                                                                           <  3:232293/149‑1 (MQ=255)
          gcgcTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGc                                                                                                                                          >  8:36194/1‑149 (MQ=255)
          gcgcTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGc                                                                                                                                          >  6:327049/1‑149 (MQ=255)
                    cTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAGCAGCCTGCTGCTTGAGCtg                                                                                                                                              >  6:319365/1‑134 (MQ=255)
                    cTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAGCAGCCTGCTGCTTGAGCtg                                                                                                                                              <  5:319365/135‑2 (MQ=255)
                          gTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGtt                                                                                                                          >  6:377375/1‑149 (MQ=255)
                          gTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGt                                                                                                                           >  7:210323/1‑148 (MQ=255)
                               cgcgTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTgg                                                                                                                     <  8:210323/149‑1 (MQ=255)
                                  gTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCa                                                                                                                  <  5:377375/149‑1 (MQ=255)
                                         tactacAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCtt                                                                                                                                        <  1:352991/120‑1 (MQ=255)
                                         tactacAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCtt                                                                                                                                        >  2:352991/1‑120 (MQ=255)
                                            tacAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGc                                                                                                           <  7:153234/146‑1 (MQ=255)
                                            tacAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGc                                                                                                           >  8:153234/1‑146 (MQ=255)
                                              cAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGATGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGAc                                                                                                       >  1:295327/1‑149 (MQ=255)
                                               caGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCGGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGa                                                                                                     <  4:3166/148‑1 (MQ=255)
                                                  gggCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGa                                                                                                  <  2:353748/149‑1 (MQ=255)
                                                    gCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATg                                                                                                   <  8:201976/146‑1 (MQ=255)
                                                    gCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATg                                                                                                   >  7:201976/1‑146 (MQ=255)
                                                      tGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGACAGGAAGACGATGACa                                                                                                >  2:100005/1‑147 (MQ=255)
                                                              cTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCgct                                                                                      >  3:172184/1‑149 (MQ=255)
                                                                    cAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTgg                                                                                 <  5:58441/148‑1 (MQ=255)
                                                                               tATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATg                                                                     >  5:251836/1‑149 (MQ=255)
                                                                                 tCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGtt                                                                        >  5:189668/1‑144 (MQ=255)
                                                                                 tCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGtt                                                                        <  6:189668/144‑1 (MQ=255)
                                                                                 tCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCa                                                                   <  4:172184/149‑1 (MQ=255)
                                                                                   aaCCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCAtt                                                                 >  7:93349/1‑149 (MQ=255)
                                                                                   aaCCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCAtt                                                                 >  1:178711/1‑149 (MQ=255)
                                                                                    aCCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTc                                                                <  8:331874/149‑1 (MQ=255)
                                                                                      cGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCg                                                               >  5:341483/1‑148 (MQ=255)
                                                                                               gTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTccatc                                                     >  5:89799/1‑149 (MQ=255)
                                                                                                     atcatcGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCt                                               <  6:89799/149‑1 (MQ=255)
                                                                                                                   gCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGctgctg                                                                                  <  2:361715/101‑1 (MQ=255)
                                                                                                                   gCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGctgctg                                                                                  >  1:361715/1‑101 (MQ=255)
                                                                                                                    cTTCCTGGAAGCGCACCTGAAGCAGCAGGGTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGGAGGCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGGTTTGTTATGCATTCGTATATCCATCCCATCTGCGTTCAGAGTCGt                                 <  1:210257/149‑1 (MQ=255)
                                                                                                                     ttCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATg                                                                     <  7:169412/112‑1 (MQ=255)
                                                                                                                     ttCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATg                                                                     >  8:169412/1‑112 (MQ=255)
                                                                                                                        cTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCgg                             >  1:233483/1‑149 (MQ=255)
                                                                                                                         tCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCTGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCgg                             >  3:234044/1‑148 (MQ=255)
                                                                                                                           gAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGAAAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCa                           <  7:36194/148‑1 (MQ=255)
                                                                                                                               ctgaGCTCAGCAGCCTGCTGCTTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTtct                      >  1:8331/4‑149 (MQ=255)
                                                                                                                                    cTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGGCGTTCGGCATTTCTctgct                 <  2:295327/149‑1 (MQ=255)
                                                                                                                                     ctcagcagcCTGCTGCTTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTctgct                 <  2:8331/146‑1 (MQ=255)
                                                                                                                                                    ctgAGCTCAGGCTTGACAATGTTCAGTTCCTGCGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATTCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTGGAAGCGGCAATGTa  <  8:93349/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTGGAAGCGGCAATGTA  >  NC_002947/2025572‑2025866

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: