Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 214,925 A→C 65.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
*Exported214,9250AC65.5% 282.1 / 128.7 156D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (23/31);  new base C (56/46);  total (79/77)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

TAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGG  >  Exported/214860‑214991
                                                                 |                                                                  
tAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCg                                                               <  1:304744/71‑1 (MQ=255)
 aaaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCg                                                               <  1:1725453/70‑1 (MQ=255)
    aCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAc                                                                 >  1:1691085/1‑65 (MQ=255)
     actcAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                  >  1:902079/4‑63 (MQ=255)
     actcAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                  >  1:2863053/4‑63 (MQ=255)
        cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                       >  1:2913314/1‑71 (MQ=255)
         aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                      >  1:1974809/1‑71 (MQ=255)
          aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     <  1:1193281/71‑1 (MQ=255)
           gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGAt                                                                   <  1:1355428/56‑1 (MQ=255)
              tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAcc                                                                >  1:1614128/1‑56 (MQ=255)
              tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAcc                                                                >  1:2901853/1‑56 (MQ=255)
              tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                  >  1:1916855/1‑70 (MQ=255)
                aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                 <  1:75156/69‑1 (MQ=255)
                 cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              <  1:2239150/71‑1 (MQ=255)
                   tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                            <  1:1644061/71‑1 (MQ=255)
                   tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                            <  1:1056531/71‑1 (MQ=255)
                   tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                            <  1:396653/71‑1 (MQ=255)
                    aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:1965873/1‑68 (MQ=255)
                    aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:1200338/1‑68 (MQ=255)
                    aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                              >  1:2452296/1‑68 (MQ=255)
                     aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:107488/1‑67 (MQ=255)
                     aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:547733/1‑67 (MQ=255)
                     aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAg                                           <  1:2014311/70‑1 (MQ=255)
                     aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                          <  1:2131395/71‑1 (MQ=255)
                     aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                          <  1:1587545/71‑1 (MQ=255)
                      aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                         <  1:2784401/71‑1 (MQ=255)
                      aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                         <  1:633628/71‑1 (MQ=255)
                      aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                         <  1:1245353/71‑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)
                              tATGAGCACTGCAATTACACGCCGGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                  >  1:2861907/1‑70 (MQ=255)
                               aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                      >  1:2439172/1‑65 (MQ=255)
                                tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                  <  1:2387130/36‑1 (MQ=255)
                                tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                       <  1:2930626/63‑1 (MQ=255)
                                 gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                               >  1:2633792/1‑70 (MQ=255)
                                 gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                               >  1:2231647/1‑70 (MQ=255)
                                  aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                     >  1:1471207/1‑63 (MQ=255)
                                  aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                              <  1:2106724/70‑1 (MQ=255)
                                  aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                              <  1:2110504/70‑1 (MQ=255)
                                  aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                              <  1:2311938/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:2790186/71‑1 (MQ=255)
                                     aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:2008292/1‑68 (MQ=255)
                                     aCTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:1583301/70‑1 (MQ=255)
                                      cTGCATTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                      <  1:1554924/58‑1 (MQ=255)
                                      cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                      <  1:1459091/58‑1 (MQ=255)
                                      cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTa                             <  1:633525/67‑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:308097/1‑71 (MQ=255)
                                        gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                       >  1:2060010/1‑71 (MQ=255)
                                        gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                       >  1:2398904/1‑71 (MQ=255)
                                         cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                               >  1:2541357/1‑46 (MQ=39)
                                          aaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        >  1:933728/1‑52 (MQ=255)
                                          aaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        >  1:1278131/1‑52 (MQ=255)
                                           aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                    >  1:2463700/1‑71 (MQ=255)
                                           aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                    <  1:929951/71‑1 (MQ=255)
                                              acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                      >  1:2474019/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:779062/71‑1 (MQ=255)
                                                 cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt              <  1:2851399/71‑1 (MQ=255)
                                                 cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt              <  1:1141546/71‑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)
                                                       aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                       <  1:1323957/40‑1 (MQ=38)
                                                       aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:1531415/52‑1 (MQ=255)
                                                       aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:1004588/1‑69 (MQ=255)
                                                       aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTg         <  1:2003566/70‑1 (MQ=255)
                                                       aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:447476/71‑1 (MQ=255)
                                                        tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:1777794/1‑68 (MQ=255)
                                                         cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                  >  1:68583/1‑43 (MQ=38)
                                                         cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                         >  1:205268/1‑52 (MQ=255)
                                                         cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:354309/1‑67 (MQ=255)
                                                         cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:2753875/1‑67 (MQ=255)
                                                         cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:2338480/1‑67 (MQ=255)
                                                         cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt           >  1:47970/1‑66 (MQ=255)
                                                          gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                    <  1:2478544/56‑1 (MQ=255)
                                                          gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt     <  1:1279488/71‑1 (MQ=255)
                                                           ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                  <  1:343965/41‑1 (MQ=38)
                                                           ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGGACTATGAAGGCCACAAGATCATTGAGAtt    <  1:2601818/71‑1 (MQ=255)
                                                           ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt    <  1:1775771/71‑1 (MQ=255)
                                                            tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcaca                  >  1:880694/1‑56 (MQ=255)
                                                            tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg   >  1:1074684/1‑71 (MQ=255)
                                                             ctcGCTACCGCAACCACCGGTATGAACCAGATTgg                                      >  1:3043456/1‑35 (MQ=37)
                                                             ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg  >  1:816950/1‑71 (MQ=255)
                                                             ctcGATACCGAACCCACCGGTACGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt              >  1:66293/1‑59 (MQ=255)
                                                                 |                                                                  
TAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGG  >  Exported/214860‑214991

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: