Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 739,195 G→C 14.9% A75G (GCG→GGG)  pilD ← leader peptidase/N‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947739,1950GC14.9% 126.6 / 15.6 47A75G (GCG→GGG) pilDleader peptidase/N‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (21/19);  new base C (4/3);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

GCCACTACCAGCGACAGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGA  >  NC_002947/739054‑739333
                                                                                                                                             |                                                                                                                                          
gCCACTACCAGCGACAGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAgtgt                                                                                                                                    >  1:1027958/1‑150 (MQ=255)
 ccACTACCAGCGACAGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTACTTACACGCACGAATCCGGTGCGCGCAGtgtg                                                                                                                                   <  1:635335/150‑1 (MQ=255)
               agcagcGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCg                                                                                                                     <  1:292054/150‑1 (MQ=255)
               agcagcGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCgcaccggattcgt                                                                                                                     >  2:708520‑M2/1‑137 (MQ=255)
                   gcgcgcTGGCCACGTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGGCCGGCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAACCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGggca                                                                                                                 <  2:37808/150‑1 (MQ=255)
                         tGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGTACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATTACCGGAATGTTCTCCAACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGcaca                                                                                                           >  2:449151/1‑150 (MQ=255)
                         tGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGcaca                                                                                                           >  2:77123/1‑150 (MQ=255)
                            ccACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCgcaccggattcgtgcgt                                                                                                                 >  2:1099569‑M2/1‑124 (MQ=255)
                            ccACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCgcaccggattcgtgcgt                                                                                                                 <  1:1099569‑M2/141‑18 (MQ=255)
                             cACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCgg                                                                                                                    >  2:1083252/1‑137 (MQ=255)
                             cACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCgg                                                                                                                    <  1:1083252/137‑1 (MQ=255)
                                   cACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAcg                                                                                                 >  2:935635/1‑150 (MQ=255)
                                        ccGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGt                                                                                            >  1:1158728/1‑150 (MQ=255)
                                                       tGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAc                                                                                                  >  2:838865/1‑129 (MQ=255)
                                                       tGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAc                                                                                                  <  1:838865/129‑1 (MQ=255)
                                                          tGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGaa                                                                                                    >  2:426206/1‑124 (MQ=255)
                                                          tGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGaa                                                                                                    <  1:426206/124‑1 (MQ=255)
                                                            cGGTTTTTGCAGGATGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCAGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGcccc                                                                        <  2:372080/150‑1 (MQ=255)
                                                                       ggACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAcgcgc                                                                                              <  1:841000/117‑1 (MQ=255)
                                                                       ggACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAcgcgc                                                                                              >  2:841000/1‑117 (MQ=255)
                                                                        gACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGcc                                                                          <  1:946872/136‑1 (MQ=255)
                                                                        gACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGcc                                                                          >  2:946872/1‑136 (MQ=255)
                                                                        gACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTgc                                                            >  2:1184642/1‑150 (MQ=255)
                                                                        gACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCgcaccggattcgt                                                                                                                     <  1:208935‑M2/93‑14 (MQ=255)
                                                                        gACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCgcaccggattcgt                                                                                                                     >  2:208935‑M2/1‑80 (MQ=255)
                                                                           gAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGcc                                                         <  2:517556/150‑1 (MQ=255)
                                                                             aCAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGAAGTGTGGGCACTGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTc                                                       <  2:887059/150‑1 (MQ=255)
                                                                                gcgcCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACg                                                    >  2:729659/1‑150 (MQ=255)
                                                                                   cccgcccAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTg                                                 >  1:60964/1‑150 (MQ=255)
                                                                                        ccAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAcgcgc                                                                                              <  1:18973/100‑1 (MQ=255)
                                                                                        ccAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAcgcgc                                                                                              >  2:18973/1‑100 (MQ=255)
                                                                                        ccAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCgcaccggattcgtgcgtggga                                                                                                             >  2:1041301‑M2/1‑64 (MQ=255)
                                                                                        ccAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCgcaccggattcgtgcgtggga                                                                                                             <  1:1041301‑M2/85‑22 (MQ=255)
                                                                                                      cTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATg                              >  1:601522/1‑150 (MQ=255)
                                                                                                         atgaCCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATc                           <  2:29778/150‑1 (MQ=255)
                                                                                                             ccGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCa                       <  2:601522/150‑1 (MQ=255)
                                                                                                             ccGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCa                       >  1:211110/1‑150 (MQ=255)
                                                                                                                gAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTg                                                 <  2:39766/121‑1 (MQ=255)
                                                                                                                gAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTg                                                 >  1:39766/1‑121 (MQ=255)
                                                                                                                gAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGc                    >  1:95129/1‑150 (MQ=255)
                                                                                                                    gTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATa                <  1:1184642/150‑1 (MQ=255)
                                                                                                                     ttCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGc                                                >  2:164997/1‑117 (MQ=255)
                                                                                                                     ttCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGc                                                <  1:164997/117‑1 (MQ=255)
                                                                                                                       ctcCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATacac             <  1:77123/150‑1 (MQ=255)
                                                                                                                                  gAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGa  >  1:882818/1‑150 (MQ=255)
                                                                                                                                        ggTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATacac             >  2:548115/1‑133 (MQ=255)
                                                                                                                                        ggTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATacac             <  1:548115/133‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                          
GCCACTACCAGCGACAGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGA  >  NC_002947/739054‑739333

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