Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,946,443 T→C 10.2% D421D (GAT→GAC bcsA → cellulose synthase, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,946,4430TC10.2% 96.3 / 3.3 39D421D (GAT→GACbcsAcellulose synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base T (14/21);  new base C (2/2);  total (16/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 = 9.77e-01

TGAACGCAACCTGGGGCGTTTCCGTAAAACGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTCCTGTGCGGTGATTCGTCGTAAGCC  >  W3110S.gb/3946374‑3946505
                                                                     |                                                              
tGAACGCAACCTGGGGCGTTTCCGTAAAACGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATg                                                               >  1:3568894/1‑71 (MQ=255)
   aCGCAACCTGGGGCGTTTCCGTAAAACGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGc                                                             >  1:1950107/1‑70 (MQ=255)
     gCAACCTGGGGCGTTTCCGTAAAACGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAAc                                                          <  1:721395/71‑1 (MQ=255)
              ggCGTTTCCGTAAAACGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTgg                                                 <  1:1244162/71‑1 (MQ=255)
              ggCGTTTCCGTAAAACGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTgg                                                 <  1:1262467/71‑1 (MQ=255)
              ggCGTTTCCGTAAAACGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTgg                                                 <  1:2631458/71‑1 (MQ=255)
              ggCGTTTCCGTAAAACGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTgg                                                 <  1:729618/71‑1 (MQ=255)
              ggCGTTTCCGTAAAACGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTgg                                                 <  1:2888407/71‑1 (MQ=255)
                  tttCCGTAAAACGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAAc                                                          <  1:3436462/58‑1 (MQ=255)
                  tttCCGTAAAACGCCGAACGAAGGCACGCTGTTCTATGGGCTGGTGCAGGATGGCAACGATATGTgg                                                 <  1:1313708/67‑1 (MQ=255)
                  tttCCGTAAAACACCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGACGGCAATGATATGTGGGa                                               <  1:3202443/69‑1 (MQ=255)
                  tttCCGTAAAACACAGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGACGGCAATGATATGTGGGa                                               <  1:3780175/69‑1 (MQ=255)
                    tCCGTAAAACGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATAtg                                                    >  1:2421328/1‑62 (MQ=255)
                    tCCGTAAAACGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATAtg                                                    >  1:3117790/1‑62 (MQ=255)
                    tCCGTAAAACGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTGGGACGc                                            >  1:1082822/1‑70 (MQ=255)
                    tCCGTAAAACGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTGGGACGc                                            >  1:1905689/1‑70 (MQ=255)
                      cGTAAAACGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTGGGACGCCa                                          <  1:349554/70‑1 (MQ=255)
                      cGTAAAACGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTGGGACGCCa                                          <  1:228993/70‑1 (MQ=255)
                      cGTAAAACGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTGGGACGCCAc                                         >  1:132231/1‑71 (MQ=255)
                       gTAAAACACCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGACGGCaat                                                          >  1:2316110/1‑52 (MQ=38)
                       gTAAAACACCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGACGGCa                                                            >  1:1767524/1‑51 (MQ=38)
                        tAAAACGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTgg                                                 <  1:1469350/61‑1 (MQ=255)
                         aaaaCGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTGGGACGCCACtt                                       >  1:2659440/1‑70 (MQ=255)
                             cGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAAc                                                          >  1:1554657/1‑47 (MQ=255)
                                cgaacgaaGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTGGGa                                               <  1:238571/55‑1 (MQ=255)
                                      aaGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTcc                         <  1:3261608/71‑1 (MQ=255)
                                       aGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTcc                         <  1:187815/70‑1 (MQ=255)
                                        ggCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTCCtg                       <  1:2003653/71‑1 (MQ=255)
                                         gCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTcc                         <  1:134968/68‑1 (MQ=255)
                                            cGCTGTTCTATGGTCTGGTGCATGATGGCAACGATACGTGGGACGCCACTttc                                     <  1:2194949/53‑1 (MQ=39)
                                            cGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTCCTGTGCg                   >  1:1071930/1‑71 (MQ=255)
                                                   cTATGGTCTGGTGCAGGATGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTCCTGTGCGGTGATtc            >  1:608429/1‑71 (MQ=255)
                                                      tGGTCTGGTGCAGGATGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTCCTGTGCGGTGATtcg           <  1:428034/69‑1 (MQ=255)
                                                      tGGTCTGGTGCAGGATGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTCCTGTGCGGTGATtcg           <  1:3727939/69‑1 (MQ=255)
                                                           tGGTGCAGGATGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTCCTGTGCGGTGATTCGTCGTAAg    >  1:872812/1‑71 (MQ=255)
                                                           tGGTGCAGGATGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTCCTGTGCGGTGATTCGTCGTAAg    >  1:1768953/1‑71 (MQ=255)
                                                             gTGCAGGATGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTCCTGTGCGGTGATTCGTCGTAAGcc  >  1:38376/1‑71 (MQ=255)
                                                               gCAGGATGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTcc                         <  1:33121/46‑1 (MQ=255)
                                                               gCAGGATGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTcc                         <  1:58881/46‑1 (MQ=255)
                                                                     tGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTCCt                        >  1:3380280/1‑41 (MQ=255)
                                                                     tGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTCCt                        >  1:3367871/1‑41 (MQ=255)
                                                                     tGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTCCt                        >  1:840656/1‑41 (MQ=255)
                                                                     tGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTCCTGTGCGGTGa               <  1:2226013/50‑1 (MQ=255)
                                                                     tGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTCCTGTGCGGTGa               <  1:1574460/50‑1 (MQ=255)
                                                                     tGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTCCTGTGCGGTGa               <  1:1350873/50‑1 (MQ=255)
                                                                     tGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTCCTGTGCGGTGa               <  1:2801423/50‑1 (MQ=255)
                                                                     tGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTCCTGTGCGGTGa               <  1:3013227/50‑1 (MQ=255)
                                                                     tGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTCCTGTGCGGTGATTCGTCGTAAGc   >  1:2982726/1‑62 (MQ=255)
                                                                     |                                                              
TGAACGCAACCTGGGGCGTTTCCGTAAAACGCCGAACGAAGGCACGCTGTTCTATGGTCTGGTGCAGGATGGCAACGATATGTGGGACGCCACTTTCTTCTGCGGTTCCTGTGCGGTGATTCGTCGTAAGCC  >  W3110S.gb/3946374‑3946505

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