Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,737,779 G→C 18.8% A42G (GCA→GGA)  hycI ← protease involved in processing C‑terminal end of HycE

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,737,7790G→C18.8% 3.9 ‑3.8 0.047125hycIprotease involved in processing C‑terminal end of HycE
Reads supporting (aligned to +/- strand):  new base (1/8):  ref base (21/15):  total (22/25)
Fisher's exact test strand distribution p-value = 2.20e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GCAGTTCACGGATGGCGACGATGTCGTTTTCCGGTGCACTACCGCCATCAATCACCACCCAGTTACCTTTC  >  REL606/2737744‑2737814
                                   |                                   
gcAGTTCGCGGATGACGACGATGTCGTTTTCCGATg                                     <  1:423325/36‑1 (MQ=25)
 cAGTTCACGGATGGCGACGATGTCGTTTTCCGGTgg                                    >  1:3215353/1‑35 (MQ=255)
  aGTTCACGGATGGCGACGATGTCGTTTTCCGGTGCa                                   >  1:1820106/1‑36 (MQ=255)
  aGTTCACGGATGGCGACGATGTCGTTTTCCGGTGCa                                   >  1:4239021/1‑36 (MQ=255)
  aGTTCACGGATGGCGACGATGTCGTTTTCCGGTGCa                                   >  1:3173442/1‑36 (MQ=255)
   gTTCACGGATTGCGACGATGTCGTTTTCCGGTGCAc                                  <  1:3730424/36‑1 (MQ=255)
   gTTCACGGATGGCGGCGATGTCGTTTTCCGGTGCAc                                  <  1:2557360/36‑1 (MQ=255)
    ttCACGGATGGCGACGATGTCGTTTTCCGGTGCACt                                 <  1:2075906/36‑1 (MQ=255)
      cACGGATGGCGACGATGTCGTTTTCCGGTGCACTAc                               >  1:1982268/1‑36 (MQ=255)
      cACGGATGGCGACGATGTCGTTTTCCGGTGCACTAc                               <  1:436684/36‑1 (MQ=255)
       gcGGATGGCGACGATGTCGTTTTCCGGTGCACTAcc                              <  1:5424957/35‑1 (MQ=255)
       aCGGATGGCGACGATGTCGTTTTCCGGTGCACTAcc                              <  1:394288/36‑1 (MQ=255)
         ggATGGCGACGATGTCGTTTTCCGGTGCACTACCGc                            <  1:5454565/36‑1 (MQ=255)
         ggATGGCGACGATGTCGTTTTCCGGTGCACTACCGc                            >  1:4898221/1‑36 (MQ=255)
         ggATGGCGACGATGTCGTTTTCCGGTGCACTACCGc                            >  1:5774485/1‑36 (MQ=255)
         ggATGGCGACGATGTCGTTTTCCGGTGCACTACCGc                            <  1:4222546/36‑1 (MQ=255)
         ggATGGCGACGATGTCGTTTTCCGGGGCACTACCGc                            >  1:24920/1‑36 (MQ=255)
         ggATGGCGACGATGTCGTCTTCCCCTCCACCACCGc                            >  1:4162149/1‑36 (MQ=16)
           aTGGCGACGATGTCGTTTTCCGGTGCACTACCGCCa                          >  1:3577578/1‑36 (MQ=255)
           aTGGCGACGATGTCGTTTTCCGGTGCACTACCGCCa                          >  1:2282481/1‑36 (MQ=255)
             ggCGACGATGTCTTTTCCCGGTGCACTACCGCCatc                        <  1:895951/36‑1 (MQ=37)
              gCGACGTTGTCTTTTCCCGGTCCACTACCGCCatca                       <  1:3060359/36‑1 (MQ=16)
              gCGACGATGTCGTTTTCCGGTGCACTACCGCCatca                       >  1:2909786/1‑36 (MQ=255)
              gCGACGATGTCGTTTTCCGGTGCACTACCGCCatca                       >  1:5116184/1‑36 (MQ=255)
                gacgaTGTCGTTTTCCGGTGCACTACCGCCatcaat                     >  1:1279349/1‑36 (MQ=255)
                 acgaTGTCGTTTTCCGGTGCACTACCGCCatcaatc                    >  1:5560648/1‑36 (MQ=255)
                   gaTGTTGTTTACCGGTGCAAAACCGCCATCAATcac                  <  1:949788/36‑1 (MQ=16)
                   gaTGTCTTTTTCCGGTGCACTACCGCCATCAATcac                  <  1:3915294/36‑1 (MQ=255)
                   gaTGTCTTTTTCCGGCACAGCGCCGCCATCAATcac                  <  1:3627077/36‑1 (MQ=2)
                   gaTGTCGTTTTCCGGTGCACTACCGCCATCAATcac                  >  1:1621159/1‑36 (MQ=255)
                   gaTGTCGTTTTCCGGTGCACTACCGCCATCAATcac                  >  1:2577043/1‑36 (MQ=255)
                   gaTGTCGTTTTCCGGTGCACTACCGCCATCAATcac                  >  1:3245071/1‑36 (MQ=255)
                    aTGTCATTTCCGGGTCCCCTCCCGCCATCAATcacc                 <  1:4217774/36‑1 (MQ=9)
                     tGTCGTTTCCCGGTGCACTACCGCCATCAATcacca                <  1:3013186/36‑1 (MQ=255)
                      gTCGTTTTCCGGTGCACTACCGCCATCAATcaccac               >  1:2371354/1‑36 (MQ=255)
                       tCGTTTCCCGGTCCACTACCGCCATCAATCaccacc              <  1:1490363/36‑1 (MQ=37)
                        ctttttCGGGCCCACAACCGCCATCAATCACCAccc             <  1:3756345/34‑1 (MQ=16)
                           tttCCGGTGCACTACCGCCATCAATCACCACCCAGt          >  1:3528608/1‑36 (MQ=255)
                            ttCCGGGCCACTACCGCCATCAATCACCACCCAGtt         <  1:5726002/36‑1 (MQ=37)
                             cccGGTCAACTACCGCCATCAATCACCACCCAGTTa        <  1:300792/35‑1 (MQ=37)
                              ccGGTACACTACCGCCATCAATCACCACCCAGTTAc       <  1:5505359/36‑1 (MQ=255)
                               cgttGCACCACCACCATCAATCACCACCCAGTTAcc      <  1:4395653/33‑1 (MQ=21)
                               cGGTCCACTACCGCCATCAATCACCACCCAGTTAcc      <  1:5044079/36‑1 (MQ=255)
                               cGGTCCACCACCGCCATCAATCACCACCCAGTTAcc      <  1:3222414/36‑1 (MQ=37)
                                ggTGCCCTACCGCCATCAATCACCACCCAGTTACCt     <  1:811058/36‑1 (MQ=255)
                                 gTGCACTACCGCCATCAATCACCACCCAGTTACCtt    >  1:4608375/1‑36 (MQ=255)
                                 gTGCACTACCGCCATCAATCACCACCCAGTTACCtt    >  1:3594038/1‑36 (MQ=255)
                                 gTGCACTACCGCCATCAATCACCACCCAGTTACCtt    <  1:5825367/36‑1 (MQ=255)
                                 gTGCACTACCGCCATCAATCACCACCCAGTTACCtt    >  1:3741025/1‑36 (MQ=255)
                                   gCACTACCGCCATCAATCACCACCCAGTTACCTTTc  <  1:579110/36‑1 (MQ=255)
                                   |                                   
GCAGTTCACGGATGGCGACGATGTCGTTTTCCGGTGCACTACCGCCATCAATCACCACCCAGTTACCTTTC  >  REL606/2737744‑2737814

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 12 ≤  ATCG  < 26 ≤  ATCG  < 37 ≤  ATCG