Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,292,843 A→G 33.9% N105N (AAT→AAC yeeO ← predicted multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,292,8430AG33.9% 37.9 / 29.2 56N105N (AAT→AACyeeOpredicted multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base A (12/25);  new base G (7/12);  total (19/37)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

CATCGCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATTCTCCATGAAGATAGGAACAGCAAGCGGGGTTATCTCGCGCCAGAACAACACTTT  >  minE/1292777‑1292899
                                                                  |                                                        
cATCGCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCAtt                                                        >  1:1340942/1‑69 (MQ=255)
   cGCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATtctc                                                     <  1:3082106/69‑1 (MQ=255)
   cGCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATtctc                                                     <  1:3395713/69‑1 (MQ=255)
    gCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATtctc                                                     >  1:3413127/1‑68 (MQ=255)
    gCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATtctc                                                     <  1:2039599/68‑1 (MQ=255)
    gCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATtctc                                                     <  1:2223871/68‑1 (MQ=255)
    gCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATtctc                                                     >  1:2845641/1‑68 (MQ=255)
    gCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATtctc                                                     <  1:3130664/68‑1 (MQ=255)
    gCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATtctc                                                     >  1:505939/1‑68 (MQ=255)
    gCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATtctc                                                     <  1:1414220/68‑1 (MQ=255)
    gCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATtctc                                                     <  1:688291/68‑1 (MQ=255)
    gCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATtctc                                                     >  1:1322292/1‑68 (MQ=255)
    gCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATtctc                                                     >  1:122370/1‑68 (MQ=255)
    gCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATtcta                                                     >  1:2761426/1‑67 (MQ=255)
    gCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATTCTcc                                                    >  1:736969/1‑69 (MQ=255)
                   cagcTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATTCTCCATGAAGATAGGAACa                                     <  1:1290474/69‑1 (MQ=255)
                   cagcTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATTCTCCATGAAGATAGGAACa                                     <  1:1926693/69‑1 (MQ=255)
                   cagcTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATTCTCCATGAAGATAGGAACa                                     <  1:1946089/69‑1 (MQ=255)
                    agcTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATtctc                                                     <  1:935427/52‑1 (MQ=255)
                     gcTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATTCTCCATGAAGATAGGAACAGc                                   <  1:1672170/69‑1 (MQ=255)
                     gcTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATTCTCCATGAAGATAGGAACAGc                                   <  1:2005163/69‑1 (MQ=255)
                     gcTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATTCTCCATGAAGATAGGAACAGc                                   <  1:1285985/69‑1 (MQ=255)
                     gcTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATTCTCCATGAAGATAGGAACAGc                                   <  1:723048/69‑1 (MQ=255)
                     gcTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATTCTCCATGAAGATAGGAACAGc                                   <  1:2375096/69‑1 (MQ=255)
                     gcTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATTCTCCATGAAGATAGGAACAGc                                   <  1:2931023/69‑1 (MQ=255)
                     gcTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATTCTCCATGAAGATAGGAACAGc                                   <  1:2887273/69‑1 (MQ=255)
                     gcTGACCAGAAAAGCGCTCAGAACCCCCATCAACAGGACACAGGCATTCTCCATGAAGATAGGAACAGc                                   <  1:204171/69‑1 (MQ=255)
                       ttaCCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATtctc                                                     <  1:2769238/47‑1 (MQ=255)
                          ccAGAAAAGTGCTCAGGACCCCCATCAACAGGACACAGGCGTTCTCCATGAAGATAGGAACAGCAAGCg                              <  1:433646/69‑1 (MQ=255)
                               aaaGTGCTCAGGACCCCCATCAACAGGACACAGGCGTTCTCCATGAAGATAGGAACAGCAAGCg                              <  1:2643761/64‑1 (MQ=255)
                               aaaGTGCTCAGGACCCCCATCAACAGGACACAGGCGTTCTCCATGAAGATAGGAACAGCAAGCg                              <  1:3242847/64‑1 (MQ=255)
                               aaaGTGCTCAGGACCCCCATCAACAGGACACAGGCGTTCTCCATGAAGATAGGAACAGCAAGCg                              <  1:3312380/64‑1 (MQ=255)
                               aaaGTGCTCAGGACCCCCATCAACAGGACACAGGCGTTCTCCATGAAGATAGGAACAGCAAGCg                              <  1:1673703/64‑1 (MQ=255)
                               aaaGTGCTCAGGACCCCCATCAACAGGACACAGGCGTTCTCCATGAAGATAGGAACAGCAAGCg                              <  1:3380068/64‑1 (MQ=255)
                               aaaGTGCTCAGGACCCCCATCAACAGGACACAGGCGTTCTCCATGAAGATAGGAACAGCAAGCg                              <  1:153266/64‑1 (MQ=255)
                               aaaGTGCTCAGGACCCCCATCAACAGGACACAGGCGTTCTCCATGAAGATAGGAACAGCAAGCg                              <  1:2769259/64‑1 (MQ=255)
                               aaaGTGCTCAGGACCCCCATCAACAGGACACAGGCGTTCTCCATGAAGATAGGAACAGCAAGCg                              <  1:652061/64‑1 (MQ=255)
                               aaaGTGCTCAGGACCCCCATCAACAGGACACAGGCGTTCTCCATGAAGATAGGAACAGCAAGCg                              <  1:2113885/64‑1 (MQ=255)
                               aaaGTGCTCAGGACCCCCATCAACAGGACACAGGCGTTCTCCATGAAGATAGGAACAGCAAGCg                              <  1:2077922/64‑1 (MQ=255)
                               aaaGTGCTCAGGACCCCCATCAACAGGACACAGGCGTTCTCCATGAAGATAGGAACAGCAAGCg                              <  1:2034303/64‑1 (MQ=255)
                                           aCCCCCATCAACAGGACACAGGCATTCTCCATGAAGATAGGAACAGCAAGCGGGGTTATCTCGCGCCAg             <  1:1341148/69‑1 (MQ=255)
                                            cccccATCAACAGGACACAGGCATTCTCCATGAAGATAGGAACAGCAAGCGGGGTTATCTCGCGCCAGa            <  1:2067688/69‑1 (MQ=255)
                                             ccccATCAACAGGACACAGGCATTCTCCATGAAGATAGGAACAGCAAGCGGGGTTATCTCGCGCCAGa            <  1:2247676/68‑1 (MQ=255)
                                                cATCAACAGGACACAGGCATTCTCCATGAAGATAGGAACAGCAAGCGGGGTTATCTCGCGCCAGaacaa        <  1:1013821/69‑1 (MQ=255)
                                                cATCAACAGGACACAGGCATTCTCCATGAAGATAGGAACAGCAAGCGGGGTTATCTCGCGCCAGaacaa        <  1:2608332/69‑1 (MQ=255)
                                                    aaCAGGACACAGGCATTCTCCATGAAGATAGGAACAGCAAGCGGGGTTATCTCGCGCCAGAACAacac     >  1:2020347/1‑68 (MQ=255)
                                                    aaCAGGACACAGGCATTCTCCATGAAGATAGGAACAGCAAGCGGGGTTATCTCGCGCCAGAACAACACt    >  1:1687313/1‑69 (MQ=255)
                                                    aaCAGGACACAGGCATTCTCCATGAAGATAGGAACAGCAAGCGGGGTTATCTCGCGCCAGAACAACACt    >  1:1946391/1‑69 (MQ=255)
                                                    aaCAGGACACAGGCATTCTCCATGAAGATAGGAACAGCAAGCGGGGTTATCTCGCGCCAGAACAACACt    >  1:2252507/1‑69 (MQ=255)
                                                              aGGCGTTCTCCATGAAGATAGGAACAGCAAGCGGGGTTATCTCGCGCCAGAACAACACttt  >  1:3280502/1‑61 (MQ=255)
                                                              aGGCGTTCTCCATGAAGATAGGAACAGCAAGCGGGGTTATCTCGCGCCAGAACAACACttt  >  1:3200621/1‑61 (MQ=255)
                                                              aGGCGTTCTCCATGAAGATAGGAACAGCAAGCGGGGTTATCTCGCGCCAGAACAACACttt  >  1:3068264/1‑61 (MQ=255)
                                                              aGGCGTTCTCCATGAAGATAGGAACAGCAAGCGGGGTTATCTCGCGCCAGAACAACACttt  >  1:1561256/1‑61 (MQ=255)
                                                              aGGCGTTCTCCATGAAGATAGGAACAGCAAGCGGGGTTATCTCGCGCCAGAACAACACttt  >  1:1557224/1‑61 (MQ=255)
                                                              aGGCGTTCTCCATGAAGATAGGAACAGCAAGCGGGGTTATCTCGCGCCAGAACAACACttt  >  1:2982774/1‑61 (MQ=255)
                                                              aGGCGTTCTCCATGAAGATAGGAACAGCAAGCGGGGTTATCTCGCGCCAGAACAACACttt  >  1:234588/1‑61 (MQ=255)
                                                                  |                                                        
CATCGCATCTTTTCCCAGCCAGCTGACCAGAAAAGTGCTCAGAACCCCCATCAACAGGACACAGGCATTCTCCATGAAGATAGGAACAGCAAGCGGGGTTATCTCGCGCCAGAACAACACTTT  >  minE/1292777‑1292899

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: