Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,558,537 C→A 34.1% P141P (CCG→CCT emrK ← EmrKY‑TolC multidrug resistance efflux pump, membrane fusion protein component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,558,5370CA34.1% 40.9 / 25.3 41P141P (CCG→CCTemrKEmrKY‑TolC multidrug resistance efflux pump, membrane fusion protein component
Reads supporting (aligned to +/- strand):  ref base C (17/10);  new base A (7/7);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

ACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATA  >  minE/1558474‑1558582
                                                               |                                             
aCTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCa                                           <  1:3420410/68‑1 (MQ=255)
aCTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCa                                           <  1:2926422/68‑1 (MQ=255)
aCTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCa                                           <  1:2385677/68‑1 (MQ=255)
   tattaACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTc                                        <  1:137985/68‑1 (MQ=255)
   tattaACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTc                                        <  1:2603351/68‑1 (MQ=255)
   tattaACGTATCTTTGGTATGCTCCAGCGCTTCTTTTGAAATAACCCCTTGCTTCGCTAAAGGCAcac                                        <  1:1329654/68‑3 (MQ=255)
   tattaACGTATCTTTGGTATGCTCCAGCGCTTCTTTTGAAATAACCCCTTGCTTCGCTAAAGGCAcac                                        <  1:593801/68‑3 (MQ=255)
   tattaACGTATCTTTGGTATGCTCCAGCGCTTCTTTTGAAATAACCCCTTGCTTCGCTAAAGGCAcac                                        <  1:2925479/68‑3 (MQ=255)
   tattaACGTATCTTTGGTATGCTCCAGCGCTTCTTTTGAAATAACCCCTTGCTTCGCTAAAGGCAcac                                        <  1:2532172/68‑3 (MQ=255)
   tattaACGTATCTTTGGTATGCTCCAGCGCTTCTTTTGAAATAACCCCTTGCTTCGCTAAAGGCAcac                                        <  1:2226633/68‑3 (MQ=255)
   tattaACGTATCTTTGGTATGCTCCAGCGCTTCTTTTGAAATAACCCCTTGCTTCGCTAAAGGCAcac                                        <  1:1861090/68‑3 (MQ=255)
   tattaACGTATCTTTGGTATGCTCCAGCGCTTCTTTTGAAATAACCCCTTGCTTCGCTAAAGGCAcac                                        <  1:1673087/68‑3 (MQ=255)
       aaCGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACgg                                             >  1:979809/1‑59 (MQ=255)
       aaCGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACgg                                             >  1:1530326/1‑59 (MQ=255)
       aaCGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACgg                                             >  1:3153439/1‑59 (MQ=255)
       aaCGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACgg                                             >  1:3087051/1‑59 (MQ=255)
       aaCGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACgg                                             >  1:3084039/1‑59 (MQ=255)
       aaCGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACgg                                             >  1:2688921/1‑59 (MQ=255)
       aaCGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACgg                                             >  1:2411394/1‑59 (MQ=255)
       aaCGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACg                                    >  1:1646379/1‑68 (MQ=255)
       aaCGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACg                                    >  1:1839692/1‑68 (MQ=255)
       aaCGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACg                                    >  1:2304429/1‑68 (MQ=255)
       aaCGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACg                                    >  1:1843721/1‑68 (MQ=255)
            aTCTTTGGTATGCTCCAGGGCTTCTTTTGAAATAACCCCTTGCTTCGCTAAAGGCACACGACGGTtata                              >  1:304207/1‑69 (MQ=255)
            aTCTTTGGTATGCTCCAGCGCTTCTTTTGAAATAACCCCTTGCTTCGCTAAAGGc                                            >  1:1863320/1‑55 (MQ=255)
            aTCTTTGGTATGCTCCAGCGCTTCTTTTGAAATAACCCCTTGCTTCGCTAAAGGCACACGATGGTtat                               >  1:1660258/1‑68 (MQ=255)
            aTCTTTGGTATGCTCCAGCGCTTCTTTTGAAATAACCCCTTGCTTCGCTAAAGGCACACGACGGTtata                              >  1:2379192/1‑69 (MQ=255)
            aTCTTTGGTATGCTCCAGCGCTTCTTTTGAAATAACCCCTTGCTTCGCTAAAGGCACACGACGGTtat                               >  1:794306/1‑68 (MQ=255)
            aTCTTTGGTATGCTCCAGCGCTTCTTTTGAAATAACCCCTTGCTTCGCTAAAGGCACACGACGGTtat                               >  1:641910/1‑68 (MQ=255)
            aTCTTTCGTATGCTCCAGCGCTTCTTTTGAAATAACCCCTTGCTTCGCTAAAGGCACACGACGGTtat                               >  1:438027/1‑68 (MQ=255)
                        cTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAg                   <  1:2246930/68‑1 (MQ=255)
                        cTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAg                   <  1:884130/68‑1 (MQ=255)
                        cTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGa                  >  1:428407/1‑69 (MQ=255)
                        cTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGa                  >  1:2415576/1‑69 (MQ=255)
                        cTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGa                  >  1:10865/1‑69 (MQ=255)
                        cTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGa                  >  1:164577/1‑69 (MQ=255)
                                  tCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATa         >  1:540173/1‑68 (MQ=255)
                                  tCTTTTGAAATAACCCCCTGCTTCGCTAACGCCACTCGACGGTTATAATCTTCTAAAGATTGTTGATa         >  1:2277848/1‑68 (MQ=255)
                                         aaaTAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATa  <  1:865928/68‑1 (MQ=255)
                                         aaaTAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATa  <  1:2566742/68‑1 (MQ=255)
                                                            tAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATa         <  1:2639341/42‑1 (MQ=255)
                                                               |                                             
ACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCCCTGCTTCGCTAACGGCACTCGACGGTTATAATCTTCTAAAGATTGTTGATACTGAATA  >  minE/1558474‑1558582

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: