Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,761,875 C→T 8.6% A408V (GCC→GTC)  PP_3325 → ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,761,8750CT 8.6% 131.7 / 2.5 48A408V (GCC→GTC) PP_3325ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base C (23/19);  new base T (2/2);  total (27/21)
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.06e-01

AGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCTGGAGCGCTTCGAGCAGACCTCGCAGAACCGCACCACGGGTCTCAGCAACA  >  NC_002947/3761743‑3762021
                                                                                                                                    |                                                                                                                                                  
aGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGa                                                                                                                                     <  5:376403/148‑1 (MQ=255)
 gACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAg                                                                                                                                    >  7:297674/1‑148 (MQ=255)
        gttgggCTTGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATTTGTACAACCCGGTTTAGGGCTAGCGGAAGCCAGCCATCAGCCGTGCCAACGACTTCTTCGAGCACACt                                                                                                                              <  4:247044/144‑1 (MQ=255)
           gggCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGaaa                                                                                                                          <  8:297674/148‑1 (MQ=255)
                    tACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGAt                                                                                                                >  5:327202/1‑149 (MQ=255)
                      cGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATc                                                                                                               <  3:50442/148‑1 (MQ=255)
                      cGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATc                                                                                                               <  3:45611/148‑1 (MQ=255)
                       gAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCaa                                                                                                             >  8:7005/1‑149 (MQ=255)
                         aaaCTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAAc                                                                                                            <  8:132514/148‑1 (MQ=255)
                         aaaCTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAAc                                                                                                            <  4:142455/148‑1 (MQ=255)
                                cGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGa                                                                                                    <  6:327202/149‑1 (MQ=255)
                                   aaCAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTAc                                                                                                                     >  8:227575/1‑129 (MQ=255)
                                   aaCAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTAc                                                                                                                     <  7:227575/129‑1 (MQ=255)
                                         cAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACACCCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCTCACTGAAAGCTACGCGAGCAACCTGCAAGACCAGATCGc                                                                                           >  8:261541/1‑149 (MQ=255)
                                                 cccGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATTTCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCAcc                                                                                    <  5:291690/148‑1 (MQ=255)
                                                 cccGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGTCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCAcc                                                                                    >  7:80434/1‑148 (MQ=255)
                                                 cccGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGTCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCAcc                                                                                    >  1:81547/1‑148 (MQ=255)
                                                  ccGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCAcc                                                                                    >  6:340610/1‑147 (MQ=255)
                                                  acGCAAAGCGCCACGTCAGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGACAACGACTACTTCGAGCACACTGAAAGCTACGCGATCAACCTGCCAGACCAGATCGCTTACAcc                                                                                    >  6:29308/2‑147 (MQ=255)
                                                        aGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCAAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCt                                                                            <  6:49239/149‑1 (MQ=255)
                                                             cACGTCCGAAGAGTATGGGCTGGATATCTACAACACGGTATAATTCAAGCCGAAGCCATCCATCACCAGAGCCAACGACTTCTTAGAGCACACTGAAAGCTACACGATCAACCTGAAAGACCAGATCCCTTTCACCGACCGCCAgcgcc                                                                       >  4:242291/1‑148 (MQ=255)
                                                              aCGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGTCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCgg                                                                      <  8:80434/149‑1 (MQ=255)
                                                              aCGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGTCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCgg                                                                      <  2:81547/149‑1 (MQ=255)
                                                                      aGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTgg                                                              <  3:242291/149‑1 (MQ=255)
                                                                               gCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGcc                                                     >  2:76812/1‑149 (MQ=255)
                                                                                cTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCt                                                    >  2:267846/1‑149 (MQ=255)
                                                                                  ggATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCTg                                                   >  1:328334/1‑148 (MQ=255)
                                                                                         tACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAAc                                                                                                            >  5:72737/1‑84 (MQ=255)
                                                                                         tACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAAc                                                                                                            <  6:72737/84‑1 (MQ=255)
                                                                                         tACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTgg                                                              >  1:278385/1‑130 (MQ=255)
                                                                                         tACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTgg                                                              <  2:278385/130‑1 (MQ=255)
                                                                                         tACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCTGGAGCGCtt                                           >  4:138363/1‑149 (MQ=255)
                                                                                         tACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGACAACGACTTCTTCGAGCACACTGAAAGCTACGAGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCTGGAGCGCtt                                           >  3:16740/1‑149 (MQ=255)
                                                                                          aCAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCgg                                                                      >  4:36123/1‑121 (MQ=255)
                                                                                          aCAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCgg                                                                      <  3:36123/121‑1 (MQ=255)
                                                                                          aCAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCTGGAGCGCTTc                                          >  5:8031/1‑149 (MQ=255)
                                                                                              cccGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGCAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCTGGAGCGCTTCGAg                                       >  5:258861/1‑148 (MQ=255)
                                                                                               ccGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCTGGAGCGCTTCGAg                                       >  6:173868/1‑147 (MQ=255)
                                                                                               ccGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCTGGAGCGCTTCGAg                                       >  4:53673/1‑147 (MQ=255)
                                                                                               ccGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCTGGAGCGCTTCGAGCa                                     <  3:367442/149‑1 (MQ=255)
                                                                                                                gAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCTGGAGCGCTTCGAGCAGACCTCGCAGAACCGca                    >  3:128375/1‑149 (MQ=255)
                                                                                                                 aaGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCTGGAGCGCTTCGAGCAGACCTCGCAGAACCGcac                   >  7:338967/1‑149 (MQ=255)
                                                                                                                   gccagccaTCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCTGGAGCGCTTCGAGCAGACCTCGCAGAACCGcacca                 >  1:221528/1‑149 (MQ=255)
                                                                                                                        ccaTCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCTGGAGCGCTTCGAGCAGACCTCGCAGAACCGCACCACGGGt            >  7:17852/1‑149 (MQ=255)
                                                                                                                           tCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCTGGAGCGCTTCGAGCAGACCTCGCAGAACCGCACCACGGGtctc         <  7:7005/149‑1 (MQ=255)
                                                                                                                             aCCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCTGGAGCGCTTCGAGCAGACCTCGCAGAACCGCACCACGGGTCTCAg       >  3:41483/1‑149 (MQ=255)
                                                                                                                              ccccttccAACGACTTCTTCTAGCACACTGAAAGCTACTCGATCAACCTTCAAGACCAGATCGGTTTGACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCTGGAGCGCTTCGAGCAGACCTCGCAGAACCGCACCACGGGTCTCAGc      <  5:29308/143‑1 (MQ=255)
                                                                                                                                cGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCTGGAGCGCTTCGAGCAGACCTCGCAGAACCGCACCACGGGTCTCAGCaa    <  5:173868/149‑1 (MQ=255)
                                                                                                                                  tGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCTGGAGCGCTTCGAGCAGACCTCGCAGAACCGCACCACGGGTCTCAGCAACa  <  5:340610/149‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                                  
AGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACCCGTGCCAACGACTTCTTCGAGCACACTGAAAGCTACGCGATCAACCTGCAAGACCAGATCGCTTTCACCGAGCGCCTGCGCGGGTTGCTGGGCGTACGCCTGGAGCGCTTCGAGCAGACCTCGCAGAACCGCACCACGGGTCTCAGCAACA  >  NC_002947/3761743‑3762021

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