Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,029,444 C→A 31.6% A467A (GCG→GCT parE ← DNA topoisomerase IV, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,029,4440CA31.6% 181.3 / 89.1 142A467A (GCG→GCTparEDNA topoisomerase IV, subunit B
Reads supporting (aligned to +/- strand):  ref base C (40/57);  new base A (27/18);  total (67/75)
Fisher's exact test for biased strand distribution p-value = 4.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

Only 100 of 145 total aligned reads displayed.

CGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAGAGACTTCCCAGG  >  minE/2029378‑2029505
                                                                  |                                                             
cgcgAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGc                                                             >  1:544189/1‑69 (MQ=255)
      gATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGa                                                       <  1:1709549/69‑1 (MQ=255)
      gATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGa                                                       <  1:3022247/69‑1 (MQ=255)
        tACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGaaa                                                     >  1:735982/1‑69 (MQ=255)
        tACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGaa                                                      >  1:290254/1‑68 (MQ=255)
        tACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGaa                                                      >  1:3390660/1‑68 (MQ=255)
                    cATAACGAAGCTGGCTCAGATCGTCGCTGTGAGGATCGATACCGATCGCTACCGAAATATCGTGCACt                                          <  1:2153203/68‑1 (MQ=255)
                    cATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACt                                          <  1:1001613/68‑1 (MQ=255)
                    cATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACt                                          <  1:1614745/68‑1 (MQ=255)
                    cATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACt                                          <  1:2854020/68‑1 (MQ=255)
                    cATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACt                                          <  1:1304065/68‑1 (MQ=255)
                    cATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACt                                          <  1:3135471/68‑1 (MQ=255)
                    cATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACt                                          <  1:80093/68‑1 (MQ=255)
                    cATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACt                                          <  1:808266/68‑1 (MQ=255)
                      tAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTcc                                       >  1:3018901/1‑69 (MQ=255)
                      tAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGAACTTcc                                       >  1:1866315/1‑69 (MQ=255)
                      tAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCAAATATCGTGCACTTCCt                                      >  1:1217184/1‑69 (MQ=255)
                        aCGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTcc                                       >  1:1324480/1‑67 (MQ=255)
                        aCGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCt                                      >  1:1286627/1‑68 (MQ=255)
                        aCGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCt                                      >  1:34069/1‑68 (MQ=255)
                        aCGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCt                                      >  1:406582/1‑68 (MQ=255)
                        aCGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCt                                      >  1:1085008/1‑68 (MQ=255)
                        aCGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCt                                      >  1:1440998/1‑68 (MQ=255)
                        aCGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCt                                      >  1:1676393/1‑68 (MQ=255)
                        aCGAAGCTGGCTCAGATCGTCGCTGACAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCt                                      >  1:1001117/1‑68 (MQ=255)
                             gCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTTCACt                                          <  1:2584349/59‑1 (MQ=255)
                                        tcgtcgCTGTCAGGATCGATACCGATAGCTACCGAAATATCGTGCAc                                           >  1:1573195/1‑47 (MQ=255)
                                        tcgtcgCTGTCAGGATCGATACCGATAGCTACCGAAATATCGTGCACt                                          >  1:3191058/1‑48 (MQ=255)
                                        tcgtcgCTGTCAGGATCGATACCGATAGCTACCGAAATATCGTGCACt                                          >  1:2610615/1‑48 (MQ=255)
                                        tcgtcgCTGTCAGGATCGATACCGATAGCTACCGAAATATCGTGCACt                                          >  1:2696254/1‑48 (MQ=255)
                                        tcgtcgCTGTCAGGATCGATACCGATAGCTACCGAAATATCGTGCACt                                          >  1:2485038/1‑48 (MQ=255)
                                        tcgtcgCTGTCAGGATCGATACCGATAGCTACCGAAATATCGTGCACt                                          >  1:493018/1‑48 (MQ=255)
                                        tcgtcgCTGTCAGGATCGATACCGATAGCTACCGAAATATCGTGCACt                                          >  1:270823/1‑48 (MQ=255)
                                        tcgtcgCTGTCAGGATCGATACCGATAGCTACCGAAATATCGTGCACt                                          >  1:43540/1‑48 (MQ=255)
                                        tcgtcgCTGTCAGGATCGATACCGATAGCTACCGAAATATCGTGCACt                                          >  1:2877564/1‑48 (MQ=255)
                                        tcgtcgCTGTCAGGATCGATACCGATAGCTACCGAAATATCGTGCACt                                          >  1:1558204/1‑48 (MQ=255)
                                        tcgtcgCTGTCAGGATCGATACCGATAGCTACCGAAATATCGTGCACt                                          >  1:1989413/1‑48 (MQ=255)
                                        tcgtcgCTGTCAGGATCGATACCGATAGCTACCGAAATATCGTGCACt                                          >  1:68396/1‑48 (MQ=255)
                                        tcgtcgCTGTCAGGATCGATACCGATAGCTACCGAAATATCGTGCACt                                          >  1:1951207/1‑48 (MQ=255)
                                        tcgtcgCTGTCAGGATCGATACCGATAGCTACCGAAATATCGTGCACt                                          >  1:1847023/1‑48 (MQ=255)
                                        tcgtcgCTGTCAGGATCGATACCGATAGCTACCGAAATATCGTGCACt                                          >  1:1793863/1‑48 (MQ=255)
                                        tcgtcgCTGTCAGGATCGATACCGATAGCTACCGAAATATCGTGCACt                                          >  1:1791456/1‑48 (MQ=255)
                                        tcgtcgCTGTCAGGATCGATACCGATAGCTACCGAAATATCGTGCACt                                          >  1:3027525/1‑48 (MQ=255)
                                        tcgtcgCAGTCAGGATCGATACCGATAGCTACAGaa                                                      >  1:2903015/1‑36 (MQ=25)
                                         cgtcgCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTg                                              <  1:2067471/43‑1 (MQ=255)
                                         cgtcgCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTg                                              <  1:1416646/43‑1 (MQ=255)
                                             gCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTtcgtc                   >  1:1885419/1‑66 (MQ=255)
                                             gCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTtcgtc                   >  1:1299527/1‑66 (MQ=255)
                                             gCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGa                <  1:261513/69‑1 (MQ=255)
                                                gTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAaga             <  1:372631/69‑1 (MQ=255)
                                                   aGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAagag            <  1:3203765/67‑1 (MQ=255)
                                                   aGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAagag            <  1:318573/67‑1 (MQ=255)
                                                   aGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAagag            <  1:929282/67‑1 (MQ=255)
                                                   aGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAagag            <  1:2462865/67‑1 (MQ=255)
                                                   aGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAagag            <  1:2354541/67‑1 (MQ=255)
                                                   aGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAagag            <  1:566892/67‑1 (MQ=255)
                                                   aGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAagag            <  1:602564/67‑1 (MQ=255)
                                                   aGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAgaga           <  1:1542099/68‑1 (MQ=255)
                                                   aGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAgaga           <  1:1030332/68‑1 (MQ=255)
                                                   aGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAgaga           <  1:2134435/68‑1 (MQ=255)
                                                   aGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAgaga           <  1:371981/68‑1 (MQ=255)
                                                   aGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAgaga           <  1:2230381/68‑1 (MQ=255)
                                                   aGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAgaga           <  1:531284/68‑1 (MQ=255)
                                                   aGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAgaga           <  1:2805502/68‑1 (MQ=255)
                                                   aGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAgaga           <  1:483732/68‑1 (MQ=255)
                                                   aGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAgaga           <  1:2591248/68‑1 (MQ=255)
                                                   aGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCAGAAgaga           <  1:1283408/68‑1 (MQ=255)
                                                    ggATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAagag            <  1:339679/66‑1 (MQ=255)
                                                        cGATACCGATAGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAGAGACTTccc     <  1:3347056/69‑1 (MQ=255)
                                                        cGATACCGATAGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAGAGACTTccc     <  1:1120414/69‑1 (MQ=255)
                                                        cGATACCGATAGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAGAGACTTccc     <  1:1611086/69‑1 (MQ=255)
                                                        cGATACCGATAGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAGAGACTTccc     <  1:1325989/69‑1 (MQ=255)
                                                        cGATACCGATAGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAGAGACTTccc     <  1:2977574/69‑1 (MQ=255)
                                                        cGATACCGATAGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAGAGACTTccc     <  1:1782881/69‑1 (MQ=255)
                                                        cGATACCGATAGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAGAGACTTccc     <  1:2915496/69‑1 (MQ=255)
                                                        cGATACCGATAGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAGAGACTTccc     <  1:1988643/69‑1 (MQ=255)
                                                        cGATACCGATAGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAGAGACTTccc     <  1:2477949/69‑1 (MQ=255)
                                                        cGATACCGATAGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAGAGACTTccc     <  1:1433346/69‑1 (MQ=255)
                                                         gATACCGATAGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAGAGACTTccc     <  1:651895/68‑1 (MQ=255)
                                                           tACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAGAGACTTCCCAgg  <  1:920118/69‑1 (MQ=255)
                                                           tACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAGAGACTTCCCAgg  <  1:3289096/69‑1 (MQ=255)
                                                           tACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAGAGACTTCCCAgg  <  1:3365564/69‑1 (MQ=255)
                                                           tACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAGAGACTTCCCAgg  <  1:1082304/69‑1 (MQ=255)
                                                            aCCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGc                               >  1:2188199/1‑39 (MQ=255)
                                                            aCCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGc                               >  1:2502967/1‑39 (MQ=255)
                                                            aCCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGc                               >  1:614579/1‑39 (MQ=255)
                                                            aCCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGc                               >  1:653458/1‑39 (MQ=255)
                                                            aCCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGc                               >  1:1716189/1‑39 (MQ=255)
                                                            aCCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGc                               >  1:1691736/1‑39 (MQ=255)
                                                            aCCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGc                               >  1:2335040/1‑39 (MQ=255)
                                                            aCCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGc                               <  1:1136458/39‑1 (MQ=255)
                                                            aCCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGc                               <  1:2343246/39‑1 (MQ=255)
                                                            aCCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGc                               >  1:3161238/1‑39 (MQ=255)
                                                            aCCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGc                               >  1:2877867/1‑39 (MQ=255)
                                                            aCCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGc                               >  1:2952991/1‑39 (MQ=255)
                                                            aCCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGc                               >  1:1413373/1‑39 (MQ=255)
                                                            aCCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGc                               >  1:3235471/1‑39 (MQ=255)
                                                            aCCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGc                               >  1:3150492/1‑39 (MQ=255)
                                                                  cGCTACCGAAATATCGTTCACTTCCTGCGAAGCCAGc                           <  1:1250870/37‑1 (MQ=255)
                                                                  cGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGc                           <  1:1216654/37‑1 (MQ=255)
                                                                  |                                                             
CGCGAGGATACAGATTTTGCCATAACGAAGCTGGCTCAGATCGTCGCTGTCAGGATCGATACCGATCGCTACCGAAATATCGTGCACTTCCTGCGAAGCCAGCACTTCGTCGGAAGAGACTTCCCAGG  >  minE/2029378‑2029505

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: