Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,089,635 C→T 6.1% R211R (CGC→CGT endA → DNA‑specific endonuclease I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,089,6350CT 6.1% 200.0 / 2.2 67R211R (CGC→CGTendADNA‑specific endonuclease I
Reads supporting (aligned to +/- strand):  ref base C (32/30);  new base T (2/2);  total (35/32)
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.09e-01

GCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAA  >  W3110S.gb/3089568‑3089703
                                                                   |                                                                    
gCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAgcgc                                                                      <  1:1113266/68‑1 (MQ=255)
  cAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATg                                                                  >  1:321727/1‑70 (MQ=255)
  cAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGa                                                                 <  1:168276/71‑1 (MQ=255)
  cAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGa                                                                 <  1:1570655/71‑1 (MQ=255)
    gCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATg                                                                  <  1:405594/68‑1 (MQ=255)
      aaaCGCAGCTGTTCAACGCATGGGACAAGATGTATCCAGTTACCGACTGGGAGTGCGAGCGTGATGAACGc                                                             <  1:1046204/71‑1 (MQ=255)
      aaaCGCAGCTGTTCAACGCATGGGACAAGATGTATCCAGTTACCGACTGGGAGTGCGAGCGTGATGAACGc                                                             <  1:996372/71‑1 (MQ=255)
      aaaCGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGc                                                             <  1:3480487/71‑1 (MQ=255)
        aCGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCAt                                                           >  1:1655161/1‑71 (MQ=255)
        aCGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCAt                                                           >  1:3535114/1‑71 (MQ=255)
        aCGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCAt                                                           <  1:3899941/71‑1 (MQ=255)
        aCGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCAt                                                           <  1:1222595/71‑1 (MQ=255)
         cGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATc                                                          >  1:1012442/1‑71 (MQ=255)
         cGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATc                                                          >  1:972432/1‑71 (MQ=255)
         cGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATc                                                          >  1:2638998/1‑71 (MQ=255)
            aGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGa                                                                    >  1:1590075/1‑58 (MQ=255)
             gCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGa                                                      <  1:1315767/71‑1 (MQ=255)
             gCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGa                                                      <  1:2056088/71‑1 (MQ=255)
             gCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGa                                                      <  1:1833156/71‑1 (MQ=255)
              cTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATcgcg                                                       >  1:1867452/1‑69 (MQ=255)
              cTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGaa                                                     <  1:3292544/71‑1 (MQ=255)
                gTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGa                                                      <  1:4378178/68‑1 (MQ=255)
                 ttCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGt                                                  >  1:521523/1‑71 (MQ=255)
                  tCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAg                                                    >  1:2395708/1‑68 (MQ=255)
                     aCGCATGGAACAAGATTTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAg                                              >  1:4212739/1‑71 (MQ=255)
                          tGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCaa                                         <  1:33496/71‑1 (MQ=255)
                           ggAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGa                                                      <  1:1664489/57‑1 (MQ=255)
                            gAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGGGATGAACGCATCGCGaa                                                     >  1:330872/1‑57 (MQ=255)
                            gAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCAt                                                           >  1:1633014/1‑51 (MQ=255)
                            gAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCAt                                                           >  1:2883269/1‑51 (MQ=255)
                            gAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGc                                                >  1:1759892/1‑62 (MQ=255)
                            gAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAAt                                        >  1:4334473/1‑70 (MQ=255)
                            gAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAAt                                        >  1:4183224/1‑70 (MQ=255)
                            gAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAAt                                        >  1:3291917/1‑70 (MQ=255)
                            gAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAAt                                        >  1:1596620/1‑70 (MQ=255)
                            gAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAAt                                        >  1:2091776/1‑70 (MQ=255)
                             aaCAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGa                                                      <  1:3199405/55‑1 (MQ=255)
                             aaCAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATc                                       >  1:824161/1‑70 (MQ=255)
                              aCAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCAt                                     <  1:3423078/71‑1 (MQ=255)
                               cAAGATGTATCCAGTTACCGACTTGGAGTGCGAGCGTGATGAACGCATCGCGAAGGTGCAGGGCAATCAt                                     >  1:3026408/1‑70 (MQ=255)
                               cAAGATGTATCCAGTTACCGACTTGGAGTGCGAGCGTGATGAACGCATCGCGAAGGTGCAGGGCAATCAt                                     >  1:3021518/1‑70 (MQ=255)
                                 aGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCAt                                     >  1:4019202/1‑68 (MQ=255)
                                     gTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCg                               >  1:2633123/1‑70 (MQ=255)
                                         ccGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAAt                                        >  1:2214662/1‑57 (MQ=255)
                                          cGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATaa                                   >  1:2785532/1‑61 (MQ=255)
                                          cGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTAtgtg                         <  1:347980/71‑1 (MQ=255)
                                           ggTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCCTATGTGc                        <  1:472272/71‑1 (MQ=255)
                                                   aCTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTg                                                 >  1:875252/1‑38 (MQ=255)
                                                   aCTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATa                                    >  1:2441117/1‑51 (MQ=255)
                                                   aCTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCt                >  1:1506353/1‑71 (MQ=255)
                                                       ggAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTg              <  1:3847276/69‑1 (MQ=255)
                                                       ggAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGcc            <  1:2874837/71‑1 (MQ=255)
                                                       ggAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGc             <  1:1532399/70‑1 (MQ=255)
                                                         aGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAg          <  1:4337556/71‑1 (MQ=255)
                                                          gTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAgg         >  1:3324654/1‑71 (MQ=255)
                                                           tGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTAt                            <  1:4236902/51‑1 (MQ=255)
                                                           tGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGgc        >  1:1835464/1‑71 (MQ=255)
                                                           tGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGgc        >  1:1328253/1‑71 (MQ=255)
                                                            gCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTAtgt                          >  1:3375619/1‑52 (MQ=255)
                                                             cGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGgcgc      >  1:4077065/1‑71 (MQ=255)
                                                               aGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGa    <  1:2925094/71‑1 (MQ=255)
                                                               aGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGa    <  1:860934/71‑1 (MQ=255)
                                                               aGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGa    <  1:855980/71‑1 (MQ=255)
                                                               aGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGa    <  1:2672927/71‑1 (MQ=255)
                                                               aGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGa    <  1:3499591/71‑1 (MQ=255)
                                                               aGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGa    <  1:402024/71‑1 (MQ=255)
                                                               aGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGa    <  1:2577950/71‑1 (MQ=255)
                                                               aGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGGACGCGCTTGCCAGGCGCGaa   <  1:4336169/71‑1 (MQ=255)
                                                                gcgcGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGa    >  1:1717738/1‑70 (MQ=255)
                                                                 cgcgATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGaaa  <  1:243929/71‑1 (MQ=255)
                                                                   |                                                                    
GCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAA  >  W3110S.gb/3089568‑3089703

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: