Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,777,889 T→C 51.4% V359A (GTG→GCG)  glvC → arbutin specific enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,777,8890TC51.4% ‑1.0 / 51.7 37V359A (GTG→GCG) glvCarbutin specific enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base T (12/6);  new base C (12/7);  total (24/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

CGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATG  >  W3110S.gb/3777820‑3777947
                                                                     |                                                          
cGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTtggtgtggtg                                                           <  1:2791310/71‑1 (MQ=255)
    gCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGCGGGc                                                        <  1:377612/70‑1 (MQ=255)
     cTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCaa                                                      <  1:2329508/71‑1 (MQ=255)
     cTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGCGGGCaa                                                      <  1:2155463/71‑1 (MQ=255)
     cTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGCGGGCaa                                                      <  1:1458205/71‑1 (MQ=255)
     cTGTTTGCAGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCaa                                                      <  1:1212027/71‑1 (MQ=255)
        tttGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGCGGGCAAc                                                     <  1:1501892/69‑1 (MQ=255)
          tGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGCGGGCAACATgg                                                 >  1:1572034/1‑71 (MQ=255)
                acacGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAg                                           >  1:569726/1‑71 (MQ=255)
                 cacGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAgg                                          >  1:248079/1‑71 (MQ=255)
                 cacGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGCGGGCAACATGGGCGGAg                                           >  1:787906/1‑70 (MQ=255)
                            ggcggcCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTggg                                                         <  1:1734477/45‑1 (MQ=255)
                                  cTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCt                                       <  1:129269/57‑1 (MQ=255)
                                    cAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTa                       >  1:2805249/1‑71 (MQ=255)
                                     aaTGTCGACCGTAATGTATCTCTTTGGTGTGGCGGGCAACa                                                    >  1:288246/1‑41 (MQ=255)
                                     aaTGTCGACCGTAATGTATCTCTTTGGTGTGGCGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTAc                      >  1:1862018/1‑71 (MQ=255)
                                     aaTGTCGACCGTAATGTATCTCTTTGGTGTGGCGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTAc                      >  1:2106959/1‑71 (MQ=255)
                                      aTGTCGACCGTAATGTATCTCTTTGGTGTGGCGGGCAACa                                                    <  1:1141843/40‑1 (MQ=255)
                                       tGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCt                                       <  1:2120667/52‑1 (MQ=255)
                                       tGTCGACCGTAATGTATCTCTTTGGTGTGGCGGGCAACATGGGCGGAGGTCTGAtt                                   <  1:2758937/56‑1 (MQ=255)
                                         tCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATggg                                                >  1:1929056/1‑41 (MQ=255)
                                         tCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGAtt                                   >  1:1540292/1‑54 (MQ=255)
                                         tCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGAtt                                   >  1:1583151/1‑54 (MQ=255)
                                          cGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCa                  >  1:2734796/1‑70 (MQ=255)
                                               gTAATGTATCTCTTTGGTGTGGCGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAAc              >  1:1675325/1‑69 (MQ=255)
                                                    gTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGt                           >  1:644548/1‑51 (MQ=255)
                                                    gTATCTCTTTGGTGTGGCGGGCAACATGGGCGGAGGTCTGATTGa                                 >  1:1019738/1‑45 (MQ=255)
                                                     tATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATcc       >  1:15876/1‑70 (MQ=255)
                                                     tATCTCTTTGGTGTGGCGGGCAACATGGGCGGAGGTc                                        >  1:2026327/1‑37 (MQ=255)
                                                     tATCTCTTTGGTGTGGCGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAAc              >  1:1568064/1‑63 (MQ=255)
                                                     tATCTCTTTGGTGTGGCGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATccc      >  1:43602/1‑71 (MQ=255)
                                                     tATCTCTTTGGTGTGGCGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATccc      >  1:1306599/1‑71 (MQ=255)
                                                      aTCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGa                                 >  1:1786578/1‑43 (MQ=255)
                                                       tctcTTTGGTGTGGTGGGCAACATGGGCGGAGGTCt                                       <  1:412151/36‑1 (MQ=255)
                                                       tctcTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGa    >  1:1668084/1‑71 (MQ=255)
                                                       tctcTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGa    >  1:1864140/1‑71 (MQ=255)
                                                         tctTTGGTGTGGCGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATg  >  1:514878/1‑71 (MQ=255)
                                                           tttGGTGTGGCGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGAt   <  1:2388684/68‑1 (MQ=255)
                                                                     |                                                          
CGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATG  >  W3110S.gb/3777820‑3777947

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