Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 236,101 A→C 24.3% 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
*AP009048236,1010AC24.3% 180.4 / 114.8 144D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (54/55);  new base C (22/13);  total (76/68)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

AACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCC  >  AP009048/236031‑236171
                                                                      |                                                                      
 aCTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGAt                                                                       <  1:840110/71‑1 (MQ=255)
  cTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATa                                                                      >  1:3183929/1‑71 (MQ=255)
  cTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATa                                                                      >  1:2019389/1‑71 (MQ=255)
  cTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATa                                                                      >  1:919177/1‑71 (MQ=255)
     tAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAcc                                                                    <  1:3709360/70‑1 (MQ=255)
      tacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACcgc                                                                  <  1:4224197‑M2/64‑2 (MQ=255)
      aaaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAcc                                                                    <  1:868003/69‑1 (MQ=255)
       aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                                 <  1:4373971/71‑1 (MQ=255)
       aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                                 <  1:1971051/71‑1 (MQ=255)
       aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                  <  1:3691367/70‑1 (MQ=255)
       aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                  <  1:1213796/70‑1 (MQ=255)
        aaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                  <  1:3531761/69‑1 (MQ=255)
         aactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                               >  1:1115160‑M2/5‑71 (MQ=255)
         aCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                               >  1:2778797/1‑71 (MQ=255)
         aCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                               >  1:1824402/1‑71 (MQ=255)
            aCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGAt                                                                       <  1:4926694/60‑1 (MQ=255)
             cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                           >  1:1340211/1‑71 (MQ=255)
                gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                         >  1:2831941/1‑70 (MQ=255)
                gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGt                                                        <  1:90282/71‑1 (MQ=255)
                gTCTGACATAAATGACCGCTATGAGCACTGAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                       <  1:123260/71‑1 (MQ=255)
                 tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                         <  1:3044233/69‑1 (MQ=255)
                  cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                       <  1:4337263/70‑1 (MQ=255)
                   tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                  <  1:1360121/58‑1 (MQ=255)
                   tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                      >  1:4797251/1‑70 (MQ=255)
                   tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                      >  1:3743582/1‑70 (MQ=255)
                   tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                      >  1:72491/1‑70 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                   >  1:1040433/1‑71 (MQ=255)
                      cATAAATGACCGCTATGAGCACTTCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                     >  1:3273029/1‑68 (MQ=255)
                      cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                      >  1:530826/1‑67 (MQ=255)
                        tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                        <  1:5090982/47‑1 (MQ=255)
                        tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                                 >  1:4169310/1‑52 (MQ=255)
                        tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                      >  1:541997/1‑65 (MQ=255)
                        tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                  >  1:3136743/1‑69 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                                 >  1:2061747/1‑51 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:4453092/1‑68 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:3582983/1‑68 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                       >  1:5100756/1‑63 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                  >  1:962658/1‑68 (MQ=255)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                               <  1:3971250/70‑1 (MQ=255)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAg                                               <  1:1896971/70‑1 (MQ=255)
                           aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                            >  1:1660693/1‑56 (MQ=255)
                           aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                             >  1:46989/1‑71 (MQ=255)
                           aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                             <  1:3546660/71‑1 (MQ=255)
                            tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                            >  1:3977921/1‑71 (MQ=255)
                              aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                           <  1:4726805/70‑1 (MQ=255)
                              aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                           <  1:3659258/70‑1 (MQ=255)
                              aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                           <  1:2033456/70‑1 (MQ=255)
                                cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                                <  1:2909846/63‑1 (MQ=255)
                                cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGt                                         >  1:1035585/1‑70 (MQ=255)
                                  cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                      >  1:3935532/1‑39 (MQ=255)
                                  cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAc                                                                     >  1:1389577/1‑40 (MQ=255)
                                  cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                  <  1:984706/43‑1 (MQ=255)
                                  cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                  <  1:3568204/43‑1 (MQ=255)
                                  cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                  <  1:2277325/43‑1 (MQ=255)
                                  cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                      >  1:1820275/1‑71 (MQ=255)
                                    aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                    >  1:3383492/1‑71 (MQ=255)
                                     tGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                  <  1:1296628/56‑1 (MQ=255)
                                      gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                   >  1:4814100/1‑54 (MQ=255)
                                      gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                   >  1:3394014/1‑70 (MQ=255)
                                       aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                   >  1:2239713/1‑53 (MQ=255)
                                          aCTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                              <  1:2692664/71‑1 (MQ=255)
                                           cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               >  1:1645950/1‑69 (MQ=255)
                                            tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:3560799/68‑1 (MQ=255)
                                             gCAATTACACGCCAGATCGTTCTCGCTACCGCAacc                                                              <  1:553995/36‑1 (MQ=38)
                                              cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                          >  1:4714963/1‑71 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                          <  1:2690919/71‑1 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                          >  1:1288945/1‑71 (MQ=255)
                                                 ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                       <  1:2711504/71‑1 (MQ=255)
                                                 ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                       <  1:1581103/71‑1 (MQ=255)
                                                  tACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                         >  1:1861067/1‑36 (MQ=255)
                                                  tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                        >  1:2916474/1‑69 (MQ=255)
                                                    cacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                    <  1:490193/71‑1 (MQ=255)
                                                    cacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                    <  1:4410177/71‑1 (MQ=255)
                                                       gCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                 <  1:346444/71‑1 (MQ=255)
                                                         cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt               >  1:4595119/1‑71 (MQ=255)
                                                          aGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              <  1:1934268/71‑1 (MQ=255)
                                                           gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              >  1:4875204/1‑70 (MQ=255)
                                                            aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt               >  1:630420/1‑68 (MQ=255)
                                                            aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga            <  1:4319253/71‑1 (MQ=255)
                                                             tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAATATGAAGGCCACAAGATCAtt              >  1:2468854/1‑68 (MQ=255)
                                                              cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              >  1:4871725/1‑67 (MQ=255)
                                                              cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:3182348/50‑1 (MQ=255)
                                                              cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                             >  1:1019230/1‑52 (MQ=255)
                                                              cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              >  1:212275/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)
                                                               gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgaga          <  1:3944323/70‑1 (MQ=255)
                                                                ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                             >  1:1500018/1‑50 (MQ=255)
                                                                ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              <  1:850199/65‑1 (MQ=255)
                                                                 tctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                              <  1:2421586/48‑1 (MQ=255)
                                                                  ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg       <  1:3329428/70‑1 (MQ=255)
                                                                  ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg       <  1:528430/70‑1 (MQ=255)
                                                                  ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg       <  1:2850074/70‑1 (MQ=255)
                                                                   tcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt     >  1:1125068/1‑71 (MQ=255)
                                                                   tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                            >  1:809580/1‑48 (MQ=255)
                                                                    cGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTg    >  1:4529244/1‑71 (MQ=255)
                                                                     gATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                       <  1:1596230/51‑1 (MQ=255)
                                                                      aTACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGcc  <  1:3810079/71‑1 (MQ=255)
                                                                      |                                                                      
AACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCC  >  AP009048/236031‑236171

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: