Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,109,964 G→A 100% C185C (TGC→TGT mdoC ← membrane protein required for modification of periplasmic glucan

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,109,9640GA82.8% 62.5 / 7.8 29C185C (TGC→TGTmdoCmembrane protein required for modification of periplasmic glucan
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base A (20/4);  total (24/5)
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.86e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACTCAGAATGGGCGGATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGGCATAAAAAAATCACCGATAGTTTTACCATCGAGAATTTTTTATTCGTTTTATCAGAATTTTCTAAAT  >  W3110S.gb/1109901‑1110031
                                                               |                                                                   
aCTCAGAATGGGCGGATACACAATAAAAATCGTTCTTCTTATTACCGCATACCCGATGCCGAGACATa                                                                 >  1:2688945/1‑68 (MQ=255)
aCTCAGAATGGGCGGATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGACAt                                                                  >  1:3565119/1‑67 (MQ=255)
aCTCAGAATGGGCGGATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATa                                                                 >  1:2748952/1‑68 (MQ=255)
      aaTGGGCGGATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGGCATAAAAAAAt                                                          <  1:2234160/69‑1 (MQ=255)
       aTGGGCGGATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGGCATAAAAAAATc                                                         >  1:3667378/1‑69 (MQ=255)
          ggCGGATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCAcc                                                      <  1:3278982/69‑1 (MQ=255)
            cGGATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCACCGa                                                    <  1:2206555/69‑1 (MQ=255)
                  cacaATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATaaaaa                                                             >  1:943814/1‑54 (MQ=255)
                  cacaATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATaaaaa                                                             >  1:3585840/1‑54 (MQ=255)
                    caATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGGCATaaaaaa                                                            >  1:1541333/1‑53 (MQ=255)
                    caATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGGCATaaaaaa                                                            >  1:377070/1‑53 (MQ=255)
                    caATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGGCATaaaaaa                                                            >  1:2727067/1‑53 (MQ=255)
                        aaaaaTCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCa                                         >  1:3328560/1‑68 (MQ=255)
                         aaaaTCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATaaaaaa                                                            >  1:2500007/1‑48 (MQ=255)
                           aaTCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAAt                                                          >  1:305976/1‑48 (MQ=255)
                           aaTCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAAt                                                          >  1:921128/1‑48 (MQ=255)
                            aTCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCgag                                    >  1:2876808/1‑69 (MQ=255)
                             tCGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCACCGa                                                    >  1:260592/1‑52 (MQ=255)
                              cGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCGAGaa                                  >  1:2915719/1‑69 (MQ=255)
                              cGTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCGAGaa                                  >  1:584131/1‑69 (MQ=255)
                               gTTCTTCTTATTACCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCGAGAAt                                 >  1:362139/1‑69 (MQ=255)
                                        attaCCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCGAGAATTTTTTATTc                        >  1:1239343/1‑69 (MQ=255)
                                        attaCCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCGAGAATTTTTTATTc                        >  1:2565354/1‑69 (MQ=255)
                                           aCCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCgaga                                   >  1:74549/1‑55 (MQ=255)
                                           aCCGCATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCGAGAATTTTTTATTCGt                      >  1:3571933/1‑68 (MQ=255)
                                              gCATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCGAGAATTTTTTATTCGtttt                   <  1:934147/68‑1 (MQ=255)
                                               cATAACCGATGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCGAGAATTTTTTATTCGtttt                   >  1:207487/1‑67 (MQ=255)
                                                       aTGCCGAGACATAAAAAAATCACCGATAGTTTTACCATCGAGAATTTTTTATTCGTTTTATCAGAAtt          <  1:1778631/68‑1 (MQ=255)
                                                         gCCGAGACATAAAAAAATCACCGATAGTTTTACCATCGAGAATTTTTTATTCGTTTTATCAGAAtttt        >  1:1648838/1‑68 (MQ=255)
                                                              ggCATAAAAAAATCACCGATAGTTTTACCATCGAGAATTTTTTATTCGTTTTATCAGAATTTTCTAAat  >  1:660219/1‑69 (MQ=255)
                                                               |                                                                   
ACTCAGAATGGGCGGATACACAATAAAAATCGTTCTTCTTATTACCGCATAACCGATGCCGAGGCATAAAAAAATCACCGATAGTTTTACCATCGAGAATTTTTTATTCGTTTTATCAGAATTTTCTAAAT  >  W3110S.gb/1109901‑1110031

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