Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,417,696 C→A 17.1% intergenic (‑130/‑286) emrK ← / → evgA EmrKY‑TolC multidrug resistance efflux pump, membrane fusion protein component/DNA‑binding response regulator in two‑component regulatory system with EvgS

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,417,6960C→A17.1% 4.4 ‑2.7 0.047intergenic (‑130/‑286)emrK/evgAEmrKY‑TolC multidrug resistance efflux pump, membrane fusion protein component/DNA‑binding response regulator in two‑component regulatory system with EvgS
Reads supporting (aligned to +/- strand):  new base (1/8):  ref base (18/20):  total (19/28)
Fisher's exact test strand distribution p-value = 6.40e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTT  >  REL606/2417662‑2417730
                                  |                                  
ccAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCa                                   <  1:6094578/36‑1 (MQ=255)
ccAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCa                                   >  1:1489182/1‑36 (MQ=255)
ccAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCa                                   >  1:5436201/1‑36 (MQ=255)
ccAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCa                                   <  1:44148/36‑1 (MQ=255)
 cAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCaa                                  >  1:555256/1‑36 (MQ=255)
 cAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCaa                                  <  1:3372337/36‑1 (MQ=255)
  aGGAACGTTGGAGGCCTAATTTAAATGAGGTGCaaa                                 <  1:5161186/36‑1 (MQ=255)
    gAACGTTGGAGGCCTAATTTAAATGAGGTGCaaaaa                               <  1:5585787/36‑1 (MQ=255)
    gAACGTTGGAGGCCTAATTTAAATGAGGTGCaaaaa                               <  1:5328736/36‑1 (MQ=255)
    gAACGTTGGAGGCCTAATTTAAATGAGGTGCaaaaa                               <  1:4400103/36‑1 (MQ=255)
     aaCGTTGGAGGCCTAATTTAAATGAGGTGCaaaaaa                              >  1:1124495/1‑36 (MQ=255)
      aCGTTGGAGGCCTAATTTAAATGAGGTGCaaaaaaa                             >  1:5044056/1‑36 (MQ=255)
        gTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAaac                           >  1:2396345/1‑36 (MQ=255)
        gTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAaac                           <  1:3171723/36‑1 (MQ=255)
         ttGGAGGCCTAATTTAAATGAGGTGCAAAAAAaaca                          >  1:4075567/1‑36 (MQ=255)
           tgAGGCCTAATTTAAATGAGGTGAAAAAAAaacaac                        <  1:2858849/35‑1 (MQ=255)
           ggAGGCCTAATTTAAATGACGTGAAAAAAAaacaac                        <  1:236244/36‑1 (MQ=38)
             aGGCCTAATTTAAATGAGGTGCAAAAAAAACAACag                      >  1:3476097/1‑36 (MQ=255)
             aGGCCTAATTTAAATGAGGTGCAAAAAAAACAACag                      >  1:7413443/1‑36 (MQ=255)
              ggCCTAATTTAAATGAGGTGCAAAAAAAACAACaga                     <  1:2657967/36‑1 (MQ=255)
               gCCTAATTTAAATGAGGTGCAAAAAAAACAACagag                    <  1:4180997/36‑1 (MQ=255)
                 cTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGa                  <  1:6321916/36‑1 (MQ=255)
                 cTAATTTAAATGAGGTGAAAAAAAAACAACAGAGGa                  <  1:2812689/36‑1 (MQ=255)
                 cTAATTTAAATGAGGTGAAAAAAAAACAACAGAGGa                  <  1:3668686/36‑1 (MQ=255)
                  tAATTTAAATGAGGTGCAAAAAAAACAACAGAGGAc                 <  1:2781896/36‑1 (MQ=255)
                  tAATTTAAATGAGGTGCAAAAAAAACAACAGAGGAc                 <  1:7902351/36‑1 (MQ=255)
                   aaTTTAAATGAGGTGCAAAAAAAACAACAGAGGACt                >  1:6026547/1‑36 (MQ=255)
                    aTTTAAATGAGGTGAAAAAAAAACAACAGAGGACTa               <  1:400327/36‑1 (MQ=255)
                     tttAAATGAGGTGAAAAAAAAACAACAGAGGACTaa              <  1:6571310/36‑1 (MQ=255)
                      ttAAATGAGGTGCAAAAAAAACAACAGAGGACTaaa             >  1:6072578/1‑36 (MQ=255)
                      ttAAATGAGGTGCAAAAAAAACAACAGAGGACTaaa             >  1:6403226/1‑36 (MQ=255)
                      ttAAATGAGGTGCAAAAAAAACAACAGAGGACTaaa             >  1:4497739/1‑36 (MQ=255)
                       tAAGTGAGGTGAAAAAAAAACAACAGAGGGCTAAAc            >  1:6324253/1‑36 (MQ=25)
                       tAAATGAGGTGCAAAAAAAACAACAGAGGACTAAAc            <  1:2156503/36‑1 (MQ=255)
                       tAAATGAGGTGAAAAAAAAACAACAGAGGACTAAAc            <  1:3811333/36‑1 (MQ=255)
                        aaaTGAGGTGCAAAAAAAACAACAGAGGACTAAAcc           >  1:2544692/1‑36 (MQ=255)
                        aaaTGAGGTGAAAAAAAAACAACAGAGGACTAAAcc           <  1:1127991/36‑1 (MQ=255)
                              ggTGCAAAAAAAACAACAGAGGACTAAACCGTGGCt     <  1:3883415/36‑1 (MQ=255)
                               gTGCAAAAAAAACAACAGAGGACTAAACCGTGGCtt    >  1:915065/1‑36 (MQ=255)
                               gTGCAAAAAAAACAACAGAGGACTAAACCGTGGCtt    >  1:95729/1‑36 (MQ=255)
                                tGCAAAAAAAACAACAGAGGACTAAGCCGTGGCttt   <  1:1447749/36‑1 (MQ=255)
                                tGCAAAAAAAACAACAGAGGACTAAACCGTGGCttt   >  1:1329553/1‑36 (MQ=255)
                                tGCAAAAAAAACAACAGAGGACTAAACCGTGGCttt   <  1:7580524/36‑1 (MQ=255)
                                 gCAAAAAAAACAACAGAGGACTAAACCGTGGCtttt  <  1:1664107/36‑1 (MQ=255)
                                 gCAAAAAAAACAACAGAGGACTAAACCGTGGCtttt  <  1:1818281/36‑1 (MQ=255)
                                 gCAAAAAAAACAACAGAGGACTAAACCGTGGCtttt  <  1:4518332/36‑1 (MQ=255)
                                 gCAAAAAAAACAACACAGGACTAAAACGTGGCtttt  >  1:2080047/1‑36 (MQ=37)
                                  |                                  
CCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTT  >  REL606/2417662‑2417730

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 37 ≤  ATCG