Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 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
*AP009048234,9560G. 5.4% 243.3 / 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

TATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGC  >  AP009048/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)
         gcgcTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAt‑tt                                                           <  1:1095747/71‑1 (MQ=255)
         gcgcTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAt‑tt                                                           <  1:4198614/71‑1 (MQ=255)
         gcgcTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAt‑tt                                                           <  1:197472/71‑1 (MQ=255)
          cgcTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAt‑tt                                                           >  1:7646016/1‑70 (MQ=255)
            cTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAg                                                        <  1:4174467/71‑1 (MQ=255)
               gAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAgcg                                                      >  1:4546391/1‑70 (MQ=255)
                   gcCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGa                                                  >  1:2602033/1‑70 (MQ=255)
                   gcCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAA‑ATTGCCAATCTTAAgca                                                      <  1:4267613/66‑2 (MQ=25)
                   gcCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAA‑ATTGCCAATCTTAAgca                                                      <  1:2630458/66‑2 (MQ=25)
                     cAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTaa                                                         <  1:4102394/61‑1 (MQ=255)
                     cAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAA‑ATTGCCAATCTTAAgcaga                                                    >  1:1477441/1‑63 (MQ=25)
                      aGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAAt                                               >  1:159085/1‑70 (MQ=255)
                      aGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAAt                                               >  1:2175729/1‑70 (MQ=255)
                      aGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATc                                              <  1:6604795/71‑1 (MQ=255)
                      aGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATc                                              <  1:3815251/71‑1 (MQ=255)
                        cTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCaa                                            <  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)
                           aaCCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAAtt                                          <  1:1257238/70‑1 (MQ=255)
                           aaCCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAA‑ATTGCCAATCTTAAg                                                        >  1:1039156/1‑56 (MQ=25)
                            aCCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGc                                        >  1:3403722/1‑71 (MQ=255)
                             cccGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGa                                                  >  1:3877471/1‑60 (MQ=255)
                               cGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAt‑tt                                                           <  1:2182094/49‑1 (MQ=255)
                               cGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAt‑tt                                                           <  1:1770795/49‑1 (MQ=255)
                                      cACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTgc                               >  1:1117493/1‑70 (MQ=255)
                                       aCTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTgcgc                             >  1:2845840/1‑71 (MQ=255)
                                       aCTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTgcg                              <  1:6573259/70‑1 (MQ=255)
                                        caaaaaTGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAgcg                                                      <  1:1556211/43‑1 (MQ=255)
                                        cTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAgcg                                                      <  1:981604/45‑1 (MQ=255)
                                        cTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGaa                                                 >  1:3315049/1‑50 (MQ=255)
                                        cTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGcgcg                            <  1:5806499/71‑1 (MQ=255)
                                         tAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGaa                                                 <  1:3231021/49‑1 (MQ=255)
                                         tAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGaa                                                 <  1:3025879/49‑1 (MQ=255)
                                         tAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAAt                                               <  1:2390648/51‑1 (MQ=255)
                                         tAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCgg                           >  1:1277035/1‑71 (MQ=255)
                                          aaaaTGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAAtt                                          >  1:5899163/1‑55 (MQ=255)
                                           aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCgg                           <  1:6920621/69‑1 (MQ=255)
                                           aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGtt                         >  1:2082765/1‑71 (MQ=255)
                                           aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGt                          <  1:7321672/70‑1 (MQ=255)
                                           aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGt                          <  1:672675/70‑1 (MQ=255)
                                           aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGt                          <  1:2418330/70‑1 (MQ=255)
                                           aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGt                          <  1:2444846/70‑1 (MQ=255)
                                           aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGt                          <  1:2495726/70‑1 (MQ=255)
                                           aaaTGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGAGAAGCGTGCGCGGtt                         >  1:5911194/1‑71 (MQ=255)
                                            aaTGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGttt                        <  1:3468013/71‑1 (MQ=255)
                                             aTGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTtc                       <  1:5544257/71‑1 (MQ=255)
                                              tGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTtct                      >  1:4208070/1‑71 (MQ=255)
                                                tATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTtctc                     <  1:3383255/70‑1 (MQ=255)
                                                  tGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCtcatc                  <  1:5336890/71‑1 (MQ=255)
                                                  tGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCtcat                   >  1:5426026/1‑70 (MQ=255)
                                                   ggTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTtctc                     >  1:5845911/1‑67 (MQ=255)
                                                      tttCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAAt                                               <  1:2160281/38‑1 (MQ=255)
                                                      tttCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGa                                      >  1:1773583/1‑47 (MQ=255)
                                                      tttCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGa                                      >  1:748859/1‑47 (MQ=255)
                                                        tCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTg             <  1:4484582/70‑1 (MQ=255)
                                                         cATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCGTCAAGTGaa           >  1:3154912/1‑71 (MQ=255)
                                                          aTCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGttt                        >  1:129530/1‑57 (MQ=255)
                                                          aTCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCtca                    >  1:4718575/1‑61 (MQ=255)
                                                           tCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGaa                                     <  1:6762619/43‑1 (MQ=255)
                                                           tCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGa                                      <  1:7375454/42‑1 (MQ=255)
                                                           tCTGATTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATg         >  1:5529879/1‑71 (MQ=255)
                                                               cTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGt        <  1:2626952/68‑1 (MQ=255)
                                                               cTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGt        <  1:1523903/68‑1 (MQ=255)
                                                               cTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTc     <  1:327880/71‑1 (MQ=255)
                                                                 tAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTg   >  1:7262553/1‑71 (MQ=255)
                                                                  aaGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTgc  >  1:35165/1‑71 (MQ=255)
                                                                    |                                                                     
TATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAAT‑TTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGC  >  AP009048/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: ‑