Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,020,505 G→T 11.7% D176Y (GAT→TAT)  ycbU → predicted fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,020,5050G→T11.7% 5.6 ‑4.1 ‑0.043526ycbUpredicted fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (15/23):  total (20/23)
Fisher's exact test strand distribution p-value = 1.61e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.76e-01

AACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACG  >  REL606/1020471‑1020540
                                  |                                   
aaCGGGCGGTAATGCGACGGCGGTTCTCTATTTTGa                                    <  1:7931210/36‑1 (MQ=255)
aaCGGGCGGGAATGCGACGGCGGGTCTCTCTTTTGa                                    >  1:148634/1‑36 (MQ=37)
 aCGGGCGGTAATGCGACGGCGGTTCTCTATTTTGAt                                   <  1:7413541/36‑1 (MQ=255)
 aCGGGCGGTAATGCGACGGCGGTTCTCTATTTTGAt                                   <  1:3954218/36‑1 (MQ=255)
 aCGGGCGGTAATGCGACGGCGGGTTTTTATTTTGAt                                   >  1:294350/1‑36 (MQ=21)
  cGGGCGGTAATGCGACGGTGGTTCTCTATTTTGATc                                  <  1:4810887/36‑1 (MQ=255)
    ggCTGTAATGCGACGGCTGTTCTCTATTTTGATCTg                                <  1:6006930/36‑1 (MQ=37)
     gCGGTAATGCGACGGCGGTTCTCTATTTTGATCTgg                               <  1:2201307/36‑1 (MQ=255)
      cGGTAATGCGACGGCGGTTCTCTATTTTTATCTTGt                              >  1:1023416/1‑36 (MQ=37)
      cGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGt                              <  1:3466627/36‑1 (MQ=255)
       ggTAATGCGACGGCGGTTCTCTATTTTTTTCTTGtt                             >  1:679892/1‑36 (MQ=25)
       ggTAATGCGACGGCGGTTCTCTATTTTTATCTTGtt                             >  1:4814291/1‑36 (MQ=37)
         tAATGCGACGGCGGGTTTCTTTTTTTTTCTTGtttt                           >  1:4992086/1‑35 (MQ=9)
           aTGCGACGGCGGTTTTCTTTTTTTATTTGTTTTTcc                         >  1:7734232/1‑36 (MQ=2)
              cGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGt                      <  1:1514817/36‑1 (MQ=255)
                    cggTTCTCTATTTTGATCTGGTTTACAATTGAAGGa                >  1:1423830/1‑36 (MQ=37)
                      gTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGt              >  1:4886005/1‑36 (MQ=255)
                      gTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGt              >  1:1997017/1‑36 (MQ=255)
                      gTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGt              <  1:2378699/36‑1 (MQ=255)
                       ttCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTg             <  1:3587135/36‑1 (MQ=255)
                       ttCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTg             >  1:3751598/1‑36 (MQ=255)
                       ttCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTg             <  1:566534/36‑1 (MQ=255)
                        tctcTATTTTGATCTGTTTTACCAGTGAAGGAGTGa            >  1:6823331/1‑36 (MQ=255)
                        tctcTATTTTGATCTGGTTTACCAGTGAAGGAGTGa            <  1:6590057/36‑1 (MQ=255)
                        tctcTATTTTGATCTGGTTTACCAGTGAAGGAGTGa            >  1:4828525/1‑36 (MQ=255)
                        tctcTATTTTGATCTGGTTTACCAGTGAAGGAGTGa            <  1:6500992/36‑1 (MQ=255)
                         ctctATTTTGATCTGGTTTACCAGTGAAGGAGTGAt           >  1:6662648/1‑36 (MQ=255)
                          tctATTTTGATCTGGTTTACCAGTGAAGGAGTGATg          >  1:334410/1‑36 (MQ=255)
                           ctATTTTGATCTGGTTTACCAGTGAAGGAGTGATGt         >  1:2599501/1‑36 (MQ=255)
                           ctATTTTGATCTGGTTTACCAGTGAAGGAGTGATGt         <  1:163460/36‑1 (MQ=255)
                           ctATTTTGATCTGGTTTACCAGTGAAGGAGTGATGt         <  1:3075404/36‑1 (MQ=255)
                            tATTTTGATCTGGTTTACCAGTGAAGGAGTGATGtt        <  1:4314143/36‑1 (MQ=255)
                            tATTTTGATCTGGTTTACCAGTGAAGGAGTGATGtt        <  1:1207763/36‑1 (MQ=255)
                              ttttGATCTGGTTTACCAGTGAAGGAGTGATGTTGa      >  1:924066/1‑36 (MQ=255)
                              ttttGATCTGGTTTACCAGTGAAGGAGTGATGTTGa      >  1:4992143/1‑36 (MQ=255)
                              ttttGATCTGGTTTACCAGTGAAGGAGTGATGTTGa      <  1:7853209/36‑1 (MQ=255)
                              ttttGATCTGGTTTACCAGTGAAGGAGTGATGTTGa      >  1:779609/1‑36 (MQ=255)
                              ttttGATCTGGTTTACCAGTGAAGGAGTGATGTTGa      <  1:4509913/36‑1 (MQ=255)
                              ttttGATCTGGTTTACCAGTGAAGGAGTGATGTTGa      <  1:2438054/36‑1 (MQ=255)
                               tttGATCTGGTTTACCAGTGAAGGAGTGATGTTGaa     >  1:1956230/1‑36 (MQ=255)
                               tttGATCTGGTTTACCAGTGAAGGAGTGATGTTGaa     <  1:5490751/36‑1 (MQ=255)
                                ttGATCTGGTTTACCAGTGAAGGAGTGATGTTGaaa    <  1:4407403/36‑1 (MQ=255)
                                ttGATCTGGTTTACCAGTGAAGGAGTGATGTTGaaa    <  1:6580497/36‑1 (MQ=255)
                                  gATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACg  <  1:7311041/36‑1 (MQ=255)
                                  gATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACg  <  1:7611950/36‑1 (MQ=255)
                                  gATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACg  <  1:2609009/36‑1 (MQ=255)
                                  |                                   
AACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACG  >  REL606/1020471‑1020540

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG