Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 739,200 G→C 13.5% H73Q (CAC→CAG pilD ← leader peptidase/N‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947739,2000GC13.5% 144.1 / 15.2 52H73Q (CAC→CAGpilDleader peptidase/N‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (23/22);  new base C (4/3);  total (27/25)
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.85e-01

GCCACTACCAGCGACAGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCA  >  NC_002947/739054‑739348
                                                                                                                                                  |                                                                                                                                                    
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                                                                                                                                    <  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                                                                                                                            <  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                                                                >  1:39766/1‑121 (MQ=255)
                                                                                                                gAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTg                                                                <  2:39766/121‑1 (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)
                                                                                                                                                 aGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCccacca  <  2:1158728/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCCACTACCAGCGACAGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCA  >  NC_002947/739054‑739348

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: