Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,465,692 A→C 22.0% A427A (GCT→GCG gyrA ← DNA gyrase (type II topoisomerase), subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,465,6920AC22.0% 72.2 / 18.4 50A427A (GCT→GCGgyrADNA gyrase (type II topoisomerase), subunit A
Reads supporting (aligned to +/- strand):  ref base A (25/14);  new base C (5/6);  total (30/20)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGTACAGACCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGATTAGCAACCAGCGCAGTTTTCGC  >  minE/1465626‑1465757
                                                                  |                                                                 
gtagtaCAGACCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTTGCCAg                                                                  <  1:2854844/68‑1 (MQ=255)
gtagtaCAGACCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAg                                                                  <  1:2382027/68‑1 (MQ=255)
gtagtaCAGACCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAg                                                                  <  1:2299469/68‑1 (MQ=255)
gtagtaCAGACCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAg                                                                  <  1:2030056/68‑1 (MQ=255)
gtagtaCAGACCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAg                                                                  <  1:1688143/68‑1 (MQ=255)
gtagtaCAGACCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAg                                                                  <  1:1970582/68‑1 (MQ=255)
 tagtaCAGACCATCACGCACGCCGATCTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAg                                                                  >  1:115673/1‑67 (MQ=255)
 tagtaCAGACCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCa                                                                <  1:1145440/69‑1 (MQ=255)
      cAGACCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGt                                                             >  1:2910475/1‑67 (MQ=255)
      cAGACCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGt                                                             >  1:130815/1‑67 (MQ=255)
      cAGACCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGt                                                             >  1:2591491/1‑67 (MQ=255)
      cAGACCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGt                                                             >  1:703152/1‑67 (MQ=255)
      cAGACCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGt                                                             >  1:2365897/1‑67 (MQ=255)
         aCCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAg                                                        <  1:881344/69‑1 (MQ=255)
         aCCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAg                                                        <  1:2185198/69‑1 (MQ=255)
          ccATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGc                                                       <  1:862461/69‑1 (MQ=255)
                     gcCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAAc                                            <  1:678414/69‑1 (MQ=255)
                     gcCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAAc                                            <  1:179873/69‑1 (MQ=255)
                     gcCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAAc                                            <  1:2527330/69‑1 (MQ=255)
                                  tCCAGCCATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGc                                  <  1:1344327/66‑1 (MQ=255)
                                  tCCAGCCATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGc                                  <  1:2790928/66‑1 (MQ=255)
                                  tCCAGCCATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGc                                  <  1:2413676/66‑1 (MQ=255)
                                  tCCAGCCATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGc                                  <  1:3415502/66‑1 (MQ=255)
                                   ccAGCCATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGc                                  <  1:256772/65‑1 (MQ=255)
                                       ccATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGCTGCCACgg                          >  1:2351103/1‑69 (MQ=255)
                                       ccATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGCTGCCACgg                          >  1:561290/1‑69 (MQ=255)
                                       ccATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGCTGCCACgg                          >  1:3168276/1‑69 (MQ=255)
                                       ccATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGCTGCCACgg                          >  1:2016210/1‑69 (MQ=255)
                                       ccATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGCTGCCACgg                          >  1:3116026/1‑69 (MQ=255)
                                          ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACgg                          >  1:487899/1‑66 (MQ=255)
                                          ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                       >  1:392579/1‑69 (MQ=255)
                                          ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                       >  1:2640429/1‑69 (MQ=255)
                                          ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                       >  1:3379064/1‑69 (MQ=255)
                                          ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                       >  1:615799/1‑69 (MQ=255)
                                          ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                       >  1:3099136/1‑69 (MQ=255)
                                          ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                       >  1:2951788/1‑69 (MQ=255)
                                          ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                       >  1:1099653/1‑69 (MQ=255)
                                          ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                       >  1:2640952/1‑69 (MQ=255)
                                          ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                       >  1:1319348/1‑69 (MQ=255)
                                          ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                       >  1:877380/1‑69 (MQ=255)
                                          ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                       >  1:1565431/1‑69 (MQ=255)
                                          ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                       >  1:1579344/1‑69 (MQ=255)
                                          ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                       >  1:2211912/1‑69 (MQ=255)
                                          ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                       >  1:1608364/1‑69 (MQ=255)
                                          ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                       >  1:1991853/1‑69 (MQ=255)
                                            ccGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGATTa                      >  1:239353/1‑68 (MQ=255)
                                            ccGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGATTa                      >  1:1503507/1‑68 (MQ=255)
                                            ccGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGATTa                      >  1:1057632/1‑68 (MQ=255)
                                                          tcgtcgCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGc                                  <  1:1020575/42‑1 (MQ=37)
                                                               gCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGATTAGCAACCAGCGCAGTTTTCGc  <  1:168240/69‑1 (MQ=255)
                                                                  |                                                                 
GTAGTACAGACCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGATTAGCAACCAGCGCAGTTTTCGC  >  minE/1465626‑1465757

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