Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,410,527 A→G 18.8% S222S (AGT→AGC bcr ← bicyclomycin/multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,410,5270AG18.8% 51.5 / 5.3 37S222S (AGT→AGCbcrbicyclomycin/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base A (10/20);  new base G (4/3);  total (14/23)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

TGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGTTACC  >  minE/1410459‑1410583
                                                                    |                                                        
tGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCa                                                          >  1:862315/1‑69 (MQ=255)
tGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCa                                                          >  1:1167716/1‑69 (MQ=255)
tGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCa                                                          >  1:2852512/1‑69 (MQ=255)
tGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCa                                                          >  1:2674248/1‑69 (MQ=255)
tGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCa                                                          >  1:1967114/1‑69 (MQ=255)
tGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCa                                                          >  1:77975/1‑69 (MQ=255)
tGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCa                                                          >  1:1959812/1‑69 (MQ=255)
  gTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACt                                                        >  1:675852/1‑69 (MQ=255)
        ttCAATGTAGACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCGCTGGCCAg                                                  <  1:165861/69‑1 (MQ=255)
        ttCAATGTAGACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCGCTGGCCAg                                                  <  1:2694811/69‑1 (MQ=255)
         tCAATGTAGACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCGCTGGCCAg                                                  <  1:2084130/68‑1 (MQ=255)
              ataAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTa                                            >  1:1376948/1‑69 (MQ=255)
              ataAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTa                                            >  1:122269/1‑69 (MQ=255)
              ataAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTa                                            >  1:2564801/1‑69 (MQ=255)
              ataAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTa                                            >  1:259921/1‑69 (MQ=255)
              ataAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTa                                            >  1:334708/1‑69 (MQ=255)
               taAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGGAAGCATGTAg                                           <  1:614709/69‑1 (MQ=255)
               taAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGc                                                 >  1:1690492/1‑63 (MQ=255)
               taAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAg                                           <  1:1437483/69‑1 (MQ=255)
                 aacaaacGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAg                                           <  1:218531/67‑1 (MQ=255)
                    aaacGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCGCTGGCCAGCATATAACTCAg                                      >  1:3055277/1‑69 (MQ=255)
                    aaacGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCGCTGGCCAGCATATAACTCAg                                      >  1:3105916/1‑69 (MQ=255)
                    aaacGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCGCTGGCCAGCATATAACTCAg                                      >  1:1292392/1‑69 (MQ=255)
                    aaacGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCGCTGGCCAGCATATAACTCAg                                      >  1:52741/1‑69 (MQ=255)
                    aaacGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAg                                      <  1:240879/69‑1 (MQ=255)
                    aaacGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAg                                      <  1:1774957/69‑1 (MQ=255)
                            cGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGtt                              <  1:866160/69‑1 (MQ=255)
                            cGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGtt                              <  1:237069/69‑1 (MQ=255)
                            aggCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGtt                              <  1:2067801/68‑1 (MQ=255)
                                gcTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGtt                              >  1:314877/1‑65 (MQ=255)
                                gcTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGtt                              >  1:3264040/1‑65 (MQ=255)
                                gcTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGtt                              >  1:3348511/1‑65 (MQ=255)
                                                     ggCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACTTTTATGGCGGAACAGCGCCGCAAAGt      <  1:2645173/68‑1 (MQ=255)
                                                     ggCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGt      <  1:3401870/68‑1 (MQ=255)
                                                     ggCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGt      <  1:1086822/68‑1 (MQ=255)
                                                     ggCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGt      <  1:1395617/68‑1 (MQ=255)
                                                     ggCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGt      <  1:2003366/68‑1 (MQ=255)
                                                     ggCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGt      <  1:440767/68‑1 (MQ=255)
                                                     ggCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGt      <  1:2938910/68‑1 (MQ=255)
                                                     ggCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGt      <  1:2854869/68‑1 (MQ=255)
                                                     ggCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGt      <  1:291539/68‑1 (MQ=255)
                                                        aaaGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGTTAcc  <  1:843184/69‑1 (MQ=255)
                                                         aaGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGTTAcc  <  1:3269000/68‑1 (MQ=255)
                                                           gCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGTAACAGCGCCGCAAAGt      <  1:1118999/62‑1 (MQ=255)
                                                                    |                                                        
TGGTTAATTTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGTTACC  >  minE/1410459‑1410583

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