Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,101 A→C 63.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,1010AC63.1% 336.2 / 159.8 192D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (35/36);  new base C (72/49);  total (107/85)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 200 total aligned reads displayed.

AACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGC  >  W3110S.gb/236031‑236170
                                                                      |                                                                     
   gcgtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                     >  1:1182459‑M2/11‑70 (MQ=255)
       aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                 <  1:2288900/70‑1 (MQ=255)
       aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCACGATACCGaa                                                                <  1:3351444/71‑1 (MQ=255)
        caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                               >  1:2084728‑M2/6‑71 (MQ=255)
            tCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                           >  1:4099162‑M2/2‑71 (MQ=255)
             cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                          >  1:4788626/1‑71 (MQ=255)
              aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                         >  1:2678533/1‑71 (MQ=255)
              aaGTCTGACATAAATGACCGCTATGAGCACTGAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                        >  1:1043925/1‑71 (MQ=255)
                gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                        <  1:4152369/70‑1 (MQ=255)
                  cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                     >  1:3025529/1‑71 (MQ=255)
                   tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                                 <  1:1518523/58‑2 (MQ=255)
                   tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                              >  1:1333407/1‑61 (MQ=255)
                   tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                     >  1:1124790/1‑70 (MQ=255)
                   tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                     >  1:4355575/1‑70 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                     <  1:3686647/68‑1 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                    <  1:1240182/69‑1 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                  >  1:2985718/1‑71 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                  <  1:24970/71‑1 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                   <  1:2446062/70‑1 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                  >  1:2909307/1‑71 (MQ=255)
                      cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                                 <  1:3790398/55‑2 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                       >  1:4583996/1‑46 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                       >  1:1024525/1‑46 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                       >  1:2432118/1‑46 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                       >  1:98915/1‑46 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                 >  1:2549623/1‑68 (MQ=255)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                 >  1:2175701/1‑67 (MQ=255)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCg                                                                  >  1:700333/1‑50 (MQ=255)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                                >  1:4617353/1‑52 (MQ=255)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                 >  1:1393894/1‑67 (MQ=255)
                            tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                           >  1:4217966/1‑71 (MQ=255)
                             gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACAGGTATGAACCAGATTg                                          >  1:4754433/1‑71 (MQ=255)
                              aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                          <  1:1676775/70‑1 (MQ=255)
                                cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                              >  1:737641/1‑48 (MQ=255)
                                  cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                               <  1:2364674/61‑1 (MQ=255)
                                  cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                    <  1:2321148/56‑1 (MQ=255)
                                  cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                     <  1:2644027/71‑1 (MQ=255)
                                    aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                             >  1:2375533/1‑61 (MQ=255)
                                     tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                   <  1:3066549/70‑1 (MQ=255)
                                     tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                  <  1:3899211/71‑1 (MQ=255)
                                       aGCACTGCAATTACACGCCAGCTCGTTCTCGCTACCGcaa                                                               >  1:4293870/1‑40 (MQ=255)
                                       aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                >  1:349807/1‑71 (MQ=255)
                                        gCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAacca                                                            >  1:3447610/1‑42 (MQ=255)
                                          aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:695002/70‑1 (MQ=255)
                                          aCTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:3849895/70‑1 (MQ=255)
                                           cTGCAATTACACGCCAGATCGTTCTCTATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                               <  1:2402584/68‑1 (MQ=255)
                                           cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                               <  1:1282706/68‑1 (MQ=255)
                                            tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                               <  1:4476240/51‑1 (MQ=255)
                                             gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                     >  1:4186514/1‑60 (MQ=255)
                                             gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                            >  1:3580077/1‑69 (MQ=255)
                                             gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                           >  1:3758437/1‑70 (MQ=255)
                                             gCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                              >  1:4572064/1‑35 (MQ=255)
                                             gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                           >  1:3757188/1‑70 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                         >  1:4687226/1‑71 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                         >  1:1636560/1‑71 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                         >  1:4245761/1‑71 (MQ=255)
                                                aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                             >  1:3323035/1‑49 (MQ=39)
                                                aTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                         >  1:4197698/1‑53 (MQ=255)
                                                aTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                     <  1:2791694/57‑1 (MQ=255)
                                                aTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                       >  1:2214345/1‑71 (MQ=255)
                                                 ttACACTCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      <  1:4020305/71‑1 (MQ=255)
                                                   acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                               >  1:1577504/1‑44 (MQ=38)
                                                   acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                               >  1:1472666/1‑44 (MQ=38)
                                                   acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                            >  1:273505/1‑63 (MQ=255)
                                                   acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      >  1:4086836/1‑69 (MQ=255)
                                                   acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      >  1:1159769/1‑69 (MQ=255)
                                                   acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      >  1:1130054/1‑69 (MQ=255)
                                                   acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      >  1:519927/1‑69 (MQ=255)
                                                   acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      >  1:154580/1‑69 (MQ=255)
                                                   acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      >  1:2138281/1‑69 (MQ=255)
                                                    cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                    <  1:3083759/70‑1 (MQ=255)
                                                       gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                 >  1:1980734/1‑70 (MQ=255)
                                                        ccAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:1881586/56‑1 (MQ=255)
                                                         cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                            >  1:2251411/1‑57 (MQ=255)
                                                          aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                           <  1:867012/41‑1 (MQ=38)
                                                          aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                            >  1:2292801/1‑56 (MQ=255)
                                                          aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:1806071/1‑71 (MQ=255)
                                                          aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:4534571/1‑71 (MQ=255)
                                                          aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:1899194/1‑71 (MQ=255)
                                                           gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                         <  1:3179388/42‑1 (MQ=255)
                                                           gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                               <  1:612733/52‑1 (MQ=255)
                                                            aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                           <  1:1485050/39‑1 (MQ=37)
                                                            aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             <  1:3212186/69‑1 (MQ=255)
                                                            aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                           <  1:588953/39‑1 (MQ=255)
                                                            aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                          <  1:4449259/40‑1 (MQ=255)
                                                            aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                          <  1:1213686/40‑1 (MQ=255)
                                                            aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga           <  1:2690386/71‑1 (MQ=255)
                                                             tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                         >  1:4035849/1‑40 (MQ=38)
                                                             tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga           <  1:3145493/70‑1 (MQ=255)
                                                             tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:3530430/1‑68 (MQ=255)
                                                               gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt        >  1:3320199/1‑71 (MQ=255)
                                                                ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      <  1:2657558/56‑1 (MQ=255)
                                                                ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga           <  1:3918338/67‑1 (MQ=255)
                                                                ttCTCGATACCGAAACCACCGGTATGAACCAGATTg                                          >  1:2225386/1‑36 (MQ=255)
                                                                  ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg     >  1:283838/1‑71 (MQ=255)
                                                                   tcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag          <  1:3734092/65‑1 (MQ=255)
                                                                   tcGATACCGAAACCCCCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa               >  1:3828323/1‑60 (MQ=255)
                                                                   tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACt                                 <  1:9340/42‑1 (MQ=255)
                                                                     gATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt    <  1:666606/69‑1 (MQ=255)
                                                                     gATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGc  >  1:3620005/1‑71 (MQ=255)
                                                                      |                                                                     
AACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGC  >  W3110S.gb/236031‑236170

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: