Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 134,149 A→T 11.6% intergenic (‑85/‑270) yacH ← / → acnB hypothetical protein/aconitate hydratase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606134,1490A→T11.6% 3.2 ‑3.6 0.032intergenic (‑85/‑270)yacH/acnBhypothetical protein/aconitate hydratase
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (10/18):  total (14/18)
Fisher's exact test strand distribution p-value = 2.78e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAG  >  REL606/134114‑134182
                                   |                                 
gTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTa                                   <  1:144496/36‑1 (MQ=255)
   cTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAAt                                <  1:6812764/36‑1 (MQ=255)
   cTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAAt                                <  1:6498573/36‑1 (MQ=255)
     gtgAATCCGAGTATAAAGAGGCGGTTGTTTTTAttt                              >  1:3826509/1‑36 (MQ=25)
     gtgAATCCGAGTATAAAGAGGCGGTAGTTTAAAttt                              <  1:6500169/36‑1 (MQ=255)
       gAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTg                            <  1:8308397/36‑1 (MQ=255)
       gAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTg                            <  1:2060414/36‑1 (MQ=255)
        aaTCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGa                           <  1:2483392/36‑1 (MQ=255)
        aaTCCGAGTATAAAGAGGCGGGAGTTTTTAtttttt                           >  1:1252433/1‑34 (MQ=21)
         aTCCGAGTATAAAGTGGCGGTGGTTTTCAttttttt                          >  1:593996/1‑33 (MQ=12)
         aTCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGAc                          <  1:7007135/36‑1 (MQ=255)
         aTCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGAc                          <  1:7006504/36‑1 (MQ=255)
            cGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTaa                       <  1:3317663/36‑1 (MQ=255)
             gAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAAt                      <  1:4520357/36‑1 (MQ=255)
                   aaaGAGGCGGTAGTTTAAATTTTGACTAATCTTggg                <  1:1552383/36‑1 (MQ=255)
                    aaGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGa               <  1:159049/36‑1 (MQ=255)
                     agagGCGGTAGTTTAAATTTTGACTAATCTTGGGAt              <  1:2613891/36‑1 (MQ=255)
                      gagGCGGTAGTTTTAATTTTTACTAATCTTTGGGtt             >  1:8008340/1‑36 (MQ=16)
                      gagGCGGTAGTTTAAATTTTGACTAATCTTGGGAtt             <  1:3733430/36‑1 (MQ=255)
                      gagGCGGTAGTTTAAATTTTGACTAATCTTGGGAtt             <  1:1939211/36‑1 (MQ=255)
                      gagGCGGTAGTTTAAATTTTGACTAATCTTGGGAtt             <  1:7814958/36‑1 (MQ=255)
                       agGCGGTAGTTTAAATTTTGACTAATCTTTGGATTc            >  1:6219617/1‑36 (MQ=255)
                          cGGTAGTTTAAATTTTGACTAATCTTGGGATTCGtt         >  1:6442479/1‑36 (MQ=255)
                          cGGTAGTTTAAATTTTGACTAATCTTGGGATTCGtt         >  1:2590786/1‑36 (MQ=255)
                           ggTAGTTTAAATTTTGACTAATCTTGGGATTCGTTg        <  1:2162691/36‑1 (MQ=255)
                            gTAGTTTAAATTTTGACTAATCTTGGGATTCGTTga       >  1:7500417/1‑36 (MQ=255)
                             tAGTTTAAATTTTGACTAATCTTGGGATTCGTTgag      >  1:384204/1‑36 (MQ=255)
                              aGTTTAAATTTTGACTAATCTTGGGATTCGTTgaga     >  1:4850204/1‑36 (MQ=255)
                              aGTTTAAATTTTGACTAATCTTGGGATTCGTTgaga     >  1:414061/1‑36 (MQ=255)
                              aGTTTAAATTTTGACTAATCTTGGGATTCGTTgaga     <  1:8057162/36‑1 (MQ=255)
                               gTTTAAATTTTGACTAATCTTGGGATTCGTTGAGaa    >  1:161652/1‑36 (MQ=255)
                                tttAAATTTTGACTAATCTTGGGATTCGTTGAGaaa   >  1:211666/1‑36 (MQ=255)
                                 ttAAATTTTGACTAATCTTGGGATTCGTTGAGAAAg  >  1:5823658/1‑36 (MQ=255)
                                   |                                 
GTGCTGTGAATCCGAGTATAAAGAGGCGGTAGTTTAAATTTTGACTAATCTTGGGATTCGTTGAGAAAG  >  REL606/134114‑134182

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG