Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 236,107 A→C 27.8% 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
*AP009048236,1070AC27.8% 114.0 / 120.0 126E14A (GAA→GCA) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (51/40);  new base C (21/14);  total (72/54)
Fisher's exact test for biased strand distribution p-value = 8.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 142 total aligned reads displayed.

AAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGA  >  AP009048/236038‑236176
                                                                     |                                                                     
aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                                      <  1:1971051/71‑1 (MQ=255)
aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                                      <  1:4373971/71‑1 (MQ=255)
aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                       <  1:2967882/70‑1 (MQ=255)
aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                       <  1:1213796/70‑1 (MQ=255)
aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                       <  1:600126/70‑1 (MQ=255)
 aaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                       <  1:3531761/69‑1 (MQ=255)
  aCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                                    >  1:1824402/1‑71 (MQ=255)
  aCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                                    >  1:2778797/1‑71 (MQ=255)
      cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                                >  1:1340211/1‑71 (MQ=255)
       aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCg                                                               >  1:4966952/1‑71 (MQ=255)
         gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                             <  1:3352391/71‑1 (MQ=255)
         gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                             <  1:4152686/71‑1 (MQ=255)
          tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                       <  1:4406794/60‑1 (MQ=255)
          tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                              <  1:3044233/69‑1 (MQ=255)
            tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                           >  1:4797251/1‑70 (MQ=255)
            tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                           >  1:72491/1‑70 (MQ=255)
            tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                           >  1:3743582/1‑70 (MQ=255)
             gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                       <  1:1365950/57‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                        >  1:1040433/1‑71 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                         <  1:919663/70‑1 (MQ=255)
               cATAAATGACCGCTATGAGCACTTCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                          >  1:3273029/1‑68 (MQ=255)
               cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                           >  1:530826/1‑67 (MQ=255)
               cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                       >  1:3202880/1‑71 (MQ=255)
                aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                       <  1:1693065/54‑1 (MQ=255)
                 tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                       >  1:3136743/1‑69 (MQ=255)
                  aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                       >  1:4453092/1‑68 (MQ=255)
                  aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                             >  1:3400781/1‑62 (MQ=255)
                  aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                            >  1:5100756/1‑63 (MQ=255)
                  aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                        >  1:5103432/1‑67 (MQ=255)
                  aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                       >  1:962658/1‑68 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                                    <  1:3971250/70‑1 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                       >  1:2067846/1‑67 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAg                                                    <  1:1896971/70‑1 (MQ=255)
                   aaTGACCGCTATGAGAACTGCAATTACACGCCAGATCGTTCTCGATAACGAAACCACCGGTATg                                                          >  1:2303565/1‑64 (MQ=255)
                    aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                                  <  1:3546660/71‑1 (MQ=255)
                       aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                                <  1:4726805/70‑1 (MQ=255)
                       aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                                <  1:3659258/70‑1 (MQ=255)
                       aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                                <  1:2033456/70‑1 (MQ=255)
                         cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGt                                              >  1:1035585/1‑70 (MQ=255)
                          gCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgc                                            >  1:3141872/1‑71 (MQ=255)
                           cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                                <  1:1536324/66‑1 (MQ=255)
                           cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                       <  1:3568204/43‑1 (MQ=255)
                           cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                       <  1:2277325/43‑1 (MQ=255)
                           cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                       <  1:3659510/43‑1 (MQ=255)
                           cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                           >  1:1820275/1‑71 (MQ=255)
                              tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcgc                                          <  1:2698835/69‑1 (MQ=255)
                               gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                        >  1:4814100/1‑54 (MQ=255)
                               gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                        >  1:3394014/1‑70 (MQ=255)
                               gAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGt                                              >  1:3903466/1‑64 (MQ=255)
                                aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                        >  1:2239713/1‑53 (MQ=255)
                                   aCTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                   <  1:2692664/71‑1 (MQ=255)
                                    cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    >  1:1645950/1‑69 (MQ=255)
                                     tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                                <  1:460367/56‑1 (MQ=255)
                                     tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACt                                       <  1:1234144/65‑1 (MQ=255)
                                     tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:3560799/68‑1 (MQ=255)
                                      gCAATTACACGCCAGATCGTTCTCGCTACCGCAacc                                                                   <  1:553995/36‑1 (MQ=38)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                               >  1:373530/1‑71 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATTAACCAGATTGGTGCGCACTATGa                                   <  1:4201377/67‑1 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                               >  1:4714963/1‑71 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                               >  1:1288945/1‑71 (MQ=255)
                                           tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                             >  1:2916474/1‑69 (MQ=255)
                                            acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                                 >  1:2165973/1‑64 (MQ=255)
                                            acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                            >  1:2594865/1‑69 (MQ=255)
                                                gCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                       >  1:1909832/1‑70 (MQ=255)
                                                gCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                      <  1:346444/71‑1 (MQ=255)
                                                gCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                      <  1:900980/71‑1 (MQ=255)
                                                  cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt                    >  1:4595119/1‑71 (MQ=255)
                                                    gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:4875204/1‑70 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                     <  1:2908946/67‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:4319253/71‑1 (MQ=255)
                                                      tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:2826866/1‑68 (MQ=255)
                                                      tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAATATGAAGGCCACAAGATCAtt                   >  1:2468854/1‑68 (MQ=255)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:4458523/1‑67 (MQ=255)
                                                       cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                  >  1:1019230/1‑52 (MQ=255)
                                                       cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:3947573/1‑67 (MQ=255)
                                                       cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:2122048/1‑67 (MQ=255)
                                                       cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt                    >  1:4162406/1‑66 (MQ=255)
                                                       cGTTCTCGATACCGAAACCACCGGTATAAACCAGATTGGTGCGCACTATGaa                                  >  1:1494794/1‑52 (MQ=255)
                                                        gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt              <  1:1404500/71‑1 (MQ=255)
                                                         ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                  >  1:1500018/1‑50 (MQ=255)
                                                         ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   <  1:850199/65‑1 (MQ=255)
                                                          tctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                   <  1:2421586/48‑1 (MQ=255)
                                                          tctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                   <  1:4895987/48‑1 (MQ=255)
                                                           ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg            <  1:3329428/70‑1 (MQ=255)
                                                           ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg            <  1:2850074/70‑1 (MQ=255)
                                                           ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg            <  1:534609/70‑1 (MQ=255)
                                                            tcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                     <  1:1336279/44‑1 (MQ=38)
                                                            tcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt          >  1:1125068/1‑71 (MQ=255)
                                                            tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                  >  1:3578452/1‑47 (MQ=255)
                                                            tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                                 >  1:809580/1‑48 (MQ=255)
                                                            tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag                <  1:4820505/65‑1 (MQ=255)
                                                             cgctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg            >  1:3819446/4‑68 (MQ=255)
                                                             cGATACCGAAACCACCGGTATGAACCAGATTGGTg                                             >  1:4334957/1‑35 (MQ=255)
                                                             cGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTg         >  1:4529244/1‑71 (MQ=255)
                                                              gctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg           <  1:1253573/66‑1 (MQ=255)
                                                              gATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                            <  1:1596230/51‑1 (MQ=255)
                                                                tACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCg      >  1:2782123/1‑71 (MQ=255)
                                                                tACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCg      >  1:3887090/1‑71 (MQ=255)
                                                                 aCCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgaga               >  1:1984894/1‑61 (MQ=255)
                                                                    gAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTTCCGTTGa  <  1:764834/71‑1 (MQ=255)
                                                                     |                                                                     
AAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGA  >  AP009048/236038‑236176

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: