Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,101 A→C 57.0% 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,1010AC57.0% 296.1 / 248.6 244D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (43/62);  new base C (75/64);  total (118/126)
Fisher's exact test for biased strand distribution p-value = 5.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

Only 100 of 245 total aligned reads displayed.

ACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTG  >  W3110S.gb/236032‑236166
                                                                     |                                                                 
aCTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGAt                                                                  <  1:2163243/71‑1 (MQ=255)
 cTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATa                                                                 <  1:5027250/71‑1 (MQ=255)
  gcgtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAc                                                                >  1:2998207‑M2/11‑71 (MQ=255)
   gTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAcc                                                               >  1:4139278/1‑71 (MQ=255)
    tAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCg                                                              <  1:4941262/71‑1 (MQ=255)
    gtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCg                                                              >  1:4040934‑M2/9‑71 (MQ=255)
     tacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCg                                                              >  1:2962896‑M2/8‑70 (MQ=255)
       caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                           >  1:1392250‑M2/6‑71 (MQ=255)
       caactCAAGTCTGACATAAATGACCGCTATGAGCACTGAAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                           >  1:2724835‑M2/6‑71 (MQ=255)
         cGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGAt                                                                  <  1:908312/62‑1 (MQ=255)
          ctCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                        >  1:876636‑M2/3‑71 (MQ=255)
           tCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                        >  1:318893‑M2/2‑70 (MQ=255)
           aCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccac                                                       >  1:2042282/1‑71 (MQ=255)
            cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccac                                                       <  1:2629142/70‑1 (MQ=255)
             aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAcc                                                               <  1:166513/61‑1 (MQ=255)
             aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                      <  1:169521/70‑1 (MQ=255)
              aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                    <  1:4562492/71‑1 (MQ=255)
              aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCg                                                     >  1:2717257/1‑70 (MQ=255)
               gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGt                                                   <  1:1349113/71‑1 (MQ=255)
               gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                   <  1:4334521/71‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                               <  1:679043/70‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                <  1:4159922/69‑1 (MQ=255)
                      aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                   >  1:4537119/1‑48 (MQ=255)
                      aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                             >  1:1321565/1‑70 (MQ=255)
                      aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAcc                                            <  1:4417343/71‑1 (MQ=255)
                       tAAATGACGGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                           <  1:3551787/71‑1 (MQ=255)
                       tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                           <  1:146267/71‑1 (MQ=255)
                        aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGATACCACCgg                                                    <  1:239558/61‑1 (MQ=255)
                        aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                             >  1:110009/1‑68 (MQ=255)
                        aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                             >  1:1140735/1‑68 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                             >  1:3880095/1‑67 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                             >  1:303510/1‑67 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                             >  1:4201602/1‑67 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                             >  1:4782320/1‑67 (MQ=255)
                           tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                       <  1:1193491/71‑1 (MQ=255)
                           tGACCGCTATGAGCACTGCAATTACACGCAAGATGGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                        <  1:1640650/70‑1 (MQ=255)
                             aCCGCTTTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACAGGTATGAACCAGATTg                                      <  1:4720175/70‑1 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                     <  1:3323001/71‑1 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                      <  1:930720/70‑1 (MQ=255)
                              ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                      <  1:1789977/69‑1 (MQ=255)
                                gCTATGAGCACTGCAATTACACGCCAGATCGTTCCCGCTACCGCAACCACCGGTATGAACCAGATTGGTgc                                  >  1:3717259/1‑71 (MQ=255)
                                 cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                 <  1:9448/71‑1 (MQ=255)
                                 cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                 >  1:74507/1‑71 (MQ=255)
                                 cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                  >  1:4806080/1‑54 (MQ=255)
                                   aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                               >  1:15875/1‑71 (MQ=255)
                                     gAGCACTGCAATTACACGCCATATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                              >  1:4909273/1‑70 (MQ=255)
                                     gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                              >  1:639908/1‑70 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                            >  1:4792187/1‑71 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                            >  1:3409663/1‑71 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                            >  1:4307440/1‑71 (MQ=255)
                                        cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                  <  1:1730603/47‑1 (MQ=39)
                                        cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                          <  1:609177/71‑1 (MQ=255)
                                         aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                          <  1:5054781/70‑1 (MQ=255)
                                          cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                               <  1:1541810/64‑1 (MQ=255)
                                          cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                          >  1:2569454/1‑69 (MQ=255)
                                           tGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                      <  1:2438512/40‑1 (MQ=255)
                                           tGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                      <  1:3456302/40‑1 (MQ=255)
                                           tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                       <  1:1519680/55‑1 (MQ=255)
                                           tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                        >  1:4434739/1‑54 (MQ=255)
                                           tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                          <  1:43517/68‑1 (MQ=255)
                                           tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACAGGTATGAACCAGATTGGTGCGCACTATGAAg                       >  1:4817561/1‑71 (MQ=255)
                                            gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                      >  1:1876237/1‑71 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAacc                                                         >  1:1021625/1‑35 (MQ=38)
                                             cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                         >  1:2772506/1‑51 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                     >  1:47519/1‑71 (MQ=255)
                                              aaTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                         <  1:2339359/66‑1 (MQ=255)
                                               aTTACACGCCACATCGTTCTCGCTACCGCAACCACAGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                   <  1:3835948/71‑1 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                           <  1:4605/46‑1 (MQ=39)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                           <  1:102354/46‑1 (MQ=39)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                       <  1:2012014/66‑1 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                   >  1:1206995/1‑70 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                  <  1:2297926/71‑1 (MQ=255)
                                                 tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                  >  1:1843768/1‑70 (MQ=255)
                                                  acacGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                  >  1:1903095/1‑37 (MQ=255)
                                                   cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                       <  1:3076654/47‑1 (MQ=39)
                                                   cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                <  1:2249354/70‑1 (MQ=255)
                                                     cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt             <  1:2201257/71‑1 (MQ=255)
                                                      gCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt             >  1:1439995/1‑70 (MQ=255)
                                                        cAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                              <  1:810879/51‑1 (MQ=255)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                              <  1:4003946/48‑1 (MQ=39)
                                                           aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                        <  1:3067525/38‑1 (MQ=255)
                                                           aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                        <  1:2938474/38‑1 (MQ=255)
                                                           aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                        <  1:5158065/38‑1 (MQ=255)
                                                           aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                <  1:4709842/62‑1 (MQ=255)
                                                           aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt         <  1:3632107/69‑1 (MQ=255)
                                                             cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                     >  1:3106754/1‑39 (MQ=37)
                                                             cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                           <  1:2504193/49‑1 (MQ=39)
                                                             cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt         >  1:2618792/1‑67 (MQ=255)
                                                             cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt         >  1:3935374/1‑67 (MQ=255)
                                                             cGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                         >  1:290853/1‑35 (MQ=255)
                                                             cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag      >  1:1056051/1‑70 (MQ=255)
                                                              gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt    <  1:3294170/71‑1 (MQ=255)
                                                               ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                <  1:4050406/58‑1 (MQ=255)
                                                               ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc            <  1:5203656/62‑1 (MQ=255)
                                                               ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt   >  1:711523/1‑71 (MQ=255)
                                                               ttCTCGATACCGAAACCACCGGTATGAACCAGAtt                                       <  1:3276218/35‑1 (MQ=255)
                                                                tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg  >  1:4554211/1‑71 (MQ=255)
                                                                tctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc            <  1:586618/61‑1 (MQ=255)
                                                                  tcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg               <  1:4240471/56‑1 (MQ=255)
                                                                  tcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga       <  1:2799263/64‑1 (MQ=255)
                                                                     |                                                                 
ACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTG  >  W3110S.gb/236032‑236166

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: