Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,492,926 T→C 25.7% L575L (TTG→CTG)  gyrB → DNA gyrase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,492,9260TC25.7% 65.7 / 19.4 58L575L (TTG→CTG) gyrBDNA gyrase, subunit B
Reads supporting (aligned to +/- strand):  ref base T (34/9);  new base C (7/8);  total (41/17)
Fisher's exact test for biased strand distribution p-value = 2.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

TCAGTACCAGATCTCTATCGCGCTGGACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAAT  >  minE/2492862‑2492981
                                                                |                                                       
tCAGTACCAGATCTCTATCGCGCTGGACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATtgg                                                      <  1:11802/68‑1 (MQ=255)
      ccAGATCTCTATCGCGCTGGACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCg                                                >  1:2064188/1‑68 (MQ=255)
      ccAGATCTCTATCGCGCTGGACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGa                                               >  1:3184477/1‑69 (MQ=255)
      ccAGATCTCTATCGCGCTGGACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGa                                               >  1:717712/1‑69 (MQ=255)
      ccAGATCTCTATCGCGCTGGACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGa                                               >  1:1632798/1‑69 (MQ=255)
      ccAGATCTCTATCGCGCTGGACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGa                                               >  1:1563910/1‑69 (MQ=255)
      ccAGATCTCTATCGCGCTGGACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGa                                               >  1:2506860/1‑69 (MQ=255)
      ccAGATCTCTATCGCGCTGGACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGa                                               >  1:2603066/1‑69 (MQ=255)
            ctctATCGCGCTGGACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGc                                            >  1:287667/1‑66 (MQ=255)
            ctctATCGCGCTGGACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGc                                            >  1:3075082/1‑66 (MQ=255)
            ctctATCGCGCTGGACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGc                                            >  1:1751075/1‑66 (MQ=255)
            ctctATCGCGCTGGACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGc                                            >  1:648736/1‑66 (MQ=255)
            ctctATCGCGCTGGACGGCGCAACGCTGCAAACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGc                                            >  1:167798/1‑66 (MQ=255)
                  cgcgCTGGACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGaa                                   >  1:2427792/1‑69 (MQ=255)
                         gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                            >  1:3122200/1‑69 (MQ=255)
                         gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                            >  1:2223759/1‑69 (MQ=255)
                         gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                            >  1:267595/1‑69 (MQ=255)
                         gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                            >  1:3189671/1‑69 (MQ=255)
                         gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                            >  1:725153/1‑69 (MQ=255)
                         gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                            >  1:739642/1‑69 (MQ=255)
                         gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                            >  1:1582637/1‑69 (MQ=255)
                         gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                            >  1:97016/1‑69 (MQ=255)
                         gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                            >  1:1404583/1‑69 (MQ=255)
                         gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                            >  1:1400622/1‑69 (MQ=255)
                         gACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTa                            >  1:1379139/1‑69 (MQ=255)
                           cGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTAt                           >  1:66919/1‑68 (MQ=255)
                           cGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTAt                           >  1:3183530/1‑68 (MQ=255)
                           cGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTAt                           >  1:2914662/1‑68 (MQ=255)
                           cGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTAt                           >  1:1631405/1‑68 (MQ=255)
                           cGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTAt                           >  1:3211778/1‑68 (MQ=255)
                           cGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTAt                           >  1:2793041/1‑68 (MQ=255)
                                       gCACACCAACGCCAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAAcgc              <  1:350869/69‑1 (MQ=255)
                                       gCACACCAACGCCAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAAcgc              <  1:1159763/69‑1 (MQ=255)
                                       gCACACCAACGCCAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAAcgc              <  1:3286607/69‑1 (MQ=255)
                                       gCACACCAACGCCAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAAcgc              <  1:2367122/69‑1 (MQ=255)
                                       gCACACCAACGCCAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAAcgc              <  1:3322203/69‑1 (MQ=255)
                                       gCACACCAACGCCAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAAcgc              <  1:374868/69‑1 (MQ=255)
                                       gCACACCAACGCCAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAAcgc              <  1:981838/69‑1 (MQ=255)
                                        cacaCCAACGCCAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAAcgc              <  1:310988/68‑1 (MQ=255)
                                              aaCGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCa        >  1:26627/1‑68 (MQ=255)
                                              aaCGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCa        >  1:2905026/1‑68 (MQ=255)
                                              aaCGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCa        >  1:1511550/1‑68 (MQ=255)
                                              aaCGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCa        >  1:2515853/1‑68 (MQ=255)
                                                cGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGaa     <  1:635604/69‑1 (MQ=255)
                                                cGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGaa     <  1:2153411/69‑1 (MQ=255)
                                                cGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGaa     <  1:2164694/69‑1 (MQ=255)
                                                cGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGaa     <  1:2219407/69‑1 (MQ=255)
                                                   cAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAACGCGACGCAGaaaa   >  1:3217499/1‑68 (MQ=255)
                                                   cAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAACGCGACGCAGAAAAt  >  1:2170665/1‑69 (MQ=255)
                                                   cAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAACGCGACGCAGAAAAt  >  1:2562448/1‑69 (MQ=255)
                                                   cAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAACGCGACGCAGAAAAt  >  1:1866227/1‑69 (MQ=255)
                                                   cAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAACGCGACGCAGAAAAt  >  1:1690240/1‑69 (MQ=255)
                                                   cAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAACGCGACGCAGAAAAt  >  1:1600848/1‑69 (MQ=255)
                                                   cAGTGCACCGGCGCTGGCTGGCGAAGCGTTAGAGAAACTGGTGTCTGAGTACAACGCGACGCAGAAAAt  >  1:2709993/1‑69 (MQ=255)
                                                       gCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGc         <  1:3024683/58‑1 (MQ=255)
                                                       gCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGc         <  1:2015828/58‑1 (MQ=255)
                                                       gCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGc         <  1:737614/58‑1 (MQ=255)
                                                       gCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGc         <  1:1177094/58‑1 (MQ=255)
                                                                |                                                       
TCAGTACCAGATCTCTATCGCGCTGGACGGCGCAACGCTGCACACCAACGCCAGTGCACCGGCATTGGCTGGCGAAGCGTTAGAGAAACTGGTATCTGAGTACAACGCGACGCAGAAAAT  >  minE/2492862‑2492981

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: