Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 236,101 A→C 36.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,1010AC36.3% 43.0 / 195.4 146D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (49/44);  new base C (31/22);  total (80/66)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 151 total aligned reads displayed.

AACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCC  >  AP009048/236031‑236171
                                                                      |                                                                      
 aCTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGAt                                                                       <  1:4427713/71‑1 (MQ=255)
 aCTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGAt                                                                       <  1:568009/71‑1 (MQ=255)
 aCTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGAt                                                                       <  1:5274709/71‑1 (MQ=255)
  cTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATa                                                                      >  1:4627545/1‑71 (MQ=255)
  cTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATa                                                                      >  1:4966632/1‑71 (MQ=255)
   gcgtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAc                                                                     >  1:4656438‑M2/11‑71 (MQ=255)
    gTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAcc                                                                    >  1:3950360/1‑71 (MQ=255)
       acaactCAAGTCTGACATAAATGCCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACcgc                                                                  >  1:1481654‑M2/7‑69 (MQ=255)
       aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                  <  1:2661479/70‑1 (MQ=255)
        caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                >  1:1484404‑M2/6‑71 (MQ=255)
            aCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccac                                                            <  1:5005206/71‑1 (MQ=255)
             cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                           >  1:1410067/1‑71 (MQ=255)
               aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                         <  1:1517955/71‑1 (MQ=255)
                   tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                      >  1:1629689/1‑70 (MQ=255)
                   tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                                 <  1:1247389/59‑1 (MQ=255)
                    gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                          >  1:914529/1‑65 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                   >  1:1886672/1‑71 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                     <  1:2115568/69‑1 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                     <  1:245/69‑1 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                   >  1:4883588/1‑71 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                   >  1:642514/1‑71 (MQ=255)
                      cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  <  1:4993218/71‑1 (MQ=255)
                       aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:1047743/1‑70 (MQ=255)
                        tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTAGATACCGAAACCACCGGTATGAACCa                                                <  1:2295115/71‑1 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                  >  1:5311774/1‑68 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                  >  1:2552045/1‑68 (MQ=255)
                         aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                  >  1:5207321/1‑68 (MQ=255)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                            >  1:2475610/1‑57 (MQ=255)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                  >  1:2371641/1‑67 (MQ=255)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAg                                               <  1:1709697/70‑1 (MQ=255)
                             gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                        >  1:4350451/1‑42 (MQ=255)
                             gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                           >  1:3333962/1‑71 (MQ=255)
                              aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                          <  1:3926194/71‑1 (MQ=255)
                              aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                           <  1:3023788/70‑1 (MQ=255)
                                  cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                      >  1:4654405/1‑71 (MQ=255)
                                  cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                      >  1:293704/1‑71 (MQ=255)
                                    aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                         >  1:4623502/1‑66 (MQ=255)
                                    aTGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                    >  1:2086393/1‑71 (MQ=255)
                                     tGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                           <  1:13514/63‑1 (MQ=255)
                                     tGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                           <  1:3439646/63‑1 (MQ=255)
                                     tGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                    <  1:2847740/70‑1 (MQ=255)
                                        gCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                          <  1:2910471/61‑1 (MQ=255)
                                        gCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                <  1:3503475/71‑1 (MQ=255)
                                         cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCTGTATGAACCAGATTGGTGCGCACTATg                               <  1:3140443/71‑1 (MQ=255)
                                          aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               >  1:4363459/1‑70 (MQ=255)
                                           cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAGTATGaa                             <  1:693717/71‑1 (MQ=255)
                                           cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                             <  1:4316437/71‑1 (MQ=255)
                                            tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                               <  1:3482468/52‑1 (MQ=255)
                                            tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                <  1:4514306/67‑1 (MQ=255)
                                             gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                       >  1:3082133/1‑43 (MQ=38)
                                             gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                         <  1:2932380/41‑1 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCTATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                          <  1:4101585/71‑1 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                          >  1:1518554/1‑71 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                       >  1:4899237/1‑42 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                          >  1:1240302/1‑55 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGt                                         >  1:1909750/1‑56 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                          >  1:335663/1‑71 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                          <  1:1773730/71‑1 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                          >  1:4170204/1‑71 (MQ=255)
                                               aaTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                           <  1:4162716/53‑1 (MQ=255)
                                               aaTTACACGCCAGATCGTTCTCGATACCGAAACCCCGGTATGAACCAGATTg                                           <  1:49064/52‑1 (MQ=255)
                                                aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                        >  1:2323195/1‑71 (MQ=255)
                                                aTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAg                                               <  1:3755736/48‑1 (MQ=255)
                                                 ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:4435891/63‑1 (MQ=255)
                                                 ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                        >  1:4517467/1‑70 (MQ=255)
                                                 ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                   <  1:4237161/43‑1 (MQ=255)
                                                  tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                           <  1:3540739/50‑1 (MQ=255)
                                                   acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                   >  1:5104055/1‑41 (MQ=255)
                                                   acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAcc                                                 >  1:1633684/1‑43 (MQ=255)
                                                   acacGCCAGATAGTTCTCGATACCGAAACCACCGGTATGaa                                                   >  1:3381816/1‑41 (MQ=255)
                                                     acGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                          <  1:2287109/48‑1 (MQ=39)
                                                       gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                 <  1:2588345/71‑1 (MQ=255)
                                                       gCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                 <  1:4865020/71‑1 (MQ=255)
                                                        ccAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCGCTATGAAGGCCACAAGATCa                <  1:2840527/71‑1 (MQ=255)
                                                         cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt               >  1:4214201/1‑71 (MQ=255)
                                                         cAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACt                                  <  1:3737436/52‑1 (MQ=255)
                                                         cAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTa                                 >  1:2607851/1‑53 (MQ=255)
                                                         cAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt               >  1:1671689/1‑71 (MQ=255)
                                                           gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              >  1:4038044/1‑70 (MQ=255)
                                                            aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                  <  1:2207662/65‑1 (MQ=255)
                                                            aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                       <  1:1915405/60‑1 (MQ=255)
                                                             tCGTTCTCGATACCGAACCCACCGGTATGAACCAGATTGGTGCGCACTAt                                >  1:4491397/1‑50 (MQ=255)
                                                             tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                     <  1:628537/61‑1 (MQ=255)
                                                             tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              >  1:1656635/1‑68 (MQ=255)
                                                             tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTg             <  1:4446272/69‑1 (MQ=255)
                                                              cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              >  1:1352024/1‑67 (MQ=255)
                                                              cGTTCTCGATACCTAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              >  1:3714269/1‑67 (MQ=255)
                                                              cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                             >  1:3790500/1‑52 (MQ=255)
                                                               gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt         >  1:519674/1‑71 (MQ=255)
                                                                ttCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                            <  1:1456185/35‑1 (MQ=37)
                                                                ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt        <  1:1010291/71‑1 (MQ=255)
                                                                 tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                 <  1:2486770/61‑1 (MQ=255)
                                                                  ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg       <  1:1704584/70‑1 (MQ=255)
                                                                  ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg      >  1:4964488/1‑71 (MQ=255)
                                                                   tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga            >  1:5081692/1‑64 (MQ=255)
                                                                   tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt     >  1:3555509/1‑71 (MQ=255)
                                                                   tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt     >  1:27838/1‑71 (MQ=255)
                                                                   tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt     >  1:1988949/1‑71 (MQ=255)
                                                                   tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt     >  1:1949552/1‑71 (MQ=255)
                                                                      aTACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGcc  >  1:4285655/1‑71 (MQ=255)
                                                                      |                                                                      
AACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCC  >  AP009048/236031‑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: