Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 234,965:1 +C 6.0% coding (150/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,9651.C 6.0% 250.5 / 13.7 67coding (150/723 nt)yafSpredicted S‑adenosyl‑L‑methionine‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (27/36);  new base C (2/2);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAAT  >  W3110S.gb/234897‑235033
                                                                     |                                                                    
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:2630458/66‑2 (MQ=25)
          gcCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAATTGCCAATCTTAAgca                                                               <  1:4267613/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)
                  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:3025879/49‑1 (MQ=255)
                                tAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGaa                                                          <  1:3231021/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:2495726/70‑1 (MQ=255)
                                  aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGt                                   <  1:7321672/70‑1 (MQ=255)
                                  aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGt                                   <  1:2418330/70‑1 (MQ=255)
                                  aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGt                                   <  1:2444846/70‑1 (MQ=255)
                                  aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGt                                   <  1:672675/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)
                                                            aTTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCa        <  1:592170/71‑1 (MQ=255)
                                                             ttGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCaa       <  1:5486/71‑1 (MQ=255)
                                                             ttGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCa        <  1:2732619/70‑1 (MQ=255)
                                                              tGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCaa       >  1:7613946/1‑70 (MQ=255)
                                                              tGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAg      >  1:6830380/1‑71 (MQ=255)
                                                                gCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCaa                         >  1:352078/1‑50 (MQ=255)
                                                                gCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGa    <  1:351047/71‑1 (MQ=255)
                                                                gCAATTTAAGCGCAGAAATCAATT‑‑GAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAAt  <  1:5767009/71‑1 (MQ=255)
                                                                  aaTTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCtca                             <  1:1246030/44‑1 (MQ=255)
                                                                    tttAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCtcatca                          <  1:2628178/45‑1 (MQ=255)
                                                                    tttAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCtcatca                          <  1:959554/45‑1 (MQ=255)
                                                                    tttAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGttt               <  1:6221959/56‑1 (MQ=255)
                                                                     |                                                                    
GCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAAT  >  W3110S.gb/234897‑235033

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: