Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 827,597 C→T 49.0% P26P (CCG→CCA ybhS ← predicted transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb827,5970CT49.0% ‑0.6 / 72.2 51P26P (CCG→CCAybhSpredicted transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (15/11);  new base T (18/7);  total (33/18)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTAATGCCGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGGATCG  >  W3110S.gb/827531‑827664
                                                                  |                                                                   
gTTAATGCCGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGAt                                                                 <  1:871257/71‑1 (MQ=255)
gTTAATGCCGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGAt                                                                 <  1:1268226/71‑1 (MQ=255)
gTTAATGCCGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGAt                                                                 <  1:1319836/71‑1 (MQ=255)
gTTAATGCCGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGAt                                                                 <  1:710354/71‑1 (MQ=255)
gTTAATGCCGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGAt                                                                 <  1:1564744/71‑1 (MQ=255)
 ttAATGTCGTAACCAAAAATAAACAGCAGTAGCAGGGGGATCACTACCGCAATCAGCCAGCTACTTGGAt                                                                 <  1:2668272/70‑1 (MQ=255)
 ttAATGCCGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGAt                                                                 <  1:2926961/70‑1 (MQ=255)
         gTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATCGCGAACGa                                                        <  1:2420056/71‑1 (MQ=255)
         gTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATCGCGAACGa                                                        <  1:2588572/71‑1 (MQ=255)
         gTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATCGCGAACGa                                                        <  1:213351/71‑1 (MQ=255)
         gTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATCGCGAACGa                                                        <  1:26208/71‑1 (MQ=255)
         gTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATCGCGAACGa                                                        <  1:3040674/71‑1 (MQ=255)
         gTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATCGCGAACGa                                                        <  1:3108108/71‑1 (MQ=255)
         gTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATCGCGAACGa                                                        <  1:1819208/71‑1 (MQ=255)
         gTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATCGCGAACGa                                                        <  1:655729/71‑1 (MQ=255)
         gTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATCGCGAACGa                                                        <  1:825814/71‑1 (MQ=255)
                    tAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATCGCGAACGATCTGCCGCGtc                                             <  1:3080105/71‑1 (MQ=255)
                    tAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATCGCGAACGATCTGCCGCGtc                                             <  1:1104803/71‑1 (MQ=255)
                                gcagcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGca                                  >  1:2224497/1‑70 (MQ=255)
                                gcagcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGca                                  >  1:89594/1‑70 (MQ=255)
                                gcagcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGca                                  >  1:134560/1‑70 (MQ=255)
                                gcagcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGca                                  >  1:1417297/1‑70 (MQ=255)
                                gcagcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGca                                  >  1:2096903/1‑70 (MQ=255)
                                gcagcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGca                                  >  1:3110600/1‑70 (MQ=255)
                                gcagcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGca                                  >  1:1893508/1‑70 (MQ=255)
                                               ccGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGa                  >  1:2262483/1‑71 (MQ=255)
                                               ccGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGa                  >  1:855953/1‑71 (MQ=255)
                                               ccGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGa                  >  1:461452/1‑71 (MQ=255)
                                               ccGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGa                  >  1:420417/1‑71 (MQ=255)
                                               ccGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGa                  >  1:2899563/1‑71 (MQ=255)
                                               ccGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGa                  >  1:2184559/1‑71 (MQ=255)
                                               ccGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGa                  >  1:2282367/1‑71 (MQ=255)
                                               ccGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGa                  >  1:2192617/1‑71 (MQ=255)
                                               ccGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGa                  >  1:263690/1‑71 (MQ=255)
                                               ccGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGa                  >  1:2718436/1‑71 (MQ=255)
                                               cagCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGa                  >  1:2921273/3‑71 (MQ=255)
                                                               aCTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGGATCg  >  1:570740/1‑71 (MQ=255)
                                                               aCTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGGATCg  >  1:1069540/1‑71 (MQ=255)
                                                               aCTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGGATCg  >  1:2422673/1‑71 (MQ=255)
                                                               aCTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGGATCg  >  1:790955/1‑71 (MQ=255)
                                                               aCTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGGATCg  >  1:2568789/1‑71 (MQ=255)
                                                               aCTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGGATCg  >  1:680203/1‑71 (MQ=255)
                                                               aCTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGGATCg  >  1:1534083/1‑71 (MQ=255)
                                                               aCTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGGATCg  >  1:1950298/1‑71 (MQ=255)
                                                               aCTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGGATCg  >  1:479964/1‑71 (MQ=255)
                                                               aCTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGGATCg  >  1:283457/1‑71 (MQ=255)
                                                               aCTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGGATCg  >  1:2074125/1‑71 (MQ=255)
                                                               aCTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGGATCg  >  1:419774/1‑71 (MQ=255)
                                                               aCTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGGATCg  >  1:2690806/1‑71 (MQ=255)
                                                               aCTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGGATCg  >  1:1933442/1‑71 (MQ=255)
                                                               aCTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCAAAGCGCCCGTACGCGACGCCAGGACAGGATCg  >  1:2093917/1‑71 (MQ=255)
                                                                  |                                                                   
GTTAATGCCGTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGGATCG  >  W3110S.gb/827531‑827664

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: