Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,246,866 T→C 62.5% 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,8660TC62.5% 26.0 / 30.0 32V81V (GTA→GTGmalGmaltose transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (3/9);  new base C (7/13);  total (10/22)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTGGAGAGCGCCACAATGCCAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGT  >  W3110S.gb/4246798‑4246929
                                                                    |                                                               
gtggAGAGCGCCACAATGCCAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGt                                                                 <  1:3050145/69‑1 (MQ=255)
gtggAGAGCGCCACAATGCCAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGt                                                                 <  1:1144848/69‑1 (MQ=255)
       gcgcCACAATGCCAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGa                                                          <  1:1663174/69‑1 (MQ=255)
            acaATGCCAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGc                                                     <  1:2797258/69‑1 (MQ=255)
             caATGCCAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCg                                                    >  1:1916305/1‑69 (MQ=255)
                 gCCAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTgg                                                <  1:3466489/69‑1 (MQ=255)
                  ccAATCGCGGAAATCACGGCGACCTTTACCGAGTTCCACAGCCACAGCAGCAc                                                               >  1:2336486/1‑53 (MQ=255)
                    aaTCGCGGAAATCCCGTCGACCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTGGCGt                                             <  1:3460974/69‑1 (MQ=255)
                    aaTCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTGGCGt                                             <  1:1459198/69‑1 (MQ=255)
                    aaTCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTGGCGt                                             <  1:193830/69‑1 (MQ=255)
                    aaTCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTGGCGt                                             <  1:3505877/69‑1 (MQ=255)
                    aaTCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTGGCGt                                             <  1:2536630/69‑1 (MQ=255)
                    aaTCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTGGCGt                                             <  1:1209004/69‑1 (MQ=255)
                    aaTCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTGGCGt                                             <  1:2412659/69‑1 (MQ=255)
                      tCGCGGAAATCCCGGCGTCCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTGGCGt                                             <  1:3019666/67‑1 (MQ=255)
                        gcgGAAATCCCGGCGACCTTTACCGAGTTCCAAAGCCACAGCAGCACCggg                                                           <  1:202993/51‑1 (MQ=39)
                          ggAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGTAATGCg                                       <  1:3631422/69‑1 (MQ=255)
                             aaTCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGTAATGCGAcc                                    >  1:2609279/1‑69 (MQ=255)
                             aaTCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGTAATGCGAcc                                    >  1:2940842/1‑69 (MQ=255)
                              aTCCCGGCGCCCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCgg                                                   >  1:2067757/1‑53 (MQ=39)
                                cccGGCGACCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTGGCGTAATGCGAc                                     <  1:2224141/65‑1 (MQ=255)
                                      gACCTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTGGCGTAATGCGAc                                     <  1:2945840/59‑1 (MQ=255)
                                         cTTTACCGAGTTCCACAGCCACAGCAGCACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTc                        >  1:2120992/1‑69 (MQ=255)
                                         cTTTACCGAGTTCCACAGACACAGCAGCACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTc                        >  1:2440800/1‑69 (MQ=255)
                                           ttACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGt                                             <  1:2509798/46‑1 (MQ=255)
                                           ttACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGt                                             <  1:321719/46‑1 (MQ=255)
                                            tACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGt                                             >  1:1286479/1‑45 (MQ=255)
                                                       acagccacagcAGTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCt           <  1:1073746/68‑1 (MQ=255)
                                                       acagccacagcAGTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCt           <  1:635007/68‑1 (MQ=255)
                                                       acagccacagcAGTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCt           <  1:1337060/68‑1 (MQ=255)
                                                       acagccacagcAGTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCt           <  1:1069388/68‑1 (MQ=255)
                                                       acagccacagcAGCACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCt           <  1:1289400/68‑1 (MQ=255)
                                                            cacagcAGCACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAAc             >  1:2890668/1‑61 (MQ=255)
                                                            cacagcAGCACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAAc             >  1:1978619/1‑61 (MQ=255)
                                                               agcagcACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGt  <  1:2471240/69‑1 (MQ=255)
                                                                gcagcaCCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGt  >  1:727855/1‑68 (MQ=255)
                                                                    |                                                               
GTGGAGAGCGCCACAATGCCAATCGCGGAAATCCCGGCGACCTTTACCGAGTTCCACAGCCACAGCAGTACCGGGAATGGCGGTGGCGTAATGCGACCATCAGCCTGTTCAACGCTAAAACCTAACGCCAGT  >  W3110S.gb/4246798‑4246929

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: