Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,101 A→C 52.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,1010AC52.5% 174.7 / 245.4 202D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (45/51);  new base C (60/46);  total (105/97)
Fisher's exact test for biased strand distribution p-value = 2.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 203 total aligned reads displayed.

TAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGC  >  W3110S.gb/236036‑236170
                                                                 |                                                                     
tAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAcc                                                                   >  1:4799769/1‑70 (MQ=255)
 tacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCg                                                                  >  1:654568‑M2/8‑70 (MQ=255)
 aaaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCg                                                                  >  1:3123107/1‑70 (MQ=255)
 aaaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCg                                                                  >  1:3411355/1‑70 (MQ=255)
  acaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACcgc                                                                 >  1:3647265‑M2/7‑69 (MQ=255)
  acaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACcgc                                                                 >  1:3719707‑M2/7‑69 (MQ=255)
  aaaCGTCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                                <  1:4822220/71‑1 (MQ=255)
   caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAccgca                                                                >  1:3203720‑M2/6‑68 (MQ=255)
   caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                               >  1:1104171‑M2/6‑71 (MQ=255)
   caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                               >  1:2740899‑M2/6‑71 (MQ=255)
     cGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                              >  1:1576920/1‑70 (MQ=255)
        cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAcc                                                                   >  1:4470113/1‑62 (MQ=255)
        cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAcc                                                                   >  1:2192584/1‑62 (MQ=255)
        cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccacc                                                          >  1:941261/1‑71 (MQ=255)
         aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                          <  1:1023851/70‑1 (MQ=255)
          aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                        <  1:4353117/71‑1 (MQ=255)
          aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                        <  1:776647/71‑1 (MQ=255)
              tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                     >  1:4756910/1‑54 (MQ=255)
              tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                     >  1:4980634/1‑54 (MQ=255)
              tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAAGCACCGGTa                                                      >  1:1171248/1‑69 (MQ=255)
              tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                     >  1:3466797/1‑70 (MQ=255)
              tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                     >  1:4914868/1‑70 (MQ=255)
                aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                   <  1:3895586/70‑1 (MQ=255)
                aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                   <  1:4263771/70‑1 (MQ=255)
                aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                   <  1:3698612/70‑1 (MQ=255)
                aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                   <  1:1338160/70‑1 (MQ=255)
                  aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                 >  1:4857816/1‑70 (MQ=255)
                  aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                 >  1:5186755/1‑70 (MQ=255)
                   tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                               <  1:4880596/71‑1 (MQ=255)
                   tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                               <  1:783994/71‑1 (MQ=255)
                   tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                               <  1:4048251/71‑1 (MQ=255)
                   tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                               <  1:4269405/71‑1 (MQ=255)
                   tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                        <  1:874379/62‑1 (MQ=255)
                     aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                 >  1:2440195/1‑67 (MQ=255)
                     aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                 >  1:1675709/1‑67 (MQ=255)
                     aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                 >  1:4399034/1‑67 (MQ=255)
                     aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                                <  1:229876/52‑1 (MQ=255)
                     aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                      >  1:103763/1‑62 (MQ=255)
                     aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                               <  1:3085338/69‑1 (MQ=255)
                     aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAg                                              <  1:384981/70‑1 (MQ=255)
                       tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                           <  1:2950930/71‑1 (MQ=255)
                        gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                          >  1:4101269/1‑71 (MQ=255)
                        gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                          >  1:5110833/1‑71 (MQ=255)
                         aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                         >  1:428507/1‑71 (MQ=255)
                         aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                         >  1:2786684/1‑71 (MQ=255)
                         aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                         >  1:4420558/1‑71 (MQ=255)
                         aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                         >  1:1953321/1‑71 (MQ=255)
                         aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                         >  1:1891866/1‑71 (MQ=255)
                             cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                       <  1:4107474/53‑1 (MQ=255)
                             cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                       <  1:5006649/53‑1 (MQ=255)
                             cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                       <  1:959100/53‑1 (MQ=255)
                             cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                       <  1:975683/53‑1 (MQ=255)
                             cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                       <  1:182089/53‑1 (MQ=255)
                             cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                       <  1:1990488/53‑1 (MQ=255)
                             cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                       <  1:2434974/53‑1 (MQ=255)
                                    cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                        <  1:109184/61‑1 (MQ=255)
                                    cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                        <  1:489776/61‑1 (MQ=255)
                                    cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:2114352/71‑1 (MQ=255)
                                    cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:2280056/71‑1 (MQ=255)
                                    cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:4250319/71‑1 (MQ=255)
                                    cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:2608336/71‑1 (MQ=255)
                                      cTGCAATTACACGTCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                               <  1:4072130/52‑1 (MQ=255)
                                       tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                             <  1:3818411/69‑1 (MQ=255)
                                       tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                             <  1:2384437/69‑1 (MQ=255)
                                        gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                   >  1:3829542/1‑62 (MQ=255)
                                        gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                   >  1:1657055/1‑62 (MQ=255)
                                        gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                          >  1:3022631/1‑71 (MQ=255)
                                        gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                          >  1:2558220/1‑71 (MQ=255)
                                         cAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                           >  1:4601291/1‑37 (MQ=38)
                                         cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                            >  1:450298/1‑52 (MQ=255)
                                         cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                            >  1:3365078/1‑52 (MQ=255)
                                         cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                         <  1:2860232/71‑1 (MQ=255)
                                         cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                         <  1:363238/71‑1 (MQ=255)
                                         cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                         <  1:4211241/71‑1 (MQ=255)
                                         cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                         >  1:1445461/1‑71 (MQ=255)
                                              acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      >  1:1613256/1‑69 (MQ=255)
                                                      gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:895581/53‑1 (MQ=255)
                                                      gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:666738/53‑1 (MQ=255)
                                                      gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:586572/53‑1 (MQ=255)
                                                      gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:531624/53‑1 (MQ=255)
                                                       aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga           <  1:1188375/71‑1 (MQ=255)
                                                       aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:1253476/52‑1 (MQ=255)
                                                       aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga           <  1:70027/71‑1 (MQ=255)
                                                        tcattCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga           <  1:474702/67‑1 (MQ=255)
                                                        tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:5384838/1‑68 (MQ=255)
                                                        tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:2762465/1‑68 (MQ=255)
                                                        tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:4009828/1‑68 (MQ=255)
                                                         cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:5217583/1‑67 (MQ=255)
                                                         cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:3761996/1‑67 (MQ=255)
                                                         cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:1540732/1‑67 (MQ=255)
                                                         cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:2909040/1‑67 (MQ=255)
                                                         cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt              >  1:3545381/1‑66 (MQ=255)
                                                         cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:1780325/1‑67 (MQ=255)
                                                         cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:1555276/1‑67 (MQ=255)
                                                         cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:490185/1‑67 (MQ=255)
                                                          gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                            >  1:1420725/1‑51 (MQ=255)
                                                           ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                   <  1:2146869/59‑1 (MQ=255)
                                                            tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:5507879/47‑1 (MQ=39)
                                                             ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg      <  1:4244804/70‑1 (MQ=255)
                                                                gATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGc  >  1:4524053/1‑71 (MQ=255)
                                                                 |                                                                     
TAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGC  >  W3110S.gb/236036‑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: