Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,021,723 C→G 13.3% A229G (GCA→GGA)  ycbF → predicted periplasmic pilini chaperone

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,021,7230C→G13.3% 2.5 ‑3.8 0.042686ycbFpredicted periplasmic pilini chaperone
Reads supporting (aligned to +/- strand):  new base (6/0):  ref base (16/19):  total (23/19)
Fisher's exact test strand distribution p-value = 2.26e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCG  >  REL606/1021688‑1021758
                                   |                                   
cTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGc                                     <  1:5027219/36‑1 (MQ=255)
cTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGc                                     <  1:1622631/36‑1 (MQ=255)
 tGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCa                                    <  1:2162231/36‑1 (MQ=255)
 tGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCa                                    <  1:715402/36‑1 (MQ=255)
    cgggcgATAAATGATTACGGCGGGTTGAGTGCAaag                                 <  1:1050713/36‑1 (MQ=255)
     gggcgATAAATGATTACGGCGGGTTGAGTGCAaaga                                <  1:3187398/36‑1 (MQ=255)
     gggcgATAAATGATTACGGCGGGGTGTGTGGaaaaa                                >  1:6420302/1‑34 (MQ=21)
       gcgATAAATGATTACGGCGGGTTTTGTGTAaaaaag                              >  1:7756739/1‑36 (MQ=19)
       gcgATAAATGATTACGGCGGGGTGGGTGGGaagaag                              >  1:4015709/1‑36 (MQ=19)
         gATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGa                            <  1:5740156/36‑1 (MQ=255)
          aTAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAg                           <  1:1909557/36‑1 (MQ=255)
          aTAAATGATTACGGCGGGTGGATTGCAAAGAAGGAg                           >  1:8166237/1‑36 (MQ=25)
          aTAAATGATTACGGCGGGGTTTGTGGAAAGAAGGAg                           >  1:3485202/1‑36 (MQ=19)
            aaaTGATTACGGCGGGTTTTGTGGAAAGAAAGAgga                         >  1:1747977/1‑34 (MQ=16)
             aaTGATTACGGCGGGGTTACTGCCTAGAAGGAGCaa                        >  1:7920191/1‑36 (MQ=16)
              aTGATTACGGCGGGTTTAGTGCCAAGAAGGAGCaaa                       >  1:8070894/1‑36 (MQ=38)
              aTGATTACGGCGGGTTTACTGGCTAGAAGGAGCaga                       >  1:424975/1‑34 (MQ=9)
               tGATTACGGCGGGTTGAGTGCAAAGAAGGAGCaaaa                      <  1:3098220/36‑1 (MQ=255)
               tGATTACGGCGGGGTGAGTGCAAAGAAGGAGGaaaa                      >  1:2402087/1‑36 (MQ=37)
                gATTACGGCGGGGTGGGTGGAAAAAAAGAGGaaaaa                     >  1:4600597/1‑35 (MQ=2)
                   tACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTg                  <  1:2596683/36‑1 (MQ=255)
                    aCGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGc                 <  1:1349981/36‑1 (MQ=255)
                     cggcggGTTGAGTGCAAAGAAGGAGCAAAATCTGcc                >  1:1481928/1‑36 (MQ=255)
                       gcggGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCt              <  1:6468501/36‑1 (MQ=255)
                        cggGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTg             >  1:1903497/1‑36 (MQ=255)
                         gggTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGa            >  1:7530408/1‑36 (MQ=255)
                         gggTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGa            >  1:6318364/1‑36 (MQ=255)
                          ggTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGaa           >  1:4396008/1‑36 (MQ=255)
                           gTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGaaa          <  1:1367880/36‑1 (MQ=255)
                           gTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGaaa          >  1:5030876/1‑36 (MQ=255)
                            ttGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAAc         <  1:848048/36‑1 (MQ=255)
                             tGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACa        >  1:7950866/1‑36 (MQ=255)
                             tGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACa        <  1:6058906/36‑1 (MQ=255)
                             tGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACa        <  1:201046/36‑1 (MQ=255)
                             tGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACa        <  1:2480791/36‑1 (MQ=255)
                              gAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAg       <  1:6197208/36‑1 (MQ=255)
                               aGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAgg      >  1:2444839/1‑36 (MQ=255)
                                 ttcAAATAAGGATCAAAATATGCCCTGAAACAGGtt    <  1:2446766/34‑1 (MQ=21)
                                 tGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGtt    >  1:2319197/1‑36 (MQ=255)
                                 tGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGtt    <  1:7742169/36‑1 (MQ=255)
                                 tGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGtt    <  1:502268/36‑1 (MQ=255)
                                  gCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTc   <  1:6158213/36‑1 (MQ=255)
                                  gCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTc   >  1:7179603/1‑36 (MQ=255)
                                  gCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTc   >  1:8087209/1‑36 (MQ=255)
                                  gCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTc   >  1:1090703/1‑36 (MQ=255)
                                   cAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCg  <  1:1898633/36‑1 (MQ=255)
                                   cAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCg  <  1:237170/36‑1 (MQ=255)
                                   |                                   
CTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCG  >  REL606/1021688‑1021758

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