Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,249,681 A→C 15.3% F111V (TTC→GTC)  malE ← maltose transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,249,6810AC15.3% 77.8 / 17.1 46F111V (TTC→GTC) malEmaltose transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (20/19);  new base C (2/5);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.58e-01

GGTAAGCAATCAGCTTGCCGTTGTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCGTAGCCACCAAAGCGGTCGTG  >  W3110S.gb/4249615‑4249744
                                                                  |                                                               
ggTAAGCAATCAGCTTGCCGTTGTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCt                                                             <  1:3986661/71‑1 (MQ=255)
   aaGCAATCAGCTTGCCGTTGTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGACCGCttt                                                           <  1:3605787/70‑1 (MQ=255)
       aaTCAGCTTGCCGTTGTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCg                                                      <  1:3365244/71‑1 (MQ=255)
        aTCAGCTTGCCGTTGTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCt                                                             <  1:3749334/63‑1 (MQ=255)
        aTCAGCTTGCCGTTGTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCt                                                             <  1:2387690/63‑1 (MQ=255)
          cAGCTTGCCGTTGTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCttt                                                           >  1:2442645/1‑63 (MQ=255)
           aGCTTGCCGTTGTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGACCGCtt                                                            <  1:786528/61‑1 (MQ=255)
           aGCTTGCCGTTGTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCgg                                                     >  1:2989492/1‑68 (MQ=255)
           aGCTTGCCGTTGTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGt                                                  >  1:3307987/1‑71 (MQ=255)
           aGCTTGCCGTTGTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGt                                                  <  1:3046257/71‑1 (MQ=255)
           aGCTTGCCGTTGTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGt                                                  <  1:283391/71‑1 (MQ=255)
           aGCTTGCCGTTGTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGt                                                  >  1:1043622/1‑71 (MQ=255)
           aGCTTGCCGTTGTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGATTTGTCCGGGGt                                                  >  1:3362649/1‑71 (MQ=255)
           aGCTTGCCGTTGTAACGTACGGCATCCAAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGt                                                  >  1:2534876/1‑71 (MQ=255)
                gCCGTTGTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGAtt                                              <  1:2086104/70‑1 (MQ=255)
                 ccGTTGTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGACCGCTTTGTCCGGGGTGAttt                                             <  1:2518369/70‑1 (MQ=255)
                   gTTGTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAg                                          >  1:2507356/1‑71 (MQ=255)
                    ttGTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTc                                                        <  1:1521743/56‑1 (MQ=255)
                      gTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCATcc                                        <  1:1576283/70‑1 (MQ=255)
                      gTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGcc                                        <  1:4244461/70‑1 (MQ=255)
                       tAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCa                                       <  1:476636/70‑1 (MQ=255)
                         aCGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGACCGCTTTGTCCGGGGTGATTTCAGCCAAc                                     <  1:3862153/70‑1 (MQ=255)
                              cGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCttt                                                           >  1:4361963/1‑43 (MQ=255)
                              cGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTcc                                                       >  1:3764504/1‑47 (MQ=255)
                                gCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGAt                                               <  1:2629143/53‑1 (MQ=255)
                                    cccAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCgg                                                     >  1:3915867/1‑43 (MQ=255)
                                     ccAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGa                             >  1:786125/1‑66 (MQ=255)
                                      cAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGc                       >  1:4505284/1‑71 (MQ=255)
                                       aGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGAt                                               <  1:3540306/46‑1 (MQ=255)
                                        ggTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACa                                    <  1:1121212/56‑1 (MQ=255)
                                        ggTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCg                      >  1:595304/1‑70 (MQ=255)
                                        ggTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCGt                     <  1:1590992/71‑1 (MQ=255)
                                          tAAACGGATACAGCTTGTCCTGGACCGCTTTGTCCggg                                                    >  1:4583625/1‑38 (MQ=255)
                                           aaaCGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTc                                            >  1:3041084/1‑45 (MQ=255)
                                               ggATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCGTAGccac               >  1:1505502/1‑70 (MQ=255)
                                                  tACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAAc                                     <  1:3915642/45‑1 (MQ=255)
                                                  tACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCGTAGCCACCaaa           <  1:4517381/71‑1 (MQ=255)
                                                  tACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCGTAGCCACCaaa           <  1:1478433/71‑1 (MQ=255)
                                                  tACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCGTAGCCACCaa            >  1:936939/1‑70 (MQ=255)
                                                   aCAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCGTa                    <  1:3921809/61‑1 (MQ=255)
                                                    cAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAAc                                     >  1:2406762/1‑43 (MQ=255)
                                                    cAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACAg                                   >  1:2853195/1‑45 (MQ=255)
                                                      gCTTGTCCTGGACCGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCGTAGCCACCAAAGCgg       >  1:549229/1‑71 (MQ=255)
                                                       cTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCGTAGCCACCAAAGCGGt      >  1:2443277/1‑71 (MQ=255)
                                                           tCCTGGACCGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCGTAGCCACCAAAGCGGTCgtg  <  1:3145597/71‑1 (MQ=255)
                                                               ggAACGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCGTAGCCACCAAAGCg        >  1:746100/1‑61 (MQ=255)
                                                                  aCGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCg                      >  1:3812993/1‑44 (MQ=255)
                                                                  |                                                               
GGTAAGCAATCAGCTTGCCGTTGTAACGTACGGCATCCCAGGTAAACGGATACAGCTTGTCCTGGAACGCTTTGTCCGGGGTGATTTCAGCCAACAGGCCAGATTGAGCGTAGCCACCAAAGCGGTCGTG  >  W3110S.gb/4249615‑4249744

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: