Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,119,132 +G 7.4% coding (117/456 nt) csgA → cryptic curlin major subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,119,1321.→G 7.4% 3.1 ‑0.7 ‑0.227coding (117/456 nt)csgAcryptic curlin major subunit
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (15/10):  total (17/10)
Fisher's exact test strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.19e-01

CGGCGGTAACCACGGTGGTGGCGGTAATAATAGC.GG.CCCAAATTCTGAGCTGAACATTTACCAGTACG  >  REL606/1119099‑1119166
                                  |                                   
cggcggTAACCACGGTGGTGGCGGTAATAATAGC.gg                                   <  1:312173/36‑1 (MQ=255)
  gcggTAACCACGGTGGTGGCGGTAATAATAGC.gggc                                 >  1:1197190/1‑34 (MQ=255)
   cggTAACCACGGTGGTGGCGGGAATAATAGC.GG.ccc                               >  1:2823194/1‑36 (MQ=255)
     gTAACCACGGTGGTGGCGGTAATAATAGC.GG.gccca                             >  1:7933847/1‑34 (MQ=38)
       aaCCACGGTGGTGGCGGTAATAATAGA.GG.GCCCAAt                           >  1:1431114/1‑36 (MQ=25)
        aCCACGGTGGTGGCGGTAATAATAGC.GG.CCCAAAtt                          <  1:4503182/36‑1 (MQ=255)
         ccACGGTGGTGGCGGTAATAATAGC.GG.CCCCAATTc                         >  1:58611/1‑36 (MQ=255)
          cACGGTGGTGGCGGTGATAATAGT.GGCCCCAAAttt                         >  1:5338334/1‑35 (MQ=16)
          cACGGTGGTGGCGGGAATAATATC.GG.CCCCAAtttt                        >  1:6373282/1‑34 (MQ=21)
          cACGGTGGTGGCGGAAATAATAGCGGG.CCCAAAtct                         >  1:8231317/1‑34 (MQ=21)
           aCGGTGGTGGCGGTAATAATAGC.GG.CCCAAATTCTg                       <  1:3597947/36‑1 (MQ=255)
           aCGGTGGTGGCGGTAATAATAACGGG.CCCAAAtttt                        >  1:3715826/1‑34 (MQ=21)
           aCGGTGGTGGCGGATATAATAGC.GG.CCCGCATTcgt                       >  1:201586/1‑34 (MQ=16)
             ggtggtGGCGGTAATAATAGC.GG.CCCAAATTCTGAg                     <  1:6362414/36‑1 (MQ=255)
               tggtggCGGTAATAATAGC.GG.CCCAAATTTTGAGCt                   >  1:8291967/1‑36 (MQ=255)
                ggtggCGGTAATAATAGC.GG.CCCCAATTCTGAGCTg                  >  1:1269216/1‑36 (MQ=255)
                  gggCGGTAATAATAGC.GG.CCCAAATTCTGAGCTGaa                <  1:4615320/35‑1 (MQ=255)
                    gCGGTAATAATAGC.GG.CCCAAATTCTGAGCTGAACa              >  1:4047429/1‑36 (MQ=255)
                    gCGGTAATAATAGC.GG.CCCAAATTCTGAGCTGAACa              <  1:2581573/36‑1 (MQ=255)
                     cGGTAATAATAGC.GG.CCCAAATTCTGAGCTGAACAt             >  1:1311964/1‑36 (MQ=255)
                      ggTAATAATAGC.GG.CCCAAATTCTGAGCTGAACAtt            <  1:6043328/36‑1 (MQ=255)
                      ggTAATAATAGC.GG.CCCAAATTCTGAGCTGAACAtt            <  1:6225061/36‑1 (MQ=255)
                        taataaTAGC.GG.CCCAAATTCTGAGCTGAACATTTa          <  1:5058543/36‑1 (MQ=255)
                        taataaTAGC.GG.CCCAAATTCTGAGCTGAACATTTa          <  1:4353484/36‑1 (MQ=255)
                          ataataGC.GG.CCCCAATTCTGAGCTGAACATTTAcc        >  1:7837103/1‑36 (MQ=255)
                            aataGC.GG.CCCAAATTCTGAGCTGAACATTTACCAg      >  1:3717395/1‑36 (MQ=255)
                              taGC.GG.CCCAAATTCTGAGCTGAACATTTACCAGTa    >  1:3266733/1‑36 (MQ=255)
                               aGC.GG.CCCAAATTCTGAGCTGAACATTTACCAGTAc   <  1:7418774/36‑1 (MQ=255)
                                gC.GG.CCCAAATTCTGAGCTGAACATTTACCAGTACg  >  1:6069622/1‑36 (MQ=255)
                                  |                                   
CGGCGGTAACCACGGTGGTGGCGGTAATAATAGC.GG.CCCAAATTCTGAGCTGAACATTTACCAGTACG  >  REL606/1119099‑1119166

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 13 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG