Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,085,470 G→A 21.0% R622C (CGT→TGT)  malT ← DNA‑binding transcriptional activator, maltotriose‑ATP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,085,4700GA21.0% 77.4 / 16.1 38R622C (CGT→TGT) malTDNA‑binding transcriptional activator, maltotriose‑ATP‑binding
Reads supporting (aligned to +/- strand):  ref base G (14/16);  new base A (5/3);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

TTGCCATTCCCCAGCAGGTTTTCCAGACGGTTCAGCTGGCTACGGGCGTTATCTAAATCACCACGGGCCAGCGAGCATTGAATCAACATTGCCAGGCACTGAAGCTGTTGCTGTGGCTGATAAGACGA  >  W3110S.gb/4085406‑4085533
                                                                |                                                               
ttGCCATTCCCCAGCAGGTTTTCCAGACGGTTCAGCTGGCTACGGGCGTTATCTAAATCACCACGGGCCAg                                                           >  1:420417/1‑71 (MQ=255)
ttGCCATTCCCCAGCAGGTTTTCCAGACGGTTCAGCTGGCTACGGGCGTTATCTAAATCACCACGGGCCAg                                                           >  1:1901708/1‑71 (MQ=255)
 tGCCATTCCCCAGCAGGTTTTCCAGACGGTTCAGCTGGCTACGGGCGTTATCTAAATCACCACGGGCCAGc                                                          >  1:1649531/1‑71 (MQ=255)
        ccccAGCAGGTTTTCCAGACGGTTCAGCTGGCTACGGGCGTTATCTAAATCACCACGGGCCAGCGAGCAtt                                                   >  1:741255/1‑71 (MQ=255)
           cagcagGTTTTCCAGACGGTTCAGCTGGCTACGGGCGTTATCTAAATCACCACAGGCCAGCGAGCATTGaa                                                >  1:2524538/1‑71 (MQ=255)
               agGTTTTCCAGACGGTTCAGCTGGCTACGGGCGTTATCTAAATCACCACGGGCCAGCGAGCATTGAATCaa                                            <  1:678942/71‑1 (MQ=255)
                   tttCCAGACGGTTCAGCTGGCTACGGGCGTTATCTAAATCACCACGGGCCAGCGAGCATTGAATCAACAtt                                        >  1:2920596/1‑71 (MQ=255)
                    ttCCAGACGGTTCAGCTGGCTACGGGCGTTATCTAAATCACCACGGGCCAGCGAGCATTGAATCAACAtt                                        <  1:289151/70‑1 (MQ=255)
                        aGACGGTTCAGCTGGCTACGGGCGTTATCTAAATCACCACAGGCCAGCGAGCATTGAATCAACATTg                                       <  1:292329/67‑1 (MQ=255)
                         gACGGTTCAGCTGGCTACGGGCGTTATCTAAATCACCACGGGCCAGCGAGCATTGAATCAACATTGCCAgg                                  <  1:2529391/71‑1 (MQ=255)
                           cGGTTCAGCTGGCTACGGGCGTTATCTAAATCACCACGGGCCAGCGAGCATTGAATCAACATTGCCAGGCa                                >  1:4044584/1‑71 (MQ=255)
                            ggTTCAGCTGGCTACGGGCGTTATCTAAATCACCACAGGCCAGCGAGCATTGAATCAACATTGCCAGGCAc                               >  1:4354946/1‑71 (MQ=255)
                              ttCAGCTGGCTACGGGCGTTATCTAAATCACCACGGGCCAGCGAGCATTGAATCAACATTGCCAGGCACt                              >  1:351499/1‑70 (MQ=255)
                                      gCTACGGGCGTTATCTAAATCACCACGGGCCAGCGAGGATTGAATCAACATTGCAGGCACTGAAGc                         >  1:4147662/1‑66 (MQ=255)
                                      gCTACGGGCGTTATCTAAATCACCACGGGCCAGCGAGCATTGAATCAACATTGCCAGGCACTGAAGCTGtt                     <  1:3941653/71‑1 (MQ=255)
                                                ttATCTAAATCACCACAGGCCAGCGAGCATTGAATCAACATTGCCAGGCACTGAAGCTGTTGCtgtg               <  1:3622267/67‑1 (MQ=255)
                                                  aTCTAAATCATCACGGGCCAGCGAGCATTGAATCAACATTGc                                      >  1:1176686/1‑42 (MQ=255)
                                                  aTCTAAATCACCACGGGCCAGCGAGCATTGAATCAACATTGCCAGGCACTGAAGCTGTTGCTGTGGCTGa          <  1:2972855/70‑1 (MQ=255)
                                                  aTCTAAATCACCACGGGCCAGCGAGCATTGAATCAACATTGCCAGGCACTGAAGCTGTTGCTGTGGCTGa          <  1:647416/70‑1 (MQ=255)
                                                  aTCTAAATCACCACGGGCCAGCGAGCATTGAATCAACATTGCCAGGCACTGAAGCTGTTGCTGTGGCTGa          >  1:3889072/1‑70 (MQ=255)
                                                  aTCTAAATCACCACGGGCCAGCGAGCATTGAATCAACATTGCCAGGCACTGAAGCTGTTGCTGTGGCTGa          <  1:79527/70‑1 (MQ=255)
                                                  aTCTAAATCACCACGGGCCAGCGAGCATTGAATCAACATTGCCAGGCACTGAAGCTGTTGCTGTGGCTGa          >  1:1465571/1‑70 (MQ=255)
                                                   tCTAAATCACCACGGGCCAGCGAGCATTGAATCAACATTGc                                      <  1:1722705/41‑1 (MQ=255)
                                                    cTAAATCACCACGGGCCAGCGAGCATTGAATCAACATTGCCAGGCACTGAAGCt                        <  1:594645/54‑1 (MQ=255)
                                                    cTAAATCACCACAGGCCAGCGAGCATTGAATCAACATTGCCAGGCACTGAAGCTGt                      >  1:857497/1‑56 (MQ=255)
                                                     tAAATCACCACGGGCCAGCGAGCATTGAATCAACATTGCCAGGCACTGAAGCTGTTGc                   <  1:4393412/58‑1 (MQ=255)
                                                      aaaTCACCACGGGCCAGCGAGCATTGAATCAACATTGCCAGGCa                                >  1:505305/1‑44 (MQ=255)
                                                       aaTCACCACGGGCCAGCGAGCATTGAATCAACAtt                                        <  1:3997016/35‑1 (MQ=255)
                                                       aaTCACCACGGGCCAGCGAGCATTGAATCAACAtt                                        <  1:1223132/35‑1 (MQ=255)
                                                       aaTCACCACGGGCCAGCGAGCATTGAATCAACATTGCCAgg                                  <  1:835274/41‑1 (MQ=255)
                                                       aaTCACCACGGGCCAGCGAGCATTGAATCAACATTGCCAGGCACTGAAGCAGTTGCTGTGGCTGATAAgac    >  1:3340992/1‑71 (MQ=255)
                                                        aTCACCACGGGCCAGCGAGCATTGAATCAACATTGCCAgg                                  <  1:4261905/40‑1 (MQ=255)
                                                         tCACCACAGGCCAGCGAGCATTGAATCAACATTGCCAGGCACTGAAGCTGTTGc                   >  1:1542658/1‑54 (MQ=255)
                                                         tCACCACAGGCCAGCGAGCATTGAATCAACATTGCCAGGCACTGAAGCTGTTGCTGTGGCTGATAAgacga  >  1:1722932/1‑71 (MQ=255)
                                                           accacAGGCCAGCGAGCATTGAATCAACATTGCCAGGCACTGAAGCTGTTGCtgt                <  1:2261374/55‑1 (MQ=255)
                                                             cacGGGCCAGCGAGCATTGAATCAACATTGCCAGGCACTGAAGCTGTTGc                   <  1:3823194/50‑1 (MQ=255)
                                                             cacGGGCCAGCGAGCATTGAATCAACATTGCCAGGCACTGAAGCTGTTGc                   <  1:2708330/50‑1 (MQ=255)
                                                              acGGGCCAGCGAGCATTGAATCAACATTGCCAGGCACTGAAGCTGTTGCTGTGGCTGATAAga     >  1:2855042/1‑63 (MQ=255)
                                                                |                                                               
TTGCCATTCCCCAGCAGGTTTTCCAGACGGTTCAGCTGGCTACGGGCGTTATCTAAATCACCACGGGCCAGCGAGCATTGAATCAACATTGCCAGGCACTGAAGCTGTTGCTGTGGCTGATAAGACGA  >  W3110S.gb/4085406‑4085533

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: