Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,130,879:1 +A 15.4% coding (275/1227 nt) yhfW → predicted mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,130,8791.A15.4% 115.1 / 21.8 39coding (275/1227 nt)yhfWpredicted mutase
Reads supporting (aligned to +/- strand):  ref base . (12/21);  new base A (4/2);  total (16/23)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGTGGCAGAGCTGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGAGCAGGCATT  >  W3110S.gb/4130821‑4130948
                                                           |                                                                     
ggCGTGGCAGAGCTGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAGGCACGCGccc                                                           <  1:3968452/71‑1 (MQ=255)
  cGTGGCAGAGCTGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAgg                                                                    >  1:1206369/1‑60 (MQ=255)
  cGTGGCAGAGCTGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAgg                                                                    >  1:724556/1‑60 (MQ=255)
      gCAGAGCTGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAAGGCACtcgc                                                             >  1:1145049/1‑64 (MQ=255)
      gCAGAGCTGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAAGGCACGCGCCCGTTAc                                                      >  1:4066532/1‑71 (MQ=255)
       cAGAGCTGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGTTAcc                                                     >  1:2149468/1‑70 (MQ=255)
          agCTGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAAGGCACGCGCCCGt                                                         >  1:2843713/1‑64 (MQ=255)
             tGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAAGGCACGCTcc                                                            <  1:3115121/58‑1 (MQ=255)
             tGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCt                                               >  1:537777/1‑70 (MQ=255)
             tGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTg                                              >  1:494077/1‑71 (MQ=255)
              gCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGc                                             >  1:4049065/1‑71 (MQ=255)
               cAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCg                                            >  1:88312/1‑71 (MQ=255)
               cAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCg                                            <  1:3043488/71‑1 (MQ=255)
               cAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCg                                            <  1:918931/71‑1 (MQ=255)
                aaCATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAAGGCACGCGCCCGTTACCGCCGCTGCg                                            <  1:3241171/71‑1 (MQ=255)
                               aTACCTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGtt                                                        <  1:713512/43‑1 (MQ=255)
                               aTACCTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGtt                                                        <  1:93026/43‑1 (MQ=255)
                                  ccTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgtg                         <  1:231225/71‑1 (MQ=255)
                                  ccTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgtg                         <  1:722196/71‑1 (MQ=255)
                                  ccTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgtg                         <  1:2666874/71‑1 (MQ=255)
                                  ccTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgtg                         <  1:2917889/71‑1 (MQ=255)
                                  ccTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgtg                         <  1:808389/71‑1 (MQ=255)
                                  ccTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgt                          <  1:868748/70‑1 (MQ=255)
                                  ccTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgt                          <  1:939539/70‑1 (MQ=255)
                                  ccTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgt                          <  1:3406845/70‑1 (MQ=255)
                                  ccTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgt                          <  1:1272248/70‑1 (MQ=255)
                                   cTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgt                          <  1:4081856/69‑1 (MQ=255)
                                          gggCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCg                 <  1:1100273/71‑1 (MQ=255)
                                           ggCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCgt                <  1:2199153/71‑1 (MQ=255)
                                           ggCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCgt                <  1:1346311/71‑1 (MQ=255)
                                                 aGGAAATTTTAAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGtt             >  1:1737887/1‑69 (MQ=255)
                                                 aGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCgtgt              >  1:835542/1‑67 (MQ=255)
                                                 aGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGa           >  1:3250214/1‑70 (MQ=255)
                                                  ggAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGAg          <  1:2393158/70‑1 (MQ=255)
                                                   gAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCgt                <  1:61442/63‑1 (MQ=255)
                                                   gAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCgt                <  1:659176/63‑1 (MQ=255)
                                                    aaaTTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGcc                                      >  1:354948/1‑40 (MQ=255)
                                                    aaaTTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGAg          >  1:3882708/1‑68 (MQ=255)
                                                    aaaTTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGAg          >  1:1514783/1‑68 (MQ=255)
                                                     attttaAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGAGCAgg      >  1:1716367/3‑71 (MQ=255)
                                                     attttaAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGAGCAgg      >  1:1984491/3‑71 (MQ=255)
                                                       ttttAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCgtg               <  1:2617571/60‑1 (MQ=255)
                                                        tttAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGAGCAGGc     <  1:1648307/69‑1 (MQ=255)
                                                          tAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGAGCAGGCAtt  <  1:1791213/70‑1 (MQ=255)
                                                           |                                                                     
GGCGTGGCAGAGCTGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGAGCAGGCATT  >  W3110S.gb/4130821‑4130948

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: