Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,557,606 G→A 12.1% intergenic (+83/‑96) PP_3141 → / → PP_3142 WecB/TagA/CpsF family glycosyl transferase/sugar transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,557,6060GA12.1% 80.1 / 3.4 33intergenic (+83/‑96)PP_3141/PP_3142WecB/TagA/CpsF family glycosyl transferase/sugar transferase
Reads supporting (aligned to +/- strand):  ref base G (15/14);  new base A (2/2);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01

GCTATTCCTGTTGCGTACCCTGTTGTACCGGCGGGCGGCGGTTTCGGCCAGGCGGGCCCTATGAATCGGCCAAAAATCGCAAAAGCGGTGATCCCGGGTACGCCAGCGCCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATGGTTTTCGAACCCAGAAGCAGTCGTTCCTTACTCCAGAGAAGAAGCAG  >  NC_002947/3557460‑3557751
                                                                                                                                                   |                                                                                                                                                  
gCTATTCCTGTTGCGTACCCTGTTGTACCGGCGGGCGGCGGTTTCGGCCAGGCGGGCCCTATGAATCGGCCAAAAATCGCAAAAGCGGTGATCCCGGGTACGCCAGCGCCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGcc                                                                                                                                                  <  7:265312/149‑1 (MQ=255)
        tGTTGCGTACCCTGTTGTACCGGCGGGCGGCGGTTTCGGCCAGGCGGGCCCTATGAATCGGCCAAAAATCGCAAAAGCGGTGATCCCGGGTACGCCAGCGCCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGcc                                                                                                                                                  <  8:47472/141‑1 (MQ=255)
        tGTTGCGTACCCTGTTGTACCGACGGGCGGCGGTTTCGGCCAGGCGGGCCCTATGAATCGGCCAAAAATCGCAAAAGCGGTGATCCCGGGTACGCCAGCGCCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGcc                                                                                                                                                  >  7:47472/1‑141 (MQ=255)
                 cccTGTTGTACCGGCGGGCGGCGGTTTCGGCCAGGCGGGCCCTATGAATCGGCCAAAAATCGCAAAAGCGGTGATCCCGGGTACGCCAGCGCCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACgg                                                                                                                                 >  7:187467/1‑149 (MQ=255)
                           ccGGCGGGCGGCGGTTTCGGCCAGGCGGGCCCTATGAATCGGCCAAAAATCGCAAAAGCGGTGATCCCGGGTACGCCAGCGCCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCtt                                                                                                                        <  6:52736/148‑1 (MQ=255)
                           ccGGCGGGCGGCGGTTTCGGCCAGGCGGGCCCTATGAATCGGCCAAAAATCGCAAAAGCGGTGATCCCGGGTACGCCAGCGCCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCtt                                                                                                                        >  1:45560/1‑148 (MQ=255)
                                             ggCCAGGCGGGCCCTATGAATCGGCCAAAAATCGCAAAAGCGGTGATCCCGGGTACGCCAGCGCCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACtt                                                                                                     >  4:31860/1‑149 (MQ=255)
                                               ccAGGCGGGCCCTATGAATCGGCCAAAAATCGCAAAAGCGGTGATCCCGGGTACGCCAGCGCCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGttttt                                                                                                                            <  4:45978/124‑1 (MQ=255)
                                               ccAGGCGGGCCCTATGAATCGGCCAAAAATCGCAAAAGCGGTGATCCCGGGTACGCCAGCGCCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGttttt                                                                                                                            >  3:45978/1‑124 (MQ=255)
                                                          cTATGAATCGGCCAAAAATCGCAAAAGCGGTGATCCCGGGTACGCCAGCGCCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGa                                                                                       <  1:302661/149‑1 (MQ=255)
                                                           tATGAATCGGCCAAAAATCGCAAAAGCGGTGATCCCAGGGTACGCCAGCGCCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTACTGTTTGATgg                                                                                        >  1:214358/1‑148 (MQ=255)
                                                                                  aaGCGGTGATCCCGGGTACGCCAGCGCCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGATGTTTGATGGAAGGCACTTGCTACGCTCATCTtga                                                               >  7:171546/1‑149 (MQ=255)
                                                                                          aTCCCGGGTACGCCAGCGCCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTa                                                                           >  8:146412/1‑129 (MQ=255)
                                                                                          aTCCCGGGTACGCCAGCGCCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTa                                                                           <  7:146412/129‑1 (MQ=255)
                                                                                              cGGGTACGCCAGCGCCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGaaa                                                   >  2:185284/1‑149 (MQ=255)
                                                                                                  gTACGCCAGCGCCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATgg                                                <  7:255581/149‑1 (MQ=255)
                                                                                                      gCCAGCGCCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATGGtttt                                            >  2:84169/1‑149 (MQ=255)
                                                                                                           cgcCGGGTGGCTGCTACAGTTATATAAATCTCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTATGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATGGTTTTCGAAc                                       <  8:187467/149‑1 (MQ=255)
                                                                                                           cgcCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATGGTTTTCGAAc                                       <  1:185284/149‑1 (MQ=255)
                                                                                                             ccGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATTGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATGGTTTTCGAAcc                                      <  1:84169/149‑1 (MQ=255)
                                                                                                                            aGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATGGTTTTCGAACCCAGAAGCAGTCGTTcc                      >  3:303968/1‑149 (MQ=255)
                                                                                                                            aGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATGGTTTTCGAACCCAGAAGCAGTCGTTcc                      <  2:45560/149‑1 (MQ=255)
                                                                                                                              tGATATCAATCGCTTGAGGGCACCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATGGTTTTCGAACCCAGAAGc                               >  7:35645/1‑138 (MQ=255)
                                                                                                                              tGATATCAATCGCTTGAGGGCACCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATGGTTTTCGAACCCAGAAGc                               >  7:40698/1‑138 (MQ=255)
                                                                                                                              tGATATCAATCGCTTGAGGGCACCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATGGTTTTCGAACCCAGAAGc                               <  8:35645/138‑1 (MQ=255)
                                                                                                                              tGATATCAATCGCTTGAGGGCACCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATGGTTTTCGAACCCAGAAGc                               <  8:40698/138‑1 (MQ=255)
                                                                                                                                 tatCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATGGTTTTCGAACCCAGAAGCAGTCGTTCCTTACt                 >  5:333211/1‑149 (MQ=255)
                                                                                                                                  atCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATGGTTTTCGAACCCAGAAGCAGTCGTTCCTTACTc                >  2:167780/1‑149 (MQ=255)
                                                                                                                                   tCAATCGCTTGAGGGCGACACGAGCCTGCTGCACGTTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATGGTTTTCGAACCCAGAAGCAGTCGTTCCTTACTc                <  8:171546/148‑1 (MQ=255)
                                                                                                                                      aTCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATGGTTTTCGAACCCAGAAGCAGTCGTTCCTTACTCCa              >  3:16625/1‑147 (MQ=255)
                                                                                                                                         gCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATGGTTTTCGAACCCAGAAGCAGTCGTTCCTTACTCCAGagaa         <  6:333211/149‑1 (MQ=255)
                                                                                                                                              aGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTACTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATGGTTTTCGAACCCAGAAGCAGTCGTTCCTTACTCCAGAGAAGAAGc    <  2:214358/149‑1 (MQ=255)
                                                                                                                                                ggCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATGGTTTTCGAACCCAGAAGCAGTCGTTCCTTACTCCAGAGAAGAAGCAg  >  1:178546/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTATTCCTGTTGCGTACCCTGTTGTACCGGCGGGCGGCGGTTTCGGCCAGGCGGGCCCTATGAATCGGCCAAAAATCGCAAAAGCGGTGATCCCGGGTACGCCAGCGCCGGGTGGCTGCTACAGTGATATCAATCGCTTGAGGGCGCCACGAGCCTGCTGCACGGTTTTTGCTTAAGTCGGTTTTCAAGACTTGCTGTTTGATGGAAGGCACTTGCTACGCTCATCTTGATGAGGTCTGAAATGGTTTTCGAACCCAGAAGCAGTCGTTCCTTACTCCAGAGAAGAAGCAG  >  NC_002947/3557460‑3557751

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: