Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,410,533 A→C 19.4% L220L (CTT→CTG bcr ← bicyclomycin/multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,410,5330AC19.4% 57.6 / 10.5 36L220L (CTT→CTGbcrbicyclomycin/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base A (9/20);  new base C (4/3);  total (13/23)
Fisher's exact test for biased strand distribution p-value = 2.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-01

TTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGTTACC  >  minE/1410467‑1410583
                                                                  |                                                  
ttCAATGTAGACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCGCTGGCCAg                                                  <  1:165861/69‑1 (MQ=255)
ttCAATGTAGACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCGCTGGCCAg                                                  <  1:2694811/69‑1 (MQ=255)
 tCAATGTAGACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCGCTGGCCAg                                                  <  1:2084130/68‑1 (MQ=255)
      ataAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTa                                            >  1:122269/1‑69 (MQ=255)
      ataAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTa                                            >  1:334708/1‑69 (MQ=255)
      ataAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTa                                            >  1:1376948/1‑69 (MQ=255)
      ataAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTa                                            >  1:259921/1‑69 (MQ=255)
      ataAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTa                                            >  1:2564801/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:3105916/1‑69 (MQ=255)
            aaacGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCGCTGGCCAGCATATAACTCAg                                      >  1:3055277/1‑69 (MQ=255)
            aaacGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCGCTGGCCAGCATATAACTCAg                                      >  1:52741/1‑69 (MQ=255)
            aaacGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCGCTGGCCAGCATATAACTCAg                                      >  1:1292392/1‑69 (MQ=255)
            aaacGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAg                                      <  1:240879/69‑1 (MQ=255)
            aaacGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAg                                      <  1:1774957/69‑1 (MQ=255)
                    cGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGtt                              <  1:237069/69‑1 (MQ=255)
                    cGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGtt                              <  1:866160/69‑1 (MQ=255)
                    aggCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGtt                              <  1:2067801/68‑1 (MQ=255)
                        gcTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGtt                              >  1:3348511/1‑65 (MQ=255)
                        gcTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGtt                              >  1:3264040/1‑65 (MQ=255)
                        gcTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGtt                              >  1:314877/1‑65 (MQ=255)
                                             ggCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACTTTTATGGCGGAACAGCGCCGCAAAGt      <  1:2645173/68‑1 (MQ=255)
                                             ggCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGt      <  1:2938910/68‑1 (MQ=255)
                                             ggCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGt      <  1:291539/68‑1 (MQ=255)
                                             ggCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGt      <  1:2854869/68‑1 (MQ=255)
                                             ggCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGt      <  1:1086822/68‑1 (MQ=255)
                                             ggCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGt      <  1:3401870/68‑1 (MQ=255)
                                             ggCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGt      <  1:440767/68‑1 (MQ=255)
                                             ggCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGt      <  1:2003366/68‑1 (MQ=255)
                                             ggCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGt      <  1:1395617/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)
                                                                  |                                                  
TTCAATATAAACAAACGGTCCGGCGCTTAAGAATGAGAACATCCCGGCAAAGCTGAAACCACTGGCAAGCATGTAGCTCAGGACACGTTTATGGCGGAACAGCGCCGCAAAGTTACC  >  minE/1410467‑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: