Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,464,999 A→G 47.5% T658T (ACT→ACC 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,464,9990AG47.5% ‑1.8 / 32.1 40T658T (ACT→ACCgyrADNA gyrase (type II topoisomerase), subunit A
Reads supporting (aligned to +/- strand):  ref base A (3/18);  new base G (10/9);  total (13/27)
Fisher's exact test for biased strand distribution p-value = 1.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCGTCAACCAGTTTGATCGCCACTTTACCGGCGGTACGCAGACGGTTGAACTCGGTGAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGCCATGAAGACTTTCACGCCTTCTTCAAACTCGGT  >  minE/1464935‑1465061
                                                                |                                                              
cGCCTTCAACCAGTTTGATCGCCACTTTACCGGCGGTACGCAGACGGTTGAACTCGGTGAGGACGGttt                                                            <  1:1913329/69‑1 (MQ=255)
cGCCTTCAACCAGTTTGATCGCCACTTTACCGGCGGTACGCAGACGGTTGAACTCGGTGAGGACGGttt                                                            <  1:1070292/69‑1 (MQ=255)
cGCCTTCAACCAGTTTGATCGCCACTTTACCGGCGGTACGCAGACGGTTGAACTCGGTGAGGACGGttt                                                            <  1:60900/69‑1 (MQ=255)
cGCCTTCAACCAGTTTGATCGCCACTTTACCGGCGGTACGCAGACGGTTGAACTCGGTGAGGACGGttt                                                            <  1:2969692/69‑1 (MQ=255)
cGCCTTCAACCAGTTTGATCGCCACTTTACCGGCGGTACGCAGACGGTTGAACTCGGTGAGGACGGttt                                                            <  1:1876844/69‑1 (MQ=255)
cGCCTTCAACCAGTTTGATCGCCACTTTACCGGCGGTACGCAGACGGTTGAACTCGGTGAGGACGGttt                                                            <  1:1712061/69‑1 (MQ=255)
cGCCTTCAACCAGTTTGATCGCCACTTTACCGGCGGTACGCAGACGGTTGAACTCGGTGAGGACGGttt                                                            <  1:1176603/69‑1 (MQ=255)
cGCCTTCAACCAGTTTGATCGCCACTTTACCGGCGGTACGCAGACGGTTGAACTCGGTGAGGACGGttt                                                            <  1:601865/69‑1 (MQ=255)
                        cTTTACCGGCGGTACGCAGACGGTTGAACTCGGTGAGGACGGttt                                                            <  1:2444241/45‑1 (MQ=255)
                           tACCGGCGGTACGCAGACGGTTGAACTCGGTTAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGcc                                  >  1:2270858/1‑68 (MQ=255)
                           tACCGGCGGTACGCAGACGGTTGAACTCGGTGAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGcc                                  >  1:984892/1‑68 (MQ=255)
                           tACCGGCGGTACGCAGACGGTTGAACTCGGTGAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGcc                                  >  1:898475/1‑68 (MQ=255)
                                gcggTACGCAGACGGTTGAACTCGGTGAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGCCATGAAg                            <  1:1689408/69‑1 (MQ=255)
                                gcggTACGCAGACGGTTGAACTCGGTGAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGCCATGAAg                            <  1:578045/69‑1 (MQ=255)
                                gcggTACGCAGACGGTTGAACTCGGTGAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGCCATGAAg                            <  1:369686/69‑1 (MQ=255)
                                gcggTACGCAGACGGTTGAACTCGGTGAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGCCATGAAg                            <  1:3351892/69‑1 (MQ=255)
                                gcggTACGCAGACGGTTGAACTCGGTGAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGCCATGAAg                            <  1:2245921/69‑1 (MQ=255)
                                gcggTACGCAGACGGTTGAACTCGGTGAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGCCATGAAg                            <  1:3072641/69‑1 (MQ=255)
                                gcggTACGCAGACGGTTGAACTCGGTGAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGCCATGAAg                            <  1:1250272/69‑1 (MQ=255)
                                gcggTACGCAGACGGTTGAACTCGGTGAGGACAGTTTTCTTCACGGTACCGTTAGAGGTCGCCATGAAg                            <  1:2284655/69‑1 (MQ=255)
                                   gTACGCAGACGGTTGAACTCGGTGAGGACGGTTTTCTTCACAGTACCGTTAGCGGTCGCCa                                 >  1:2715020/1‑61 (MQ=255)
                                   gTACGCAGACGGTTGAACTCGGTGAGGACGGTTTTCTTCACAGTACCGTTAGCGGTCGCCATGAAGACt                         >  1:2868687/1‑69 (MQ=255)
                                   gTACGCAGACGGTTGAACTCGGTGAGGACGGTTTTCTTCACAGTACCGTTAGCGGTCGCCATGAAGACt                         >  1:2899928/1‑69 (MQ=255)
                                   gTACGCAGACGGTTGAACTCGGTGAGGACGGTTTTCTTCACAGTACCGTTAGCGGTCGCCATGAAGACt                         >  1:637687/1‑69 (MQ=255)
                                   gTACGCAGACGGTTGAACTCGGTGAGGACGGTTTTCTTCACAGTACCGTTAGCGGTCGCCATGAAGACt                         >  1:2985555/1‑69 (MQ=255)
                                   gTACGCAGACGGTTGAACTCGGTGAGGACGGTTTTCTTCACAGTACCGTTAGCGGTCGCCATGAAGACt                         >  1:219729/1‑69 (MQ=255)
                                   gTACGCAGACGGTTGAACTCGGTGAGGACGGTTTTCTTCACAGTACCGTTAGCGGTCGCCATGAAGACt                         >  1:755809/1‑69 (MQ=255)
                                   gTACGCAGACGGTTGAACTCGGTGAGGACGGTTTTCTTCACAGTACCGTTAGCGGTCGCCATGAAGACt                         >  1:1395707/1‑69 (MQ=255)
                                   gTACGCAGACGGTTGAACTCGGTGAGGACGGTTTTCTTCACAGTACCGTTAGCGGTCGCCATGAAGACt                         >  1:1116656/1‑69 (MQ=255)
                                   gTAAGCAGACGGTTGAACTCGGTGAGGACGGTTTTCTTCACAGTACCGTTAGCGGTCGCCATGAAGAc                          >  1:3052784/1‑68 (MQ=255)
                                     aCGCAGACGGTTGAACTCGGTGAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGCCATGAAGACttt                       <  1:2865075/69‑1 (MQ=255)
                                     aCGCAGACGGTTGAACTCGGTGAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGCCATGAAGACttt                       <  1:2635235/69‑1 (MQ=255)
                                     aCGCAGACGGTTGAACTCGGTGAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGCCATGAAGACttt                       <  1:2630176/69‑1 (MQ=255)
                                     aCGCAGACGGTTGAACTCGGTGAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGCCATGAAGACttt                       <  1:151345/69‑1 (MQ=255)
                                                          gAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGCCATGAAGACTTTCACGCCTTCTTCAAACTCGGt  <  1:689657/69‑1 (MQ=255)
                                                          gAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGCCATGAAGACTTTCACGCCTTCTTCAAACTCGGt  <  1:2023053/69‑1 (MQ=255)
                                                          gAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGCCATGAAGACTTTCACGCCTTCTTCAAACTCGGt  <  1:858885/69‑1 (MQ=255)
                                                          gAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGCCATGAAGACTTTCACGCCTTCTTCAAACTCGGt  <  1:1729975/69‑1 (MQ=255)
                                                          gAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGCCATGAAGACTTTCACGCCTTCTTCAAACTCGGt  <  1:99244/69‑1 (MQ=255)
                                                           aGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGCCATGAAGACTTTCACGCCTTCTTCAAACTCGGt  <  1:997902/68‑1 (MQ=255)
                                                                |                                                              
CGCCGTCAACCAGTTTGATCGCCACTTTACCGGCGGTACGCAGACGGTTGAACTCGGTGAGGACAGTTTTCTTCACGGTACCGTTAGCGGTCGCCATGAAGACTTTCACGCCTTCTTCAAACTCGGT  >  minE/1464935‑1465061

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: