Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,560,351 G→A 16.4% P62L (CCC→CTC)  rbn ← tRNA processing exoribonuclease BN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,560,3510GA16.4% 136.8 / 18.0 55P62L (CCC→CTC) rbntRNA processing exoribonuclease BN
Reads supporting (aligned to +/- strand):  ref base G (30/16);  new base A (4/5);  total (34/21)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

AGGCAGAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAGGCAAG  >  W3110S.gb/3560293‑3560421
                                                          |                                                                      
aggcagAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAgcggcg                                                             >  1:2744496/1‑70 (MQ=255)
  gcagAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATgg                                                                        <  1:2263244/57‑1 (MQ=255)
   cagAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGaa                                                           >  1:682476/1‑69 (MQ=255)
    agAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAg                                                                  >  1:2674891/1‑61 (MQ=255)
    agAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAg                                                        >  1:745864/1‑71 (MQ=255)
       aaGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAgcg                                                      >  1:1096630/1‑70 (MQ=255)
       aaGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAgcg                                                      >  1:1456290/1‑70 (MQ=255)
        aGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAgcgc                                                     >  1:2269431/1‑70 (MQ=255)
        aGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGAGAAAAGCGGCGaa                                                           <  1:2292525/64‑1 (MQ=255)
        aGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGAGAAAAGCGGCGaa                                                           <  1:789279/64‑1 (MQ=255)
         gTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGa                                                            <  1:1358166/62‑1 (MQ=255)
         gTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCaa                                                   <  1:725551/71‑1 (MQ=255)
          ttGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCaaa                                                  <  1:2710399/71‑1 (MQ=255)
          ttGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCaaa                                                  <  1:2099988/71‑1 (MQ=255)
           tGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGAGAAAAGCGGCGAAAAGCGCaa                                                   >  1:2513794/1‑69 (MQ=255)
            ggCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGAGAAAAGCGGCGAAAAGCGCaaaaa                                                >  1:368831/1‑71 (MQ=255)
             gCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaac                                               <  1:1634367/71‑1 (MQ=255)
             gCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaac                                               <  1:2753517/71‑1 (MQ=255)
               aaaaaTAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaacaa                                             >  1:1706211/1‑71 (MQ=255)
               aaaaaTAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaacaa                                             >  1:2129523/1‑71 (MQ=255)
                           aCGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCg                                 >  1:1853007/1‑71 (MQ=255)
                           aCGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCg                                 >  1:356874/1‑71 (MQ=255)
                             gTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaac                                               <  1:1705655/55‑1 (MQ=255)
                             gTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCgg                                >  1:2280596/1‑70 (MQ=255)
                              tAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCa                              >  1:1350157/1‑71 (MQ=255)
                                aCTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAAcc                                     >  1:1717166/1‑62 (MQ=255)
                                aCTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCAc                             >  1:2090506/1‑70 (MQ=255)
                                aCTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCAc                             >  1:2145394/1‑70 (MQ=255)
                                   gAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTaa                          >  1:285433/1‑70 (MQ=255)
                                   gAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAAt                         >  1:776476/1‑71 (MQ=255)
                                    aaTGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAAt                         <  1:1681570/70‑1 (MQ=255)
                                     aTGCTGACGTCGGAAAACATGAGAAAAGCGGCGAAAAGCGCAAAAACAACGGc                                         >  1:905995/1‑53 (MQ=255)
                                      tGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaac                                               <  1:1884592/46‑1 (MQ=255)
                                        cTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACt                            <  1:209107/63‑1 (MQ=255)
                                        cTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGagag                    >  1:1703781/1‑71 (MQ=255)
                                        cTGACGTCGGAAAACATGAGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGagag                    >  1:1267559/1‑71 (MQ=255)
                                        cTGACGTCGGAAAACATAGGAAAAGCGGCGAAAAGCGCAAAAACAACTGCAACCAGCGGCACt                            <  1:1738807/63‑1 (MQ=255)
                                           aCGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCaa                 >  1:2560804/1‑71 (MQ=255)
                                           aCGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCaa                 >  1:1104393/1‑71 (MQ=255)
                                              tCGGAAAACATGAGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGa              <  1:2675560/71‑1 (MQ=255)
                                                ggAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaacaa                                             >  1:2494378/1‑38 (MQ=255)
                                                ggAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGaca            <  1:1413701/71‑1 (MQ=255)
                                                ggAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGaca            >  1:1529627/1‑71 (MQ=255)
                                                ggAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGaca            >  1:16184/1‑71 (MQ=255)
                                                ggAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGaca            >  1:23782/1‑71 (MQ=255)
                                                ggAAAACATGAGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGaca            <  1:2275268/71‑1 (MQ=255)
                                                 gAAAACATGAGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGaca            <  1:474669/70‑1 (MQ=255)
                                                  aaaaCATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGAcaca          >  1:192130/1‑71 (MQ=255)
                                                  aaaaCATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGAcaca          <  1:224413/71‑1 (MQ=255)
                                                   aaaCATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGAcaca          >  1:2360864/1‑70 (MQ=255)
                                                   aaaCATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGAcaca          >  1:1318837/1‑70 (MQ=255)
                                                   aaaCATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGAcaca          >  1:793212/1‑70 (MQ=255)
                                                    aaCATGGGAAAAGCGTCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGaca            <  1:2355430/67‑1 (MQ=255)
                                                      cATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAgg      >  1:2238602/1‑71 (MQ=255)
                                                        tGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCa                  <  1:2029042/57‑1 (MQ=255)
                                                        tGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCa                  <  1:1473685/57‑1 (MQ=255)
                                                          ggAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAGGCaa   <  1:2291234/70‑1 (MQ=255)
                                                          ggAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAGGCaa   <  1:1498023/70‑1 (MQ=255)
                                                          ggAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAGGCaa   <  1:1154352/70‑1 (MQ=255)
                                                          ggAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAGGCAAg  <  1:900264/71‑1 (MQ=255)
                                                          |                                                                      
AGGCAGAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAGGCAAG  >  W3110S.gb/3560293‑3560421

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