Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,101 A→C 88.1% 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,1010AC88.1% 1911.3 / 122.4 601D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (37/35);  new base C (299/230);  total (336/265)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 601 total aligned reads displayed.

ACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGG  >  W3110S.gb/236032‑236167
                                                                     |                                                                  
     tacaactCAAGTCTGACATAAATGACCTCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACcgc                                                              <  1:15603‑M2/64‑2 (MQ=255)
     aaaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCg                                                               >  1:21507/1‑70 (MQ=255)
      acaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACcgc                                                              <  1:3020743‑M2/64‑2 (MQ=255)
       caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:1951874‑M2/6‑71 (MQ=255)
       caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:2769128‑M2/6‑71 (MQ=255)
       caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:3777601‑M2/6‑71 (MQ=255)
       caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:423391‑M2/6‑71 (MQ=255)
        aactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                           <  1:1063378‑M2/67‑1 (MQ=255)
         actCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                           >  1:2668172‑M2/4‑70 (MQ=255)
            cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                       >  1:788474/1‑71 (MQ=255)
            cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                       >  1:2989540/1‑71 (MQ=255)
            cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                       >  1:484525/1‑71 (MQ=255)
            cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAacc                                                          >  1:3473104/1‑68 (MQ=255)
             aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                      >  1:1767881/1‑71 (MQ=255)
             aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                      >  1:1207097/1‑71 (MQ=255)
              aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     <  1:1335247/71‑1 (MQ=255)
              aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     <  1:1231532/71‑1 (MQ=255)
                  tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                   >  1:2278217/1‑69 (MQ=255)
                  tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                  >  1:2004576/1‑70 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                 <  1:1973103/69‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                <  1:3938422/70‑1 (MQ=255)
                     cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                               >  1:463558/1‑70 (MQ=255)
                       tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     <  1:304912/62‑1 (MQ=255)
                       tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                            <  1:2370378/71‑1 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACAGGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:1110019/1‑53 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGt                                                    >  1:682343/1‑61 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                  >  1:3176095/1‑63 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:3207454/1‑67 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:1169937/1‑67 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:2228753/1‑67 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                           <  1:2817151/70‑1 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAGGAAc                                              >  1:3973600/1‑67 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTAGCTACCGCAACCACCGGTATGAAc                                              >  1:1411487/1‑67 (MQ=255)
                          aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                         <  1:95088/71‑1 (MQ=255)
                           tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        <  1:3920359/71‑1 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                           <  1:1135639/66‑1 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                          >  1:653594/1‑67 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                      <  1:1352884/71‑1 (MQ=255)
                               cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTg                                    <  1:2044776/71‑1 (MQ=255)
                                 cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                  <  1:3253341/71‑1 (MQ=255)
                                   aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                >  1:52012/1‑71 (MQ=255)
                                    tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                              <  1:1004718/40‑2 (MQ=255)
                                    tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                >  1:2470122/1‑70 (MQ=255)
                                    tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                               <  1:3729591/71‑1 (MQ=255)
                                    tGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                               <  1:3410196/71‑1 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                              <  1:2861209/70‑1 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:2505958/1‑71 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:1825515/1‑71 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                      >  1:560158/1‑62 (MQ=255)
                                       gCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                              <  1:3435389/37‑1 (MQ=255)
                                        cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     >  1:3605046/1‑45 (MQ=255)
                                        cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                 >  1:3539625/1‑49 (MQ=39)
                                          cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           >  1:179626/1‑69 (MQ=255)
                                           tGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                           <  1:1901403/36‑1 (MQ=38)
                                           tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     <  1:2914811/42‑1 (MQ=255)
                                           tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        <  1:2599510/55‑1 (MQ=255)
                                            gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                >  1:2084213/1‑62 (MQ=255)
                                            gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                       >  1:3853850/1‑71 (MQ=255)
                                            gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                        >  1:1702514/1‑70 (MQ=255)
                                            gCAATTACACGACAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                         >  1:1740044/1‑69 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCACCCACCGGTATGaa                                               >  1:233186/1‑46 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                         >  1:2702519/1‑68 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:3638627/1‑71 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:2653319/1‑71 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      <  1:59354/68‑1 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:2556656/63‑1 (MQ=255)
                                                 tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                   >  1:3491775/1‑38 (MQ=37)
                                                  acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:3227092/1‑69 (MQ=255)
                                                  acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:1299239/1‑69 (MQ=255)
                                                  acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:2190151/1‑59 (MQ=255)
                                                   cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                <  1:715444/71‑1 (MQ=255)
                                                    acGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                           <  1:4058536/43‑1 (MQ=38)
                                                    acGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:3331296/1‑67 (MQ=255)
                                                     cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:1533025/1‑66 (MQ=255)
                                                     cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGac              <  1:3276929/71‑2 (MQ=255)
                                                         aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          <  1:85316/71‑1 (MQ=255)
                                                          gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt              <  1:4021549/66‑1 (MQ=255)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcaaa                  <  1:3125428/61‑3 (MQ=255)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:761444/71‑1 (MQ=255)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:3358006/71‑1 (MQ=255)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:1617565/71‑1 (MQ=255)
                                                           aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:2740212/1‑69 (MQ=255)
                                                            tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:883762/1‑68 (MQ=255)
                                                            tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:1863194/1‑68 (MQ=255)
                                                            tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:3819849/70‑1 (MQ=255)
                                                            tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:1594728/1‑68 (MQ=255)
                                                             cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgga                                     >  1:1628012/1‑39 (MQ=37)
                                                             cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                      >  1:3083438/1‑39 (MQ=37)
                                                             cGTTCTCGCTACCGCAACCACCGGTATGAACAAGATTGGTGCGCACTATGAAGGCCACAAGAt              >  1:3316889/1‑63 (MQ=255)
                                                               ttCTCTATACCGAAACCACCGGTATGAACCAGAtt                                        <  1:3690417/35‑1 (MQ=255)
                                                               ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt              <  1:2876623/61‑1 (MQ=255)
                                                               ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt              <  1:1500109/61‑1 (MQ=255)
                                                               ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt    <  1:1663651/71‑1 (MQ=255)
                                                               ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt    <  1:2021299/71‑1 (MQ=255)
                                                               ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt    <  1:2173566/71‑1 (MQ=255)
                                                               ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt     <  1:2428989/70‑1 (MQ=255)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                         >  1:125934/1‑48 (MQ=39)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt    <  1:3884650/69‑1 (MQ=255)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg  >  1:1449669/1‑71 (MQ=255)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg  >  1:1024417/1‑71 (MQ=255)
                                                                     |                                                                  
ACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGG  >  W3110S.gb/236032‑236167

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: