Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 214,925 A→C 52.9% 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,9250AC52.9% 222.1 / inf 278D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (58/73);  new base C (84/63);  total (142/136)
Fisher's exact test for biased strand distribution p-value = 4.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 279 total aligned reads displayed.

ACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGT  >  Exported/214856‑214992
                                                                     |                                                                   
      aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                               <  1:3097199/70‑1 (MQ=255)
         cGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                            >  1:2482950/1‑70 (MQ=255)
           aCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccac                                                         <  1:3392258/71‑1 (MQ=255)
             aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                       >  1:3644375/1‑71 (MQ=255)
               gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                      >  1:3071653/1‑70 (MQ=255)
               gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                      >  1:1668619/1‑70 (MQ=255)
               gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                     <  1:1129891/71‑1 (MQ=255)
               gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                     <  1:124178/71‑1 (MQ=255)
                tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                    >  1:3284769/1‑71 (MQ=255)
                 cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                    >  1:4472036/1‑70 (MQ=255)
                  tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                              <  1:1804055/59‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                 <  1:1547752/70‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                >  1:188524/1‑71 (MQ=255)
                      aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                               >  1:35274/1‑70 (MQ=255)
                        aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                               >  1:1061314/1‑68 (MQ=255)
                        aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                               >  1:3537983/1‑68 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                >  1:1260779/1‑66 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                >  1:3872717/1‑66 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                           <  1:2724982/71‑1 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccacc                                                        >  1:2339750/1‑58 (MQ=255)
                           tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                     >  1:2741717/1‑43 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                       >  1:2798762/1‑71 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                       >  1:3452540/1‑71 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                       >  1:1868003/1‑71 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                       <  1:4385137/71‑1 (MQ=255)
                                  tATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                   >  1:3082786/1‑70 (MQ=255)
                                   aTGAGCACTGCAATTACACGCCAGATAGTTCTCGATAcc                                                                 <  1:1990463/39‑1 (MQ=255)
                                    tGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                                >  1:1303076/1‑71 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                              >  1:4400338/1‑71 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                              >  1:4600587/1‑71 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                    >  1:2273052/1‑49 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                    >  1:3569971/1‑49 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACt                               <  1:4001730/70‑1 (MQ=255)
                                        cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                             <  1:3850539/54‑1 (MQ=255)
                                        cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                             <  1:782337/54‑1 (MQ=255)
                                          cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                          <  1:2508769/71‑1 (MQ=255)
                                          cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                          <  1:622726/71‑1 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                           >  1:3889757/1‑67 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                         <  1:4303246/69‑1 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                         <  1:872682/69‑1 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                       >  1:2154318/1‑71 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                       >  1:3206419/1‑71 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                           <  1:3801427/48‑1 (MQ=39)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                         <  1:4630093/50‑1 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            <  1:1846324/63‑1 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACt                               <  1:169074/60‑1 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCtac                    <  1:90299/71‑1 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                    <  1:1410018/71‑1 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                    <  1:2440193/71‑1 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                    <  1:4252025/71‑1 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                    <  1:1277137/71‑1 (MQ=255)
                                                 tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            <  1:396571/62‑1 (MQ=255)
                                                  acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                              >  1:1761177/1‑59 (MQ=255)
                                                  acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgc                                    >  1:3759372/1‑53 (MQ=255)
                                                  acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgc                                    >  1:3623848/1‑53 (MQ=255)
                                                  acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                    >  1:800361/1‑69 (MQ=255)
                                                  acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                    >  1:3355878/1‑69 (MQ=255)
                                                     cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                    >  1:73192/1‑66 (MQ=255)
                                                     cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                 <  1:1398073/53‑1 (MQ=255)
                                                     cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                    >  1:2124053/1‑66 (MQ=255)
                                                      gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa                >  1:4012141/1‑69 (MQ=255)
                                                      gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt               >  1:4540186/1‑70 (MQ=255)
                                                      gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt               >  1:1145855/1‑70 (MQ=255)
                                                      gCCAGATCGTTCTAGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt               >  1:4517568/1‑70 (MQ=255)
                                                       ccAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            <  1:1713032/56‑1 (MQ=255)
                                                       ccAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                        <  1:318052/44‑1 (MQ=255)
                                                        cAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt            >  1:2761058/1‑71 (MQ=255)
                                                        cAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt            >  1:939753/1‑71 (MQ=255)
                                                          gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt           >  1:302036/1‑70 (MQ=255)
                                                          gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt           >  1:2702549/1‑70 (MQ=255)
                                                          gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                 <  1:182306/48‑1 (MQ=255)
                                                          gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                 <  1:1449201/48‑1 (MQ=255)
                                                           aTCGTTCTCGCTACCGGAACCACCGGTATGAACCa                                             <  1:2229441/35‑1 (MQ=37)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                             <  1:1244178/51‑1 (MQ=255)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                  >  1:1045023/1‑62 (MQ=255)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                  >  1:2689814/1‑62 (MQ=255)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                  >  1:4732583/1‑62 (MQ=255)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga         <  1:2944849/71‑1 (MQ=255)
                                                           aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            <  1:1939943/52‑1 (MQ=255)
                                                           aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                            <  1:4611009/52‑1 (MQ=255)
                                                           aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga         <  1:2643765/71‑1 (MQ=255)
                                                            tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTg                                     >  1:2995009/1‑42 (MQ=38)
                                                            tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt           >  1:586337/1‑68 (MQ=255)
                                                            tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt           >  1:4190544/1‑68 (MQ=255)
                                                            tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga         <  1:4489747/70‑1 (MQ=255)
                                                            tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                   >  1:4136350/1‑44 (MQ=255)
                                                            tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                   >  1:538208/1‑44 (MQ=255)
                                                             cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt           >  1:1499561/1‑67 (MQ=255)
                                                             cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt           >  1:685851/1‑67 (MQ=255)
                                                             cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt           >  1:405431/1‑67 (MQ=255)
                                                             cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt           >  1:489704/1‑67 (MQ=255)
                                                              gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt      >  1:2092974/1‑71 (MQ=255)
                                                              gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt      >  1:142712/1‑71 (MQ=255)
                                                               ttCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                         <  1:3498178/35‑1 (MQ=37)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt               >  1:457599/1‑59 (MQ=255)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg   >  1:59002/1‑71 (MQ=255)
                                                                 ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                         >  1:2208173/1‑49 (MQ=255)
                                                                  tcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt  >  1:1560310/1‑71 (MQ=255)
                                                                  tcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt  >  1:4208024/1‑71 (MQ=255)
                                                                  tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt               <  1:2890651/58‑1 (MQ=255)
                                                                     |                                                                   
ACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGT  >  Exported/214856‑214992

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: