Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 739,197 G→A 14.9% C74C (TGC→TGT pilD ← leader peptidase/N‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947739,1970GA14.9% 108.6 / 12.1 47C74C (TGC→TGTpilDleader peptidase/N‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (21/19);  new base A (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 = 4.75e-01

GCCACTACCAGCGACAGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCAC  >  NC_002947/739054‑739346
                                                                                                                                               |                                                                                                                                                     
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                                                                                                                              <  1:1099569‑M2/141‑18 (MQ=255)
                            ccACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCgcaccggattcgtgcgt                                                                                                                              >  2:1099569‑M2/1‑124 (MQ=255)
                             cACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCgg                                                                                                                                 <  1:1083252/137‑1 (MQ=255)
                             cACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCgg                                                                                                                                 >  2:1083252/1‑137 (MQ=255)
                                   cACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAcg                                                                                                              >  2:935635/1‑150 (MQ=255)
                                        ccGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGt                                                                                                         >  1:1158728/1‑150 (MQ=255)
                                                       tGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAc                                                                                                               <  1:838865/129‑1 (MQ=255)
                                                       tGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAc                                                                                                               >  2:838865/1‑129 (MQ=255)
                                                          tGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGaa                                                                                                                 >  2:426206/1‑124 (MQ=255)
                                                          tGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGaa                                                                                                                 <  1:426206/124‑1 (MQ=255)
                                                            cGGTTTTTGCAGGATGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCAGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGcccc                                                                                     <  2:372080/150‑1 (MQ=255)
                                                                       ggACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAcgcgc                                                                                                           >  2:841000/1‑117 (MQ=255)
                                                                       ggACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAcgcgc                                                                                                           <  1:841000/117‑1 (MQ=255)
                                                                        gACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGcc                                                                                       >  2:946872/1‑136 (MQ=255)
                                                                        gACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGcc                                                                                       <  1:946872/136‑1 (MQ=255)
                                                                        gACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTgc                                                                         >  2:1184642/1‑150 (MQ=255)
                                                                        gACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCgcaccggattcgt                                                                                                                                  >  2:208935‑M2/1‑80 (MQ=255)
                                                                        gACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCgcaccggattcgt                                                                                                                                  <  1:208935‑M2/93‑14 (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                                                                                                                          <  1:1041301‑M2/85‑22 (MQ=255)
                                                                                        ccAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCgcaccggattcgtgcgtggga                                                                                                                          >  2:1041301‑M2/1‑64 (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)
                                                                                                                                              cgcAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTc                                                  <  2:158911/103‑1 (MQ=255)
                                                                                                                                              cgcAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTc                                                  >  1:158911/1‑103 (MQ=255)
                                                                                                                                               gcAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGa                               <  1:377184/121‑1 (MQ=255)
                                                                                                                                               gcAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGa                               >  2:377184/1‑121 (MQ=255)
                                                                                                                                               gcAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCccac  >  2:696450/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GCCACTACCAGCGACAGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCAC  >  NC_002947/739054‑739346

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: