Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,941,124 C→T 8.0% V26I (GTC→ATC)  PP_3475 ← type II secretion system pseudopilin PulG‑like signal peptide protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,941,1240CT 8.0% 141.9 / 2.2 50V26I (GTC→ATC) PP_3475type II secretion system pseudopilin PulG‑like signal peptide protein
Reads supporting (aligned to +/- strand):  ref base C (20/26);  new base T (2/2);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.66e-01

CCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCGGCCACGCCGCCCGCCACACCCTCGGGCGGCGGCACCAGCTG  >  NC_002947/3940986‑3941272
                                                                                                                                          |                                                                                                                                                    
ccAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagc                                                                                                                                            <  5:366/149‑1 (MQ=255)
ccAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcag                                                                                                                                             >  2:207228/1‑148 (MQ=255)
  aGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGATTGCGGTcagcag                                                                                                                                          <  8:228187/149‑1 (MQ=255)
    cccTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGAGTTCTCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagcag                                                                                                                                          >  2:353701/1‑147 (MQ=255)
     ccTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagcagcag                                                                                                                                       <  3:126001/149‑1 (MQ=255)
         gCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGGCCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCAcc                                                                                                                                   <  5:383326/149‑1 (MQ=255)
               cgcgGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAgg                                                                                                                             <  2:294882/149‑1 (MQ=255)
                gcgGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGt                                                                                                                            <  1:207228/149‑1 (MQ=255)
                    tagtagCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGt                                                                                                                        <  2:374581/149‑1 (MQ=255)
                         gCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAgc                                                                                                                   >  6:114138/1‑149 (MQ=255)
                               cAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGt                                                                                                             >  3:321990/1‑149 (MQ=255)
                                    cTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGcc                                                                                                        <  1:318767/149‑1 (MQ=255)
                                     ttGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGccg                                                                                                       >  8:13221/1‑148 (MQ=255)
                                     ttGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCt                                                                                                       >  5:212953/1‑149 (MQ=255)
                                        ggCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGGCGCCAGGGCCTCGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTAc                                                                                                    <  3:311075/149‑1 (MQ=255)
                                          cATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGt                                                                                                                            <  2:93785/123‑1 (MQ=255)
                                          cATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGt                                                                                                                            >  1:93785/1‑123 (MQ=255)
                                                   tGCCGACACACAGCAACTGGCGTTCTCGGTCGCGCTGGGCGTCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGaa                                                                                         <  8:130553/149‑1 (MQ=255)
                                                     ccGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAag                                                                                        <  2:347962/148‑1 (MQ=255)
                                                          ccACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGt                                                                                                                            <  7:336785/107‑1 (MQ=255)
                                                          ccACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGt                                                                                                                            >  8:336785/1‑107 (MQ=255)
                                                                   cTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCCGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCAGTCAGCGAAGAGGCTGCCATCAcgc                                                                         >  7:265300/1‑149 (MQ=255)
                                                                    tGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCCCCGTTGGCCGTACCATTCAGCGAAGCGGCTGCCATCAcgcg                                                                        >  4:10954/1‑149 (MQ=255)
                                                                      gCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGcc                                                                      >  5:117407/1‑149 (MQ=255)
                                                                        gTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGccct                                                                    >  5:102188/1‑149 (MQ=255)
                                                                           cgcgGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGccctcc                                                                  <  7:43041/148‑1 (MQ=255)
                                                                           cgcgGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGccctcc                                                                  <  6:212953/148‑1 (MQ=255)
                                                                            gcgGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGccctccc                                                                 <  2:293205/148‑1 (MQ=255)
                                                                                  gcgcTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGc                                                           <  7:192856/148‑1 (MQ=255)
                                                                                      tGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGc                                                                                    <  8:191495/119‑1 (MQ=255)
                                                                                      tGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGc                                                                                    >  7:191495/1‑119 (MQ=255)
                                                                                        ggcggcGCTCGCCCAGATCGTGCCGGTGGCCGCCCGGGCCCCGCTGCTGACCGCAATCCGCAGCAGCACCCCCAGGGATGTAAAGCCGCCGTTGGCCGTCCCATTCAGCGAAGCGGCTGCCATCACGCGCCCTCCCGGTGGCACCGctc                                                    >  3:200273/1‑149 (MQ=255)
                                                                                         gcggcgCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGctct                                                   <  6:117407/149‑1 (MQ=255)
                                                                                           ggcgCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTg                                                 <  5:114138/149‑1 (MQ=255)
                                                                                                 ggCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCg                                           <  7:262712/148‑1 (MQ=255)
                                                                                                  gCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCgg                                          <  7:206703/149‑1 (MQ=255)
                                                                                                         tCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCCGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCGGCCAcgcc                                   >  8:177173/1‑149 (MQ=255)
                                                                                                          cGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCGCCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCGGCCAcgccg                                  >  3:145038/1‑149 (MQ=255)
                                                                                                            tGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCGGCCACgccgcc                                >  3:17915/1‑149 (MQ=255)
                                                                                                             gCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCg                                           >  2:358183/1‑137 (MQ=255)
                                                                                                             gCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCg                                           <  1:358183/137‑1 (MQ=255)
                                                                                                               cGGTGGCCGCCAGTGCCACGCTGCTGACCGCCATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCGTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCGGCCACGccgcccgc                             <  4:200273/149‑1 (MQ=255)
                                                                                                                          aGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCGGCCACGCCGCCCGCCACACCCTCgg                  <  4:321990/149‑1 (MQ=255)
                                                                                                                                 tgctgctGATTGCGGTCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCGGCCACGCCGCCCGCCACACCCTCGggcggcgg           >  2:378377/2‑149 (MQ=255)
                                                                                                                                  gctgctGACCGCAATCAGCAGCAGCGCCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCGGCCACGCCGCCCGCCACACCCTCGggcggcggc          <  4:145038/149‑1 (MQ=255)
                                                                                                                                   ctgctgATTGCGGTCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCtgt                                                              >  3:131728/1‑96 (MQ=255)
                                                                                                                                   ctgctgATTGCGGTCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCtgt                                                              <  4:131728/96‑1 (MQ=255)
                                                                                                                                    tgctgACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCGGCCACGCCGCCCGCCACACCCTCGGGCGGCGGCAc        >  7:43487/1‑149 (MQ=255)
                                                                                                                                      ctgACCGCAATCAGCAGCAGCACCCCCACGTATCTAAAGCCGCCGGTCGCCGTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCGGCCACGCCGCCCGCCACACCCTCGGGCGGCGGCAc        <  4:17915/147‑1 (MQ=255)
                                                                                                                                         aCCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCAATCACCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCGGCCACGCCGCCCGCCACACCCACCGGCGGCGGCACCAGCt   >  3:341460/1‑149 (MQ=255)
                                                                                                                                          ccGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCGGCCACGCCGCCCGCCACACCCTCGGGCGGCGGCACCAGCTg  <  8:262226/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
CCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCGGCCACGCCGCCCGCCACACCCTCGGGCGGCGGCACCAGCTG  >  NC_002947/3940986‑3941272

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: ‑