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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,454,3740C→A13.0% 2.5 ‑4.5 ‑0.0231926hrpAATP‑dependent helicase
Reads supporting (aligned to +/- strand):  new base (0/3):  ref base (17/3):  total (17/6)
Fisher's exact test strand distribution p-value = 1.13e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.49e-01

GCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCGAAATGGGTAATGGTGGCGGAACT  >  REL606/1454339‑1454409
                                   |                                   
gcgcgTTTCTCCATCTTCCCCGGTTCTGGTTTATTc                                     >  1:5685322/1‑36 (MQ=255)
gcgcgTTTCTCCATCTTCCCCGGTTCTGGTTTATTc                                     <  1:425959/36‑1 (MQ=255)
  gcgTTTCTCCATCTTCCCCGGTTCTGGTTTATTCaa                                   >  1:7254076/1‑36 (MQ=255)
     tttCTCCATCTTCCCCGGTTCTGGTTTATTCaaaaa                                >  1:1164602/1‑36 (MQ=255)
     tttCTCCATCTTCCCCGGTTCTGGTTTATTCaaaaa                                >  1:6961683/1‑36 (MQ=255)
       tctcCATCTTCCCCGGTTCTGGTTTATTCAAAAAAc                              >  1:3680787/1‑36 (MQ=255)
        ctcCATCTTCCCCGGTTCTGGTTTATTCAAAAAAcc                             <  1:5277791/36‑1 (MQ=255)
        ctcCATCTTCCCCGGTTCTGGTTTATTCAAAAAAcc                             <  1:4388429/36‑1 (MQ=255)
          ccATCTTCCCCGGTTCTGGTTTATTCAAAAAAccgc                           <  1:5520892/36‑1 (MQ=255)
            aTCTTCCCCGGTTCTGGTTTATTCAAAAAAccgccg                         >  1:6037/1‑36 (MQ=255)
             tCTTCCCCGTTTCTGTTTTATTAAAAAAACCGCCGa                        <  1:1031290/36‑1 (MQ=25)
              cTTCCCCGGTTCTGGTTTATTCAAAAAAccgccaga                       >  1:3453354/1‑33 (MQ=255)
                  cccGGTTCTGGTTTATTCAAAAAACCGCCGAAATgg                   >  1:4011540/1‑36 (MQ=255)
                  cccGGTTCTGGTTTATTCAAAAAACCGCCGAAATgg                   >  1:3316956/1‑36 (MQ=255)
                   cgggTTCTGGTTTATAAAAAAAACCGCCGAAATggg                  <  1:4022442/34‑1 (MQ=25)
                    gggTTTTGTTTTATAAAAAAAACCGCCGAAATGGGt                 <  1:4124000/35‑1 (MQ=16)
                        tCTGGTTTATTCAAAAAACCGCCGAAATGGGTAAtg             >  1:4656519/1‑36 (MQ=255)
                           ggTTTATTCAAAAAACCGCCGAAATGGGTAAtggtg          >  1:2233108/1‑36 (MQ=255)
                            gTTTATTCAAAAAACCGCCGAAATGGGTAAtggtgg         >  1:2662570/1‑36 (MQ=255)
                             tttATTCAAAAAACCGCCGAAATGGGTAATGGTGGc        >  1:1584267/1‑36 (MQ=255)
                              ttATTCAAAAAACCGCCGAAATGGGTAATGGTGGCg       >  1:3923305/1‑36 (MQ=255)
                               tATTCAAAAAACCGCCGAAATGGGTAATGGTGGCgg      >  1:4896418/1‑36 (MQ=255)
                                aTTCAAAAAACCGCCGAAATGGGTAATGGTGGCGGa     >  1:7238602/1‑36 (MQ=255)
                                 ttCAAAAAACCGCCGAAATGGGTAATGGTGGCGGaa    >  1:4097550/1‑36 (MQ=255)
                                 ttCAAAAAACCGCCGAAATGGGTAATGGTGGCGGaa    >  1:3166344/1‑36 (MQ=255)
                                   cAAAAAACCGCCGAAATGGGTAATGGTGGCGGAACt  >  1:3856622/1‑36 (MQ=255)
                                   |                                   
GCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCGAAATGGGTAATGGTGGCGGAACT  >  REL606/1454339‑1454409

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