Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,246,990 C→T 26.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,9900CT26.7% 48.6 / 18.2 30R40H (CGT→CAT) malGmaltose transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (8/14);  new base T (2/6);  total (10/20)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

AACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAAGTAGCAGTAGC  >  W3110S.gb/4246921‑4247055
                                                                     |                                                                 
aaCGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACg                                                                  >  1:468587/1‑71 (MQ=255)
      aGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCga                                                            <  1:2518867/71‑1 (MQ=255)
        tttCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCgaga                                                          <  1:1265696/71‑1 (MQ=255)
          tCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGAt                                                         <  1:1209827/70‑1 (MQ=255)
          tCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGAt                                                         <  1:2321284/70‑1 (MQ=255)
           ccAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCgaga                                                          >  1:2054819/1‑68 (MQ=255)
             aGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGcg                                                     <  1:384882/71‑1 (MQ=255)
                    ccAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCGAGATAGCGACGAcc                                               <  1:1037010/70‑1 (MQ=255)
                         agaTTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATcag                                          <  1:1162956/70‑1 (MQ=255)
                           aTTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGc                                                      <  1:208865/56‑1 (MQ=255)
                             ttGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCg                                     >  1:982866/1‑71 (MQ=255)
                                 tCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGaa                                 >  1:1392053/1‑71 (MQ=255)
                                          aGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGc                        <  1:247197/71‑1 (MQ=255)
                                          aGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGcga                                                    >  1:1230080/1‑43 (MQ=255)
                                                 cGGTCGCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATa                 <  1:1744130/71‑1 (MQ=255)
                                                        aaaGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagt           >  1:190692/1‑70 (MQ=255)
                                                        aaaGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagt           >  1:2474882/1‑70 (MQ=255)
                                                        aaaGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAag            >  1:1264952/1‑69 (MQ=255)
                                                         aaGTTTCCCTGATGCAGCAAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAag            <  1:1243608/68‑1 (MQ=255)
                                                          aGTTTCCCTGATGCAGCGAGATAGCGACGACCATc                                            <  1:729341/35‑1 (MQ=255)
                                                          aGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagtagc        <  1:1962800/71‑1 (MQ=255)
                                                          aGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagtag         >  1:1653771/1‑70 (MQ=255)
                                                          aGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagtagc        <  1:2013413/71‑1 (MQ=255)
                                                           gTTTCCCTGACGCAGCGAGATAGCGACGACCATcagcag                                       <  1:2232245/39‑1 (MQ=255)
                                                           gTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagtag         >  1:1792545/1‑69 (MQ=255)
                                                            tttCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACAt                              <  1:2275005/47‑1 (MQ=255)
                                                            tttCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATaaaa              <  1:2016869/63‑1 (MQ=255)
                                                              tCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACAt                              <  1:157237/45‑1 (MQ=255)
                                                                ccTGACGCAGCGAGATAGCGACGACATCAGCAGCGGGAACATAATCGCTGCGATAAAAagtagcagtagc  <  1:2004544/70‑1 (MQ=255)
                                                                 cTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATaaa               <  1:1479973/57‑1 (MQ=255)
                                                                     |                                                                 
AACGCCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAAGTAGCAGTAGC  >  W3110S.gb/4246921‑4247055

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: