Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 214,931 A→C 54.0% 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
*minE214,9310AC54.0% 253.7 / inf 290E14A (GAA→GCA) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (60/73);  new base C (90/66);  total (150/140)
Fisher's exact test for biased strand distribution p-value = 3.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 306 total aligned reads displayed.

AAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGA  >  minE/214862‑215000
                                                                     |                                                                     
    ctcAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                                  >  1:2269417/3‑71 (MQ=255)
          tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                             >  1:3075560/1‑70 (MQ=255)
           cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                            >  1:4472036/1‑70 (MQ=255)
            tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                                      <  1:3549345/59‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                          <  1:1682569/69‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                         <  1:1547752/70‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaaa                                                                     >  1:3518301/1‑58 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGACAGATCGTTCTCGATACCGaaa                                                                     >  1:2768778/1‑58 (MQ=255)
                 tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAACGCAACCACCGGTATGAACCa                                                     <  1:3853146/71‑1 (MQ=255)
                 tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                                     <  1:88284/71‑1 (MQ=255)
                  aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                               <  1:3605741/60‑1 (MQ=255)
                  aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                       >  1:1061314/1‑68 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                        >  1:1260779/1‑66 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                                   <  1:3441056/71‑1 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                                   <  1:2724982/71‑1 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccacc                                                                >  1:2502317/1‑58 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                       >  1:1800533/1‑67 (MQ=255)
                       aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                               <  1:4385137/71‑1 (MQ=255)
                              tGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                                        >  1:1303076/1‑71 (MQ=255)
                              tGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                                        >  1:4545965/1‑71 (MQ=255)
                                aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                      >  1:4400338/1‑71 (MQ=255)
                                aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                            >  1:3266374/1‑49 (MQ=255)
                                aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                            >  1:1843728/1‑49 (MQ=255)
                                  cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:1410363/71‑1 (MQ=255)
                                  cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                                     <  1:782337/54‑1 (MQ=255)
                                  cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                                     <  1:604813/54‑1 (MQ=255)
                                  cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCCGTATGAACCa                                                     <  1:3529917/54‑1 (MQ=255)
                                   aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:2480391/70‑1 (MQ=255)
                                   aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                   <  1:3629662/71‑1 (MQ=255)
                                    cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    >  1:3649915/1‑69 (MQ=255)
                                    cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    >  1:3389041/1‑69 (MQ=255)
                                    cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                  <  1:3484710/71‑1 (MQ=255)
                                     tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                                     <  1:587276/51‑1 (MQ=255)
                                     tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                                     <  1:3880029/51‑1 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                          <  1:2369952/44‑1 (MQ=38)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                               >  1:2154318/1‑71 (MQ=255)
                                          ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                                 <  1:4630093/50‑1 (MQ=255)
                                          ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                                 <  1:114360/50‑1 (MQ=255)
                                          ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGt                                              <  1:4010577/53‑1 (MQ=255)
                                          ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                            <  1:1277137/71‑1 (MQ=255)
                                          ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                            <  1:4252025/71‑1 (MQ=255)
                                           tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:396571/62‑1 (MQ=255)
                                           tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:1437303/62‑1 (MQ=255)
                                               cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                            >  1:73192/1‑66 (MQ=255)
                                               cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                            >  1:4299690/1‑66 (MQ=255)
                                               cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                         <  1:4026742/53‑1 (MQ=255)
                                               cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                         <  1:1398073/53‑1 (MQ=255)
                                               cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                            >  1:2124053/1‑66 (MQ=255)
                                               cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                       <  1:3784037/71‑1 (MQ=255)
                                                gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                       >  1:1106382/1‑70 (MQ=255)
                                                gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                       >  1:2746663/1‑70 (MQ=255)
                                                gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                       >  1:115725/1‑70 (MQ=255)
                                                gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                      <  1:821815/71‑1 (MQ=255)
                                                 ccAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                     <  1:4716835/71‑1 (MQ=255)
                                                  cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                     <  1:3818897/70‑1 (MQ=255)
                                                   aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   <  1:1711942/71‑1 (MQ=255)
                                                   aGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:3347400/1‑71 (MQ=255)
                                                    gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:302036/1‑70 (MQ=255)
                                                    gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:2702549/1‑70 (MQ=255)
                                                    gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTg                                             <  1:960582/44‑1 (MQ=255)
                                                    gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTg                                             <  1:1822897/44‑1 (MQ=255)
                                                    gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                      >  1:552022/1‑67 (MQ=255)
                                                     aTCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                     <  1:30936/35‑1 (MQ=37)
                                                     aTCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                     <  1:1932786/35‑1 (MQ=37)
                                                     aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                     <  1:1244178/51‑1 (MQ=255)
                                                     aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                          >  1:1045023/1‑62 (MQ=255)
                                                     aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:2944849/71‑1 (MQ=255)
                                                     aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:4518012/71‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:2888843/71‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:2223393/71‑1 (MQ=255)
                                                      tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTg                                             >  1:2995009/1‑42 (MQ=38)
                                                      tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:4190544/1‑68 (MQ=255)
                                                      tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:4489747/70‑1 (MQ=255)
                                                      tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:4610914/70‑1 (MQ=255)
                                                      tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                           >  1:4136350/1‑44 (MQ=255)
                                                      tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:3994627/1‑68 (MQ=255)
                                                      tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTg                  <  1:1878793/69‑1 (MQ=255)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:1499561/1‑67 (MQ=255)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:685851/1‑67 (MQ=255)
                                                        gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt              >  1:2092974/1‑71 (MQ=255)
                                                        gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt              >  1:1964335/1‑71 (MQ=255)
                                                           ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                       >  1:2674113/1‑59 (MQ=255)
                                                           ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                       >  1:457599/1‑59 (MQ=255)
                                                           ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                                 >  1:2208173/1‑49 (MQ=255)
                                                           ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                                 >  1:910389/1‑49 (MQ=255)
                                                            tcGCTACCGCAACCACCGGTATTAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt          >  1:3186523/1‑71 (MQ=255)
                                                            tcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt          >  1:4208024/1‑71 (MQ=255)
                                                            tcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt          >  1:1848480/1‑71 (MQ=255)
                                                              gctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGc        >  1:2533951/3‑71 (MQ=255)
                                                              gctACCGCAACCACCGGTATGAACCAGATTGCTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:2624818/3‑60 (MQ=255)
                                                               ctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGcc       >  1:1464029/2‑71 (MQ=255)
                                                                  ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCg      <  1:1731585/69‑1 (MQ=255)
                                                                  ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGtt    >  1:504640/1‑71 (MQ=255)
                                                                  ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGtt    <  1:2068242/71‑1 (MQ=255)
                                                                  ccGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGtt    <  1:1657209/71‑1 (MQ=255)
                                                                  ccGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGtt    <  1:4262048/71‑1 (MQ=255)
                                                                    gAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGTCACAAGATCATTGAGATTGGTGCCGTTGa  <  1:3172012/71‑1 (MQ=255)
                                                                    gAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgaga               >  1:495181/1‑58 (MQ=255)
                                                                    gAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTg   >  1:1604516/1‑70 (MQ=255)
                                                                     aaaCCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGtt    <  1:4746039/68‑1 (MQ=255)
                                                                     |                                                                     
AAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGA  >  minE/214862‑215000

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: