Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,057,524 G→A 7.8% S46S (AGC→AGT)  PP_3575 ← ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,057,5240GA 7.8% 146.9 / 2.5 51S46S (AGC→AGT) PP_3575ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (21/26);  new base A (2/2);  total (23/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 = 9.87e-01

GCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCGGGTGCACACAACAGCAGGCTCATCGCGAACAATGGATGTGCTGCGCGCGCCAGTGGGCGCAGGCGGCAAGATTTTACTGCAGGCATTGAACA  >  NC_002947/4057377‑4057667
                                                                                                                                                   |                                                                                                                                               
gCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAggc                                                                                                                                                >  6:50574/1‑149 (MQ=255)
       cTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGgctggctg                                                                                                                                          <  5:107327/148‑1 (MQ=255)
        tGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCtggctggc                                                                                                                                        >  7:236170/1‑149 (MQ=255)
               caACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGCTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCt                                                                                                                                 >  5:298451/1‑149 (MQ=255)
                   tGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGaa                                                                                                                             >  6:241576/1‑149 (MQ=255)
                     gTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACg                                                                                                                           >  3:191657/1‑149 (MQ=255)
                      ttGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGt                                                                                                                          <  8:137625/149‑1 (MQ=255)
                        gAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCAc                                                                                                                                   >  6:4632/1‑138 (MQ=255)
                        gAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCAc                                                                                                                                   <  5:4632/138‑1 (MQ=255)
                          aGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGCACGTACgc                                                                                                                      >  8:194890/1‑149 (MQ=255)
                          aGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACgc                                                                                                                      <  5:29545/149‑1 (MQ=255)
                           gCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACgcc                                                                                                                     <  8:99217/149‑1 (MQ=255)
                             cTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCtggctggc                                                                                                                                        <  6:218773/128‑1 (MQ=255)
                             cTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCtggctggc                                                                                                                                        >  5:218773/1‑128 (MQ=255)
                                  tCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACgccgggccgg                                                                                                              <  1:248453/149‑1 (MQ=255)
                                  tCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACgccgggccgg                                                                                                              <  5:75662/149‑1 (MQ=255)
                                   cGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGa                                                                                                                              <  8:39509/132‑1 (MQ=255)
                                   cGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGa                                                                                                                              >  7:39509/1‑132 (MQ=255)
                                   cGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGACTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGc                                                                                                             >  6:174688/1‑149 (MQ=255)
                                   cGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGACTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGc                                                                                                             >  6:170485/1‑149 (MQ=255)
                                       acTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCATAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCAt                                                                                                         >  1:251896/1‑149 (MQ=255)
                                         tGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCg                                                                                                       >  5:166069/1‑149 (MQ=255)
                                          gTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCgg                                                                                                      >  3:301285/1‑149 (MQ=255)
                                           tAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCgg                                                                                                      <  7:114175/148‑1 (MQ=255)
                                           tAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCAGGAACGTACGCCGGGCCGGCTCATCg                                                                                                       <  5:12747/147‑1 (MQ=255)
                                             aCGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGccc                                                                                                   <  5:241576/149‑1 (MQ=255)
                                             aCGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGACTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGccc                                                                                                   <  5:170485/149‑1 (MQ=255)
                                             aCGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGACTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGccc                                                                                                   <  5:174688/149‑1 (MQ=255)
                                               gACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAg                                                                                                 >  1:144335/1‑149 (MQ=255)
                                               gACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAg                                                                                                 >  5:311744/1‑149 (MQ=255)
                                                   ggACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCgg                                                                                             <  4:78884/149‑1 (MQ=255)
                                                      cAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCGGGTg                                                                                          <  2:283589/149‑1 (MQ=255)
                                                       aaGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCt                                                                                                                                 >  3:3256/1‑109 (MQ=255)
                                                       aaGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCt                                                                                                                                 <  4:3256/109‑1 (MQ=255)
                                                       aaGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCGGGTGc                                                                                         <  4:191657/149‑1 (MQ=255)
                                                        aggcGTTACTTTGCCGGCCCTCTACCAATGCCTGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCGGGTGca                                                                                        <  6:298451/146‑1 (MQ=255)
                                                         gCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCGGGTGcac                                                                                       <  5:50574/149‑1 (MQ=255)
                                                             ttACTTTGCCGGGCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCGGGTGCACacaa                                                                                   <  4:182033/149‑1 (MQ=255)
                                                                cTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCGGGTGCACacaaca                                                                                 <  5:291621/148‑1 (MQ=255)
                                                                 tttGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCGGGTGCACACAAcag                                                                                >  4:252740/1‑148 (MQ=255)
                                                                  ttGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACgccgggccgg                                                                                                              >  8:289421/1‑117 (MQ=255)
                                                                  ttGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACgccgggccgg                                                                                                              <  7:289421/117‑1 (MQ=255)
                                                                          ccTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCGGGTGCACAAAACAGCAGGCTCATc                                                                      >  5:192396/1‑149 (MQ=255)
                                                                                            ccAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCGGGTGCACACAACAGCAGGCTCATCGCGAACAATGGATGtgct                                                    <  6:311744/149‑1 (MQ=255)
                                                                                                          acacCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCGGGTGCACACAACAGCAGGCTCATCGCGAACAATGGATGTGCTGCGCGCGCCAGt                                        <  6:166069/147‑1 (MQ=255)
                                                                                                          acacCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCGGGTGCACACAACAGCAGGCTCATCGCGAACAATGGATGTGCTGCGCGCGCCAGTg                                       <  7:148880/148‑1 (MQ=255)
                                                                                                          acacCGGCCTGGGCGGCGAAATGGGGCAAAGCGGCCCCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCGCATCGGCCCAGGCGGGTGCACACAACAGCAGGCTCATCGCGAACAACCGATGTGCTGCGCGCGCCAGTgg                                      <  7:194890/149‑1 (MQ=255)
                                                                                                                        ggCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCGGGTGCACACAACAGCAGGCTCATCGCGAACAATGGATGTGCTGCGCGCGCCAGTGGGCGCAggcggc                           >  2:112136/1‑146 (MQ=255)
                                                                                                                         gCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCGGGTGCACACAACAGCAGGCTCATCGCGAACAATGGATGTGCTGCGCGCGCCAGTGGGCGCAGGCGGCaa                         >  4:281378/1‑147 (MQ=255)
                                                                                                                         gCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCGGGTGCACACAACAGCAGGCTCATCGCGAACAATGGATGTGCTGCGCGCGCCAGTGGGCGCAGGCGGCaa                         <  3:281378/147‑1 (MQ=255)
                                                                                                                                gggTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCGGGTGCACACAACAGCAGGCTCATCGCGAACAATGGATGTGCTGCGCGCGCCAGTGGGCGCAGGCGGCAAGATTTTACt                >  7:117301/1‑149 (MQ=255)
                                                                                                                                        gCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCGGGTGCACACAACAGCAGGCTCATCGCGAACAATGGATGTGCTGCGCGCGCCAGTGGGCGCAGGCGGCAAGATTTTACTGCAGGCAt        >  5:34620/1‑149 (MQ=255)
                                                                                                                                              cccAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCGGGTGCACACAACAGCAGGCTCATCGCGAACAATGGATGTGCTGCGCGCGCCAGTGGGCGCAGGCGGCAAGATTTTACTGCAGGCATTGAACa  <  7:112847/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GCCGCTACTGCGCAGCAACTGGTTGAAGCCTTCATCGACACTGTAACGACCGGACAAGCCGTTACTTTGCCGGCCCTCTACCAATGCCGGGTCCAGCGAAAGGCTGACACCGGCCTGGTCGGCGAAATGGGTCAAAGCGGCACCCAGGCTGGCTGGCGCCACCTGGAACGTACGCCGGGCCGGCTCATCGGCCCAGGCGGGTGCACACAACAGCAGGCTCATCGCGAACAATGGATGTGCTGCGCGCGCCAGTGGGCGCAGGCGGCAAGATTTTACTGCAGGCATTGAACA  >  NC_002947/4057377‑4057667

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