Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 267,502 G→T 12.4% 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,5020GT12.4% 223.7 / 10.2 72P228H (CCT→CAT) ykfApredicted GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (24/38);  new base T (4/5);  total (29/43)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GAGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCAC  >  W3110S.gb/267432‑267572
                                                                      |                                                                      
gAGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAaaggaag                                                                        >  1:4504498/1‑71 (MQ=255)
gAGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAaaggaag                                                                        <  1:3230950/71‑1 (MQ=255)
 aGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAaaggaagg                                                                       <  1:1741564/71‑1 (MQ=255)
   tACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAATGAt                                                                     <  1:4371674/71‑1 (MQ=255)
     cGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATgg                                                                   <  1:727933/71‑1 (MQ=255)
      ggCCCCGACGGTTTCACCAAAATCGCTACGAGCCAGTTCCCGGACTACCGTCGTGCGAAAGGAatgatg                                                                    >  1:3296758/1‑69 (MQ=255)
       gCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAg                                                                 <  1:3405861/71‑1 (MQ=255)
         cccGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGtt                                                               <  1:3806413/71‑1 (MQ=255)
         cccGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGtt                                                               <  1:772842/71‑1 (MQ=255)
         cccGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGtt                                                               <  1:284716/71‑1 (MQ=255)
          ccGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGtt                                                               <  1:252281/70‑1 (MQ=255)
           cGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGc                                                             >  1:2441432/1‑71 (MQ=255)
           cGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGc                                                             <  1:536189/71‑1 (MQ=255)
            gACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGc                                                             >  1:1637414/1‑70 (MQ=255)
             aCGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGa                                                           >  1:4373833/1‑71 (MQ=255)
             aCGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGa                                                           >  1:1597666/1‑71 (MQ=255)
              cGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTg                                                              >  1:362370/1‑67 (MQ=255)
                 tttCACCAAAATCTCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGacac                                                        >  1:1763707/1‑70 (MQ=255)
                 tttCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGacac                                                        >  1:1202850/1‑70 (MQ=255)
                   tCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGAcaccac                                                     >  1:777257/1‑71 (MQ=255)
                   tCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGAcaccac                                                     >  1:1315008/1‑71 (MQ=255)
                   tCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAACGATg                                                                    >  1:4050276/1‑56 (MQ=255)
                     aCCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGacac                                                        <  1:629174/66‑1 (MQ=255)
                     aCCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGAcacca                                                      <  1:3448190/68‑1 (MQ=255)
                     aCCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACaccacc                                                    >  1:3087983/1‑70 (MQ=255)
                     aCCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCg                                                   <  1:4177810/71‑1 (MQ=255)
                      ccAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATg                                                                    <  1:3802589/53‑1 (MQ=255)
                       cAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCggg                                                 >  1:57445/1‑71 (MQ=255)
                       cAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCgg                                                  >  1:2177816/1‑70 (MQ=255)
                       cAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCgg                                                  >  1:3012892/1‑70 (MQ=255)
                       cAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCAACggg                                                 >  1:1688052/1‑71 (MQ=255)
                         aaaTCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAATGATGGAGTTGCGACACCACCGGGCt                                               >  1:1250712/1‑71 (MQ=255)
                          aaTCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCtg                                              <  1:4462824/71‑1 (MQ=255)
                           aTCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAGAGGAAGGATGGAGTTGCGACACCACCGGGCtgg                                             >  1:4079141/1‑71 (MQ=255)
                            tCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCtggt                                            <  1:604730/71‑1 (MQ=255)
                            tCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCtggt                                            <  1:3119421/71‑1 (MQ=255)
                              gCTACGAGCCTTTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCtggtgg                                          <  1:1118696/71‑1 (MQ=255)
                                tACGAGCCTGTTCCCGGACTACCGTCGTGCGAAaggaagga                                                                      >  1:928767/1‑41 (MQ=255)
                                tACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGt                                                                >  1:2714390/1‑47 (MQ=255)
                                tACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGt                                                                >  1:3914524/1‑47 (MQ=255)
                                 aCGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCt                                       >  1:1451723/1‑71 (MQ=255)
                                  cGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAATGATGGAGTTGCGACACCACCGGGCtggt                                            >  1:2512992/1‑65 (MQ=255)
                                  cGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAATGATGGAGTTGCGACACCACCGGGCTGGTGGCCt                                       >  1:4128932/1‑70 (MQ=255)
                                  cGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCt                                       >  1:4246320/1‑70 (MQ=255)
                                  cGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCt                                       <  1:20320/70‑1 (MQ=255)
                                   gAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACaccacc                                                    <  1:3058817/56‑1 (MQ=255)
                                     gCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACg                                   >  1:4507633/1‑71 (MQ=255)
                                      ccTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGc                                         <  1:1607888/64‑1 (MQ=255)
                                      ccTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGt                                  >  1:1381316/1‑71 (MQ=255)
                                       cTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTg                                 <  1:510299/71‑1 (MQ=255)
                                         gTTCCCGGACTACCGTCGTGCGAAAGGAATGATGGAGtt                                                               <  1:1713897/39‑1 (MQ=255)
                                         gTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCg                                                            >  1:1149417/1‑42 (MQ=255)
                                          ttCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTggc                               <  1:3973818/70‑1 (MQ=255)
                                          ttCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTgg                                <  1:2255104/69‑1 (MQ=255)
                                          ttCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTgg                                <  1:691478/69‑1 (MQ=255)
                                            cccGGACTACCGTCGTGCGAAAGGAATGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTggcag                             <  1:3353535/70‑1 (MQ=255)
                                            cccGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCtggtg                                           >  1:2838516/1‑56 (MQ=255)
                                              cGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTggcaggca                          <  1:1308812/71‑1 (MQ=255)
                                                gACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTggc                               <  1:2015575/64‑1 (MQ=255)
                                                gACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCAc                         >  1:3632771/1‑70 (MQ=255)
                                                 aCTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACtt                       <  1:3777657/71‑1 (MQ=255)
                                                  cTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTa                      <  1:2037990/71‑1 (MQ=255)
                                                   tACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCtgg                                             <  1:4430126/47‑1 (MQ=255)
                                                    aCCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTggca                              <  1:231681/61‑1 (MQ=255)
                                                     ccGTCGTGCGAAAGGAATGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTaa                     >  1:56073/1‑69 (MQ=255)
                                                     ccGTCGTGCGAAAGGAATGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTaa                     <  1:4372228/69‑1 (MQ=255)
                                                     ccGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTgg                                <  1:4444906/58‑1 (MQ=255)
                                                     ccGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTaa                     <  1:3381265/69‑1 (MQ=255)
                                                      cgtcgtGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAatca                  <  1:665051/71‑1 (MQ=255)
                                                      cgtcgtGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAatca                  <  1:4424954/71‑1 (MQ=255)
                                                       gtcgtGCGAAAGGAAGGATGGAGTTGCGACACCACCgg                                                  <  1:3113211/38‑1 (MQ=255)
                                                        tcgtGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTggcaggc                           <  1:2817057/60‑1 (MQ=255)
                                                          gtGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCg              <  1:344593/71‑1 (MQ=255)
                                                             cGAAAGGAATGATGGAGTTGCGACACCACCGGGCtg                                              <  1:4535219/36‑1 (MQ=255)
                                                              gAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAat                    <  1:2414326/61‑1 (MQ=255)
                                                              gAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCg          <  1:3781760/71‑1 (MQ=255)
                                                               aaaGGAAGGATGGAGTTGCGACACCACCGGGCtggt                                            <  1:4121698/36‑1 (MQ=255)
                                                                 aggaatgaTGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGc        >  1:2995573/1‑70 (MQ=255)
                                                                 aggaaggaTGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTg                                 >  1:3204133/1‑45 (MQ=255)
                                                                 aggaaggaTGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGcc       <  1:3474773/71‑1 (MQ=255)
                                                                 aggaaggaTGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGcc       <  1:1654165/71‑1 (MQ=255)
                                                                 aggaaggaTGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGc        >  1:383588/1‑70 (MQ=255)
                                                                 aggaaggaTGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGc        >  1:3795748/1‑70 (MQ=255)
                                                                   gaatgaTGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTg                                 <  1:3999684/43‑1 (MQ=255)
                                                                      ggaTGGAGTTGCGACACCACCGGGCTGGTGGCCTCAc                                    <  1:1082633/37‑1 (MQ=255)
                                                                      ggaTGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCAc  <  1:946531/71‑1 (MQ=255)
                                                                      |                                                                      
GAGTACGGCCCCGACGGTTTCACCAAAATCGCTACGAGCCTGTTCCCGGACTACCGTCGTGCGAAAGGAAGGATGGAGTTGCGACACCACCGGGCTGGTGGCCTCACGTGGCAGGCACTTAATCATCCGCTCGGCCATCAC  >  W3110S.gb/267432‑267572

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: