Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 218,776 +G 6.2% coding (666/711 nt) cutF → lipoprotein involved with copper homeostasis and adhesion

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606218,7761.→G 6.2% 2.8 ‑0.7 ‑0.032coding (666/711 nt)cutFlipoprotein involved with copper homeostasis and adhesion
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (18/12):  total (20/12)
Fisher's exact test strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.63e-01

GTGCGCCGACTAAAGTATTGGCACCCGATACGGCA.GGGAAATTTTACCCCAACCAGGATTGCAGTAGTTT  >  REL606/218742‑218811
                                   |                                   
gTGCGCCGACTAAAGTATTGGCACCCGATACGGCA.g                                    >  1:2833096/1‑36 (MQ=255)
  gcgcCGACTAAAGTATTGGCACCCGATACGGCA.ggg                                  <  1:255199/36‑1 (MQ=255)
   cgcCGACTAAAGTATTGGCACCCGATACGGCA.GGGa                                 <  1:4667493/36‑1 (MQ=255)
    gcCGACTAAAGTATTGGCACCCGATACGGCA.gggga                                >  1:2183557/1‑34 (MQ=255)
    gcCGACTAAAGTATTGGCACCCGATACGGCA.gggga                                >  1:4342964/1‑34 (MQ=255)
    gcCGACTAAAGTATTGGCACCCGATACGGCA.GGGaa                                >  1:3986524/1‑36 (MQ=255)
      cGACTAAAGTATTGGCACCCGATACGGCA.GGGgaaa                              >  1:3772234/1‑35 (MQ=255)
       gACTAAAGTATTGGCAACCGATACGGCA.GGGcaaat                             >  1:3165752/1‑34 (MQ=37)
        aCTAAAGTATTGGCACCCGATACGGCA.GGGGAAttt                            >  1:1632875/1‑36 (MQ=255)
        aCTAAAGTATTGGCACCCGATACGGCA.GGGGAAttt                            >  1:4046982/1‑36 (MQ=255)
        aCTAAAGTATTGGCACCCGATACGGCA.GGGGAAttt                            >  1:2178286/1‑36 (MQ=255)
        aCTAAAGTATTGGCACCCGATACGGCA.GGGGAAatt                            >  1:4113021/1‑36 (MQ=38)
         cTAAAGTATTGGCACCCGATACGGCA.GGGGAAattt                           >  1:3011468/1‑36 (MQ=38)
          tAAAGTATTGGCACCCGAGACGGCAGGGGAAAtttt                           >  1:3005088/1‑36 (MQ=21)
           aaaGTATTGGCAGCCGATACGGCA.GGGAAATTTTAc                         <  1:5699893/36‑1 (MQ=255)
           aaaGTATTGGCACCCGATACCGCA.GGGGAAtttttc                         >  1:3684657/1‑34 (MQ=37)
            aaGTATTGGCACCCGATACGGCA.GGGGAATTTTatt                        >  1:1408382/1‑34 (MQ=37)
             aGTATTGGCACCCGATACGGCA.GGGAAAATTTAccc                       >  1:1264630/1‑36 (MQ=255)
             aGTATTGGCACCCGATACGGAA.GGGGAATTTTTccc                       >  1:4892215/1‑36 (MQ=25)
              gTATTGGCACCCGATACGGCA.GGGGAATTTTAcccc                      >  1:1442922/1‑36 (MQ=255)
               tATTGGCACCCGATACGGCA.GGGAAATTTTACCCCa                     <  1:1704544/36‑1 (MQ=255)
               tATTGGCACCCGATACGGCA.GGGAAATTTTACCCCa                     <  1:1361200/36‑1 (MQ=255)
                aTTGGCACCCGATACGGCA.GGGAAATTTTAccccca                    >  1:5714967/1‑34 (MQ=255)
                aTTGGCACCCGATACGGCA.GGGAAAATTTACCCCaa                    >  1:5050643/1‑36 (MQ=255)
                 ttGGCACCCGATACGGCA.GGGAAATTTTACCCCAAc                   <  1:2306984/36‑1 (MQ=255)
                 ttGGCACCCGATACGGCA.GGGAAATTTTACCCCAAc                   <  1:4693600/36‑1 (MQ=255)
                  tGGCTCCCGATACGGCA.GGGAAATTTTACCCCAAcc                  <  1:3735381/36‑1 (MQ=255)
                   ggCACCCGATACGGCAGGGGAAATTTTACCCCACcc                  >  1:3028797/1‑36 (MQ=25)
                    gCACCCGATACGGCA.GGGAAATTTTACCCCAACCAg                <  1:4207191/36‑1 (MQ=255)
                        ccGATACGGCA.GGGAAATTTTACCCCAACCAGGAtt            <  1:3414473/36‑1 (MQ=255)
                         cGATACGGCA.GGGGAATTTTACCCCCACCAGGATTg           >  1:3004221/1‑36 (MQ=37)
                          gATACGGCA.GGGGAATTTTACCCCCACCACGATTGc          >  1:1551858/1‑36 (MQ=25)
                             aCGGCA.GGGAAATTTTACCCCAACCAGGATTGCagt       <  1:1150676/36‑1 (MQ=255)
                              cGGCA.GGGAAATTTTACCCCAACCAGGATTGCagta      <  1:2792449/36‑1 (MQ=255)
                               ggCA.GGGAAATTTTACCCCAACCAGGATTGCagtag     >  1:284016/1‑36 (MQ=255)
                               ggCA.GGGAAATTTTACCCCAACCAGGATTGCagtag     <  1:1299258/36‑1 (MQ=255)
                                  a.GGGAAATTTTACCCCAACCAGGATTGCAGTAGttt  >  1:2819710/1‑36 (MQ=255)
                                   |                                   
GTGCGCCGACTAAAGTATTGGCACCCGATACGGCA.GGGAAATTTTACCCCAACCAGGATTGCAGTAGTTT  >  REL606/218742‑218811

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