Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,246,990 C→T 16.7% R40H (CGT→CAT)  malG ← maltose transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,246,9900CT16.7% 90.4 / 10.5 36R40H (CGT→CAT) malGmaltose transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (18/12);  new base T (3/3);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

AGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAAGTAGCAGTAGCA  >  W3110S.gb/4246927‑4247056
                                                               |                                                                  
aGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTATCCCTGACGCAGCga                                                             <  1:4024903/71‑1 (MQ=255)
  tttCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCgaga                                                           <  1:4178738/71‑1 (MQ=255)
    tCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCgaga                                                           >  1:644086/1‑69 (MQ=255)
    tCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGAt                                                          <  1:2088378/70‑1 (MQ=255)
          gATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGcgac                                                    >  1:1234767/1‑70 (MQ=255)
           aTCCCATGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGcgacga                                                  >  1:1048170/1‑71 (MQ=255)
                 ggAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCAt                                              >  1:3090201/1‑69 (MQ=255)
                       ttGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCg                                      >  1:1507000/1‑71 (MQ=255)
                        tGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCgg                                     >  1:1091960/1‑71 (MQ=255)
                        tGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCgg                                     >  1:531950/1‑71 (MQ=255)
                         gCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCg                                      >  1:4678859/1‑69 (MQ=255)
                         gCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCg                                      >  1:21589/1‑69 (MQ=255)
                             cGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATc                                             <  1:2314917/58‑1 (MQ=255)
                              gggATCAGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACAt                               <  1:1367596/71‑1 (MQ=255)
                                    aGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGc                         >  1:1422448/1‑71 (MQ=255)
                                          ccGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGa                    >  1:373395/1‑70 (MQ=255)
                                          ccGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGa                    >  1:2036188/1‑70 (MQ=255)
                                            ggTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATa                  >  1:4644129/1‑70 (MQ=255)
                                            ggTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATa                  >  1:4312835/1‑70 (MQ=255)
                                             gTCGCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATaaa                <  1:1894338/71‑1 (MQ=255)
                                              tCGCAAAGTTTCCCTGACGCAGCGAGATAGCgacgac                                                 <  1:2448880/37‑1 (MQ=255)
                                              tCGCAAAGTTTCCCTGACGCAGCGAGATAGCgacgac                                                 <  1:1501470/37‑1 (MQ=255)
                                                gCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATaaa                >  1:3150459/1‑68 (MQ=255)
                                                 cAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAag             >  1:2061951/1‑70 (MQ=255)
                                                   aaGTTTCCCTGATGCAGCGAGATAGCGACGACCATTAGCAGCGGGAACATAATCGCTGCGAGAAAAagtag          >  1:91318/1‑71 (MQ=255)
                                                   aaGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagtag          <  1:1853432/71‑1 (MQ=255)
                                                    aGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagtag          >  1:1851730/1‑70 (MQ=255)
                                                      tttCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagtag          <  1:2799199/68‑1 (MQ=255)
                                                        tCCCTGACGCAGGGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagtagca        >  1:4403414/1‑68 (MQ=255)
                                                        tCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagtagca        >  1:3588939/1‑68 (MQ=255)
                                                          ccTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagtagcagtag    <  1:802072/70‑1 (MQ=255)
                                                           cTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagtagcagtag    >  1:1208180/1‑69 (MQ=255)
                                                           cTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGAt                   <  1:4751476/54‑1 (MQ=255)
                                                           cTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAAgtagcagtagca  <  1:338187/71‑1 (MQ=255)
                                                              aCGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAag             <  1:3056667/57‑1 (MQ=255)
                                                              aCGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAag             <  1:1111706/57‑1 (MQ=255)
                                                               |                                                                  
AGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAAGTAGCAGTAGCA  >  W3110S.gb/4246927‑4247056

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: