Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 234,956 Δ1 bp 5.4% coding (141/723 nt) yafS → predicted S‑adenosyl‑L‑methionine‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb234,9560G. 5.4% 243.4 / 4.9 74coding (141/723 nt)yafSpredicted S‑adenosyl‑L‑methionine‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (27/43);  new base . (2/2);  total (29/45)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

TATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGC  >  W3110S.gb/234888‑235024
                                                                    |                                                                     
tATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAg                                                                       <  1:7699491/69‑1 (MQ=255)
tATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGAt                                                                     >  1:7533452/1‑71 (MQ=255)
 aTCGCAAGGCGCTGTAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGAt                                                                     <  1:2180261/70‑1 (MQ=255)
 aTCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGAt                                                                     <  1:3499427/70‑1 (MQ=255)
  tCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGa                                                                      >  1:3428717/1‑68 (MQ=255)
  tCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGAtt                                                                    <  1:7362855/70‑1 (MQ=255)
         gcgcTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAttt                                                           <  1:1095747/71‑1 (MQ=255)
         gcgcTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAttt                                                           <  1:4198614/71‑1 (MQ=255)
         gcgcTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAttt                                                           <  1:197472/71‑1 (MQ=255)
          cgcTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAttt                                                           >  1:7646016/1‑70 (MQ=255)
            cTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAg                                                        <  1:4174467/71‑1 (MQ=255)
               gAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAgcg                                                      >  1:4546391/1‑70 (MQ=255)
                   gcCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGa                                                  >  1:2602033/1‑70 (MQ=255)
                   gcCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAATTGCCAATCTTAAgca                                                      <  1:4267613/66‑2 (MQ=25)
                   gcCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAATTGCCAATCTTAAgca                                                      <  1:2630458/66‑2 (MQ=25)
                     cAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTaa                                                         <  1:4102394/61‑1 (MQ=255)
                     cAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAATTGCCAATCTTAAgcaga                                                    >  1:1477441/1‑63 (MQ=25)
                      aGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAAt                                               >  1:159085/1‑70 (MQ=255)
                      aGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAAt                                               >  1:2175729/1‑70 (MQ=255)
                      aGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATc                                              <  1:6604795/71‑1 (MQ=255)
                      aGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATc                                              <  1:3815251/71‑1 (MQ=255)
                        cTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCaa                                            <  1:6084786/71‑1 (MQ=255)
                         tCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGAtt                                                                    <  1:7363820/47‑1 (MQ=255)
                         tCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGAtt                                                                    <  1:3599552/47‑1 (MQ=255)
                         tCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGAtt                                                                    <  1:5170112/47‑1 (MQ=255)
                          cAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGAtt                                                                    <  1:6707085/46‑1 (MQ=255)
                          cAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGAtt                                                                    <  1:2575983/46‑1 (MQ=255)
                           aaCCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAAtt                                          <  1:1257238/70‑1 (MQ=255)
                           aaCCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAATTGCCAATCTTAAg                                                        >  1:1039156/1‑56 (MQ=25)
                            aCCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGc                                        >  1:3403722/1‑71 (MQ=255)
                             cccGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGa                                                  >  1:3877471/1‑60 (MQ=255)
                               cGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAttt                                                           <  1:2182094/49‑1 (MQ=255)
                               cGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAttt                                                           <  1:1770795/49‑1 (MQ=255)
                                      cACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTgc                               >  1:1117493/1‑70 (MQ=255)
                                       aCTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTgcgc                             >  1:2845840/1‑71 (MQ=255)
                                       aCTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTgcg                              <  1:6573259/70‑1 (MQ=255)
                                        caaaaaTGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAgcg                                                      <  1:1556211/43‑1 (MQ=255)
                                        cTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAgcg                                                      <  1:981604/45‑1 (MQ=255)
                                        cTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGaa                                                 >  1:3315049/1‑50 (MQ=255)
                                        cTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGcgcg                            <  1:5806499/71‑1 (MQ=255)
                                         tAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGaa                                                 <  1:3231021/49‑1 (MQ=255)
                                         tAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGaa                                                 <  1:3025879/49‑1 (MQ=255)
                                         tAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAAt                                               <  1:2390648/51‑1 (MQ=255)
                                         tAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCgg                           >  1:1277035/1‑71 (MQ=255)
                                          aaaaTGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAAtt                                          >  1:5899163/1‑55 (MQ=255)
                                           aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCgg                           <  1:6920621/69‑1 (MQ=255)
                                           aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGtt                         >  1:2082765/1‑71 (MQ=255)
                                           aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGt                          <  1:7321672/70‑1 (MQ=255)
                                           aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGt                          <  1:672675/70‑1 (MQ=255)
                                           aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGt                          <  1:2418330/70‑1 (MQ=255)
                                           aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGt                          <  1:2444846/70‑1 (MQ=255)
                                           aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGt                          <  1:2495726/70‑1 (MQ=255)
                                           aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGAGAAGCGTGCGCGGtt                         >  1:5911194/1‑71 (MQ=255)
                                            aaTGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGttt                        <  1:3468013/71‑1 (MQ=255)
                                             aTGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTtc                       <  1:5544257/71‑1 (MQ=255)
                                              tGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTtct                      >  1:4208070/1‑71 (MQ=255)
                                                tATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTtctc                     <  1:3383255/70‑1 (MQ=255)
                                                  tGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCtcatc                  <  1:5336890/71‑1 (MQ=255)
                                                  tGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCtcat                   >  1:5426026/1‑70 (MQ=255)
                                                   ggTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTtctc                     >  1:5845911/1‑67 (MQ=255)
                                                      tttCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAAt                                               <  1:2160281/38‑1 (MQ=255)
                                                      tttCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGa                                      >  1:1773583/1‑47 (MQ=255)
                                                      tttCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGa                                      >  1:748859/1‑47 (MQ=255)
                                                        tCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTg             <  1:4484582/70‑1 (MQ=255)
                                                         cATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCGTCAAGTGaa           >  1:3154912/1‑71 (MQ=255)
                                                          aTCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGttt                        >  1:129530/1‑57 (MQ=255)
                                                          aTCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCtca                    >  1:4718575/1‑61 (MQ=255)
                                                           tCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGaa                                     <  1:6762619/43‑1 (MQ=255)
                                                           tCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGa                                      <  1:7375454/42‑1 (MQ=255)
                                                           tCTGATTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATg         >  1:5529879/1‑71 (MQ=255)
                                                               cTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGt        <  1:2626952/68‑1 (MQ=255)
                                                               cTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGt        <  1:1523903/68‑1 (MQ=255)
                                                               cTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTc     <  1:327880/71‑1 (MQ=255)
                                                                 tAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTg   >  1:7262553/1‑71 (MQ=255)
                                                                  aaGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTgc  >  1:35165/1‑71 (MQ=255)
                                                                    |                                                                     
TATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGC  >  W3110S.gb/234888‑235024

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: