Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,487,979 G→A 16.9% K267K (AAG→AAA yjhC → predicted oxidoreductase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,487,9790G→A16.9% 5.5 ‑4.0 0.052801yjhCpredicted oxidoreductase
Reads supporting (aligned to +/- strand):  new base (9/0):  ref base (23/19):  total (32/20)
Fisher's exact test strand distribution p-value = 1.83e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ACAGCAGGGTCACTTAGGATTGGCGGTCAGACAAAGCATTTTTTGGTCCATGAAACACAAGAAGAAGATG  >  REL606/4487944‑4488013
                                   |                                  
aCAGCAGGGTCACTTAGGATTGGCGGTCAGACAAAg                                    <  1:389928/36‑1 (MQ=255)
 cagcagGGTCACTTAGGATTGGCGGTCAGACAAAGc                                   <  1:7265343/36‑1 (MQ=255)
  agcagGGTCACTTAGGATTGGCGGTCCGGCAAAACa                                  >  1:2280295/1‑36 (MQ=25)
    cagGGTCACTTAGGATTGGCGGTCAGACAAATCAtt                                <  1:3586366/36‑1 (MQ=255)
    cagGGTCACTTAGGATTGGCGGTCAGACAAAGCAtt                                <  1:5723191/36‑1 (MQ=255)
     agGGTCACTTAGGATTGGCGGTCAGACAAAGCAttt                               <  1:4514608/36‑1 (MQ=255)
     agGGTCACTTAGGATTGGCGGTCAGACAAAAAAttt                               >  1:6053128/1‑36 (MQ=37)
      gggTCACTTAGGATTGGCGGTCCTACACAACAtttt                              >  1:3274085/1‑36 (MQ=19)
      gggTCACTTAGGATTGGCGGTCAGACAAAGCAtttt                              <  1:759252/36‑1 (MQ=255)
      gggTCACTTAGGATTGGCGGTCAGACAAAGCAtttt                              <  1:2990753/36‑1 (MQ=255)
      gggTCACTTAGGATTGGCGGGCAGACAAAGAAattt                              >  1:7163076/1‑36 (MQ=25)
       ggTCACTTAGGATTGGCGGTCTGGCAAAACAttttt                             >  1:4788557/1‑36 (MQ=25)
       ggTCACTTAGGATTGGCGGTCAGAAAAAACAttttt                             >  1:1340516/1‑36 (MQ=37)
        gTCACTTAGGATTGGCGGTCGGAGAAAACAtttttt                            >  1:228145/1‑36 (MQ=37)
        gTCACTTAGGATTGGCGGTCAGACAAAGCAtttttt                            <  1:4283447/36‑1 (MQ=255)
         tcccTTAGGATTGGCGGTCAGACAAAGCATTTTTTg                           <  1:4839178/33‑1 (MQ=255)
         tCACTTAGGATTGGCGGTCAGACAAAGAAttttttt                           >  1:3523532/1‑35 (MQ=255)
         tCACTTAGGATTGGCGGTCAGAAAAAAAAttttttt                           >  1:5554251/1‑35 (MQ=21)
           aCTTAGGATTGGCGGTCCGACAAAACAttttttttt                         >  1:1683584/1‑33 (MQ=37)
           aCTTAGGATTGGCGGTCAGACAAAGAAttttttttt                         >  1:5893442/1‑33 (MQ=37)
           aCTTAGGATTGGCGGTCAGACAAAACAttttttttt                         >  1:778718/1‑33 (MQ=37)
               aGGATTGGCGGTCAGACAAAGCATTTTTTGGTCCAt                     <  1:1391666/36‑1 (MQ=255)
                 gATTGGCGGTCAGACAAAGCATTTTTTGGTCCATGa                   <  1:4584675/36‑1 (MQ=255)
                 gATTGGCGGTCAGACAAAGCATTTTTTGGTCCATGa                   <  1:1817047/36‑1 (MQ=255)
                 gATTAGCGGTCAGACAAAGCACTTTTTGTTCCATGa                   >  1:4233391/1‑36 (MQ=25)
                    tGGCGGTCAGACAAAGCATTTTTTGGTCCATGAAac                <  1:4622430/36‑1 (MQ=255)
                     ggCGGTCAGACAAAGCATTTTTTTGTCCATTAAaca               >  1:729689/1‑36 (MQ=37)
                      gCGGTCAGACAAAGCATTTTTTGGTCCATGAAacac              <  1:7083476/36‑1 (MQ=255)
                       cGGTCAGACAAAGCATTTTTTGGTCCATGAAAcaca             >  1:4979108/1‑36 (MQ=255)
                       cGGTCAGACAAAGCATTTTTTGGTCCATGAAAcaca             >  1:5228770/1‑36 (MQ=255)
                         gTCAGACAAAGCATTTTTTGGTCCATGAAACACaag           >  1:3284705/1‑36 (MQ=255)
                         gTCAGACAAAGCATTTTTTGGTCCATGAAACACaag           >  1:5746271/1‑36 (MQ=255)
                         gTCAGACAAAGCATTTTTTGGTCCATGAAACACaag           >  1:2433261/1‑36 (MQ=255)
                          tCAGACAAAGCATTTTTTGGTCCATGAAACACaaga          <  1:5765591/36‑1 (MQ=255)
                          tCAGACAAAGCATTTTTTGGTCCATGAAACACaaga          >  1:1151109/1‑36 (MQ=255)
                          tCAGACAAAGCATTTTTTGGTCCATGAAACACaaga          >  1:54607/1‑36 (MQ=255)
                          tCAGACAAAGCATTTTTTGGTCCATGAAACACaaga          <  1:1522508/36‑1 (MQ=255)
                          tCAGACAAAGCATTTTTTGGTCCATGAAACACaaga          <  1:1874541/36‑1 (MQ=255)
                           cAGACAAAGCATTTTTTGGTCCATGAAACACaagaa         <  1:1290880/36‑1 (MQ=255)
                            aGACAAAGCATTTTTTGGTCCATGAAACACaagaag        <  1:5550766/36‑1 (MQ=255)
                            aGACAAAGCATTTTTTGGTCCATGAAACACaagaag        <  1:3445724/36‑1 (MQ=255)
                            aGACAAAGCATTTTTTGGTCCATGAAACACaagaag        >  1:271391/1‑36 (MQ=255)
                              aCAAAGCATTTTTTGGTCCATGAAACACaagaagaa      >  1:6754166/1‑36 (MQ=255)
                              aCAAAGCATTTTTTGGTCCATGAAACACaagaagaa      >  1:6753496/1‑36 (MQ=255)
                               cAAAGCATTTTTTGGTCCATGAAACACaagaagaag     >  1:3939208/1‑36 (MQ=255)
                               cAAAGCATTTTTTGGCCCATGAAACACaagaagaag     >  1:7349097/1‑36 (MQ=255)
                                aaaGCATTTTTTTGTCCTTTAAACACCataagaaga    >  1:1809169/1‑36 (MQ=9)
                                aaaGCATTTTTTGGTCCATGAAACACAagaagaaga    >  1:6533476/1‑36 (MQ=255)
                                aaaGCATTTTTTGGTCCATGAAACACAagaagaaga    >  1:4365752/1‑36 (MQ=255)
                                 aaGCATTTTTTGGTCCATGAAACACAAGAAGAAgat   <  1:6605174/36‑1 (MQ=255)
                                 aaGCATTTTTTGGTCCATGAAACACAAGAAGAAgat   >  1:2491068/1‑36 (MQ=255)
                                  aGCATTTTTTGGTCCATGAAACACAAGAAGAAgatg  >  1:5641307/1‑36 (MQ=255)
                                  aGCATTTTTTGGTCCATGAAACACAAGAAGAAgatg  >  1:895734/1‑36 (MQ=255)
                                   |                                  
ACAGCAGGGTCACTTAGGATTGGCGGTCAGACAAAGCATTTTTTGGTCCATGAAACACAAGAAGAAGATG  >  REL606/4487944‑4488013

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 20 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG