Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,246,990 C→T 77.3% 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,9900CT77.3% 81.7 / 24.2 44R40H (CGT→CAT) malGmaltose transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (3/7);  new base T (16/18);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

CCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAAGTAGCAGTAGCA  >  W3110S.gb/4246925‑4247056
                                                                 |                                                                  
ccAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGc                                                               >  1:8796/1‑71 (MQ=255)
   gTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCgag                                                            <  1:2666897/71‑1 (MQ=255)
     ttCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCgaga                                                           >  1:539712/1‑70 (MQ=255)
      tCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCgaga                                                           >  1:1222983/1‑69 (MQ=255)
      tCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCgaga                                                           >  1:1226273/1‑69 (MQ=255)
      tCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCgaga                                                           >  1:2069457/1‑69 (MQ=255)
           tGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGGGAGATAGcgac                                                    >  1:262404/1‑71 (MQ=255)
             aTCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCGAGATAGcgacga                                                  >  1:821609/1‑71 (MQ=255)
              tCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCgacgac                                                 <  1:2142185/71‑1 (MQ=255)
                      gaTTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATcag                                           >  1:2095669/1‑69 (MQ=255)
                      gaTTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATcag                                           >  1:2513743/1‑69 (MQ=255)
                            cTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCggg                                    <  1:2672617/70‑1 (MQ=255)
                            cTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGa                                   <  1:954956/71‑1 (MQ=255)
                                gggATCAGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGaa                                  <  1:1903320/68‑1 (MQ=255)
                                  gATCAGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATaa                             <  1:1471228/71‑1 (MQ=255)
                                  gATCAGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATaa                             <  1:2451907/71‑1 (MQ=255)
                                   aTCAGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAAt                            >  1:2318130/1‑71 (MQ=255)
                                      aGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCg                          <  1:2299738/70‑1 (MQ=255)
                                      aGGCTGCCGGTCGCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCg                          <  1:22050/70‑1 (MQ=255)
                                        gCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAAt                            >  1:447134/1‑66 (MQ=255)
                                         cTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGc                      <  1:249472/71‑1 (MQ=255)
                                         cTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGc                      <  1:1879015/71‑1 (MQ=255)
                                                tCGCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATaaa                <  1:905750/70‑1 (MQ=255)
                                                tCGCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATaaa                <  1:137456/70‑1 (MQ=255)
                                                tCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGc                      <  1:563706/64‑1 (MQ=255)
                                                  gCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATaaa                >  1:1299127/1‑68 (MQ=255)
                                                  gCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATaaa                >  1:702959/1‑68 (MQ=255)
                                                  gCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATaaa                >  1:1598664/1‑68 (MQ=255)
                                                  gCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATaaa                >  1:1329317/1‑68 (MQ=255)
                                                  gCAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATaaa                >  1:1068595/1‑68 (MQ=255)
                                                   cAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATa                              <  1:246574/53‑1 (MQ=255)
                                                   cAAAGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATa                              <  1:306366/53‑1 (MQ=255)
                                                   cAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATaaa                <  1:1940843/67‑1 (MQ=255)
                                                   cAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATaaa                <  1:1376507/67‑1 (MQ=255)
                                                      aGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATaa                             >  1:1439755/1‑51 (MQ=255)
                                                      aGTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATaa                             >  1:1024404/1‑51 (MQ=255)
                                                       gTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGc                         <  1:435479/54‑1 (MQ=255)
                                                       gTTTCCCTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGc                         <  1:2058467/54‑1 (MQ=255)
                                                          tCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAagtagcag       >  1:2144951/1‑69 (MQ=255)
                                                             cTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAAgtagcagtagca  <  1:1721557/71‑1 (MQ=255)
                                                             cTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAAgtagcagtagca  <  1:2417358/71‑1 (MQ=255)
                                                             cTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAAgtagcagtagca  <  1:1024174/71‑1 (MQ=255)
                                                             cTGATGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAAgtagcagtagca  <  1:2187942/71‑1 (MQ=255)
                                                             cTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTTCGATAAAAAgtagcagtagca  <  1:2426953/71‑1 (MQ=255)
                                                                 |                                                                  
CCAGTTTCCAGTGATCCCAGGAGATTTGCTCCGGGATCAGGCTGCCGGTCGCAAAGTTTCCCTGACGCAGCGAGATAGCGACGACCATCAGCAGCGGGAACATAATCGCTGCGATAAAAAGTAGCAGTAGCA  >  W3110S.gb/4246925‑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: