Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,107 A→C 78.0% E14A (GAA→GCA)  dnaQ → DNA polymerase III epsilon subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb236,1070AC78.0% 497.4 / 82.5 186E14A (GAA→GCA) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (21/20);  new base C (74/71);  total (95/91)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

Only 100 of 189 total aligned reads displayed.

CAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTT  >  W3110S.gb/236044‑236174
                                                               |                                                                   
 aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                             <  1:138897/71‑1 (MQ=255)
 aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAGCGAAACCACCg                                                             <  1:359241/71‑1 (MQ=255)
 aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCg                                                             <  1:375300/71‑1 (MQ=255)
   gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGt                                                           <  1:1877129/71‑1 (MQ=255)
    tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGt                                                           >  1:1579273/1‑70 (MQ=255)
      tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                          <  1:93999/69‑1 (MQ=255)
      tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                         >  1:1638290/1‑70 (MQ=255)
        aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                              <  1:1451202/63‑1 (MQ=255)
        aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                        <  1:1764904/69‑1 (MQ=255)
        aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                      >  1:1937978/1‑71 (MQ=255)
         cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                      >  1:1532905/1‑70 (MQ=255)
         cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                     <  1:617122/71‑1 (MQ=255)
         cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                     <  1:1122217/71‑1 (MQ=255)
         cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                     >  1:424707/1‑71 (MQ=255)
          aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAcc                                                    <  1:2089554/71‑1 (MQ=255)
          aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                        <  1:880138/67‑1 (MQ=255)
            aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                     >  1:1063558/1‑68 (MQ=255)
            aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                     >  1:1043842/1‑68 (MQ=255)
            aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                            >  1:1111476/1‑61 (MQ=255)
             aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                     >  1:1836812/1‑67 (MQ=255)
             aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                     >  1:1021700/1‑67 (MQ=255)
             aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                                  <  1:1757787/70‑1 (MQ=255)
              aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                                >  1:228549/1‑71 (MQ=255)
              aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                                >  1:1320326/1‑71 (MQ=255)
                 aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                             <  1:1186257/71‑1 (MQ=255)
                  ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaaa                                                                   >  1:268727/1‑48 (MQ=255)
                    gCTATTAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgc                                          >  1:349353/1‑71 (MQ=255)
                    gCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgc                                          >  1:704859/1‑71 (MQ=255)
                     cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                                 >  1:8644/1‑63 (MQ=255)
                      tATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                                  <  1:2024350/61‑1 (MQ=255)
                        tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                      <  1:519784/71‑1 (MQ=255)
                         gAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAcc                                                    <  1:1900510/56‑1 (MQ=255)
                          aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                            >  1:789585/1‑63 (MQ=255)
                          aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                                     <  1:585960/70‑1 (MQ=255)
                          aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACt                                     <  1:1054267/70‑1 (MQ=255)
                            cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                  <  1:2041441/71‑1 (MQ=255)
                            cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                   >  1:286174/1‑70 (MQ=255)
                             aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                  <  1:129356/70‑1 (MQ=255)
                              cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                  <  1:823890/69‑1 (MQ=255)
                              cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                  >  1:659434/1‑69 (MQ=255)
                               tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                               >  1:129658/1‑71 (MQ=255)
                               tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                               >  1:1978970/1‑71 (MQ=255)
                               tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                               >  1:900744/1‑71 (MQ=255)
                               tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                       <  1:14603/47‑1 (MQ=255)
                               tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                  <  1:1088053/68‑1 (MQ=255)
                                 cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                      >  1:194415/1‑46 (MQ=255)
                                 cAATAACACGCCAGATCGTTCTCGCTACAGCAacca                                                                >  1:1829584/1‑36 (MQ=16)
                                  aaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                  <  1:2077993/65‑1 (MQ=255)
                                  aaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                  <  1:1279332/65‑1 (MQ=255)
                                  aaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                              <  1:472892/69‑1 (MQ=255)
                                   aTTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                              <  1:1152686/36‑1 (MQ=38)
                                   aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                >  1:1948865/1‑66 (MQ=255)
                                   aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                 <  1:791176/65‑1 (MQ=255)
                                    ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                          <  1:316924/71‑1 (MQ=255)
                                    ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                          <  1:1572765/71‑1 (MQ=255)
                                      acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                          >  1:1504924/1‑69 (MQ=255)
                                      acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                          >  1:437251/1‑69 (MQ=255)
                                       cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                       <  1:612708/71‑1 (MQ=255)
                                       cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAACAGATTGGTGCGCACTATGAAGGCcac                          >  1:246362/1‑68 (MQ=255)
                                        acGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa                      <  1:1240947/71‑1 (MQ=255)
                                         cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                     <  1:428637/71‑1 (MQ=255)
                                         cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                     <  1:1263929/71‑1 (MQ=255)
                                         cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTTCGCACTATGAAGGCCACAAGa                      >  1:643426/1‑70 (MQ=255)
                                          gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                      <  1:1392823/53‑1 (MQ=255)
                                            cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                   <  1:1032249/70‑1 (MQ=255)
                                             aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                               <  1:1838504/41‑1 (MQ=38)
                                             aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                 <  1:1200929/71‑1 (MQ=255)
                                             aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                 <  1:1309064/71‑1 (MQ=255)
                                             aGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                 <  1:458408/71‑1 (MQ=255)
                                              gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                 >  1:1170472/1‑70 (MQ=255)
                                               aTCGTTCTCGCTACCGCAACCACCTGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                          <  1:1598243/60‑1 (MQ=255)
                                               aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAATATCATTga               <  1:221830/71‑1 (MQ=255)
                                               aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                 >  1:127474/1‑69 (MQ=255)
                                               aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                     <  1:1667104/65‑1 (MQ=255)
                                                tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                       >  1:979851/1‑46 (MQ=39)
                                                  gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt            <  1:1831978/71‑1 (MQ=255)
                                                  gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt            <  1:685870/71‑1 (MQ=255)
                                                   ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt           >  1:1510811/1‑71 (MQ=255)
                                                    tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                 >  1:1624570/1‑64 (MQ=255)
                                                    tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt            <  1:142055/69‑1 (MQ=255)
                                                    tctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                >  1:1550265/1‑49 (MQ=255)
                                                     ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt            >  1:1377655/1‑68 (MQ=255)
                                                     ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg         >  1:54771/1‑71 (MQ=255)
                                                     ctcGATACCGAAACCACCGGTATGAACCAGATTgg                                             >  1:924695/1‑35 (MQ=255)
                                                       cGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTg       >  1:814885/1‑71 (MQ=255)
                                                        gctACCGCAACCACCGGTATGATCCAGATTGGTGCGCACTATGAAg                               >  1:1498104/3‑46 (MQ=255)
                                                        gctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgaga             <  1:1702801/62‑1 (MQ=255)
                                                        gctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt        <  1:564747/67‑1 (MQ=255)
                                                        gctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGc      >  1:749391/3‑71 (MQ=255)
                                                        gctACCGCAACCACCGGTATGAACCAGATTGGTCCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGc      >  1:1872231/3‑71 (MQ=255)
                                                        gATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt        <  1:379219/69‑1 (MQ=255)
                                                         ctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                <  1:1151531/43‑1 (MQ=255)
                                                         ctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGcc     <  1:1412788/70‑1 (MQ=255)
                                                          tACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCg    >  1:1935492/1‑71 (MQ=255)
                                                          tACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCg    >  1:1076720/1‑71 (MQ=255)
                                                          tACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTCAGATTGGTGcc     <  1:1718471/70‑1 (MQ=255)
                                                           aCCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGt   >  1:2102887/1‑71 (MQ=255)
                                                            ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg          >  1:1681502/1‑63 (MQ=255)
                                                            ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCg    <  1:666649/69‑1 (MQ=255)
                                                            ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGtt  <  1:898883/71‑1 (MQ=255)
                                                               |                                                                   
CAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTT  >  W3110S.gb/236044‑236174

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