Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,101 A→C 75.5% 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,1010AC75.5% 836.2 / 164.5 318D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (45/33);  new base C (157/82);  total (202/116)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

Only 100 of 318 total aligned reads displayed.

CTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGG  >  W3110S.gb/236033‑236167
                                                                    |                                                                  
    tacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACcgc                                                              <  1:3353352‑M2/64‑2 (MQ=255)
      caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:941619‑M2/6‑71 (MQ=255)
      caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:846831‑M2/6‑71 (MQ=255)
      caacgCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:3536826‑M2/6‑71 (MQ=255)
        actCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:764150‑M2/4‑69 (MQ=255)
        actCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:3584108‑M2/4‑69 (MQ=255)
         ctCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                         >  1:2397605‑M2/3‑71 (MQ=255)
            aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                      >  1:1654330/1‑71 (MQ=255)
               tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     <  1:1928921/69‑1 (MQ=255)
               tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                   >  1:1666829/1‑71 (MQ=255)
                 tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                  >  1:2435564/1‑70 (MQ=255)
                   aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                 <  1:1163781/69‑1 (MQ=255)
                    cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:2163127/1‑71 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                           >  1:1795326/1‑54 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                              >  1:2788374/1‑67 (MQ=255)
                         aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                         <  1:2937351/71‑1 (MQ=255)
                            aCCGCTATGAGCACTGCAATTTCACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                      <  1:3445641/71‑1 (MQ=255)
                            aCCGCTATGAGCACTGCAATTTCACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                      <  1:3087485/71‑1 (MQ=255)
                            aCCGCTATGAGCACTGCAATTTCACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                      <  1:2852468/71‑1 (MQ=255)
                            aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAACGAAACCACCGGTATGAACCAGATTgg                                      >  1:3285067/1‑71 (MQ=255)
                            aCCGCTATAAGCACTGCAATTTCACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                      <  1:2220966/71‑1 (MQ=255)
                                 tATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                   >  1:1369924/1‑53 (MQ=255)
                                   taaGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                  <  1:1239401/66‑1 (MQ=255)
                                   tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                               <  1:2661727/71‑1 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     <  1:3290087/47‑1 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     <  1:3420673/47‑1 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:1207695/1‑71 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:371634/1‑71 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:3200941/1‑71 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:2143100/1‑71 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:2104170/1‑71 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:1257481/1‑71 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:347933/1‑71 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:1101205/1‑71 (MQ=255)
                                      gCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     <  1:566233/46‑1 (MQ=255)
                                      gCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                            <  1:3614273/71‑1 (MQ=255)
                                      gCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                            <  1:1604607/71‑1 (MQ=255)
                                       cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                            >  1:2913970/1‑70 (MQ=255)
                                         cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                      >  1:1422101/1‑58 (MQ=255)
                                         cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           >  1:1891199/1‑69 (MQ=255)
                                          tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:3663645/68‑1 (MQ=255)
                                           gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                         >  1:1457431/1‑69 (MQ=255)
                                           gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                       >  1:2970282/1‑71 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                         >  1:1003215/1‑68 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                         >  1:837833/1‑68 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:638408/1‑71 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:1048058/1‑71 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:2893604/1‑71 (MQ=255)
                                               ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                  <  1:1352108/56‑1 (MQ=255)
                                               ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAtt                           <  1:105857/63‑2 (MQ=255)
                                               ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:17509/63‑1 (MQ=255)
                                                tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:2821616/1‑43 (MQ=38)
                                                tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:2462330/1‑70 (MQ=255)
                                                tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:2499693/1‑70 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCAGCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:2997038/1‑69 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                >  1:916755/1‑56 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                >  1:520077/1‑56 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATTAAGGCcac                   >  1:2254663/1‑69 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:445695/1‑69 (MQ=255)
                                                  cacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                      <  1:460439/49‑1 (MQ=255)
                                                     gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc             <  1:2181659/71‑1 (MQ=255)
                                                     gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc             <  1:2712066/71‑1 (MQ=255)
                                                      ccAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa            <  1:1311553/71‑1 (MQ=255)
                                                        aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                           <  1:182847/38‑1 (MQ=37)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                            <  1:2628541/35‑1 (MQ=37)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcgc                                 <  1:1593092/46‑1 (MQ=39)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:3556832/52‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:737657/1‑69 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:54640/1‑69 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:802828/71‑1 (MQ=255)
                                                          aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                               <  1:1984382/48‑1 (MQ=255)
                                                          aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                               <  1:1950678/48‑1 (MQ=255)
                                                          aTCGTCCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                               <  1:3370953/48‑1 (MQ=255)
                                                           tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt              >  1:888270/1‑64 (MQ=255)
                                                           tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt              >  1:1523183/1‑64 (MQ=255)
                                                           tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt              >  1:3499362/1‑64 (MQ=255)
                                                           tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:169475/1‑68 (MQ=255)
                                                           tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag       <  1:231638/71‑1 (MQ=255)
                                                           tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag       <  1:2556511/71‑1 (MQ=255)
                                                           tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag       <  1:259411/71‑1 (MQ=255)
                                                           tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag       <  1:1294965/71‑1 (MQ=255)
                                                           tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag       <  1:1578811/71‑1 (MQ=255)
                                                            cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:710806/1‑67 (MQ=255)
                                                            cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag       <  1:2361939/70‑1 (MQ=255)
                                                            cGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                          >  1:2018236/1‑35 (MQ=255)
                                                            cGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                          >  1:1565470/1‑35 (MQ=255)
                                                            cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                               >  1:3245686/1‑46 (MQ=255)
                                                            cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:1636953/1‑67 (MQ=255)
                                                            cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:2345376/1‑67 (MQ=255)
                                                            cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:2684028/1‑67 (MQ=255)
                                                            cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGAGCACTATGAAGGCCACAAGATCAtt          >  1:686959/1‑67 (MQ=255)
                                                            cGTTCTCGATACCGAAACCAACGGTATGAACCAGATTGGTGCGCa                                >  1:2812538/1‑45 (MQ=255)
                                                              ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTTCGCACTATGAAGGCCACAAGATCATTGAGAtt    <  1:3637518/71‑1 (MQ=255)
                                                              ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                               <  1:1314997/44‑1 (MQ=38)
                                                              ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt    <  1:3508426/71‑1 (MQ=255)
                                                              ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt    <  1:3169738/71‑1 (MQ=255)
                                                              ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt    <  1:146322/71‑1 (MQ=255)
                                                              ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt    <  1:1092266/71‑1 (MQ=255)
                                                                ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg  >  1:3294920/1‑71 (MQ=255)
                                                                ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg  >  1:3060976/1‑71 (MQ=255)
                                                                    |                                                                  
CTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGG  >  W3110S.gb/236033‑236167

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: