Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,673,551 G→C 11.6% G108A (GGG→GCG) ‡ PP_4137 → outer membrane siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,673,5510GC11.6% 113.3 / 7.6 43G108A (GGG→GCG) ‡PP_4137outer membrane siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (20/18);  new base C (3/2);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

CTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCCTG  >  NC_002947/4673406‑4673696
                                                                                                                                                 |                                                                                                                                                   
cTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGgc                                                                                                                                                  >  1:129169/1‑149 (MQ=255)
   gtggACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGGGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGgcg                                                                                                                                                 <  5:191779/147‑1 (MQ=255)
     ggACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACg                                                                                                                                             >  4:101129/1‑149 (MQ=255)
      gACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACgg                                                                                                                                            <  6:263964/149‑1 (MQ=255)
       aCATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCgagggtcgccc                                                                                                                                             >  1:110958/1‑141 (MQ=255)
        cATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGtt                                                                                                                                          <  3:101129/149‑1 (MQ=255)
           gCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑g                                                                                                                                      >  2:18912/1‑149 (MQ=255)
           gCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑g                                                                                                                                      >  8:25159/1‑149 (MQ=255)
               cgcgGGGGGAAAGGCGGAGTGGCAAGCAACCAGGACGTGCTGGCCTACCTGCCAGGGGTGTTCGCCAAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑Ggtgt                                                                                                                                  <  2:110958/149‑1 (MQ=255)
                                             cAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTTCCGTTGAc                                                                                                   >  3:18153/1‑149 (MQ=255)
                                               ggACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCg                                                                                                 <  7:25159/149‑1 (MQ=255)
                                                   gTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGccc                                                                                             >  4:129555/1‑149 (MQ=255)
                                                      ccggCCTACCCGCCAGGGGGGTCCGCCCAATCCGCAGGTAACGATGGCAGCACGCTGTCGATACGCGGCGCGGGCCTCAACCGCGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCggg                                                                                          <  3:54276/147‑1 (MQ=255)
                                                           cTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCa                                                                                     >  5:184210/1‑149 (MQ=255)
                                                            tACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGtt                                                                                                                                          >  6:169068/1‑97 (MQ=255)
                                                            tACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGtt                                                                                                                                          <  5:169068/97‑1 (MQ=255)
                                                            tACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGATCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACAGGCCCGGGCGGCAc                                                                                    >  5:29092/1‑149 (MQ=255)
                                                               cagccagGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATGCTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGc                                                                                  <  1:219950/149‑1 (MQ=255)
                                                                               cccAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGaa                                                                 <  3:227106/149‑1 (MQ=255)
                                                                                ccAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAAc                                                                >  4:282305/1‑149 (MQ=255)
                                                                                 cAATCGGCAGTTAACGTTGTCATCATGGTGGCTTTACGCGGGTGGGGCCTCAACGGTGGGCCGGGGGGGCTCGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAAcc                                                               <  6:29092/149‑1 (MQ=255)
                                                                                       gCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCtgtg                                                         <  3:282305/149‑1 (MQ=255)
                                                                                         aGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGc                                                       >  6:106619/1‑149 (MQ=255)
                                                                                                  tGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTg                                              >  7:67559/1‑149 (MQ=255)
                                                                                                   ggCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGc                                             >  4:29328/1‑149 (MQ=255)
                                                                                                    gCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGcc                                            <  3:29328/149‑1 (MQ=255)
                                                                                                       tCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGaa                                         <  5:106619/149‑1 (MQ=255)
                                                                                                              gTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGtgctgc                                  >  7:2428/1‑149 (MQ=255)
                                                                                                                  aTACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTgg                              <  7:13535/149‑1 (MQ=255)
                                                                                                                   tACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGgc                             >  1:159669/1‑149 (MQ=255)
                                                                                                                    aCGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGaa                                                                 <  4:9255/112‑1 (MQ=255)
                                                                                                                    aCGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGaa                                                                 >  3:9255/1‑112 (MQ=255)
                                                                                                                     cgcgGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGgcgc                           <  1:18912/149‑1 (MQ=255)
                                                                                                                      gcgGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGCCCGGCCCGGGCGGCCCGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGtggtgc                                  >  1:224099/1‑142 (MQ=255)
                                                                                                                      gcgGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGtgctgc                                  <  2:224099/141‑1 (MQ=255)
                                                                                                                         gCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAAt                       <  8:67559/149‑1 (MQ=255)
                                                                                                                             gggCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATgg                     <  3:129555/147‑1 (MQ=255)
                                                                                                                             gggCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCt                   >  4:88036/1‑149 (MQ=255)
                                                                                                                               gCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCttt                 <  7:5742/149‑1 (MQ=255)
                                                                                                                                   cAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGAc              >  6:78127/1‑148 (MQ=255)
                                                                                                                                    aaCCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTCCGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGAcc             >  7:192084/1‑148 (MQ=255)
                                                                                                                                       cGTGCGCCCCCGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAgg         >  7:250722/1‑149 (MQ=255)
                                                                                                                                        gTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCATGCCGAAGTGATGCGTGGCGCCAATGGCTTTGACCAAgg         >  3:82742/1‑148 (MQ=255)
                                                                                                                                              cGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCCTg  >  6:260738/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCC‑GGTGTCTACACCATG‑TTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCCTG  >  NC_002947/4673406‑4673696

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