Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,803,831 A→G 8.8% F167S (TTT→TCT)  pilN ← type IV pili biogenesis protein PilN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,803,8310AG 8.8% 123.6 / 2.2 45F167S (TTT→TCT) pilNtype IV pili biogenesis protein PilN
Reads supporting (aligned to +/- strand):  ref base A (17/24);  new base G (2/2);  total (19/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 = 1.00e+00

CATCGCGCAGCCCAGGCTGAACAAC‑ACAACAGCGATCAGCGCTGGCGCCAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCGCGCATGAACTGCGCC‑ACCACTGCACCGGACGCAGCCAACCCTGTCACCTGCAAGCGGCCACTATCAAGGCTC  >  NC_002947/5803684‑5803977
                                                                                                                                                    |                                                                                                                                                   
cagggcgcAGCCCAGGCTGAAGAAG‑TAATGAGCGTTGAGCGCTGGCGCCAAGCATCGGAAGCGGCCTTATGCCTCGAGCACCGACAGCCAAGCGAGCATCTGCGCTGGACTCACGACGAGAACGCCGCAACACGAGAGGTGAGCAGaat                                                                                                                                                    <  4:312918/145‑2 (MQ=255)
      gcAGCCCAGGCTGAACAACAACAACAGCGATCAGCGCTGGCGCCAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCAt                                                                                                                                               <  5:377554/149‑1 (MQ=255)
               gCTGAACAAC‑ACAACAGCGATCAGCGCTGGCGCCAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCgg                                                                                                                                     >  4:310960/1‑149 (MQ=255)
                   aacaac‑ACAACAGCGATCAGCGCTGGCGCCAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCa                                                                                                                                   <  7:54258/147‑1 (MQ=255)
                       ac‑acaacaGCGATCAGCGCTGGCGCCAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGctct                                                                                                                             >  2:24473/1‑149 (MQ=255)
                        c‑acaacaGCGATCAGCGCTGGCGCCAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCtt                                                                                                                            <  3:11476/149‑1 (MQ=255)
                        c‑acaacaGCGATCAGCGCTGGCGCCAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCtt                                                                                                                            <  4:37907/149‑1 (MQ=255)
                            aacaGCGATCAGCGCTGGCGCCAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAGACTCATCGCCGTCCGGCAgg                                                                                                                                 <  4:214700/141‑1 (MQ=255)
                            aacaGCGATCAGCGCTGGCGCCAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAGACTCATCGCCGTCCGGCAgg                                                                                                                                 <  4:209974/141‑1 (MQ=255)
                            aacaGCGATCAGCGCTGGCGCCAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAGACTCATCGCCGTCCGGCAgg                                                                                                                                 >  3:209974/1‑141 (MQ=255)
                            aacaGCGATCAGCGCTGGCGCCAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAGACTCATCGCCGTCCGGCAgg                                                                                                                                 >  3:214700/1‑141 (MQ=255)
                                 cGATCAGCGCTGGCGCCAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCt                                                                                                                    >  1:121416/1‑149 (MQ=255)
                                  gATCAGCGCTGGCGCCAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCtt                                                                                                                   <  4:312667/149‑1 (MQ=255)
                                  gATCAGCGCTGGCGCCAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCtt                                                                                                                   <  8:315866/149‑1 (MQ=255)
                                  gATCAGCGCTGGCGCCAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCtt                                                                                                                   >  3:81914/1‑149 (MQ=255)
                                   aTCAGCGCTGGCGCCAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTg                                                                                                                  <  2:121416/149‑1 (MQ=255)
                                              cgcCAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTCGATCtggagctc                                                                                                       <  7:323278/149‑8 (MQ=255)
                                                ccAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATc                                                                                                       >  2:286825/1‑147 (MQ=255)
                                                ccAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGa                                                                                                     <  3:310960/149‑1 (MQ=255)
                                                ccAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGa                                                                                                     <  7:328205/149‑1 (MQ=255)
                                                    ccATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCa                                                                                                  >  5:279373/1‑148 (MQ=255)
                                                       tCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAAc                                                                                              >  8:217188/1‑149 (MQ=255)
                                                       tCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAAc                                                                                              <  1:163940/149‑1 (MQ=255)
                                                           aaGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAATTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTc                                                                                          <  1:286825/149‑1 (MQ=255)
                                                                     aGCCCTAGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTc                                                                                 >  4:299697/1‑148 (MQ=255)
                                                                            gACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGa                                                                                                     <  5:205320/121‑1 (MQ=255)
                                                                            gACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGa                                                                                                     >  6:205320/1‑121 (MQ=255)
                                                                            gACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTcgcg                                                                         >  3:156015/1‑149 (MQ=255)
                                                                                  cGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAg                                                                                       <  7:59041/129‑1 (MQ=255)
                                                                                  cGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAg                                                                                       >  8:59041/1‑129 (MQ=255)
                                                                                      agccaAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCgcgc                                                                        <  8:325209/140‑1 (MQ=255)
                                                                                      agccaAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCgcgc                                                                        >  7:325209/1‑140 (MQ=255)
                                                                                          aaGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCGCGCATGAACTGCGcc‑a                                                          >  4:186524/1‑149 (MQ=255)
                                                                                           aGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCGCGCATGAACTGCGcc‑ac                                                         <  6:279373/149‑1 (MQ=255)
                                                                                                       cgcTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCGCGCATGAACTGCGCC‑ACCACTGCAcc                                                <  3:299697/146‑1 (MQ=255)
                                                                                                         cTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCGCGCATGAACTGCGCC‑ACCACTGCACCGGACg                                           >  2:69855/1‑149 (MQ=255)
                                                                                                            gACTCACGACGAGAAGGCCGCATCACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCGCGCATGAACTGCGCC‑ACCACTGCACCGGACGCa                                         <  5:375277/148‑1 (MQ=255)
                                                                                                                  cGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCGCGCATGAACTGCGCC‑ACCACTGCACCGGACGCAGCCAAcc                                  <  4:81914/149‑1 (MQ=255)
                                                                                                                      cAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCGCGCATGAACTGCGCC‑ACCACTGCACCGGACGCAGCCAACCCTGt                              >  7:236591/1‑149 (MQ=255)
                                                                                                                          aCGCCGCAACACGAGCGGTTAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCGCGCATGAACTGCGCC‑ACCACTGCACCGGACGCAGCCAACCCTGTCAcc                          <  5:80140/149‑1 (MQ=255)
                                                                                                                              cgcAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCGCGCATGAACTGCGCC‑ACCACTGCAcc                                                >  8:298636/1‑123 (MQ=255)
                                                                                                                              cgcAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCGCGCATGAACTGCGCC‑ACCACTGCAcc                                                <  7:298636/123‑1 (MQ=255)
                                                                                                                                       aGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCGCGCATGAACTGCGCCTACCACTGCACCGGACGCAGCCAACCCTGTCACCTGCAAGCGGCCa              <  1:260984/149‑1 (MQ=255)
                                                                                                                                        gCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCGCGCATGAACTGCGCC‑ACCACTGCACCGGACGCAGCCAACCCTGTCACCTGCAAGCGGCCACt            <  3:186524/149‑1 (MQ=255)
                                                                                                                                              aGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCGCGCATGAACTGCGCC‑ACCACTGCACCGGACGCAGCCAACCCTGTCACCTGCAAGCGGCCACTATCAAg      >  7:326750/1‑149 (MQ=255)
                                                                                                                                                  gAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCGCGCATGAACTGCGCC‑ACCACTGCACCGGACGCAGCCAACCCTGTCACCTGCAAGCGGCCACTATCAAGGCTc  <  8:236591/149‑1 (MQ=255)
                                                                                                                                                    aaCTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCGCGCATGAACTGCGCC‑ACCACTGCACCGGa                                             >  3:201423/1‑104 (MQ=255)
                                                                                                                                                    aaCTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCGCGCATGAACTGCGCC‑ACCACTGCACCGGa                                             <  4:201423/104‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CATCGCGCAGCCCAGGCTGAACAAC‑ACAACAGCGATCAGCGCTGGCGCCAACCATCGGAAGCGGCCTGAGCCCTCGACCACCGCCAGCCAAGCGAGCATCTGCGCTGGACTCACGACCAGAACGCCGCAACACGAGCGGTGAGCAGAAACTCATCGCCGTCCGGCAGGCTCTTGATGCGCTTGAGCTCCAGATCGAGCAAAACGTTCGAGCGGTTCAAGTCGCGCATGAACTGCGCC‑ACCACTGCACCGGACGCAGCCAACCCTGTCACCTGCAAGCGGCCACTATCAAGGCTC  >  NC_002947/5803684‑5803977

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