Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,465,701 G→C 26.8% L424L (CTC→CTG 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,7010GC26.8% 56.2 / 20.2 41L424L (CTC→CTGgyrADNA gyrase (type II topoisomerase), subunit A
Reads supporting (aligned to +/- strand):  ref base G (19/11);  new base C (5/6);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGATTAGCAACCAGCGCAGTTTTCGCTTCTGCA  >  minE/1465635‑1465764
                                                                  |                                                               
aCCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAg                                                               <  1:881344/69‑1 (MQ=255)
aCCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAg                                                               <  1:2185198/69‑1 (MQ=255)
 ccATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGc                                                              <  1:862461/69‑1 (MQ=255)
            gcCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAAc                                                   <  1:179873/69‑1 (MQ=255)
            gcCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAAc                                                   <  1:2527330/69‑1 (MQ=255)
            gcCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAAc                                                   <  1:678414/69‑1 (MQ=255)
                         tCCAGCCATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGc                                         <  1:3415502/66‑1 (MQ=255)
                         tCCAGCCATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGc                                         <  1:1344327/66‑1 (MQ=255)
                         tCCAGCCATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGc                                         <  1:2790928/66‑1 (MQ=255)
                         tCCAGCCATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGc                                         <  1:2413676/66‑1 (MQ=255)
                          ccAGCCATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGc                                         <  1:256772/65‑1 (MQ=255)
                              ccATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGCTGCCACgg                                 >  1:3116026/1‑69 (MQ=255)
                              ccATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGCTGCCACgg                                 >  1:561290/1‑69 (MQ=255)
                              ccATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGCTGCCACgg                                 >  1:2016210/1‑69 (MQ=255)
                              ccATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGCTGCCACgg                                 >  1:3168276/1‑69 (MQ=255)
                              ccATTCCGGACGCGCAGCATCGTCGCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGCTGCCACgg                                 >  1:2351103/1‑69 (MQ=255)
                                 ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACgg                                 >  1:487899/1‑66 (MQ=255)
                                 ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                              >  1:2211912/1‑69 (MQ=255)
                                 ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                              >  1:1991853/1‑69 (MQ=255)
                                 ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                              >  1:615799/1‑69 (MQ=255)
                                 ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                              >  1:1608364/1‑69 (MQ=255)
                                 ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                              >  1:1579344/1‑69 (MQ=255)
                                 ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                              >  1:2640429/1‑69 (MQ=255)
                                 ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                              >  1:2640952/1‑69 (MQ=255)
                                 ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                              >  1:1565431/1‑69 (MQ=255)
                                 ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                              >  1:2951788/1‑69 (MQ=255)
                                 ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                              >  1:3099136/1‑69 (MQ=255)
                                 ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                              >  1:877380/1‑69 (MQ=255)
                                 ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                              >  1:1319348/1‑69 (MQ=255)
                                 ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                              >  1:3379064/1‑69 (MQ=255)
                                 ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                              >  1:1099653/1‑69 (MQ=255)
                                 ttCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGAtt                              >  1:392579/1‑69 (MQ=255)
                                   ccGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGATTa                             >  1:1057632/1‑68 (MQ=255)
                                   ccGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGATTa                             >  1:239353/1‑68 (MQ=255)
                                   ccGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGATTa                             >  1:1503507/1‑68 (MQ=255)
                                                 tcgtcgCCCGCACGTTCCAGCATCGCAGCAACGTTGCCCAGc                                         <  1:1020575/42‑1 (MQ=37)
                                                      gCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGATTAGCAACCAGCGCAGTTTTCGc         <  1:168240/69‑1 (MQ=255)
                                                                tCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGATTAGCAACCAGCGCAGTTTTCGCTTCTGCa  <  1:1792059/66‑1 (MQ=255)
                                                                tCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGATTAGCAACCAGCGCAGTTTTCGCTTCTGCa  <  1:69815/66‑1 (MQ=255)
                                                                tCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGATTAGCAACCAGCGCAGTTTTCGCTTCTGCa  <  1:1380744/66‑1 (MQ=255)
                                                                tCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGATTAGCAACCAGCGCAGTTTTCGCTTCTGCa  <  1:928590/66‑1 (MQ=255)
                                                                  |                                                               
ACCATCACGCACGCCGAACTCTGGCTCCAGCCATTCCGGACGCGCAGCATCGTCGCCAGCACGTTCGAGCATCGCGGCAACGTTGCCCAGCTGCCACGGATTAGCAACCAGCGCAGTTTTCGCTTCTGCA  >  minE/1465635‑1465764

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