Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,417,696 C→A 17.3% 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.3% 2.4 ‑3.9 0.043intergenic (‑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 (20/13):  total (22/21)
Fisher's exact test strand distribution p-value = 2.05e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTT  >  REL606/2417661‑2417730
                                   |                                  
cccAGGAACGTTGGAGGCCTAATTTAAATGAGGTGc                                    <  1:1622347/36‑1 (MQ=255)
 ccAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCa                                   <  1:1119911/36‑1 (MQ=255)
 ccAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCa                                   >  1:2105789/1‑36 (MQ=255)
  cAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCaa                                  <  1:2539507/36‑1 (MQ=255)
  cAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCaa                                  <  1:642189/36‑1 (MQ=255)
   aGGAACGTTGGAGGCCTAATTTAAATGAGGTGCaaa                                 >  1:5079918/1‑36 (MQ=255)
   aGGAACGTTGGAGGCCTAATTTAAATGAGGTGCaaa                                 <  1:1890088/36‑1 (MQ=255)
    ggAACGTTGGAGGCCTAATTTAAATGAGGTGCaaaa                                <  1:4451904/36‑1 (MQ=255)
    ggAACGTTGGAGGCCTAATTTAAATGAGGTGCaaaa                                >  1:2050530/1‑36 (MQ=255)
    ggAACGTTGGAGGCCTAATTTAAATGAGGTGCaaaa                                >  1:2961002/1‑36 (MQ=255)
     gAACGTTGGAGGCCTAATTTAAATGAGGTGCaaaaa                               <  1:5677536/36‑1 (MQ=255)
      aaCGTTGGAGGCCTAATTTAAATGAGGTGCaaaaaa                              <  1:2353005/36‑1 (MQ=255)
        cGTTGGAGGCCTAATTTAAATGAGGTGCaaaaaaaa                            <  1:5241325/36‑1 (MQ=255)
         gTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAaac                           <  1:2380866/36‑1 (MQ=255)
         gTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAaac                           >  1:2694271/1‑36 (MQ=255)
           tGGAGGCCTAATTTAAATGAGGTGCAAAAAAaacaa                         >  1:5626442/1‑36 (MQ=255)
           tGGAGGCCTAATTTAAATGAGGTGCAAAAAAaacaa                         >  1:3513509/1‑36 (MQ=255)
            ggAGGCCTAATTTAAATGAGGTGTAAAAAAaacaac                        >  1:2755648/1‑36 (MQ=255)
            ggAGGCCTAATTTAAATGAGGTGCAAAAAAaacaaa                        >  1:4522012/1‑35 (MQ=255)
              aGGCCTAATTTAAATGAGGTGCAAAAAAAACAACag                      <  1:1008125/36‑1 (MQ=255)
              aGGCCTAATTTAAATGAGGTGCAAAAAAAACAACag                      >  1:5290288/1‑36 (MQ=255)
                gCCTAATTTAAATGAGGTGCAAAAAAAACAACagag                    <  1:3335298/36‑1 (MQ=255)
                  ttAATTAAAAGGAGGGAAAAAAAAAACAACAGAGGa                  <  1:5340143/35‑1 (MQ=12)
                   tAATTTAAATGAGGTGCAAAAAAAACAACAGAGGAc                 >  1:4540312/1‑36 (MQ=255)
                   tAATTTAAATGAGGTGC.AAAAAAACAACAGAGGACt                >  1:2283357/1‑36 (MQ=37)
                   tAATTTAAATGAGGTGAAAAAAAAACAACAGAGGAc                 >  1:1461720/1‑36 (MQ=255)
                   tAATTTAAATGAGGTGAAAAAAAAACAACAGAGGAc                 <  1:2521654/36‑1 (MQ=255)
                    aaTTTAAATGAGGTGCAAAAAAAACTACAGAGGACt                >  1:5850106/1‑36 (MQ=255)
                    aaTTTAAATGAGGTGAAAAAAAAACAACAGAGGACt                <  1:5507202/36‑1 (MQ=255)
                     aTTTAAATGAGGTGCAAAAAAAACAACAGAGGACTa               >  1:4800458/1‑36 (MQ=255)
                     aTTTAAATGAGGTGCAAAAAAAACAACAGAGGACTa               <  1:4894087/36‑1 (MQ=255)
                      tttAAATGAGGTGAAAAAAAAACAACAGAGGACTaa              <  1:1841224/36‑1 (MQ=255)
                          aaTGAGGTGCAAAAAAAACAACAGAGGACTAAACCg          >  1:4361919/1‑36 (MQ=255)
                           aTGAGGTGCAAAAAAAACAACAGAGGACTAAACCgg         >  1:4186637/1‑35 (MQ=255)
                           aTGAGGTGCAAAAAAAACAACAGAGGACTAAACCGt         >  1:119428/1‑36 (MQ=255)
                             gAGGTGCAAAAAAAACAACAGAGGACTAAACCGTgg       >  1:5596780/1‑36 (MQ=255)
                             gAGGTGAAAAAAAAACAACAGAGGACTAAACCGTgg       <  1:1669997/36‑1 (MQ=255)
                              aGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGc      >  1:4113315/1‑36 (MQ=255)
                              aGGTGAAAAAAAAACAACAGAGGACTAAACCGTGGc      <  1:120282/36‑1 (MQ=255)
                               ggTGCAAAAAAAACAACAGAGGACTAAACCGTGGCt     >  1:387667/1‑36 (MQ=255)
                               ggTGAAAAAAAAACAACAGAGGTCTAAACCGTGGCt     <  1:1736360/36‑1 (MQ=37)
                               ggTGAAAAAAAAACAACAGAGGACTAAACCGTGGCt     <  1:3507770/36‑1 (MQ=255)
                                 tGCAAAAAAAACACCAGCGGACTAAACCGTGGCttt   <  1:5236218/36‑1 (MQ=37)
                                  gCAAAAAAAACAACAGAGGACTAAACCGTGGCtttt  >  1:554973/1‑36 (MQ=255)
                                   |                                  
CCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTT  >  REL606/2417661‑2417730

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG