Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,673,550 G→C 8.7% G108R (GGG→CGG)  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,5500GC 8.7% 166.9 / 4.4 55G108R (GGG→CGG) PP_4137outer membrane siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (25/25);  new base C (2/3);  total (27/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 = 7.92e-01

GCAACCTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCC  >  NC_002947/4673401‑4673694
                                                                                                                                                     |                                                                                                                                                 
gCAACCTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCA‑GGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCgg                                                                                                                                                   >  2:283711/1‑149 (MQ=255)
      tggtggACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGgcg                                                                                                                                              <  4:238775/149‑1 (MQ=255)
         tggACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCAc                                                                                                                                           <  7:101781/149‑1 (MQ=255)
            aCATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGt                                                                                                                                        <  7:30825/149‑1 (MQ=255)
            aCATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCA‑GGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGtt                                                                                                                                       <  1:283711/149‑1 (MQ=255)
               tGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGttga                                                                                                                                     <  2:145068/149‑3 (MQ=255)
                gCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCg                                                                                                                                    >  8:45275/1‑149 (MQ=255)
                gCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCg                                                                                                                                    >  3:65383/1‑149 (MQ=255)
                     gcgTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTc                                                                                                                               <  1:96090/149‑1 (MQ=255)
                         ggggCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTaca                                                                                                                           >  4:307664/1‑149 (MQ=255)
                          gggCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTCACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTacac                                                                                                                          <  5:245248/149‑1 (MQ=255)
                            tcAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCACAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGCGGGCGCACGGTTCCGGTGTCTACACCa                                                                                                                        >  1:44429/2‑149 (MQ=255)
                                  ccGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGAAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTc                                                                                                                   >  5:147266/1‑148 (MQ=255)
                                   cGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGa                                                                                                                 <  3:307664/149‑1 (MQ=255)
                                            aGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGtt                                                                                                                                       >  2:113127/1‑118 (MQ=255)
                                            aGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGtt                                                                                                                                       <  1:113127/118‑1 (MQ=255)
                                                 ccAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGa                                                                                                  >  2:66995/1‑149 (MQ=255)
                                                 ccAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGa                                                                                                  >  8:165288/1‑149 (MQ=255)
                                                           cTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCggg                                                                                        >  2:73426/1‑149 (MQ=255)
                                                             ggCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCACGCTGCCGAACCGCGGCGCGGGCCTCACCCGGGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGggcg                                                                                      <  5:201866/149‑1 (MQ=255)
                                                                cTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCa                                                                                   >  5:310686/1‑149 (MQ=255)
                                                                   ccagccagGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGc                                                                                >  1:178578/1‑149 (MQ=255)
                                                                   ccagccagGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGc                                                                                >  5:59218/1‑149 (MQ=255)
                                                                   ccagccagGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGc                                                                                <  6:59218/149‑1 (MQ=255)
                                                                    cagccagGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGtt                                                                                                                                       <  8:43704/94‑1 (MQ=255)
                                                                    cagccagGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGtt                                                                                                                                       >  7:43704/1‑94 (MQ=255)
                                                                        cagGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTAc                                                                           <  6:310686/149‑1 (MQ=255)
                                                                         atgggTGTTCGCCCAATCGGCATGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACg                                                                          <  4:190196/147‑1 (MQ=255)
                                                                           gggTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGTGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGaa                                                                        <  1:66995/149‑1 (MQ=255)
                                                                                    cccAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGc                                                                                             >  1:328471/1‑119 (MQ=255)
                                                                                     ccAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGaa                                                               <  1:73426/149‑1 (MQ=255)
                                                                                             cAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTgg                                                      <  2:178578/149‑1 (MQ=255)
                                                                                               ggTAACGTTGGCATCAAGCTGTCTATACGCGGCGCGGGCCTCAACCGGGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCt                                                    <  6:23676/149‑1 (MQ=255)
                                                                                               ggTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCt                                                    <  7:45275/149‑1 (MQ=255)
                                                                                                gTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCt                                                    >  7:209049/1‑148 (MQ=255)
                                                                                                gTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCt                                                    >  8:101845/1‑148 (MQ=255)
                                                                                                           aTCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAAcc                                                             >  6:197126/1‑128 (MQ=255)
                                                                                                           aTCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAAcc                                                             <  5:197126/128‑1 (MQ=255)
                                                                                                           aTCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCg                                         >  3:266737/1‑148 (MQ=255)
                                                                                                           aTCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCg                                         >  3:136109/1‑148 (MQ=255)
                                                                                                                  tGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGtgctg                                 >  2:111271/1‑149 (MQ=255)
                                                                                                                  tGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CATGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGtgctg                                 >  8:299571/1‑149 (MQ=255)
                                                                                                                        tACGCGGCTCGGGCCTGAACCGGGCGGCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGgc                           <  4:65383/149‑1 (MQ=255)
                                                                                                                        tACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCGCGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTTCTGCGTGgc                           >  3:210004/1‑149 (MQ=255)
                                                                                                                          cgcgGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGGAGTGCTGCGTGgcg                          <  8:209049/148‑1 (MQ=255)
                                                                                                                              gCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCaa                      >  1:280899/1‑148 (MQ=255)
                                                                                                                                    gCCTCAACGGGGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCttt               <  6:42097/149‑1 (MQ=255)
                                                                                                                                    gCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCttt               <  7:101845/149‑1 (MQ=255)
                                                                                                                                    gCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCttt               <  2:280899/149‑1 (MQ=255)
                                                                                                                                    gCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCttt               <  4:266737/149‑1 (MQ=255)
                                                                                                                                          tccGTGCGCCGGGGGGGGACGGTTCCGCTTTCTACACCATGTTCGACTG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGATGCGGTGTGGGTGAGCCGTGCCGAAGTTCTGCGTGGGGCCATTGGCTTTGAGCaa         <  4:210004/148‑1 (MQ=255)
                                                                                                                                           ccGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAg        >  7:187179/1‑149 (MQ=255)
                                                                                                                                           ccGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAg        <  7:165288/149‑1 (MQ=255)
                                                                                                                                            cGTGCGCCCCCGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTACCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGATGCGTGGCGCCAATGGCTTTGACCAAgg       >  1:246860/1‑149 (MQ=255)
                                                                                                                                              tGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCt     >  1:57397/1‑149 (MQ=255)
                                                                                                                                                 gcCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTccc  >  3:222994/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GCAACCTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGG‑CCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCC  >  NC_002947/4673401‑4673694

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