Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 166,265 C→A 10.8% A437E (GCA→GAA)  hrpB → predicted ATP‑dependent helicase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606166,2650C→A10.8% 2.5 ‑4.2 ‑0.1371310hrpBpredicted ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  new base (0/4):  ref base (23/10):  total (23/14)
Fisher's exact test strand distribution p-value = 1.52e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.90e-01

TAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCA  >  REL606/166230‑166296
                                   |                               
tAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGc                                 >  1:1526924/1‑36 (MQ=255)
  gcgcGAAGAACGACGACGAAGCTGCTACCGCGGCaa                               <  1:6975496/36‑1 (MQ=255)
  gcgcGAAGAACGACGACGAAGCTGCTACCGCGGCaa                               <  1:4330276/36‑1 (MQ=255)
    gcgAAGAACGACGACGAAGCTGCTACCGCGGCaaaa                             >  1:6474231/1‑36 (MQ=255)
    gcgAAGAACGAAGACGAAGCTGCTACCGCGGCaaaa                             >  1:1649298/1‑36 (MQ=255)
    gagaagaaCGACGACGAAGCTGCTACCGCGGCaaaa                             <  1:674025/34‑1 (MQ=255)
     cgAAGAACGACGACGAAGCTGCTACCGCGGCaaaaa                            <  1:6886520/36‑1 (MQ=255)
     cgAAGAACGACGACGAAGCTGCTACCGCGGCaaaaa                            >  1:2409755/1‑36 (MQ=255)
     cgAAGAACGACGACGAAGCTGCTACCGCGGCaaaaa                            >  1:4178626/1‑36 (MQ=255)
      gaagaaCGACGACGAAGCTGCTACCGCGGCAAAAAt                           <  1:6217225/36‑1 (MQ=255)
      gaagaaCGACGACGAAGCTGCTACCGCGGCAAAAAt                           <  1:5736391/36‑1 (MQ=255)
       aagaaCGACGACGAAGCTGCTACCGCGGCAAAAAtt                          <  1:4696068/36‑1 (MQ=255)
       aagaaCGACGACGAAGCTGCTACCGCGGCAAAAAtt                          >  1:3056454/1‑36 (MQ=255)
        agaaCGACGACGAAGCTGCTACCGCGGCAAAAATTg                         >  1:97417/1‑36 (MQ=255)
         gaaCGACGACGAAGCTGCTACCGCGGCAAAAATTgc                        <  1:7403629/36‑1 (MQ=255)
         gaaCGACGACGAAGCTGCTACCGCGGCAAAAATTgc                        >  1:2724560/1‑36 (MQ=255)
          aaCGACGACGAAGCTGCTACCGCGGCAAAAATTgcc                       >  1:7078150/1‑36 (MQ=255)
          aaCGACGACGAAGCTGCTACCGCGGCAAAAATTgcc                       <  1:2512970/36‑1 (MQ=255)
           acgacgacgAAGCTGCTACCGCGGCAAAAATTgccg                      >  1:2287674/1‑36 (MQ=255)
           acgacgacgAAGCTGCTACCGCGGCAAAAATTgccg                      >  1:5095587/1‑36 (MQ=255)
           acgacgacgAAGCTGCTACCGCCGCTAAAATTgccg                      >  1:5661233/1‑36 (MQ=37)
             gacgacgaAGCTGCTACCGCGGCAAAAATTgccgcc                    >  1:2371941/1‑36 (MQ=255)
             gacgacgaAGCTGCTACCGCGGCAAAAATTgccgcc                    >  1:3006000/1‑36 (MQ=255)
              acgacgaAGCTGCTACCGCGGCAAAAATTGCCGCCa                   >  1:4364946/1‑36 (MQ=255)
              acgacgaAGCTGCTACCGCGGCAAAAATTGCCGCCa                   <  1:542780/36‑1 (MQ=255)
                aaagaAGCTGCTACCGGGGAAAAAATTGCCGCCAtt                 <  1:2428334/33‑1 (MQ=21)
                 acgaAGCTGCTACCGCGGCAAAAATTGCCGCCATtc                >  1:1687846/1‑36 (MQ=255)
                   gaAGCTGCTACCGCGGCAAAAATTGCCGCCATtctc              >  1:2790317/1‑36 (MQ=255)
                      gctgctACCGCGGCAAAAATTGCCGCCATTCTCgta           >  1:5205584/1‑34 (MQ=255)
                       ctgctACCGCGGCAAAAATTGCCGCCATTCTCGAag          >  1:1970918/1‑36 (MQ=255)
                           tACCGCGGCAAAAATTGCCGCCATTCTCGAAGAgcc      >  1:1239840/1‑36 (MQ=255)
                           tACCGCGGCAAAAATTGCCGCCATTCTCGAAGAgcc      >  1:2814249/1‑36 (MQ=255)
                            aCCGCGGCAAAAATTGCCGCCATTCTCGAAGAgccg     >  1:1110900/1‑36 (MQ=255)
                             gcgcgGAAAAAATTGCCGCCATTCTCGAAGAgccgc    <  1:3491030/35‑1 (MQ=255)
                              cgcgGCAAAAATTGCCGCCATTCTCGAAGAgccgcc   >  1:2334177/1‑36 (MQ=255)
                               gcgGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCa  >  1:6099551/1‑36 (MQ=255)
                               gcgGAAAAAATTGCCGCCATTCTCGAAGAGCCGCCa  <  1:537525/36‑1 (MQ=255)
                               gcgGAAAAAATTGCCGCCATTCTCGAAGAGCCGCCa  <  1:3439960/36‑1 (MQ=255)
                                   |                               
TAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCA  >  REL606/166230‑166296

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG