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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606166,2650C→A28.4% 13.6 ‑4.4 0.0381310hrpBpredicted ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  new base (2/9):  ref base (17/9):  total (20/18)
Fisher's exact test strand distribution p-value = 1.28e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGA  >  REL606/166230‑166300
                                   |                                   
tAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGc                                     <  1:4532830/36‑1 (MQ=255)
tAGCGCGAAGAACGACGACGAAGCTGCTACAGCGGc                                     >  1:1405334/1‑36 (MQ=255)
 aGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCa                                    >  1:7667788/1‑36 (MQ=255)
 aGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCa                                    <  1:1800806/36‑1 (MQ=255)
 aGCGCGAAGAACGACGACGAAGCTGCTACAGCGGCa                                    >  1:7639718/1‑36 (MQ=255)
   cgcgAAGAACGACGACGAAGCTGCTACCGCGGCaaa                                  <  1:5297081/36‑1 (MQ=255)
    gcgAAGAACGACGACGCCGCTGCTACCGCGGTaaaa                                 >  1:694541/1‑36 (MQ=25)
    gcgAAGAACGACGACGAAGCTGCTACCGCGGCaaaa                                 >  1:4687896/1‑36 (MQ=255)
    gcgAAGAACGACGACGAAGCTGCTACCGCGGCaaaa                                 >  1:480786/1‑36 (MQ=255)
    gcgAAGAACGACGACGAAGCTGCTACCGCGGAaaaa                                 >  1:1005013/1‑36 (MQ=255)
     cgAAGAACGACGACGAAGCTGCTACCGCGGCaaaaa                                >  1:7703911/1‑36 (MQ=255)
     cgAAGAACGACGACGAAGCTGCTACCGCGGCaaaaa                                >  1:6040725/1‑36 (MQ=255)
       aagaaCGACGACGAAGCTGCTACCGCGGCAAAAAtt                              >  1:6719919/1‑36 (MQ=255)
        agaaCGAAGACGAAGCTGCTACCGCGGCAAAAATTg                             <  1:2209538/36‑1 (MQ=255)
          aaCGACGACGAAGCTGCTACCGCGGCAAAAATTgcc                           <  1:3766228/36‑1 (MQ=255)
           acgacgacgAAGCTGCTAACGCGGCAAAAATTgccg                          >  1:73638/1‑36 (MQ=255)
            cgacgaggaAGCTGCTACCGCGGCAAAAATTgccgc                         <  1:6010025/36‑1 (MQ=255)
            cgacgacgaAGCTGCTACCGCGGCAAAAATTgccgc                         <  1:1786929/36‑1 (MQ=255)
             gacgacgaAGCTGCTACCGCGGCAAAAATTgccgcc                        >  1:6554132/1‑36 (MQ=255)
             gacgacgaAGCTGCTACCGCGGCAAAAATTgccgcc                        <  1:1150805/36‑1 (MQ=255)
             gacgacgaAGCTGCTACCGCGGAAAAAATTgccgcc                        <  1:4116673/36‑1 (MQ=255)
              acgacgaAGCTGCTACCGCGGCAAAAATTGCCGCCa                       >  1:3574043/1‑36 (MQ=255)
                gacgaAGCTGCTACCGCGGAAAAAATTGCCGCCAtt                     >  1:319408/1‑36 (MQ=255)
                 acgaAGCTGCTACCGCGGCAAAAATTGCCGCCATtc                    <  1:921323/36‑1 (MQ=255)
                  cgaAGCTGCTACCGCGGCAAAAATTGCCGCCATtct                   >  1:2846673/1‑36 (MQ=255)
                  caaaGCTGCTACCGCGGAAAAAATTGCCGCCATtct                   <  1:1052397/34‑1 (MQ=37)
                   gaAGCTGCTACCGCGGCAAAAATTGCCGCCATtctc                  <  1:312639/36‑1 (MQ=255)
                    aaGCTGCAACCGCGGAAAAAATTGCCGCCATTCTCg                 <  1:571452/36‑1 (MQ=25)
                    aaGCTCCTACCGCGGAAAAAATTGCCGCCATTCTCg                 <  1:3714408/36‑1 (MQ=25)
                      gctactACCGCGGAAAAAATTGCCGCCATTCTCGaa               <  1:5515183/36‑1 (MQ=25)
                       ctgctACCGCGGAAAAAATTGCCGCCATTCTCGAag              <  1:138866/36‑1 (MQ=255)
                        gcctACCGCGGAAAAAATTGCCGCCATTCTCGAaga             <  1:7065856/34‑1 (MQ=37)
                          ctACCGCGGCAAAAATTGCCGCCATTCTCGCAGAgc           >  1:3639853/1‑36 (MQ=255)
                           acccGCGGAAAAACTTGCCGCCATTCTCGAAGAgcc          <  1:3691362/34‑1 (MQ=37)
                           aaCCGCGAAAAAAATTGCCGCCATTCTCGAAGAgcc          <  1:6814039/35‑1 (MQ=37)
                            aCCGCGGCAAAAATTGCCGCCATTCTCGAAGAgccg         >  1:1595632/1‑36 (MQ=255)
                              cgcgGCAAAAATTGCCGCCATTCTCGAAGAgccgcc       >  1:6445544/1‑36 (MQ=255)
                               gcgGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCa      >  1:2103055/1‑36 (MQ=255)
                                cgGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCAc     >  1:3632144/1‑36 (MQ=255)
                                 ggCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACg    >  1:6894041/1‑36 (MQ=255)
                                   cAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGa  >  1:2201127/1‑36 (MQ=255)
                                   |                                   
TAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGA  >  REL606/166230‑166300

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 17 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG