Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 267,502 G→T 32.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,5020GT32.7% 106.7 / 37.7 58P228H (CCT→CAT) ykfApredicted GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (20/19);  new base T (7/12);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

GGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCA  >  W3110S.gb/267438‑267568
                                                                |                                                                  
ggCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAATGATGGa                                                              <  1:2248932/71‑1 (MQ=255)
ggCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATg                                                                >  1:4955495/1‑69 (MQ=255)
     cGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGc                                                         <  1:2715526/71‑1 (MQ=255)
      gACGGTTTCCCCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAATGATGGAg                                                             >  1:668444/1‑66 (MQ=255)
      gACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAATGGAAGGATGGAGTTGc                                                         >  1:891279/1‑70 (MQ=255)
      gACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGc                                                         >  1:3207248/1‑70 (MQ=255)
      gACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGc                                                         >  1:3977844/1‑70 (MQ=255)
      gACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGc                                                         <  1:4497244/70‑1 (MQ=255)
       aCGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAATGATGGAGTTGCGa                                                       <  1:4321338/71‑1 (MQ=255)
       aCGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGa                                                       >  1:1610013/1‑71 (MQ=255)
       aCGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGa                                                       >  1:3871045/1‑71 (MQ=255)
          gTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAATGATgg                                                               <  1:3997252/60‑1 (MQ=255)
          gTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAATGATgg                                                               <  1:565995/60‑1 (MQ=255)
          gTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAATGATGGAGTTGCGacac                                                    <  1:95475/71‑1 (MQ=255)
          gTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGaca                                                     <  1:720735/70‑1 (MQ=255)
           tttCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAATGATGGAGtt                                                           <  1:3204326/63‑1 (MQ=255)
           tttCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAATGATGGAGtt                                                           <  1:2627451/63‑1 (MQ=255)
            ttCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAaggaagga                                                                  >  1:4656710/1‑55 (MQ=255)
            ttCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAATGATGGAg                                                             >  1:2305251/1‑60 (MQ=255)
             tCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGa                                                       >  1:519822/1‑65 (MQ=255)
                 cAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCTAAAGGAAGGATGGAGTTGCGACACCACCgg                                              >  1:4007472/1‑70 (MQ=255)
                  aaaaTCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGa                                                              >  1:3261391/1‑53 (MQ=255)
                    aaTCGCTACGAGCCTGTTGCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGtt                                                           <  1:312514/54‑1 (MQ=255)
                     aTCGCTACGAGCCTGTTCCCGGACTTCCGTCGTGCGAAAGGAATGATGGAGTTGCGACACCACCGGGCtgg                                         >  1:2782622/1‑71 (MQ=255)
                     aTCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAATGATGGAGTTGCGACACCACCGGGCtgg                                         >  1:1835898/1‑71 (MQ=255)
                     aTCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCtgg                                         >  1:1907354/1‑71 (MQ=255)
                      tCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGc                                                         >  1:606428/1‑54 (MQ=255)
                      tCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGc                                                         >  1:3784186/1‑54 (MQ=255)
                         cTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGc                                     <  1:4484120/71‑1 (MQ=255)
                          tACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAATGATGGAGTTGCGCCACCACCGGGCTGGTGGc                                     >  1:3962296/1‑70 (MQ=255)
                          tACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAATGATGGAGTTGCGACACCACCGGGCTGGTGGc                                     >  1:4566097/1‑70 (MQ=255)
                           aCGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCt                                   >  1:4923040/1‑71 (MQ=255)
                           aCGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCt                                   >  1:3503090/1‑71 (MQ=255)
                            cGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACCCCACCGGGCTGGTGGCCt                                   >  1:164839/1‑70 (MQ=255)
                            cGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGcc                                    <  1:1260165/69‑1 (MQ=255)
                                ccTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGt                              <  1:1010747/71‑1 (MQ=255)
                                   gTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCAc                                <  1:2954201/66‑1 (MQ=255)
                                   gTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCAc                                <  1:1038634/66‑1 (MQ=255)
                                    ttCCCGGACTACCGTCGTGCGAAAGGAAGGATGGa                                                              <  1:3866301/35‑1 (MQ=255)
                                       ccGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTggcagg                        <  1:1759901/70‑1 (MQ=255)
                                         tgACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCAc                     >  1:4789710/2‑71 (MQ=255)
                                         ggACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCg                                               >  1:2314956/1‑45 (MQ=255)
                                           aCTACCGTTGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCtg                                          <  1:2806602/48‑1 (MQ=255)
                                            cTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACtt                   <  1:4255877/70‑1 (MQ=255)
                                                 gtcgtGCGAAAGGAATGATGGAGTTGCGACACCACCGGGCt                                           <  1:703408/41‑1 (MQ=255)
                                                 gtcgtGCGAAAGGAATGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAatcat             <  1:2790506/71‑1 (MQ=255)
                                                  tcgtGCGAAAGGAATGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTg                             <  1:136857/54‑1 (MQ=255)
                                                  tcgtGCGAAAGGAATGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAatcatc            <  1:1864672/71‑1 (MQ=255)
                                                  tcgtGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCtg                                          <  1:4413569/41‑1 (MQ=255)
                                                  tcgtGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCtg                                          <  1:4151192/41‑1 (MQ=255)
                                                  tcgtGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCtg                                          <  1:857627/41‑1 (MQ=255)
                                                  tcgtGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCtg                                          <  1:1359003/41‑1 (MQ=255)
                                                  tcgtGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAatcat             <  1:2043141/70‑1 (MQ=255)
                                                    gtGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCg          <  1:1594507/71‑1 (MQ=255)
                                                       cGAAAGGAATGATGGAGTTGCGACACCCCCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTc       >  1:45368/1‑71 (MQ=255)
                                                        gAAAGGAATGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAatc               <  1:2700178/62‑1 (MQ=255)
                                                         aaaGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCgg     >  1:2019602/1‑71 (MQ=255)
                                                         aaaGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCgg     >  1:4367334/1‑71 (MQ=255)
                                                         aaaGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCgg     >  1:3079119/1‑71 (MQ=255)
                                                          aaggaaggATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGc    <  1:2183750/71‑1 (MQ=255)
                                                           aggaaggaTGGAGTTGCGACACCACCGGGCTGGTGGCCTc                                  >  1:1380622/1‑40 (MQ=255)
                                                            ggaatgaTGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCa  <  1:1148116/71‑1 (MQ=255)
                                                            ggaatgaTGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCa  <  1:935399/71‑1 (MQ=255)
                                                                |                                                                  
GGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCA  >  W3110S.gb/267438‑267568

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: