Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 236,101 A→C 79.8% 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,1010AC79.8% 1885.2 / 278.4 681D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (70/68);  new base C (293/249);  total (364/317)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

Only 100 of 689 total aligned reads displayed.

ACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCC  >  AP009048/236032‑236171
                                                                     |                                                                      
aCTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                        >  1:2622251/1‑70 (MQ=255)
  gcgtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAc                                                                     >  1:2255740‑M2/11‑71 (MQ=255)
    gtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCg                                                                   >  1:6949643‑M2/9‑71 (MQ=255)
    gtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCg                                                                   >  1:6164143‑M2/9‑71 (MQ=255)
       caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                >  1:3007664‑M2/6‑71 (MQ=255)
       caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                >  1:4753763‑M2/6‑71 (MQ=255)
       caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                >  1:5468908‑M2/6‑71 (MQ=255)
         cGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAacc                                                              <  1:6610053/71‑1 (MQ=255)
          ctCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                      <  1:2369407‑M2/60‑1 (MQ=255)
          ctCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                             >  1:1007987‑M2/3‑71 (MQ=255)
          ctCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                             >  1:6059205‑M2/3‑71 (MQ=255)
           tCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAccgca                                                                 <  1:94161‑M2/65‑3 (MQ=255)
            cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccacc                                                           <  1:4450826/71‑1 (MQ=255)
             aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                          <  1:4168394/71‑1 (MQ=255)
             aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                          >  1:3562732/1‑71 (MQ=255)
             aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                          >  1:2724213/1‑71 (MQ=255)
             aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccac                                                            >  1:5850072/1‑69 (MQ=255)
             aaGTCTGACATAAATGACCGCTATGAGCACTCCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                           <  1:5887597/70‑1 (MQ=255)
              aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                         <  1:5155581/71‑1 (MQ=255)
               gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                                 <  1:6837388/62‑3 (MQ=255)
               gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                         >  1:7107839/1‑70 (MQ=255)
                  tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                      >  1:7019119/1‑70 (MQ=255)
                  tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                      >  1:1098055/1‑70 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                   >  1:1927370/1‑71 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                   >  1:4283961/1‑71 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                    <  1:7627486/70‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                    <  1:146536/70‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                    <  1:4910447/70‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                   >  1:6453366/1‑71 (MQ=255)
                      aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                               <  1:6685854/57‑1 (MQ=255)
                      aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                               <  1:3984915/57‑1 (MQ=255)
                       tAAATTACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                <  1:6086221/71‑1 (MQ=255)
                       tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                <  1:2061799/71‑1 (MQ=255)
                        aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCg                                                                   <  1:4982398/51‑1 (MQ=255)
                        aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                             >  1:3867051/1‑57 (MQ=255)
                        aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:6262299/1‑68 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCGGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:5228035/1‑67 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                      >  1:640929/1‑63 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:2437997/1‑67 (MQ=255)
                          aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                         <  1:2925599/59‑1 (MQ=255)
                          aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                             <  1:2330768/71‑1 (MQ=255)
                            gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                           <  1:850088/55‑1 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                                 <  1:1612060/48‑3 (MQ=255)
                               cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCg                                                                   >  1:4531475/1‑44 (MQ=255)
                                gCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgc                                       >  1:3435216/1‑71 (MQ=255)
                                 cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                      >  1:7891808/1‑71 (MQ=255)
                                   aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                       >  1:4702941/1‑52 (MQ=255)
                                   aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                    >  1:1984573/1‑71 (MQ=255)
                                    tGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                                   <  1:1756076/71‑1 (MQ=255)
                                     gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                >  1:4095601/1‑57 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                >  1:1147802/1‑40 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                >  1:5042650/1‑40 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAcc                                                 <  1:1560773/55‑1 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                 <  1:7275392/71‑1 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                 <  1:1806257/71‑1 (MQ=255)
                                        cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:203284/71‑1 (MQ=255)
                                        cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:7060756/71‑1 (MQ=255)
                                        cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:5308080/71‑1 (MQ=255)
                                          cTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                >  1:2151342/1‑36 (MQ=38)
                                          cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                       <  1:4660062/45‑1 (MQ=38)
                                           tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                          >  1:4600539/1‑41 (MQ=255)
                                           tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                     <  1:393394/46‑1 (MQ=39)
                                           tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                <  1:6224793/67‑1 (MQ=255)
                                            gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                           >  1:3093526/1‑71 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCTCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                          <  1:274497/71‑1 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                           >  1:899808/1‑70 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                          >  1:5773138/1‑71 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTa                                 >  1:691255/1‑64 (MQ=255)
                                             aaaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                     >  1:5664377/2‑44 (MQ=38)
                                               aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:1253268/1‑45 (MQ=38)
                                               aTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                        <  1:745732/71‑1 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                <  1:3169465/62‑1 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               >  1:1317819/1‑63 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               >  1:6742845/1‑63 (MQ=255)
                                                  acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                    <  1:5103860/40‑1 (MQ=38)
                                                  acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:556902/1‑42 (MQ=38)
                                                  acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                       >  1:7207730/1‑69 (MQ=255)
                                                   cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                >  1:6004346/1‑43 (MQ=38)
                                                   cacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGTTTGGTGCGCACTATGAAGGCCACaa                     <  1:4365580/70‑1 (MQ=255)
                                                   cacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTa                                 <  1:1838047/58‑1 (MQ=255)
                                                    acGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                      >  1:3714299/1‑52 (MQ=255)
                                                      gCCAGATCGTTTTCGCTACCGCAACCACCGGTATGAACCAGAtt                                            <  1:4217774/44‑1 (MQ=25)
                                                         aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAATATCAtt              <  1:3354295/71‑1 (MQ=255)
                                                         aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              <  1:7567058/71‑1 (MQ=255)
                                                          gATCGTTCTCGCTACCGCAACCACCTGTATGAACCAGATTGGTGCGCACTATg                               <  1:3231219/53‑1 (MQ=255)
                                                          gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTTCGCACTATGaa                             >  1:3768270/1‑55 (MQ=255)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                <  1:7765007/35‑1 (MQ=37)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:2976309/52‑1 (MQ=255)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga            <  1:2846552/71‑1 (MQ=255)
                                                            tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              >  1:2528684/1‑68 (MQ=255)
                                                             cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              >  1:3636719/1‑67 (MQ=255)
                                                              gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt         <  1:3527767/71‑1 (MQ=255)
                                                              gTTATCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgaga          <  1:2671277/70‑1 (MQ=255)
                                                               ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTg             <  1:635297/66‑1 (MQ=255)
                                                               ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt        >  1:228117/1‑71 (MQ=255)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg      >  1:3467980/1‑71 (MQ=255)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg       <  1:7331915/70‑1 (MQ=255)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg       <  1:1894453/70‑1 (MQ=255)
                                                                  tcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt               <  1:1689959/61‑1 (MQ=255)
                                                                  tcGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                      <  1:5376142/38‑1 (MQ=255)
                                                                     |                                                                      
ACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCC  >  AP009048/236032‑236171

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: