Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 214,931 A→C 62.6% 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,9310AC62.6% 246.5 / 153.1 163E14A (GAA→GCA) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (26/35);  new base C (55/47);  total (81/82)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

Only 100 of 164 total aligned reads displayed.

CAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTG  >  Exported/214868‑214999
                                                               |                                                                    
cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                               >  1:2913314/1‑71 (MQ=255)
 aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                              >  1:1974809/1‑71 (MQ=255)
  aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                             <  1:1193281/71‑1 (MQ=255)
        aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                         <  1:75156/69‑1 (MQ=255)
         cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                      >  1:47988/1‑71 (MQ=255)
           tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                    <  1:1056531/71‑1 (MQ=255)
           tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                    <  1:1644061/71‑1 (MQ=255)
            aaaTGACCGCTATGAGCACTGCAATTACCCGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                      >  1:2859164/1‑68 (MQ=255)
            aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                      >  1:2452296/1‑68 (MQ=255)
             aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                      >  1:547733/1‑67 (MQ=255)
             aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAg                                                   <  1:2014311/70‑1 (MQ=255)
                gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                               >  1:1508055/1‑71 (MQ=255)
                 aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                              >  1:1244745/1‑71 (MQ=255)
                  ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                    <  1:798223/64‑1 (MQ=255)
                  ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGt                                             <  1:484524/71‑1 (MQ=255)
                     cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                          <  1:1005332/71‑1 (MQ=255)
                     cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgc                                           >  1:2475783/1‑70 (MQ=255)
                       aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                              >  1:2786437/1‑65 (MQ=255)
                        tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                               <  1:2930626/63‑1 (MQ=255)
                         gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                       >  1:2231647/1‑70 (MQ=255)
                         gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                       >  1:2633792/1‑70 (MQ=255)
                          aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                             >  1:1471207/1‑63 (MQ=255)
                          aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                                      <  1:2106724/70‑1 (MQ=255)
                          aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGGGCACt                                      >  1:2961324/1‑70 (MQ=255)
                            cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                    >  1:733020/1‑70 (MQ=255)
                            cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   <  1:163837/71‑1 (MQ=255)
                            cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   <  1:2064412/71‑1 (MQ=255)
                             aCTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   <  1:1583301/70‑1 (MQ=255)
                              cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                              <  1:1459091/58‑1 (MQ=255)
                              cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTa                                     <  1:633525/67‑1 (MQ=255)
                              cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                 <  1:2997230/71‑1 (MQ=255)
                              cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                 <  1:328287/71‑1 (MQ=255)
                               tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                              <  1:977758/57‑1 (MQ=255)
                                gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                                  <  1:1525340/52‑1 (MQ=255)
                                gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                          >  1:1525667/1‑60 (MQ=255)
                                gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                          >  1:1398741/1‑60 (MQ=255)
                                gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                               >  1:2398904/1‑71 (MQ=255)
                                gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                               >  1:2060010/1‑71 (MQ=255)
                                 cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                       >  1:1573836/1‑46 (MQ=39)
                                 cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                              >  1:1591736/1‑71 (MQ=255)
                                 cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                              >  1:2492969/1‑71 (MQ=255)
                                 cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                           >  1:2429607/1‑42 (MQ=255)
                                  aaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                                >  1:933728/1‑52 (MQ=255)
                                   aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                            <  1:929951/71‑1 (MQ=255)
                                   aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                            >  1:2463700/1‑71 (MQ=255)
                                      acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                              >  1:1276727/1‑50 (MQ=255)
                                      acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                 >  1:751821/1‑63 (MQ=255)
                                      acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                 >  1:1813492/1‑63 (MQ=255)
                                        acGCCATATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           >  1:390562/1‑67 (MQ=255)
                                        acGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           >  1:2568038/1‑67 (MQ=255)
                                         cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                      <  1:650227/71‑1 (MQ=255)
                                         cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                      <  1:2851399/71‑1 (MQ=255)
                                         cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                      <  1:1141546/71‑1 (MQ=255)
                                          gCCAGATCGTTCTCTCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                     <  1:212159/71‑1 (MQ=255)
                                          gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                  <  1:2132076/58‑1 (MQ=255)
                                          gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                  <  1:1033580/58‑1 (MQ=255)
                                          gCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                 >  1:569970/1‑59 (MQ=255)
                                          gCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                 >  1:2756034/1‑59 (MQ=255)
                                             aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  <  1:1893813/71‑1 (MQ=255)
                                             aGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACt                                      >  1:298287/1‑51 (MQ=255)
                                              gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:236864/1‑70 (MQ=255)
                                              gATCGTTCTCGATACCGAAACCACCGGTATGAACCAg                                                   <  1:1809079/37‑1 (MQ=255)
                                               aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                               <  1:1323957/40‑1 (MQ=38)
                                               aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                               <  1:1946933/40‑1 (MQ=38)
                                               aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   <  1:1531415/52‑1 (MQ=255)
                                               aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:356461/1‑69 (MQ=255)
                                               aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:1004588/1‑69 (MQ=255)
                                               aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTg                 <  1:2003566/70‑1 (MQ=255)
                                               aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTg                 <  1:1215868/70‑1 (MQ=255)
                                               aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                <  1:1342709/71‑1 (MQ=255)
                                                tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag               <  1:1140421/71‑1 (MQ=255)
                                                tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag               <  1:815260/71‑1 (MQ=255)
                                                 cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                          >  1:68583/1‑43 (MQ=38)
                                                 cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                          >  1:437167/1‑43 (MQ=38)
                                                 cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  >  1:2338480/1‑67 (MQ=255)
                                                 cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt                   >  1:47970/1‑66 (MQ=255)
                                                  gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt             <  1:1279488/71‑1 (MQ=255)
                                                   ttCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                               <  1:2759997/36‑1 (MQ=37)
                                                   ttCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                               <  1:2006393/36‑1 (MQ=37)
                                                   ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                          <  1:343965/41‑1 (MQ=38)
                                                   ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGGACTATGAAGGCCACAAGATCATTGAGAtt            <  1:2601818/71‑1 (MQ=255)
                                                   ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt            <  1:1775771/71‑1 (MQ=255)
                                                   ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt            <  1:3055432/71‑1 (MQ=255)
                                                    tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcaca                          >  1:880694/1‑56 (MQ=255)
                                                    tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg           >  1:1074684/1‑71 (MQ=255)
                                                     ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg          >  1:2890682/1‑71 (MQ=255)
                                                     ctcGATACCGAACCCACCGGTACGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                      >  1:66293/1‑59 (MQ=255)
                                                       cgctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTg        >  1:1766332/4‑71 (MQ=255)
                                                        gctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGc       >  1:1706159/3‑71 (MQ=255)
                                                        gctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGc       >  1:2017828/3‑71 (MQ=255)
                                                           aCCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  <  1:234917/57‑1 (MQ=255)
                                                           aCCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                  <  1:129772/57‑1 (MQ=255)
                                                           aCCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGt    >  1:2740526/1‑71 (MQ=255)
                                                            ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           <  1:2521753/47‑1 (MQ=255)
                                                            ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                           <  1:1951250/47‑1 (MQ=255)
                                                            ccGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGtt   <  1:2659893/71‑1 (MQ=255)
                                                            ccGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGt    <  1:672703/70‑1 (MQ=255)
                                                             cGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                   <  1:2752051/38‑1 (MQ=255)
                                                              gAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTg  >  1:2282502/1‑70 (MQ=255)
                                                              gAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTg  >  1:763448/1‑70 (MQ=255)
                                                               |                                                                    
CAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTG  >  Exported/214868‑214999

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