Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,112,202 T→C 27.3% F431L (TTT→CTT)  mdoG → glucan biosynthesis protein, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,112,2020TC27.3% 68.3 / 29.7 44F431L (TTT→CTT) mdoGglucan biosynthesis protein, periplasmic
Reads supporting (aligned to +/- strand):  ref base T (18/14);  new base C (3/9);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 9.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACGGGGGATGTGAAGCAGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGTCACAGCGCAAACCAGCATTGGTGA  >  W3110S.gb/1112134‑1112269
                                                                    |                                                                   
cAACGGGGGATGTGAAGCAGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGAttt                                                                    <  1:3389774/70‑1 (MQ=255)
cAACGGGGGATGTGAAGCAGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGAttt                                                                    <  1:641376/70‑1 (MQ=255)
cAACGGGGGATGTGAAGCAGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGAttt                                                                    <  1:2432232/70‑1 (MQ=255)
 aaCGGGGGATGTGAAGCAGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGAtctt                                                                   <  1:467208/70‑1 (MQ=255)
            tgAAGCAGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCaga                                                        <  1:1437196/70‑1 (MQ=255)
                  aGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGaa                                                  >  1:2467598/1‑70 (MQ=255)
                   gTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGAt                                                     >  1:169750/1‑66 (MQ=255)
                   gTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGaaaa                                                <  1:4049640/71‑1 (MQ=255)
                   gTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGaaaa                                                <  1:353601/71‑1 (MQ=255)
                    tCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGaaaa                                                <  1:2684834/70‑1 (MQ=255)
                      gAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGaaaaaa                                              <  1:2847826/70‑1 (MQ=255)
                       aaCCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATCTTACCGGCGCAGAGATGAAAAAACt                                            <  1:1034143/71‑1 (MQ=255)
                        aCCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCgg                                                              >  1:2983977/1‑52 (MQ=255)
                        aCCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTg                                           >  1:1215244/1‑71 (MQ=255)
                        aCCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTg                                           >  1:5076162/1‑71 (MQ=255)
                        aCCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACAGGCGCAGAGATGAAAAAAc                                             >  1:84940/1‑69 (MQ=255)
                        aCCTGATTCGCAAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCgg                                                              >  1:4401525/1‑52 (MQ=255)
                         ccTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGc                                          >  1:3110849/1‑71 (MQ=255)
                         ccTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATCTTACCGGCGCAGAGATGaaaa                                                <  1:3474241/65‑1 (MQ=255)
                         ccTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATCTTACCGGCGCAGAGATGaaaa                                                <  1:2584547/65‑1 (MQ=255)
                         ccTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATCTTACCGGCGCAGAGATGaaaa                                                <  1:3396266/65‑1 (MQ=255)
                            gATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGAttt                                                                    <  1:2916193/42‑1 (MQ=255)
                            gATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGAttt                                                                    <  1:678754/42‑1 (MQ=255)
                            gATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGAttt                                                                    <  1:3345850/42‑1 (MQ=255)
                            gATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGAttt                                                                    <  1:4370902/42‑1 (MQ=255)
                              ttCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCagag                                     <  1:4953905/71‑1 (MQ=255)
                              ttCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATCTTACCGGCGCAGAGATGAAAAAACTGCCagag                                     <  1:1830468/71‑1 (MQ=255)
                                cGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATCTTACCGGCGCAGAGATGAAAAAACTGCCagag                                     <  1:2718676/69‑1 (MQ=255)
                                   cAGCCTGACGGTACTATCGCCTTTGTGTTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATAc                                <  1:2471712/71‑1 (MQ=255)
                                    aGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATAcc                               >  1:2114527/1‑71 (MQ=255)
                                       cTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCg                            >  1:1672584/1‑71 (MQ=255)
                                          aCGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCa                                                          <  1:4282054/38‑1 (MQ=255)
                                          aCGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGt                          <  1:3274247/70‑1 (MQ=255)
                                          aCGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGt                          <  1:657983/70‑1 (MQ=255)
                                            ggTACTATCGCCTTTGTGGTCGATCTTACCGGCGCAGAGATGaaa                                                 >  1:3638247/1‑45 (MQ=255)
                                            ggTACTATCGCCTTTGTGGTCGATCTTACCGGCGCAGAGATGaaa                                                 >  1:4588419/1‑45 (MQ=255)
                                                 tATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGc                                          <  1:2185853/47‑1 (MQ=255)
                                                 tATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGc                                          <  1:4761383/47‑1 (MQ=255)
                                                  aTCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGc                                          <  1:2825155/46‑1 (MQ=255)
                                                   tCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGaaaaaa                                              >  1:4897951/1‑41 (MQ=255)
                                                   tCGCCTTTGTGGTCGATCTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGTCACAGCGCa                 >  1:868746/1‑70 (MQ=255)
                                                     gCCTTTGTGGTCGATTTTACCGGCGCAGAGATGaaaaaa                                              >  1:3327073/1‑39 (MQ=255)
                                                     gCCTTTGTGGTCGATTTTACCGGCGCAGAGATGaaaaa                                               >  1:1522742/1‑38 (MQ=255)
                                                     gCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGTCACAGCGCAAAc              >  1:1676507/1‑71 (MQ=255)
                                                       cTTTGTGGTCGATCTTACCGGCGCAGAGATGAAAAAACTGCCagag                                     <  1:145093/46‑1 (MQ=255)
                                                         ttGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGTCACAGCGCAAACCAg           >  1:3557469/1‑70 (MQ=255)
                                                          tgtgGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGTCACAGCGCAAACCAg           <  1:4960193/69‑1 (MQ=255)
                                                            tgGTCGATCTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGTCACAGCGCAAACCAGCAtt       <  1:2161783/71‑1 (MQ=255)
                                                              gTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGTCACAGCGCAAACCAGCATTgg     >  1:515034/1‑71 (MQ=255)
                                                                cGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGTCACAGCGCAAACCAGCATTg      >  1:1835259/1‑68 (MQ=255)
                                                                 gATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGTCACAGCGCAAACCAGCATTGGTGa  >  1:236554/1‑71 (MQ=255)
                                                                 gATCTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGt                          <  1:2236997/47‑1 (MQ=255)
                                                                    |                                                                   
CAACGGGGGATGTGAAGCAGTCGAACCTGATTCGCCAGCCTGACGGTACTATCGCCTTTGTGGTCGATTTTACCGGCGCAGAGATGAAAAAACTGCCAGAGGATACCCCGGTCACAGCGCAAACCAGCATTGGTGA  >  W3110S.gb/1112134‑1112269

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: