Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,019,082 T→G 18.1% V62G (GTA→GGA)  ycbT → predicted fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,019,0820T→G18.1% 3.0 ‑3.4 0.033185ycbTpredicted fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (14/12):  total (21/12)
Fisher's exact test strand distribution p-value = 3.19e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGA  >  REL606/1019047‑1019117
                                   |                                   
tAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGt                                     <  1:5698095/36‑1 (MQ=255)
  aaaCCCTGGAGGAAGAAGCAAATACGGCAGGGGGAg                                   >  1:4877257/1‑36 (MQ=255)
    aCCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGtt                                 <  1:5626310/36‑1 (MQ=255)
        tGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAg                             >  1:6812811/1‑36 (MQ=255)
        tGGAGGAAGAAGCAAATACGGCAGGGGGAGCTTAAg                             >  1:340223/1‑36 (MQ=25)
            ggaAGAAGCAAATACGGCAGGGGTAGTTAAAGACaa                         <  1:549576/36‑1 (MQ=255)
             gaagaaGCAAATACGGCAGGGGTAGTTAAAGACaaa                        <  1:623920/36‑1 (MQ=255)
              aagaagCAAATACGGCAGGGGTAGTTAAAGACAAAt                       <  1:4961005/36‑1 (MQ=255)
              aagaagCAAATACGGCAGGGGTAGTTAAAGACAAAt                       >  1:291949/1‑36 (MQ=255)
              aagaagCAAATACGGCAGGGGTAGTGAAAGACAAAt                       >  1:374261/1‑36 (MQ=255)
              aagaagCAAATACGGCAGGGGGAGTTAAAAACaaaa                       >  1:1116798/1‑35 (MQ=37)
                gaagAAAATACGGCAGGGGTAGTTAAAGACAAAttt                     <  1:3011971/36‑1 (MQ=255)
                 aatcAAATACGGCAGGGGTAGTTAAAGACAAAtttt                    <  1:7027804/33‑1 (MQ=255)
                 aagCAAATACGGCAGGGGTAGTTAAAGACAAAtttt                    >  1:1395/1‑36 (MQ=255)
                 aagCAAATACGGCAGGGGGAGGTAAAAAAAAAattt                    >  1:5777489/1‑36 (MQ=16)
                    cAAATTCGTCAGGGGTAGTTAAAGACAAATTTTATc                 <  1:6467296/36‑1 (MQ=37)
                    cAAATACGGCAGGGGGGGTTTAAGAGAAAatttttt                 >  1:2603607/1‑33 (MQ=255)
                    cAAATACGGCAGGGGGACTTTAAGACAAATTTTTTc                 >  1:6177539/1‑36 (MQ=255)
                      aaTACGGCAGGGGTAGTTAAAGACAAATTTTATCAg               >  1:1114447/1‑36 (MQ=255)
                        tACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTg             <  1:6745046/36‑1 (MQ=255)
                           ggCAGGGGTAGTTAAAGAAAAATTTTATCAGTGGaa          >  1:4236797/1‑36 (MQ=255)
                            gCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAc         >  1:978642/1‑36 (MQ=255)
                            gCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAc         >  1:7113529/1‑36 (MQ=255)
                            gCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAc         >  1:2580317/1‑36 (MQ=255)
                            gCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAc         >  1:3644836/1‑36 (MQ=255)
                            gCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAc         >  1:4558247/1‑36 (MQ=255)
                              aGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcga       >  1:2931343/1‑36 (MQ=255)
                              aGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcga       <  1:6588957/36‑1 (MQ=255)
                              aGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcga       >  1:7265365/1‑36 (MQ=255)
                               ggggTAGTTAAAGACAAATTTTATCAGTGGAAcgaa      >  1:4052915/1‑36 (MQ=255)
                                gggTAGTTAAAGACAAATTTTATCAGTGGAAcgaat     >  1:2078010/1‑36 (MQ=255)
                                 ggTAGTTAAAGACAAATTTTATCAGTGGAAcgaatc    <  1:5824494/36‑1 (MQ=255)
                                 ggTAGTTAAAGACAAATTTTATCAGTGGAAcgaatc    <  1:1768790/36‑1 (MQ=255)
                                 ggTAGTTAAAGACAAATTTTATCAGTGGAAcgaatc    <  1:4029200/36‑1 (MQ=255)
                                   tAGTTAAAGACAAATTTTATCTGTGGAAcgaatcga  <  1:6332456/36‑1 (MQ=255)
                                   |                                   
TAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGA  >  REL606/1019047‑1019117

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 22 ≤  ATCG  < 31 ≤  ATCG  < 37 ≤  ATCG