Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 218,780 A→G 14.4% K224E (AAA→GAA)  cutF → lipoprotein involved with copper homeostasis and adhesion

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606218,7800A→G14.4% 2.1 ‑2.6 0.043670cutFlipoprotein involved with copper homeostasis and adhesion
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (21/15):  total (28/15)
Fisher's exact test strand distribution p-value = 7.65e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CGCCGACTAAAGTATTGGCACCCGATACGGCA.GGGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGG  >  REL606/218745‑218814
                                    |                                  
cgcCGACTAAAGTATTGGCACCCGATACGGCA.GGGa                                    <  1:1874289/36‑1 (MQ=255)
 gcCGACTAAAGTATTGGCACCCGATACGGCA.GGGaa                                   <  1:944490/36‑1 (MQ=255)
   cGACTAAAGTATTGGCACCCGATACGGCA.GGGAAAt                                 <  1:373860/36‑1 (MQ=255)
    gACTAAAGTATTGGCACCCGATACGGCA.GGGgaaat                                >  1:4793780/1‑34 (MQ=38)
     aCTAAAGTATTGGCACCCGATACGGCA.GGGAAAttt                               >  1:711481/1‑36 (MQ=255)
      cTAAAGTATTGTCACCCCTTCCGGCA.GGGAAAtttt                              <  1:6211964/36‑1 (MQ=19)
      cTAAAGTATTGGCACCCGATACGGCA.GGGAAAtttt                              <  1:3060820/36‑1 (MQ=255)
       tAAAGTATTGGCACCCGATACGGCAGGGGAAAtttt                              >  1:5376397/1‑36 (MQ=37)
       tAAAGTATTGGCACCCGATACGGCA.GGGAAATTTTa                             >  1:937759/1‑36 (MQ=255)
        aaaGTATTGGCACCCGATACGGCA.GGGAAATTTTAc                            >  1:7144745/1‑36 (MQ=255)
        aaaGTATTGGCACCCGATACGGCA.GGGAAATTTTAc                            <  1:3948575/36‑1 (MQ=255)
         aaGTATTGGCACCCGATACGGCA.GGGAAATTTTAcc                           <  1:644396/36‑1 (MQ=255)
         aaGTATTGGCACCCGATACGGCA.GGGAAATTTTAcc                           <  1:2361860/36‑1 (MQ=255)
           gTATTGGCACCCGATACGGCA.GGGGAATTTTAcccc                         >  1:2078261/1‑36 (MQ=255)
           gTATTGGCACCCGATACGGCA.GGGAAATTTTAcccc                         >  1:2069066/1‑36 (MQ=255)
           gTATTGGCACCCGATACGGCA.GGGAAATTTTAcccc                         >  1:6388528/1‑36 (MQ=255)
           gTATTGGCACCCGATACGGAA.GGGAAATTTTAcccc                         >  1:573892/1‑36 (MQ=255)
            tATTGGCACCCGATACGGCA.TGTAAATTTTACCCCa                        <  1:1166996/36‑1 (MQ=38)
            tATTGGCACCCGATACGGCA.GGGGNATTTTAccccc                        >  1:3913807/1‑35 (MQ=38)
             aTTGGCACCCGATACGGCA.GGGGAATTTTTCCCCaa                       >  1:2810246/1‑36 (MQ=37)
             aTTGGCACCCGATACGGCA.GGGAAATTTTAccccca                       >  1:3517045/1‑34 (MQ=255)
             aTTGGCACCCGATACGGCA.GGGAAATTTTAccccca                       >  1:2504212/1‑34 (MQ=255)
             aTTGGCACCCGATACGGCA.GGGAAATTTTACCCCaa                       >  1:7302616/1‑36 (MQ=255)
             aTTGGCACCCGATACGGCA.GGCGAATTTTActccca                       >  1:4245390/1‑34 (MQ=18)
              ttGGCACCCGATACGGCA.GGGGAATTTTTCCCCCAc                      >  1:2086333/1‑36 (MQ=25)
              ttGGCACCCGATACGGCA.GGGAAATTTTACCCCAAc                      <  1:434926/36‑1 (MQ=255)
                  cACCCGATACGGCA.GGGAAATTTTACCCCAACCAgg                  >  1:6661431/1‑36 (MQ=255)
                    cccGATACGGCA.GGGGAATTTTACCCCAACCAGGAt                >  1:1657872/1‑36 (MQ=255)
                    cccGATACGGCA.GGGAAATTTTACCCCAACCAGGAt                >  1:3580346/1‑36 (MQ=255)
                    accGATCCGGCA.GGGAAATGTTACCCCAACCAGGAt                >  1:3602229/2‑36 (MQ=37)
                     ccGATACGGCA.GGNGAATTTTACCCCAACCAGGAtt               >  1:3250531/1‑36 (MQ=255)
                     ccGATACGGCA.GGGAAATTTTACCCCAACCAGGAtt               <  1:6612462/36‑1 (MQ=255)
                        aTACGGCA.GGGAAATTTTACCCCAACCAGGATTGCa            >  1:7290665/1‑36 (MQ=255)
                        aTACGGCA.GGGAAATTTTACCCCAACCAGGATTGCa            >  1:3118741/1‑36 (MQ=255)
                         tACGGCA.GGGAAATTTTACCCCAACCAGGATTGCag           <  1:3175570/36‑1 (MQ=255)
                          aCGGCA.GGGAAATTTTACCCCAACCAGGATTGCagt          >  1:6773472/1‑36 (MQ=255)
                           ctgCA.GGGAAATTTTGCCCCCGCCAGGATTGCcgta         <  1:1012508/34‑1 (MQ=16)
                           cGGCA.GGGAAATTTTACCCCACCCAGGATTGCagta         <  1:3757902/36‑1 (MQ=255)
                           cGGCA.GGGAAATTTTACCCCAACCAGGATTGCagta         >  1:6352132/1‑36 (MQ=255)
                             gCA.GGGAAATTTTACCCCAACCAGGATTGCagtagt       <  1:2788818/36‑1 (MQ=255)
                               a.GGGAAATTTTACCCCAACCAGGATTGCAGTAGttt     >  1:1193330/1‑36 (MQ=255)
                                 gggAAATTTTACCCCAACCAGGATTGCAGTAGTTTg    >  1:7165312/1‑36 (MQ=255)
                                 gggAAATTTTACCCCAACCAGGATTGCAGTAGTTTg    >  1:3352026/1‑36 (MQ=255)
                                  ggAAATTTTACCCCAACCAGGATTGCAGTAGTTTgg   <  1:4999704/36‑1 (MQ=255)
                                   gAAATTTTTCCCCAACCAGGATTGCAATAGTTTggg  <  1:128388/36‑1 (MQ=37)
                                   gAAATTTTACCCCAACCAGGATTGCAGTAGTTTggg  >  1:7014216/1‑36 (MQ=255)
                                    |                                  
CGCCGACTAAAGTATTGGCACCCGATACGGCA.GGGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGG  >  REL606/218745‑218814

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