Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 236,101 A→C 59.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
*AP009048236,1010AC59.5% 238.3 / 158.1 170D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (34/35);  new base C (63/38);  total (97/73)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 175 total aligned reads displayed.

TGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCC  >  AP009048/236034‑236171
                                                                   |                                                                      
tGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAc                                                                     >  1:2647893/1‑71 (MQ=255)
gcgtacaactCAAGTCTGACATAAATGCCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                      >  1:2171781‑M2/11‑70 (MQ=255)
gcgtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAc                                                                     <  1:171720‑M2/61‑1 (MQ=255)
    acaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACcgc                                                                  >  1:804068‑M2/7‑69 (MQ=255)
    aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                  >  1:1808352/1‑70 (MQ=255)
    aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                  >  1:1102077/1‑70 (MQ=255)
    aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                  <  1:2882324/70‑1 (MQ=255)
    aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                  <  1:1985856/70‑1 (MQ=255)
       actCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                      <  1:2107917‑M2/60‑1 (MQ=255)
        gACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                        <  1:523075/60‑1 (MQ=255)
        ctCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                             >  1:1631163‑M2/3‑71 (MQ=255)
          cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                           >  1:1018888/1‑71 (MQ=255)
           aaGTCTGACATAAATGACCGCTATGAGCACTGCATTTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                          >  1:776802/1‑71 (MQ=255)
           aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                          <  1:56254/71‑1 (MQ=255)
            aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                         <  1:1758449/71‑1 (MQ=255)
             gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                                 <  1:2035251/62‑1 (MQ=255)
               cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                      >  1:3016513/1‑71 (MQ=255)
                tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                        >  1:3536296/1‑52 (MQ=255)
                tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                      >  1:2072413/1‑70 (MQ=255)
                 gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                     >  1:735874/1‑70 (MQ=255)
                  aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                   >  1:3612071/1‑71 (MQ=255)
                  aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                   >  1:1923054/1‑71 (MQ=255)
                  aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                   <  1:2906242/71‑1 (MQ=255)
                   cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                   <  1:3974689/70‑1 (MQ=255)
                   cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                                 <  1:1879911/56‑1 (MQ=255)
                    aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:128386/1‑70 (MQ=255)
                     tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                <  1:3308436/71‑1 (MQ=255)
                      aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:2949788/1‑68 (MQ=255)
                      aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                  >  1:687152/1‑68 (MQ=255)
                       aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:645049/1‑67 (MQ=255)
                       aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:4468703/1‑67 (MQ=255)
                       aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                  >  1:1523601/1‑67 (MQ=255)
                       aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                  >  1:2532543/1‑67 (MQ=255)
                        aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                        >  1:4139563/1‑44 (MQ=255)
                        aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAc                                                                     >  1:3097270/1‑47 (MQ=255)
                        aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                             >  1:4556471/1‑71 (MQ=255)
                        aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccac                                                            >  1:4548287/1‑56 (MQ=255)
                         tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:4455319/1‑65 (MQ=255)
                         tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                            <  1:775128/71‑1 (MQ=255)
                          gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                      >  1:1934633/1‑60 (MQ=255)
                          gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                                <  1:4076712/66‑1 (MQ=255)
                            ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                  >  1:33364/1‑62 (MQ=255)
                             cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                                 <  1:4001402/46‑3 (MQ=255)
                             cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                  <  1:4120824/45‑1 (MQ=255)
                              gCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgc                                       <  1:2727459/71‑1 (MQ=255)
                               cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                >  1:1163104/1‑45 (MQ=255)
                               cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                      <  1:892945/71‑1 (MQ=255)
                                tATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                      >  1:1031431/1‑70 (MQ=255)
                                 aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                    >  1:4412412/1‑71 (MQ=255)
                                 aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                    >  1:1754475/1‑71 (MQ=255)
                                  tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                            >  1:4287420/1‑46 (MQ=255)
                                  tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                               <  1:3578573/59‑1 (MQ=255)
                                  tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                   >  1:473913/1‑71 (MQ=255)
                                  tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                   <  1:612194/71‑1 (MQ=255)
                                   gAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                                   <  1:21560/70‑1 (MQ=255)
                                    aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                      <  1:320037/66‑1 (MQ=255)
                                    aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTg                                        >  1:2591869/1‑64 (MQ=255)
                                    aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                 >  1:2930062/1‑71 (MQ=255)
                                      cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:3181781/71‑1 (MQ=255)
                                      cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:3704378/71‑1 (MQ=255)
                                       aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                              <  1:2629249/71‑1 (MQ=255)
                                       aCTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                              <  1:3058244/71‑1 (MQ=255)
                                         tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                             <  1:2449769/54‑1 (MQ=255)
                                          gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                            >  1:2517657/1‑70 (MQ=255)
                                           cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                   >  1:314312/1‑46 (MQ=255)
                                           cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                          >  1:3677408/1‑71 (MQ=255)
                                           cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                          >  1:1080139/1‑71 (MQ=255)
                                            aaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                               <  1:1226935/49‑1 (MQ=39)
                                            aaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                         <  1:390638/71‑1 (MQ=255)
                                            aaTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                <  1:433407/64‑1 (MQ=255)
                                             aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                        >  1:426152/1‑71 (MQ=255)
                                             aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                        >  1:2247696/1‑71 (MQ=255)
                                             aTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                        <  1:2650592/71‑1 (MQ=255)
                                              ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGAGCACTATGAAGGCcac                       <  1:517851/71‑1 (MQ=255)
                                              ttACACGCCAGATCGTTCTCGCTACCGCAACCACAGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                        >  1:771361/1‑70 (MQ=255)
                                              ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                                   <  1:3025376/59‑1 (MQ=255)
                                                acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                              >  1:150265/1‑62 (MQ=255)
                                                acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                       >  1:2816492/1‑69 (MQ=255)
                                                acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                       >  1:4495341/1‑69 (MQ=255)
                                                acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                   >  1:1148637/1‑41 (MQ=255)
                                                acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                        >  1:3851087/1‑68 (MQ=255)
                                                 cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                           <  1:3632635/48‑1 (MQ=39)
                                                 cacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                    <  1:3287371/71‑1 (MQ=255)
                                                      cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                      <  1:495818/48‑1 (MQ=39)
                                                      cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                      <  1:3465049/48‑1 (MQ=39)
                                                      cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt               >  1:2090051/1‑71 (MQ=255)
                                                      cAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                        >  1:4536423/1‑62 (MQ=255)
                                                        gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                               <  1:4315388/37‑1 (MQ=37)
                                                           cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              >  1:1620452/1‑67 (MQ=255)
                                                            gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt         <  1:882124/71‑1 (MQ=255)
                                                            gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt         <  1:1314902/71‑1 (MQ=255)
                                                               ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg      >  1:1957147/1‑71 (MQ=255)
                                                               ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg      >  1:925885/1‑71 (MQ=255)
                                                               ctcGATACCGAAACCACCGGTATGAACCAGATTgg                                          >  1:395255/1‑35 (MQ=255)
                                                               ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg       <  1:1301796/70‑1 (MQ=255)
                                                               ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg       <  1:817471/70‑1 (MQ=255)
                                                                tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt     >  1:1357378/1‑71 (MQ=255)
                                                                  gATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg      <  1:4148110/68‑1 (MQ=255)
                                                                  gATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt     <  1:288025/69‑1 (MQ=255)
                                                                   aTACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGcc  >  1:941017/1‑71 (MQ=255)
                                                                   |                                                                      
TGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCC  >  AP009048/236034‑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: