Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,963,488 T→C intergenic (+51/‑5) wzzE → / → wecB ECA polysaccharide chain length modulation protein/UDP‑N‑acetylglucosamine 2‑epimerase (non‑hydrolyzing)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,963,4880TC95.0% 62.4 / ‑5.7 20intergenic (+51/‑5)wzzE/wecBECA polysaccharide chain length modulation protein/UDP‑N‑acetylglucosamine 2‑epimerase (non‑hydrolyzing)
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (11/8);  total (12/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAAAAGATCCTTTTTTTGAGGCTAAAGTTTGCGTCACTGCGCAGCAT  >  NZ_CP009273/3963355‑3963609
                                                                                                                                     |                                                                                                                         
gtGCCTTCCTGATGATTATGTGGGTCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAACCGAt                                                                                                                       <  1:326367/138‑1 (MQ=255)
 tGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAACCGAtgt                                                                                                                     >  1:18379/1‑139 (MQ=255)
   ccTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAACCGAtgtg                                                                                                                    >  1:22699/1‑138 (MQ=255)
                       ggCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAACCGATGTGAAAGTACTGACTGTATTTgg                                                                                                >  2:354595/1‑138 (MQ=255)
                                  ggggCTGATCGGGGCGGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAACCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGa                                                                                    <  2:310708/139‑1 (MQ=255)
                                   gggCTGATCGGGGTGGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGCCGCTCGCCGCTTATTCGAAGAGAACCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGaa                                                                                   <  2:370399/139‑1 (MQ=255)
                                        gATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAACCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCAt                                                                              >  1:334092/1‑139 (MQ=255)
                                        gATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGGGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAACCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCAt                                                                              >  1:13403/1‑139 (MQ=255)
                                                   gtgtCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAACCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGaa                                                                                   >  2:421301/1‑123 (MQ=255)
                                                   gtgtCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAACCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGaa                                                                                   <  1:421301/123‑1 (MQ=255)
                                                      tCGCATTACCCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAACCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGt                                                                <  2:202819/139‑1 (MQ=255)
                                                      tCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAACCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGt                                                                <  1:354595/139‑1 (MQ=255)
                                                                                     acacTGCTGCGGTGAGCGCAAAGGCGCTCTCCGCTTATTCGAAGAGAACCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATTGCGCCGTTGGTGCATGCGTTGTCAAAAGATCCtttttt                                 <  1:251590/139‑1 (MQ=255)
                                                                                           ctgcGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAAAAGATCCTTTTTTTGAGGc                           >  1:209588/1‑139 (MQ=255)
                                                                                           ctgcGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAACCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAAAAGATCCTTTTTTTGAGGc                           >  1:94623/1‑139 (MQ=255)
                                                                                           ctgcGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAACCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAAAAGATCCTTTTTTTGAGGc                           >  1:309668/1‑139 (MQ=255)
                                                                                                gTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAACCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAAAAGATCCTTTTTTTGAGGCTAAAg                      >  2:97391/1‑139 (MQ=255)
                                                                                                   aGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAACCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAAAAGATCCTTTTTTTGAGGCTAAAGttt                   >  2:366745/1‑139 (MQ=255)
                                                                                                             cgcTCGCCGCTTATTCGAAGAGAACCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAAAAGATCCTTTTTTTGAGGCTAAAGTTTGCGTCACTgc         >  2:119317/1‑139 (MQ=255)
                                                                                                                    cgcTTATTCGAAGAGAACCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAAAAGATCCTTTTTTTGAGGCTAAAGTTTGCGTCACTGCGCAGCAt  <  2:22699/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                         
GTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAAAAGATCCTTTTTTTGAGGCTAAAGTTTGCGTCACTGCGCAGCAT  >  NZ_CP009273/3963355‑3963609

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