Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 267,502 G→T 7.7% P228H (CCT→CAT)  ykfA ← predicted GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb267,5020GT 7.7% 205.2 / 2.4 64P228H (CCT→CAT) ykfApredicted GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base T (3/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

GAGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCC  >  W3110S.gb/267432‑267567
                                                                      |                                                                 
gAGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGACTACCGTCGTGCGAaaggaag                                                                   <  1:1605796/70‑1 (MQ=255)
     cGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATg                                                               <  1:2909200/70‑1 (MQ=255)
     cGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATg                                                               >  1:2683482/1‑70 (MQ=255)
      ggCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATgg                                                              <  1:2280587/70‑1 (MQ=255)
         cccGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGtt                                                          <  1:2443442/71‑1 (MQ=255)
         cccGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGt                                                           >  1:1017009/1‑70 (MQ=255)
           cGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGc                                                        <  1:1065632/71‑1 (MQ=255)
            gACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAATGATGGAGTTGCg                                                       <  1:2914986/71‑1 (MQ=255)
             aCGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGt                                                           >  1:1814010/1‑66 (MQ=255)
             aCGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGa                                                      >  1:698285/1‑71 (MQ=255)
             aCGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGa                                                      >  1:646287/1‑71 (MQ=255)
               ggTTTCACAAAAATCGCTACGAGCCTGTTCCCGGACTACAGTCGTGCGAAAGGAAGGATGGAGTTGCg                                                       >  1:182826/1‑68 (MQ=255)
                gTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGacac                                                   <  1:3632019/71‑1 (MQ=255)
                  ttCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGAcacca                                                 >  1:4053378/1‑71 (MQ=255)
                   tCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATg                                                               <  1:515549/56‑1 (MQ=255)
                   tCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGa                                                      >  1:2690516/1‑65 (MQ=255)
                   tCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGa                                                      >  1:405412/1‑65 (MQ=255)
                   tCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGAcaccac                                                >  1:3527152/1‑71 (MQ=255)
                   tCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGAcaccac                                                >  1:942932/1‑71 (MQ=255)
                    cACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGAcaccac                                                >  1:3635463/1‑70 (MQ=255)
                     aCCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGa                                                             <  1:1346177/56‑1 (MQ=255)
                      ccAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCg                                              >  1:1494102/1‑70 (MQ=255)
                       cAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAg                                                            <  1:1495723/55‑1 (MQ=255)
                       cAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCggg                                            >  1:2768509/1‑71 (MQ=255)
                        aaaaTCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAaacgaagg                                                                  >  1:3104077/1‑48 (MQ=255)
                         aaaTCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAATGATGGAg                                                            >  1:2361905/1‑53 (MQ=255)
                         aaaTCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAATGATGGAGTTGCGACACCACCGGGCt                                          >  1:1485415/1‑71 (MQ=255)
                         aaaTCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCt                                          >  1:2624306/1‑71 (MQ=255)
                         aaaTCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCt                                          <  1:3241527/71‑1 (MQ=255)
                         aaaTCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCt                                          >  1:3408327/1‑71 (MQ=255)
                          aaTCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCtg                                         <  1:3968995/71‑1 (MQ=255)
                           aTCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCggg                                            >  1:2028930/1‑67 (MQ=255)
                           aTCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCtgg                                        <  1:3220723/71‑1 (MQ=255)
                               cTACGAGCCTGTTCCCGGACTACCGTCGTGCGAaaggaagg                                                                  <  1:3867329/41‑1 (MQ=255)
                               cTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATg                                                               <  1:2725677/44‑1 (MQ=255)
                                  cGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCg                                                       >  1:506166/1‑49 (MQ=255)
                                  cGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCt                                  <  1:417358/70‑1 (MQ=255)
                                    aGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCAc                               <  1:3433648/71‑1 (MQ=255)
                                        tGTTCCCGGACTACCGTCGTGCGAAAGGAAGTATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTgg                           >  1:758466/1‑71 (MQ=255)
                                         gTTCCCGGACTACCGTCGTGCGAAAGGAATGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCAc                               >  1:2874163/1‑66 (MQ=255)
                                         gTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCtggtg                                      <  1:989677/59‑1 (MQ=255)
                                         gTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCtggtg                                      <  1:2338295/59‑1 (MQ=255)
                                          ttCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTggca                         <  1:3486187/71‑1 (MQ=255)
                                          ttCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTggc                          <  1:3351901/70‑1 (MQ=255)
                                          ttCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTggc                          <  1:2000545/70‑1 (MQ=255)
                                          ttCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTggc                          <  1:2879343/70‑1 (MQ=255)
                                           tCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCt                                  >  1:690878/1‑61 (MQ=255)
                                             ccGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTggcag                        <  1:471849/69‑1 (MQ=255)
                                              cGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTggcaggca                     >  1:2739415/1‑71 (MQ=255)
                                               ggACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTggcaggca                     >  1:126650/1‑70 (MQ=255)
                                               ggACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCAc                    >  1:4205238/1‑71 (MQ=255)
                                                gACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTggc                          <  1:2119007/64‑1 (MQ=255)
                                                gACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACt                   <  1:295788/71‑1 (MQ=255)
                                                 aCTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACtt                  <  1:300689/71‑1 (MQ=255)
                                                  cTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTa                 <  1:1749605/71‑1 (MQ=255)
                                                   tACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACg                              >  1:3739390/1‑57 (MQ=255)
                                                   tACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACtt                  <  1:2627816/69‑1 (MQ=255)
                                                     ccGTCGTGCGAAAGGAATGATGGAGTTGCGACACCACCGGGCTGGTGGCCt                                  <  1:1874479/51‑1 (MQ=255)
                                                     ccGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCAc                    >  1:2470061/1‑65 (MQ=255)
                                                     ccGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTaa                <  1:1202146/69‑1 (MQ=255)
                                                     ccGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAat               <  1:2829296/70‑1 (MQ=255)
                                                     ccGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAat               <  1:3657460/70‑1 (MQ=255)
                                                        tcgtGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCtggt                                       >  1:77676/1‑43 (MQ=255)
                                                          gtGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCtg                                         >  1:4017324/1‑39 (MQ=255)
                                                             cGAAAGGAAGGATGGAGTTGCGACACCACCGGGCtg                                         <  1:963877/36‑1 (MQ=255)
                                                              gAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCg     >  1:3280822/1‑71 (MQ=255)
                                                               aaaGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTggc                          >  1:2527744/1‑49 (MQ=255)
                                                                aaggaatgATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGc   >  1:2401154/1‑71 (MQ=255)
                                                                 aggaaggaTGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGcc  <  1:2763402/71‑1 (MQ=255)
                                                                 aggaaggaTGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGcc  <  1:3900597/71‑1 (MQ=255)
                                                                      |                                                                 
GAGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCC  >  W3110S.gb/267432‑267567

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: