Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 480,987 A→G 8.3% Y36H (TAC→CAC)  ycgB ← type IV piliation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947480,9870AG 8.3% 131.6 / 2.8 48Y36H (TAC→CAC) ycgBtype IV piliation protein
Reads supporting (aligned to +/- strand):  ref base A (20/24);  new base G (2/2);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CGCGGCTGTAGGACTTTTCCGTACTGAGGAACTGCTTGCCGTAGGACCAGTGGTGGTAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCGTCTCTGTGCTCTGGCGGTCATGTGGCTAACCTGCGCTGGAAGAGTTCACGGAAGACCGGGTA  >  NC_002947/480840‑481133
                                                                                                                                                   |                                                                                                                                                  
cgcgGCTGTAGGACTTTTCCGTACTGAGGAACTGCTTGCCGTAGGACCAGTGGTGGTAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTaa                                                                                                                                                   <  8:68976/149‑1 (MQ=255)
  cgGCTGTAGGACTTTTGCGTACTGAGTTACTGCTTGCCGTAGGACCAGTGGTGGTAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCGAGGGCGTAACg                                                                                                                                                 <  4:59627/149‑1 (MQ=255)
    gCTGTAGGACTTTTCCGTACTGAGGAACTGCTTGCCGTAGGACCAGTGGTGGTAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCt                                                                                                                                               <  1:10150/149‑1 (MQ=255)
       gTAGGACTTTTCCGTACTGAGGAACTGCTTGCCGTAGGACCAGTGGTGGTAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCgg                                                                                                                                            >  6:28630/1‑149 (MQ=255)
        tAGGACTTTTCCGTACTGAGGAACTGCTTGCCGTAGGACCAGTGGTGGTAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACGTCGATCTGGTTGGGGTAGGTGTACAGGGCGTAACGCTCGGc                                                                                                                                           >  1:135846/1‑149 (MQ=255)
        tAGGACTTTTCCGTACTGAGGAACTGCTTGCCGTAGGACCAGTGGTGGTAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGc                                                                                                                                           >  1:130470/1‑149 (MQ=255)
              ttttCCGTACTGAGGAACTGCTTGCCGTAGGACCAGTGGTGGTAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAg                                                                                                                                       >  1:183546/1‑147 (MQ=255)
              ttttCCGTACTGAGGAACTGCTTGCCGTAGGACCAGTGGTGGTAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAg                                                                                                                                       <  2:183546/147‑1 (MQ=255)
                 tCCGTACTGAGGAACTGCTTGCCGTAGGACCAGTGGTGGTAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGc                                                                                                                                   <  6:11321/148‑1 (MQ=255)
                   cGTACTGAGGAACTGCTTGCCGTAGGACCAGTGGTGGTAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGa                                                                                                                                <  5:28630/149‑1 (MQ=255)
                       cTGAGGAACTGCTTGCCGTAGGACCAGTGGTGGTAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGTAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTc                                                                                                                            <  1:290088/149‑1 (MQ=255)
                                     gCCGTAGGACCAGTGGTGGTAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTc                                                                                                                                              <  8:122060/117‑1 (MQ=255)
                                     gCCGTAGGACCAGTGGTGGTAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTc                                                                                                                                              >  7:122060/1‑117 (MQ=255)
                                      ccGTAGGACCAGTGGTGGTAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCt                                                                                                              <  1:119205/148‑1 (MQ=255)
                                       cGTAGGACCAGTGGTGGTAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCgg                                                                                                                                    <  7:114426/125‑1 (MQ=255)
                                       cGTAGGACCAGTGGTGGTAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCgg                                                                                                                                    >  8:114426/1‑125 (MQ=255)
                                         acggcccAGGGGTGGTAGCCCAGCGGCAGGCCGACCGCGGCGTACGCGTCCATCACCTGCTCCGCAGTGATGACCTCGAGCGGGTTGGGGGAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGCAGCTCTGGAt                                                                                                          <  2:173975/144‑1 (MQ=255)
                                             aCCAGTGGTGGTAGCCCAGCGGCATGCCGACGGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGc                                                                                                      <  3:206335/149‑1 (MQ=255)
                                                      ggtAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCgg                                                                                                                                    <  4:178062/110‑1 (MQ=255)
                                                      ggtAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGt                                                                                              >  3:234746/1‑148 (MQ=255)
                                                       gtAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTGACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTcc                                                                                            >  4:91618/1‑149 (MQ=255)
                                                       gtAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTGACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTc                                                                                             >  4:96566/1‑148 (MQ=255)
                                                          gCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCAc                                                                                          <  4:234746/148‑1 (MQ=255)
                                                             cAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCgg                                                                                      >  5:238371/1‑149 (MQ=255)
                                                                      gCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTg                                                                              >  2:111337/1‑148 (MQ=255)
                                                                        cgaccgaGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGa                                                                            >  5:353117/1‑148 (MQ=255)
                                                                                   taggcgtCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTGACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCg                                                                <  3:96566/149‑1 (MQ=255)
                                                                                   taggcgtCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTGACGCGCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCg                                                                <  3:91618/149‑1 (MQ=255)
                                                                                     ggcgtCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCGt                                                               <  6:353117/148‑1 (MQ=255)
                                                                                        gtCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCGTctct                                                           >  3:96156/1‑149 (MQ=255)
                                                                                            atcatcTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAg                                                                                >  1:295818/1‑124 (MQ=255)
                                                                                            atcatcTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAg                                                                                <  2:295818/124‑1 (MQ=255)
                                                                                            atcatcTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCGTCTCTGTGc                                                       >  3:29218/1‑149 (MQ=255)
                                                                                               atcTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGTGTTGGCGTCTCTGTGctct                                                    <  4:29218/149‑1 (MQ=255)
                                                                                                tcTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCg                                                                                                                     <  8:172726/83‑1 (MQ=255)
                                                                                                tcTGTTCGGCAGTGATGACCTCGATCTGCTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGACAAGCGGCTGATTTCCCGGTCg                                                                                                                     >  7:172726/1‑83 (MQ=255)
                                                                                                                gaCCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCGTCTCTGTGCTCTGGCGGTCATGTGGCTaa                                   >  8:302228/1‑149 (MQ=255)
                                                                                                                   cTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCGTCTCTGTGCTCTGGCGGTCATGTGGCTAACCt                                >  5:270260/1‑149 (MQ=255)
                                                                                                                   cTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCGTCTCTGTGCTCTGGCGGTCATGTGGCTAACCt                                >  5:231547/1‑149 (MQ=255)
                                                                                                                            gTTTGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCGTCTCTGTGCTCTGGCGGTCATGTGGCTAACCTGCGCTGGaa                       <  2:135846/149‑1 (MQ=255)
                                                                                                                            gTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCGTCTCTGTGCTCTGGCGGTCATGTGGCTAACCTgcgc                            <  4:96156/144‑1 (MQ=255)
                                                                                                                            gTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCGTCTCTGTGCTCTGGCGGTCATGTGGCTAACCTGCGCTGGaa                       <  2:130470/149‑1 (MQ=255)
                                                                                                                                gggTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCGTCTCTGTGCTCTGGCGGTCATGTGGCTAACCTGCGCTGGAAGAGt                   >  4:256410/1‑149 (MQ=255)
                                                                                                                                 ggTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATgg                                                                       >  4:82465/1‑96 (MQ=255)
                                                                                                                                 ggTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATgg                                                                       <  3:82465/96‑1 (MQ=255)
                                                                                                                                     ggTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCGTCTCTGTGCTCTGGCGGTCATGTGGCTAACCTGCGCTGGAAGAGTTCACg              <  7:302228/149‑1 (MQ=255)
                                                                                                                                            aGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGc                                                                                                      <  6:284820/54‑1 (MQ=255)
                                                                                                                                            aGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGc                                                                                                      >  5:284820/1‑54 (MQ=255)
                                                                                                                                                 gTAACGCTCGGCCAAGCGGCTGATGTCGCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCGTCTCTGTGCTCTGGCGGTCATGTGGCTAACCTGCGCTGGAAGAGTTCACGGAAGACCGGGTa  <  8:16420/149‑1 (MQ=255)
                                                                                                                                                   aaCGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCGTCTCTGTGCTCTGGCGGTCATGTGGCTAACCTGCGCTg                          <  6:52770/123‑1 (MQ=255)
                                                                                                                                                   aaCGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCGTCTCTGTGCTCTGGCGGTCATGTGGCTAACCTGCGCTg                          >  5:52770/1‑123 (MQ=255)
                                                                                                                                                   aaCGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCGTCTCTGTGCTCTGGCGGTCATGTGGCTAACCTGCGCTGGAAGAGtt                  <  7:252616/131‑1 (MQ=255)
                                                                                                                                                   aaCGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCGTCTCTGTGCTCTGGCGGTCATGTGGCTAACCTGCGCTGGAAGAGtt                  >  8:252616/1‑131 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCGGCTGTAGGACTTTTCCGTACTGAGGAACTGCTTGCCGTAGGACCAGTGGTGGTAGCCCAGCGGCATGCCGACCGAGGCGTAGGCGTCCATCATCTGTTCGGCAGTGATGACCTCGATCTGGTTGGGGTAGGTGTCCAGGGCGTAACGCTCGGCCAAGCGGCTGATTTCCCGGTCGTAGGTCTGGATCAGCTCGAACGTCCACTCGGACCCAGTGGAAATGGGTTGGCGTCTCTGTGCTCTGGCGGTCATGTGGCTAACCTGCGCTGGAAGAGTTCACGGAAGACCGGGTA  >  NC_002947/480840‑481133

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑