Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,454,374 C→A 26.1% F642L (TTC→TTA hrpA → ATP‑dependent helicase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,454,3740C→A26.1% 11.4 ‑4.3 0.0231926hrpAATP‑dependent helicase
Reads supporting (aligned to +/- strand):  new base (0/6):  ref base (11/6):  total (11/12)
Fisher's exact test strand distribution p-value = 1.37e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCGAAATGGGTAATGGTGGCGGAA  >  REL606/1454339‑1454407
                                   |                                 
gcgcgTTTCTCCATCTTCCCCGGTTCTGGTTTATTc                                   <  1:1773284/36‑1 (MQ=255)
gcgcaTTTCTCCATCTTCCCCGGTTCTGGTTTATTc                                   <  1:5919334/36‑1 (MQ=255)
 cgcgTTTCTCCATCTTCCCCGGTTCTGGTTTATTCa                                  <  1:1484080/36‑1 (MQ=255)
 cgcgTTTCTCCATCTTCCCCGGTTCTGGTTTATTCa                                  <  1:2788535/36‑1 (MQ=255)
 cgcgTTTCTCAATCTTCCCCGGTTCTGGTTTATTCa                                  <  1:5368725/36‑1 (MQ=255)
       tctcCATCTTCCCCGGTTCTGGTTTATTCAAAAAAc                            <  1:7801862/36‑1 (MQ=255)
       tctcCATCTTCCCCGGTTCTGGTTTATTCAAAAAAc                            <  1:5988177/36‑1 (MQ=255)
        ctcCATCTTCCCCGGTTCTGGTTTATTCAAAAAAcc                           >  1:4896063/1‑36 (MQ=255)
          ccATCTTCCCCGGTTTTGGTTTATTCAAAAAAccgc                         >  1:2995663/1‑36 (MQ=255)
          ccATCTTCCCCGGTTCTGGTTTATTCAAAAAAccgc                         <  1:7731548/36‑1 (MQ=255)
               ttCCCCGGTTCTGGTTTATTCAAAAAACCGCCGaaa                    >  1:4704242/1‑36 (MQ=255)
                   cgggTTCTGGTTTATAAAAAAAACCGCCGAAATggg                <  1:7309133/34‑1 (MQ=25)
                   ccGGTTTTGTTTTATTAAAAAACCCGCCGAAATggg                <  1:7007065/36‑1 (MQ=16)
                   ccGGTTTTGTTTTATAAAAAAAACCGCCGAAATggg                <  1:6681818/36‑1 (MQ=16)
                    cgattCTGCTTTATAAAAAAAGCCGCCGAAATGGGt               <  1:292401/33‑1 (MQ=12)
                      tttCTGGTATATAAAAAAACCCGCCGAAATGGGTaa             <  1:641764/35‑1 (MQ=9)
                           gttttATAAAAAAAACCGCCGAAATGGGTAAtggtg        <  1:7165254/34‑1 (MQ=25)
                           ggTTTATTCAAAAAACCGCCGTAATGGGGAAtggtg        >  1:706464/1‑36 (MQ=37)
                           ggTTTATTCAAAAAACCGCCGAAATGGGTAAtggtg        >  1:4196363/1‑36 (MQ=255)
                           ggTTTATTCAAAAAACCGCCGAAATGGGTAAtggtg        >  1:3048656/1‑36 (MQ=255)
                            gTTTATTCAAAAAACCGCCGAGATGGGAAAtggtgg       >  1:4263177/1‑36 (MQ=37)
                               tATTCAAAAAACCGCCGAAATGGGTAATGGGGGCgg    >  1:2278096/1‑36 (MQ=255)
                                aTTCAAAAAACCGCCGAAATGGGTAATGGTGGCGGa   >  1:86008/1‑36 (MQ=255)
                                aTTCAAAAAACCGCCGAAATGGGTAATGGTGGCGGa   >  1:3752905/1‑36 (MQ=255)
                                 ttCAAAAAACCGCCGAAATGGGGAATGGTGTCGGaa  >  1:2672996/1‑36 (MQ=25)
                                   |                                 
GCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCGAAATGGGTAATGGTGGCGGAA  >  REL606/1454339‑1454407

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