Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,464,254 T→C 22.8% intergenic (‑58/+91) yfaA ← / ← gyrA hypothetical protein/DNA gyrase (type II topoisomerase), subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,464,2540TC22.8% 74.7 / 17.3 61intergenic (‑58/+91)yfaA/gyrAhypothetical protein/DNA gyrase (type II topoisomerase), subunit A
Reads supporting (aligned to +/- strand):  ref base T (12/35);  new base C (7/7);  total (19/42)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-01

CACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTC  >  minE/1464187‑1464321
                                                                   |                                                                   
caccacTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGt                                                                     <  1:2719660/68‑1 (MQ=255)
caccacTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGt                                                                     <  1:1799593/68‑1 (MQ=255)
  ccacTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGc                                                                   <  1:1143175/68‑1 (MQ=255)
  ccacTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGc                                                                   <  1:2981554/68‑1 (MQ=255)
  ccacTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGc                                                                   <  1:1611902/68‑1 (MQ=255)
    acTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCt                                                                <  1:1398596/69‑1 (MQ=255)
    acTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCt                                                                <  1:2328078/69‑1 (MQ=255)
     cTCATTGCCACATTCCTTGTGTATAGCCATCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCt                                                                <  1:3281895/68‑1 (MQ=255)
     cTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCt                                                                <  1:1370398/68‑1 (MQ=255)
     cTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTATACCGTGCCCt                                                                <  1:1619997/68‑1 (MQ=255)
       cATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCt                                                                >  1:2169162/1‑66 (MQ=255)
       cATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCt                                                                >  1:1137956/1‑66 (MQ=255)
       cATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCt                                                                <  1:211451/66‑1 (MQ=255)
       cATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCt                                                                <  1:2091262/66‑1 (MQ=255)
       cATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCt                                                                >  1:3238830/1‑66 (MQ=255)
       cATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCt                                                                >  1:2751859/1‑66 (MQ=255)
       cATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCt                                                                >  1:1742676/1‑66 (MQ=255)
       cATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCt                                                                <  1:3429133/66‑1 (MQ=255)
       cATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCt                                                                >  1:535172/1‑66 (MQ=255)
       cATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCt                                                                <  1:579342/66‑1 (MQ=255)
       cATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCt                                                                >  1:968390/1‑66 (MQ=255)
                aTTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGc                                                                   <  1:1677426/54‑1 (MQ=255)
                          tataGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACtt                                          >  1:2650131/1‑69 (MQ=255)
                          tataGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACtt                                          >  1:2071303/1‑69 (MQ=255)
                                cagccaTTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGc                                    <  1:2892266/69‑1 (MQ=255)
                                cagccaTTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGc                                    <  1:2562650/69‑1 (MQ=255)
                                cagccaTTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGc                                    <  1:3121568/69‑1 (MQ=255)
                                cagccaTTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGc                                    <  1:182401/69‑1 (MQ=255)
                                cagccaTTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGc                                    <  1:1366566/69‑1 (MQ=255)
                                cagccaTTTTTTACGGGCACAGCCAAACTTTAACGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGc                                    <  1:3350315/69‑1 (MQ=255)
                                     aTTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGa                               >  1:1480186/1‑69 (MQ=255)
                                      ttttttACGGGCACAGCCAAACTTTACCGCGCCCTAATACGACAAAAGCCCAGACTTTACAGCCTGGAc                              >  1:836913/1‑69 (MQ=255)
                                      ttttttACGGGCACAGCCAAACTTTACCGCGCCCTAATACGACAAAAGCCCAGACTTTACAGCCTGGAc                              >  1:278296/1‑69 (MQ=255)
                                      ttttttACGGGCACAGCCAAACTTTACCGCGCCCTAATACGACAAAAGCCCAGACTTTACAGCCTGGAc                              >  1:373180/1‑69 (MQ=255)
                                      ttttttACGGGCACAGCCAAACTTTACCGCGCCCTAATACGACAAAAGCCCAGACTTTACAGCCTGGAc                              >  1:1666704/1‑69 (MQ=255)
                                      ttttttACGGGCACAGCCAAACTTTACCGCGCCCTAATACGACAAAAGCCCAGACTTTACAGCCTGGAc                              >  1:2161157/1‑69 (MQ=255)
                                      ttttttACGGGCACAGCCAAACTTTACCGCGCCCTAATACGACAAAAGCCCAGACTTTACAGCCTGGAc                              >  1:331009/1‑69 (MQ=255)
                                      ttttttACGGGCACAGCCAAACTTTACCGCGCCCTAATACGACAAAAGCCCAGACTTTACAGCCTGGAc                              >  1:3266599/1‑69 (MQ=255)
                                         tttACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACtt                            >  1:1341946/1‑68 (MQ=255)
                                          ttACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACttt                           <  1:153002/68‑1 (MQ=255)
                                                  acaGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTttcaattca                   >  1:352058/1‑68 (MQ=255)
                                                      ccAAACTTTACCGCGCCCTAATACGACAAAAGCCGAGACtt                                          <  1:1571428/41‑1 (MQ=38)
                                                      ccAAACTTTACCGCGCCCTAATACGACAAAAGCCCAGACtt                                          <  1:1925792/41‑1 (MQ=255)
                                                      ccAAACTTTACCGCGCCCTAATACGACAAAAGCCCAGACtt                                          <  1:3190985/41‑1 (MQ=255)
                                                      ccAAACTTTACCGCGCCCTAATACGACAAAAGCCCAGACtt                                          <  1:186655/41‑1 (MQ=255)
                                                      ccAAACTTTACCGCGCCCTAATACGACAAAAGCCCAGACtt                                          <  1:3350335/41‑1 (MQ=255)
                                                      ccAAACTTTACCGCGCCCTAATACGACAAAAGCCCAGACtt                                          <  1:480197/41‑1 (MQ=255)
                                                      ccAAACTTTACCGCGCCCTAATACGACAAAAGCCCAGACtt                                          <  1:3364817/41‑1 (MQ=255)
                                                          aCTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGga          <  1:2551462/69‑1 (MQ=255)
                                                          aCTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGga          <  1:3018490/69‑1 (MQ=255)
                                                          aCTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGga          <  1:9659/69‑1 (MQ=255)
                                                          aCTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGga          <  1:2537180/69‑1 (MQ=255)
                                                          aCTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGga          <  1:55015/69‑1 (MQ=255)
                                                          aCTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGga          <  1:1612786/69‑1 (MQ=255)
                                                          aCTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGga          <  1:1821179/69‑1 (MQ=255)
                                                           cTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGga          <  1:3356454/68‑1 (MQ=255)
                                                           cTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGga          <  1:2989434/68‑1 (MQ=255)
                                                                  gTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTc  <  1:2924688/69‑1 (MQ=255)
                                                                  gTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTc  <  1:2510582/69‑1 (MQ=255)
                                                                  gTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTc  <  1:1710306/69‑1 (MQ=255)
                                                                  gTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTc  <  1:965154/69‑1 (MQ=255)
                                                                  gTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTc  <  1:3033637/69‑1 (MQ=255)
                                                                  gTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTc  <  1:3025165/69‑1 (MQ=255)
                                                                   tGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTc  <  1:2994010/68‑1 (MQ=255)
                                                                   |                                                                   
CACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTC  >  minE/1464187‑1464321

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: