Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,107 A→C 78.5% 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.5% 1762.2 / 282.1 650E14A (GAA→GCA) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (66/74);  new base C (266/243);  total (333/317)
Fisher's exact test for biased strand distribution p-value = 2.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 678 total aligned reads displayed.

AACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTG  >  W3110S.gb/236039‑236175
                                                                    |                                                                    
caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                    >  1:5468908‑M2/6‑71 (MQ=255)
  cGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAacc                                                                  <  1:6610053/71‑1 (MQ=255)
   ctCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                                 >  1:2049019‑M2/3‑71 (MQ=255)
   ctCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                                 >  1:6059205‑M2/3‑71 (MQ=255)
    tCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                                   >  1:300007‑M2/2‑68 (MQ=255)
     cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                               >  1:7026189/1‑71 (MQ=255)
      aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccac                                                                >  1:5850072/1‑69 (MQ=255)
      aaGTCTGACATAAATGACCGCTATGAGCACTCCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                               <  1:5887597/70‑1 (MQ=255)
        gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                             >  1:7107839/1‑70 (MQ=255)
          cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                          >  1:2539827/1‑71 (MQ=255)
             aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                        <  1:1535168/70‑1 (MQ=255)
             aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                        <  1:7627486/70‑1 (MQ=255)
             aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                        <  1:6349245/70‑1 (MQ=255)
               aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                                   <  1:6685854/57‑1 (MQ=255)
                tAAATTACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                    <  1:6086221/71‑1 (MQ=255)
                tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAcc                                                     <  1:2257000/70‑1 (MQ=255)
                 aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                                 >  1:3867051/1‑57 (MQ=255)
                 aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                      >  1:4133361/1‑68 (MQ=255)
                 aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                      >  1:6262299/1‑68 (MQ=255)
                 aaaTGACCGCTATGAGCACTGCAATTACACGCCAGAACGTTCTCGCTACCGCAACCACCGGTATGAAc                                                      >  1:1764831/1‑68 (MQ=255)
                   aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                                 >  1:1901778/1‑71 (MQ=255)
                   aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCg                                                              >  1:6771451/1‑58 (MQ=255)
                    tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                                >  1:4748996/1‑71 (MQ=255)
                     gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                               >  1:7282953/1‑71 (MQ=255)
                      aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                               <  1:280832/70‑1 (MQ=255)
                       ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                             <  1:2956636/71‑1 (MQ=255)
                        cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                      <  1:5159107/61‑1 (MQ=255)
                        cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                      <  1:6900451/45‑1 (MQ=255)
                        cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTg                                            <  1:6968637/71‑1 (MQ=255)
                          cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                          >  1:7891808/1‑71 (MQ=255)
                          cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                          <  1:142722/71‑1 (MQ=255)
                            aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                           >  1:4702941/1‑52 (MQ=255)
                            aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                        >  1:5648926/1‑71 (MQ=255)
                              gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                                 >  1:287350/1‑44 (MQ=255)
                              gAGCACTGCAATTACACGCCAGATCGTTCTCG‑‑CCCGCAACCACCgg                                                             >  1:4050282/1‑46 (MQ=255)
                               aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                          >  1:6445459/1‑50 (MQ=255)
                                 cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   <  1:5308080/71‑1 (MQ=255)
                                 cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   <  1:7060756/71‑1 (MQ=255)
                                  aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   >  1:939987/1‑70 (MQ=255)
                                  aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                  <  1:881490/71‑1 (MQ=255)
                                  aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                  <  1:842764/71‑1 (MQ=255)
                                  aCTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCTCACTATGa                                  <  1:3023026/71‑1 (MQ=255)
                                   cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   >  1:673104/1‑69 (MQ=255)
                                   cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                              >  1:4211894/1‑58 (MQ=255)
                                   cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                 <  1:2009642/71‑1 (MQ=255)
                                    tGCTATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                              >  1:4655486/1‑41 (MQ=255)
                                    tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                           <  1:6208072/44‑1 (MQ=38)
                                    tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                         <  1:393394/46‑1 (MQ=39)
                                    tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                        <  1:1814706/63‑1 (MQ=255)
                                     gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                               >  1:6373281/1‑71 (MQ=255)
                                     gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                               >  1:3531505/1‑71 (MQ=255)
                                      cAATTACACGCGAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                              >  1:398432/1‑71 (MQ=255)
                                        aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                      >  1:1253268/1‑45 (MQ=38)
                                        aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                             >  1:4280058/1‑70 (MQ=255)
                                        aTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                            <  1:745732/71‑1 (MQ=255)
                                         ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                        <  1:4971218/42‑1 (MQ=38)
                                         ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                    <  1:3169465/62‑1 (MQ=255)
                                         ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   >  1:1317819/1‑63 (MQ=255)
                                           acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                        <  1:5103860/40‑1 (MQ=38)
                                           acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcgc                                         >  1:2678465/1‑55 (MQ=255)
                                           acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           >  1:7207730/1‑69 (MQ=255)
                                           acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           >  1:3673080/1‑69 (MQ=255)
                                              cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                      <  1:4405167/39‑1 (MQ=37)
                                               gCCAGATCGTTTTCGCTACCGCAACCACCGGTATGAACCAGAtt                                                <  1:1840905/44‑1 (MQ=25)
                                               gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAcc                                                     >  1:1083816/1‑39 (MQ=37)
                                               gCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                      >  1:3440712/1‑70 (MQ=255)
                                                  aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAATATCAtt                  <  1:3354295/71‑1 (MQ=255)
                                                  aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  <  1:7567058/71‑1 (MQ=255)
                                                   gATCGTTCTCGCTACCGCAACCACCTGTATGAACCAGATTGGTGCGCACTATg                                   <  1:3231219/53‑1 (MQ=255)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                    <  1:7765007/35‑1 (MQ=37)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                    <  1:5754617/35‑1 (MQ=37)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                          <  1:2442682/45‑1 (MQ=38)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                      <  1:4486571/65‑1 (MQ=255)
                                                      cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:34939/1‑67 (MQ=255)
                                                      cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                  >  1:1134130/1‑51 (MQ=255)
                                                      cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:4563005/1‑67 (MQ=255)
                                                       gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgaga              <  1:2729840/70‑1 (MQ=255)
                                                       gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt             >  1:5228301/1‑71 (MQ=255)
                                                        ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           <  1:5036041/56‑1 (MQ=255)
                                                         tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                        >  1:2402570/1‑42 (MQ=38)
                                                         tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg           >  1:602601/1‑71 (MQ=255)
                                                          ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg          >  1:2628446/1‑71 (MQ=255)
                                                          ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg           <  1:7331915/70‑1 (MQ=255)
                                                          ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg           <  1:3586549/70‑1 (MQ=255)
                                                          ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg           <  1:1607871/70‑1 (MQ=255)
                                                           tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt             <  1:3102823/67‑1 (MQ=255)
                                                            cgctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTg        >  1:1953594/4‑71 (MQ=255)
                                                             gctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGc       >  1:4316518/3‑71 (MQ=255)
                                                             gctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGc       >  1:6139908/3‑71 (MQ=255)
                                                              ctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGcc      >  1:20995/2‑71 (MQ=255)
                                                              aTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   <  1:2291982/42‑1 (MQ=255)
                                                              aTACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt             <  1:3727679/64‑1 (MQ=255)
                                                               tACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                 <  1:2169160/43‑1 (MQ=255)
                                                                aCCGCAACCACCGGTATGAACGAGATTGGTGCGCACTATGa                                  <  1:2338826/41‑1 (MQ=38)
                                                                aCCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt            >  1:5376795/1‑63 (MQ=255)
                                                                 ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                  >  1:100279/1‑40 (MQ=255)
                                                                 ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag               <  1:3781144/59‑1 (MQ=255)
                                                                  cGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                 <  1:4851924/40‑1 (MQ=255)
                                                                  cGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTg  <  1:1706442/71‑1 (MQ=255)
                                                                  cGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCAACAAGATCAtt                  >  1:5562867/1‑55 (MQ=255)
                                                                    |                                                                    
AACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTG  >  W3110S.gb/236039‑236175

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: