Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 214,925 A→C 73.1% 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
*minE214,9250AC73.1% 519.5 / 139.5 245D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (30/36);  new base C (88/91);  total (118/127)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 247 total aligned reads displayed.

ACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTG  >  minE/214856‑214993
                                                                     |                                                                    
aCTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                      <  1:1528231/70‑1 (MQ=255)
            cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                         >  1:753669/1‑71 (MQ=255)
             aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                          <  1:3057111/69‑1 (MQ=255)
             aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                        >  1:2510338/1‑71 (MQ=255)
             aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                        >  1:1385613/1‑71 (MQ=255)
              aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                       <  1:2644889/71‑1 (MQ=255)
              aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                       <  1:1108308/71‑1 (MQ=255)
                 cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                     <  1:1956974/70‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                                <  1:195432/56‑2 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                    <  1:3768357/68‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                    <  1:325098/68‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                   <  1:2879524/69‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                   <  1:2054010/69‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                 <  1:2415973/71‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                      >  1:3254268/1‑66 (MQ=255)
                     cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                <  1:3395824/71‑1 (MQ=255)
                     cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                >  1:2140393/1‑71 (MQ=255)
                       tAAATGACCGCTATGAGCACTGCAATTACAGGCCAGATCGTTCTCGCTACCgca                                                               <  1:1468016/54‑3 (MQ=255)
                       tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                              <  1:3622388/71‑1 (MQ=255)
                        aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                >  1:3490156/1‑68 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                             >  1:2974369/1‑54 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                >  1:3192951/1‑67 (MQ=255)
                          aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAc                                                                   <  1:168064/47‑1 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                        <  1:2023949/71‑1 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                         <  1:2691329/70‑1 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                         <  1:2462283/70‑1 (MQ=255)
                              ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                       <  1:3210722/71‑1 (MQ=255)
                               cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTg                                      >  1:316936/1‑71 (MQ=255)
                               cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                            <  1:2932599/65‑1 (MQ=255)
                                gCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgc                                     >  1:883302/1‑71 (MQ=255)
                                 cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                            <  1:925363/63‑1 (MQ=255)
                                   aTGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                               >  1:1696785/1‑42 (MQ=255)
                                    tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                 <  1:1720855/71‑1 (MQ=255)
                                    tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                 <  1:3678894/71‑1 (MQ=255)
                                    tGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                         <  1:2446184/63‑1 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGAGCACTa                               >  1:2178470/1‑71 (MQ=255)
                                       gCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTa                               <  1:2112001/70‑1 (MQ=255)
                                       gCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                              <  1:1187691/71‑1 (MQ=255)
                                        cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                  <  1:487751/66‑1 (MQ=255)
                                        cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                             <  1:2835364/71‑1 (MQ=255)
                                        cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                             <  1:3842046/71‑1 (MQ=255)
                                        cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                             <  1:3860999/71‑1 (MQ=255)
                                         aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                            <  1:9608/71‑1 (MQ=255)
                                          cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                              >  1:1093074/1‑68 (MQ=255)
                                          cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                           <  1:650709/71‑1 (MQ=255)
                                          cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                           <  1:131509/71‑1 (MQ=255)
                                           tGCAATTACACGCCAGATCGTTCTCGCTACCGCAacc                                                            <  1:3093705/37‑1 (MQ=38)
                                           tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                              <  1:2758983/51‑1 (MQ=255)
                                           tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                          >  1:1011366/1‑71 (MQ=255)
                                           tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                          >  1:824426/1‑71 (MQ=255)
                                           tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAg                                             <  1:2634089/52‑1 (MQ=255)
                                            gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                          >  1:1456727/1‑70 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                  >  1:2385892/1‑45 (MQ=38)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                         <  1:313796/70‑1 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                        >  1:2790228/1‑71 (MQ=255)
                                              aaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                           <  1:344971/51‑1 (MQ=255)
                                              aaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                        >  1:3534980/1‑70 (MQ=255)
                                               aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                      >  1:3367004/1‑71 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                     <  1:1471416/71‑1 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                       <  1:1931240/37‑1 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                     <  1:60248/71‑1 (MQ=255)
                                                  acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                     >  1:178512/1‑69 (MQ=255)
                                                    acGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa                 <  1:1631316/71‑1 (MQ=255)
                                                    acGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa                 <  1:3830718/71‑1 (MQ=255)
                                                       ccAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa              <  1:2538356/71‑1 (MQ=255)
                                                       ccAGATCGTTCTCGATACCGAAACCACGGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa              >  1:2077664/1‑71 (MQ=255)
                                                       ccAGATCGTTCTCGATACCAAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa              <  1:3513246/71‑1 (MQ=255)
                                                        cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt             <  1:396849/71‑1 (MQ=255)
                                                        cAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                  >  1:52779/1‑50 (MQ=255)
                                                        cAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt             >  1:2295879/1‑71 (MQ=255)
                                                         aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt             >  1:358095/1‑70 (MQ=255)
                                                         aGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAACATCAtt            >  1:1799099/1‑71 (MQ=255)
                                                          gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                  >  1:2210791/1‑48 (MQ=39)
                                                          gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                  >  1:2670892/1‑48 (MQ=39)
                                                          gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTg           <  1:2355266/71‑1 (MQ=255)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                          <  1:1211607/39‑1 (MQ=37)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga          <  1:2273476/71‑1 (MQ=255)
                                                           aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa                 >  1:288683/1‑64 (MQ=255)
                                                            tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag         <  1:3041685/71‑1 (MQ=255)
                                                             cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt            >  1:1577954/1‑67 (MQ=255)
                                                              gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt       >  1:894219/1‑71 (MQ=255)
                                                              gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt       <  1:1730875/71‑1 (MQ=255)
                                                              gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt       <  1:1238740/71‑1 (MQ=255)
                                                               ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt      >  1:3871930/1‑71 (MQ=255)
                                                               ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGAGCACTATGAAGGCCACAAGATCATTGAGAtt      <  1:1439737/71‑1 (MQ=255)
                                                               ttCTCGCTACCGAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg     >  1:3009120/1‑71 (MQ=255)
                                                                tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg     >  1:360409/1‑71 (MQ=255)
                                                                tctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                         <  1:1604813/51‑1 (MQ=255)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                                >  1:2250271/1‑43 (MQ=38)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                                >  1:3703089/1‑43 (MQ=38)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                                >  1:1411889/1‑43 (MQ=38)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg    >  1:3926281/1‑71 (MQ=255)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg    >  1:1710883/1‑71 (MQ=255)
                                                                 ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg    >  1:3836418/1‑71 (MQ=255)
                                                                 atcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg    >  1:1865187/2‑71 (MQ=255)
                                                                  tcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt      >  1:1335995/1‑68 (MQ=255)
                                                                  tcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                       <  1:1388290/51‑1 (MQ=255)
                                                                   cGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa              <  1:1742086/59‑1 (MQ=255)
                                                                    gATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTg  >  1:1258388/1‑70 (MQ=255)
                                                                     aTACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt   >  1:323201/1‑68 (MQ=255)
                                                                     |                                                                    
ACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTG  >  minE/214856‑214993

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: