Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,246,866 T→C 31.4% V81V (GTA→GTG malG ← maltose transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,246,8660TC31.4% 41.6 / 28.1 35V81V (GTA→GTGmalGmaltose transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (16/8);  new base C (2/9);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 1.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCACAATGCCAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGC  >  W3110S.gb/4246806‑4246925
                                                            |                                                           
cgcCACAATGCCAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATg                                                   >  1:503995/1‑71 (MQ=255)
cgcCACAATGCCAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATg                                                   >  1:3068927/1‑71 (MQ=255)
       aTGCCAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTgg                                            >  1:956366/1‑71 (MQ=255)
       aTGCCAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTgg                                            >  1:3540583/1‑71 (MQ=255)
       aTGCCAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTg                                             >  1:2154186/1‑70 (MQ=255)
          ccAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGt                                         <  1:963418/71‑1 (MQ=255)
          ccAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGt                                         <  1:133735/71‑1 (MQ=255)
          ccAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTGGCGt                                         <  1:981021/71‑1 (MQ=255)
          ccAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTGGCGt                                         <  1:3625228/71‑1 (MQ=255)
           cAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTGGCGt                                         <  1:59275/70‑1 (MQ=255)
               cgcgGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGTaa                                       >  1:518651/1‑68 (MQ=255)
               cgcgGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGTaa                                       >  1:2946116/1‑68 (MQ=255)
                gcgGAAATCCCGGGGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGTaa                                       >  1:2792957/1‑67 (MQ=255)
                gcgGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGTaa                                       >  1:988804/1‑67 (MQ=255)
                gcgGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGTaa                                       >  1:3208333/1‑67 (MQ=255)
                gcgGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGTaa                                       >  1:29138/1‑67 (MQ=255)
                gcgGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGTaa                                       >  1:2675657/1‑67 (MQ=255)
                gcgGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATg                                                   >  1:3694845/1‑55 (MQ=255)
                gcgGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTGGCGTAATGc                                    <  1:707298/70‑1 (MQ=255)
                         ccGGCGACCTTTACCAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTGGCGTAATGCGACCATCAGcc                         <  1:1268769/71‑1 (MQ=255)
                         ccGGCGACCTTTACCAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTGGCGTAATGCGACCATCAGcc                         <  1:1272959/71‑1 (MQ=255)
                          cGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGt                                         <  1:3687288/55‑1 (MQ=255)
                          cGGCGACCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTGGCGTAATGCGACCATCAGcc                         <  1:3073960/71‑1 (MQ=255)
                          cGGCGACCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTGGCGTAATGCGACCATCAGcc                         <  1:1641760/71‑1 (MQ=255)
                                    tACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCaa                  >  1:1089339/1‑68 (MQ=255)
                                    tACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCaa                  >  1:230362/1‑68 (MQ=255)
                                    tACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCaa                  >  1:240755/1‑68 (MQ=255)
                                    tACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCaa                  >  1:265486/1‑68 (MQ=255)
                                             ccacagccacagCAGTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCt       <  1:90267/70‑1 (MQ=255)
                                             ccacagccacagCAGTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCt       <  1:135844/70‑1 (MQ=255)
                                             ccacagccacagCAGTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCt       <  1:2736496/70‑1 (MQ=255)
                                             ccacagccacagCAGTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCt       <  1:985881/70‑1 (MQ=255)
                                             ccacagccacagCAGTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCt       <  1:2902493/70‑1 (MQ=255)
                                             ccacagccacagCAGCACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCt       <  1:2250917/70‑1 (MQ=255)
                                                  gccacagcAGCACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGc  >  1:2367612/1‑70 (MQ=255)
                                                            tACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGtcc                    <  1:393415/42‑3 (MQ=255)
                                                            tACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTc                    <  1:3164643/42‑1 (MQ=255)
                                                            tACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTaaa           <  1:3662169/51‑1 (MQ=255)
                                                            tACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTaaa           <  1:1375082/51‑1 (MQ=255)
                                                            |                                                           
CGCCACAATGCCAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGC  >  W3110S.gb/4246806‑4246925

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: