Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 214,931 A→C 59.3% 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
*Exported214,9310AC59.3% 232.7 / 177.0 179E14A (GAA→GCA) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (40/33);  new base C (61/45);  total (101/78)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 186 total aligned reads displayed.

AAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTG  >  Exported/214862‑214999
                                                                     |                                                                    
aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                                     <  1:651895/71‑1 (MQ=255)
aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                      >  1:1111766/1‑70 (MQ=255)
aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                      >  1:2284414/1‑70 (MQ=255)
   actcAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                                   >  1:323615/4‑70 (MQ=255)
   actcAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                                   >  1:3754532/4‑70 (MQ=255)
       aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                              >  1:3856553/1‑71 (MQ=255)
        aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                             <  1:2046901/71‑1 (MQ=255)
        aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                             <  1:1604579/71‑1 (MQ=255)
          tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                      <  1:2535304/60‑1 (MQ=255)
           cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                          >  1:2394724/1‑71 (MQ=255)
           cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                          >  1:2104254/1‑71 (MQ=255)
            ttaCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                          >  1:3620847/3‑70 (MQ=255)
            tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                          >  1:1357862/1‑70 (MQ=255)
            tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                          >  1:280384/1‑70 (MQ=255)
            tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                          >  1:2839764/1‑70 (MQ=255)
            tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                          >  1:926086/1‑70 (MQ=255)
            tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                          >  1:951067/1‑70 (MQ=255)
             gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                         >  1:1463658/1‑70 (MQ=255)
             gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                         >  1:3400368/1‑70 (MQ=255)
             gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                         >  1:3164205/1‑70 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                         <  1:961907/69‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                       >  1:3376900/1‑71 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                        <  1:2048648/70‑1 (MQ=255)
               cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                      >  1:2267009/1‑71 (MQ=255)
                 tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                                    <  1:629233/71‑1 (MQ=255)
                  aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                      >  1:1586860/1‑68 (MQ=255)
                  aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                      >  1:2366385/1‑68 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                                   <  1:861723/70‑1 (MQ=255)
                        ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                         <  1:1091951/59‑1 (MQ=255)
                        ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                         <  1:1570897/59‑1 (MQ=255)
                         cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCTGTATGAACCAGa                                                  <  1:1319208/65‑1 (MQ=255)
                           cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                    <  1:3478426/61‑1 (MQ=255)
                           cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                          >  1:476613/1‑71 (MQ=255)
                           cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                          >  1:2430405/1‑71 (MQ=255)
                              tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                       >  1:808593/1‑71 (MQ=255)
                               gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                        <  1:1811304/53‑1 (MQ=255)
                                  cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   <  1:2118704/71‑1 (MQ=255)
                                   aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                  <  1:3176102/71‑1 (MQ=255)
                                    cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   >  1:1901996/1‑69 (MQ=255)
                                    cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   >  1:302549/1‑69 (MQ=255)
                                    cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                 <  1:2710634/71‑1 (MQ=255)
                                      gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                               >  1:1271650/1‑71 (MQ=255)
                                      gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                               >  1:26844/1‑71 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                               >  1:2062236/1‑70 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                              >  1:2805829/1‑71 (MQ=255)
                                          ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                         <  1:3561550/41‑1 (MQ=38)
                                          ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                         <  1:1216136/41‑1 (MQ=38)
                                          ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           <  1:746207/71‑1 (MQ=255)
                                          ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                            >  1:2628962/1‑70 (MQ=255)
                                          ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                                    <  1:2597436/46‑1 (MQ=255)
                                          ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                                    <  1:3694489/46‑1 (MQ=255)
                                           tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           >  1:2138201/1‑70 (MQ=255)
                                           tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           >  1:1680840/1‑70 (MQ=255)
                                           tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTCGTGCGCACTATGAAGGCcac                           >  1:3072485/1‑70 (MQ=255)
                                            acacGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                             >  1:3771099/1‑35 (MQ=38)
                                            acacGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                             >  1:1664280/1‑35 (MQ=38)
                                            acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                                  >  1:3631421/1‑46 (MQ=39)
                                            acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                                  >  1:1106066/1‑46 (MQ=39)
                                            acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                               >  1:2387973/1‑49 (MQ=39)
                                            acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                               >  1:2123754/1‑49 (MQ=39)
                                            acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                                  >  1:387369/1‑46 (MQ=255)
                                            acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           >  1:281241/1‑69 (MQ=255)
                                             cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           >  1:981899/1‑68 (MQ=255)
                                             cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           >  1:311043/1‑68 (MQ=255)
                                             cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           >  1:2557333/1‑68 (MQ=255)
                                             cacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           >  1:2214715/1‑68 (MQ=255)
                                              acGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa                       <  1:2948754/71‑1 (MQ=255)
                                               cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                      <  1:3084819/71‑1 (MQ=255)
                                                gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   <  1:3336593/57‑1 (MQ=255)
                                                  cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                  >  1:1954293/1‑56 (MQ=255)
                                                   aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt                   >  1:2495641/1‑70 (MQ=255)
                                                    gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                       <  1:303488/49‑1 (MQ=39)
                                                    gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                       <  1:2979531/49‑1 (MQ=39)
                                                    gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:1741491/1‑70 (MQ=255)
                                                    gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                              <  1:973439/42‑1 (MQ=255)
                                                     aTCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                    <  1:3404282/35‑1 (MQ=37)
                                                     aTCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                    <  1:3512670/35‑1 (MQ=37)
                                                     aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                <  1:3820466/71‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                <  1:1720890/71‑1 (MQ=255)
                                                      tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                     >  1:1048496/1‑49 (MQ=39)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                                  >  1:2779928/1‑35 (MQ=37)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:1657426/1‑67 (MQ=255)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:100565/1‑67 (MQ=255)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:3165839/1‑67 (MQ=255)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:2963051/1‑67 (MQ=255)
                                                        gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt             <  1:2105629/71‑1 (MQ=255)
                                                        gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt             >  1:459079/1‑71 (MQ=255)
                                                        gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgaga              <  1:3058477/70‑1 (MQ=255)
                                                        gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt             <  1:1139455/71‑1 (MQ=255)
                                                         ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                    <  1:3841156/47‑1 (MQ=39)
                                                         ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt            <  1:3779756/71‑1 (MQ=255)
                                                           ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                  >  1:3040722/1‑47 (MQ=39)
                                                             cGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt                   >  1:1170778/1‑60 (MQ=255)
                                                              gctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg          <  1:1153219/66‑1 (MQ=255)
                                                              gATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                          >  1:2724824/1‑36 (MQ=255)
                                                               ctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGcc      >  1:326204/2‑71 (MQ=255)
                                                                tACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                     <  1:3729291/39‑1 (MQ=255)
                                                                  ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt                   >  1:183887/1‑55 (MQ=255)
                                                                  ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCg     <  1:670723/69‑1 (MQ=255)
                                                                   cGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTg  <  1:2309903/71‑1 (MQ=255)
                                                                     |                                                                    
AAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTG  >  Exported/214862‑214999

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: