Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,101 A→C 93.7% 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,1010AC93.7% 3543.6 / 95.0 1058D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (37/31);  new base C (534/456);  total (571/487)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.36e-01

Only 100 of 1061 total aligned reads displayed.

CTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGG  >  W3110S.gb/236033‑236167
                                                                    |                                                                  
 gcgtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAc                                                                 >  1:3343480‑M2/11‑71 (MQ=255)
     acaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACcgc                                                              >  1:1608658‑M2/7‑69 (MQ=255)
     acaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACCCCAGATCGTTCTCGCTACCGCAac                                                           >  1:2515841‑M2/7‑71 (MQ=255)
       aactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:3550189‑M2/5‑70 (MQ=255)
        actCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacc                                                          <  1:3751927‑M2/68‑1 (MQ=255)
          tCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                  <  1:842700‑M2/60‑1 (MQ=255)
            aaGTCTGACATAAATGCCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                      >  1:786204/1‑71 (MQ=255)
             aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     <  1:1273241/71‑1 (MQ=255)
             aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     >  1:2091394/1‑71 (MQ=255)
              gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     >  1:3832397/1‑70 (MQ=255)
               tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGt                                                    <  1:24543/70‑1 (MQ=255)
               tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                   >  1:1337426/1‑71 (MQ=255)
                   aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                               >  1:1904062/1‑71 (MQ=255)
                    cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:3592444/1‑71 (MQ=255)
                     aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:3951146/1‑70 (MQ=255)
                      tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                 <  1:4287636/66‑1 (MQ=255)
                      tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                              >  1:2575442/1‑69 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:1005307/1‑67 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:928838/1‑67 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:1697558/1‑67 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                              >  1:435321/1‑67 (MQ=255)
                         aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACAACCGGTATGAACCAGAt                                         >  1:888449/1‑71 (MQ=255)
                          tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                           >  1:2817967/1‑68 (MQ=255)
                          tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        <  1:580189/71‑1 (MQ=255)
                          tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        >  1:3705956/1‑71 (MQ=255)
                          tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        >  1:3730529/1‑71 (MQ=255)
                          tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        <  1:2059875/71‑1 (MQ=255)
                            aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                         >  1:2209428/1‑52 (MQ=255)
                            aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                      <  1:2775637/71‑1 (MQ=255)
                            aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                       <  1:4623322/70‑1 (MQ=255)
                            aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                       <  1:4314327/70‑1 (MQ=255)
                                cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                  <  1:4014085/71‑1 (MQ=255)
                                cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                  <  1:3066327/71‑1 (MQ=255)
                                cTATGAGCACTGCAATTAAACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                  <  1:2611118/71‑1 (MQ=255)
                                  aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                       <  1:2413461/48‑1 (MQ=255)
                                  aTGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                >  1:3314254/1‑71 (MQ=255)
                                   tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCTCAACCACCGGTATGAACCAGATTGGTgcgc                                 <  1:507630/69‑1 (MQ=255)
                                      gCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTg                                    <  1:3468891/63‑1 (MQ=255)
                                      gCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                              >  1:1417533/1‑69 (MQ=255)
                                      gCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                            <  1:1778631/71‑1 (MQ=255)
                                      gCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                            <  1:2745061/71‑1 (MQ=255)
                                      gCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                            <  1:4152399/71‑1 (MQ=255)
                                       cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:2532659/71‑1 (MQ=255)
                                        aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           >  1:3002131/1‑70 (MQ=255)
                                         cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           >  1:1814142/1‑69 (MQ=255)
                                          tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        <  1:4416237/55‑1 (MQ=255)
                                          tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTg                                    <  1:291365/59‑1 (MQ=255)
                                          tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:3168748/68‑1 (MQ=255)
                                           gCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                           >  1:2843271/1‑35 (MQ=38)
                                           gCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                           >  1:550213/1‑35 (MQ=38)
                                           gCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                           >  1:4120904/1‑35 (MQ=38)
                                           gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                       >  1:1453808/1‑71 (MQ=255)
                                           gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                       >  1:4562929/1‑71 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                         >  1:688048/1‑36 (MQ=38)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                          >  1:956344/1‑51 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:4179747/1‑71 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:4040550/1‑71 (MQ=255)
                                             aaTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:110809/1‑70 (MQ=255)
                                             aaTTAAACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                       >  1:3279438/1‑69 (MQ=255)
                                              aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                    >  1:3419440/1‑71 (MQ=255)
                                               ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                     <  1:2943831/53‑1 (MQ=255)
                                               ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:2703753/63‑1 (MQ=255)
                                                tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:3638057/1‑70 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                         >  1:1949477/1‑63 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:1671667/1‑69 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:3395699/1‑69 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:2000814/1‑69 (MQ=255)
                                                  cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        <  1:3856153/47‑1 (MQ=39)
                                                    cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                <  1:215393/37‑1 (MQ=37)
                                                    cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt              <  1:1058842/71‑1 (MQ=255)
                                                    cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt              <  1:1141132/71‑1 (MQ=255)
                                                     gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt              >  1:4255857/1‑70 (MQ=255)
                                                     gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt              >  1:47488/1‑70 (MQ=255)
                                                     gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc             <  1:1520483/71‑1 (MQ=255)
                                                      ccAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa            <  1:2676913/71‑1 (MQ=255)
                                                      ccAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa            >  1:1568838/1‑71 (MQ=255)
                                                       cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGAGCACTATGAAGGCCACAAGATCAt           >  1:4459862/1‑71 (MQ=255)
                                                        aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                    >  1:607023/1‑61 (MQ=255)
                                                        aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          <  1:4072167/71‑1 (MQ=255)
                                                         gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:2330342/53‑1 (MQ=255)
                                                         gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:3234221/53‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                           <  1:3897569/36‑1 (MQ=37)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                         <  1:3663603/38‑1 (MQ=37)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                               <  1:229518/48‑1 (MQ=39)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                <  1:737097/63‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:1848713/71‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:1744552/71‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:3112044/71‑1 (MQ=255)
                                                           tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:1192877/1‑68 (MQ=255)
                                                           tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:1375180/70‑1 (MQ=255)
                                                            cGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                          >  1:2370282/1‑35 (MQ=37)
                                                             gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                            >  1:2639698/1‑48 (MQ=39)
                                                             gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgaga      <  1:3506164/70‑1 (MQ=255)
                                                             gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt     <  1:2891650/71‑1 (MQ=255)
                                                              ttCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        <  1:4519116/35‑1 (MQ=37)
                                                              ttCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        <  1:2248947/35‑1 (MQ=37)
                                                              ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt    <  1:1239110/71‑1 (MQ=255)
                                                               tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa            >  1:4224066/1‑62 (MQ=255)
                                                               tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg   >  1:633452/1‑71 (MQ=255)
                                                                ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa               >  1:4378199/1‑58 (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: