Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,101 A→C 91.4% D12A (GAT→GCT)  dnaQ → DNA polymerase III epsilon subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb236,1010AC91.4% 2696.1 / 104.4 851D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (43/31);  new base C (434/342);  total (477/374)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

Only 100 of 853 total aligned reads displayed.

AACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGC  >  W3110S.gb/236031‑236170
                                                                      |                                                                     
  cgcgtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                     >  1:3795848‑M2/12‑71 (MQ=255)
    gTAAAACGCCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAcc                                                                   >  1:3243706/1‑71 (MQ=255)
       aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                 <  1:697835/70‑1 (MQ=255)
        caactCAAGTCTGACATAAATGACCGCTATGGGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                               >  1:2557175‑M2/6‑71 (MQ=255)
        caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAccgca                                                                >  1:2943515‑M2/6‑68 (MQ=255)
        caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                               >  1:1359744‑M2/6‑71 (MQ=255)
         aCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                              >  1:4318873/1‑71 (MQ=255)
          actCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                              >  1:593045‑M2/4‑70 (MQ=255)
           ctCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                            >  1:752109‑M2/3‑71 (MQ=255)
            tCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacc                                                             <  1:806208‑M2/68‑1 (MQ=255)
            tCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                              <  1:1403242‑M2/67‑1 (MQ=255)
             cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                          >  1:447682/1‑71 (MQ=255)
             cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                          >  1:3034730/1‑71 (MQ=255)
             cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                          >  1:3905548/1‑71 (MQ=255)
              aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                         <  1:183187/71‑1 (MQ=255)
              aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                         >  1:200848/1‑71 (MQ=255)
               aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                        <  1:3389124/71‑1 (MQ=255)
                gTCTGACATAAATGGCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAACGCAACCACCgg                                                        >  1:4017488/1‑70 (MQ=255)
                gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGt                                                       <  1:2092922/71‑1 (MQ=255)
                  cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                     >  1:1878984/1‑71 (MQ=255)
                   tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                     >  1:3707590/1‑70 (MQ=255)
                   tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                     >  1:335247/1‑70 (MQ=255)
                    gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                    >  1:2119970/1‑70 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                    <  1:1783895/69‑1 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                    <  1:2392119/69‑1 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                    <  1:4233710/69‑1 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                   <  1:1446583/70‑1 (MQ=255)
                        tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                        <  1:2987399/62‑1 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                 >  1:1908260/1‑68 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                 >  1:3457299/1‑68 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                 >  1:4268316/1‑68 (MQ=255)
                          aaTGCCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                     >  1:3184634/1‑47 (MQ=39)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                 >  1:2652075/1‑67 (MQ=255)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                 >  1:129217/1‑67 (MQ=255)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                 >  1:3972051/1‑67 (MQ=255)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                 >  1:3624505/1‑67 (MQ=255)
                           aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                            >  1:2277603/1‑71 (MQ=255)
                            tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACAGGTATGAACCAGAtt                                           >  1:224197/1‑71 (MQ=255)
                              aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                         <  1:1938431/71‑1 (MQ=255)
                                  cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAc                                                                    >  1:2907190/1‑40 (MQ=255)
                                  cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                              >  1:3551667/1‑46 (MQ=255)
                                   tATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                        <  1:958940/51‑1 (MQ=255)
                                    aTGAGCACTGCAATTCCACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                      >  1:486144/1‑52 (MQ=255)
                                    aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                      >  1:363903/1‑52 (MQ=255)
                                    aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                   >  1:3867117/1‑71 (MQ=255)
                                     tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcgc                                    <  1:1977401/69‑1 (MQ=255)
                                       aGCACTGCAATTACACGCGAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                            >  1:2434618/1‑59 (MQ=255)
                                       aGCACTGCAATTACACGCCAGATCGTTCTCGCTAc                                                                    >  1:2878386/1‑35 (MQ=255)
                                       aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTa                                >  1:3756181/1‑71 (MQ=255)
                                       aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTa                                >  1:3216362/1‑71 (MQ=255)
                                        gCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                        >  1:351435/1‑46 (MQ=255)
                                         cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAATCACCgg                                                        <  1:5378/45‑1 (MQ=255)
                                         cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:2351800/71‑1 (MQ=255)
                                         cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:83480/71‑1 (MQ=255)
                                         cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:2719698/71‑1 (MQ=255)
                                           cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                   <  1:2199474/64‑1 (MQ=255)
                                            tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                   <  1:3282329/47‑1 (MQ=39)
                                            tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                           <  1:1626760/55‑1 (MQ=255)
                                             gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                           >  1:406717/1‑54 (MQ=255)
                                             gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                        >  1:1257273/1‑57 (MQ=255)
                                             gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                               >  1:4092993/1‑66 (MQ=255)
                                             gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                          >  1:3085919/1‑71 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                         >  1:3832311/1‑71 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                         >  1:2528580/1‑71 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                         >  1:1568356/1‑71 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                         >  1:1310837/1‑71 (MQ=255)
                                                 ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                         <  1:1136340/52‑1 (MQ=255)
                                                 ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      <  1:4409488/71‑1 (MQ=255)
                                                  tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                             >  1:1494238/1‑47 (MQ=39)
                                                   acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                      >  1:1518779/1‑37 (MQ=37)
                                                   acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                   >  1:2592729/1‑56 (MQ=255)
                                                   acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      >  1:2825882/1‑69 (MQ=255)
                                                      cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                 <  1:3431645/71‑1 (MQ=255)
                                                      cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                 <  1:1209321/71‑1 (MQ=255)
                                                       gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                           <  1:2477987/44‑1 (MQ=38)
                                                       gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                 <  1:3684781/70‑1 (MQ=255)
                                                       gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                 >  1:918591/1‑70 (MQ=255)
                                                       gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                <  1:2623142/71‑1 (MQ=255)
                                                         cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:3309696/55‑1 (MQ=255)
                                                         cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      <  1:2757313/63‑1 (MQ=255)
                                                         cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt              >  1:666914/1‑71 (MQ=255)
                                                           gATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                               <  1:2680789/36‑1 (MQ=37)
                                                           gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:1115010/1‑70 (MQ=255)
                                                            aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                <  1:3927335/66‑1 (MQ=255)
                                                            aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga           <  1:4215188/71‑1 (MQ=255)
                                                            aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga           <  1:2167409/71‑1 (MQ=255)
                                                             tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:4130880/1‑68 (MQ=255)
                                                             tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:1660115/1‑68 (MQ=255)
                                                             tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:357447/1‑68 (MQ=255)
                                                             tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:1167514/1‑67 (MQ=255)
                                                             tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                                  <  1:1697250/47‑1 (MQ=255)
                                                              cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa               >  1:1022279/1‑65 (MQ=255)
                                                              cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:4166141/1‑67 (MQ=255)
                                                              cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                         >  1:315177/1‑55 (MQ=255)
                                                               gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt        <  1:1002559/71‑1 (MQ=255)
                                                                ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt       <  1:876242/71‑1 (MQ=255)
                                                                  ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgaga         >  1:175397/1‑67 (MQ=255)
                                                                  ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg     >  1:1066896/1‑71 (MQ=255)
                                                                  ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg      >  1:4353533/1‑70 (MQ=255)
                                                                  ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                            >  1:2037194/1‑48 (MQ=255)
                                                                      |                                                                     
AACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGC  >  W3110S.gb/236031‑236170

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: