Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,229,243 C→G 8.7% intergenic (‑137/+39) PP_4606 ← / ← PP_4607 ferric siderophore receptor/transmembrane sensor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,229,2430CG 8.7% 137.8 / 4.4 46intergenic (‑137/+39)PP_4606/PP_4607ferric siderophore receptor/transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base C (21/21);  new base G (2/2);  total (23/23)
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.02e-01

GGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGAC  >  NC_002947/5229095‑5229392
                                                                                                                                                    |                                                                                                                                                      
ggTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCt                                                                                                                                                       >  1:1168961/1‑150 (MQ=255)
       ttCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggggg                                                                                                                                               >  1:1502/1‑150 (MQ=255)
              aCGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGTGTAAACCCGCCCCTACGGGGGGTTTAcc                                                                                                                                        >  1:689447/1‑150 (MQ=255)
               cGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGGACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTAccc                                                                                                                                       <  1:576063/150‑1 (MQ=255)
               cGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGTGTAAACCCGCCCCTACGGGGGGTTTAccc                                                                                                                                       <  2:689447/150‑1 (MQ=255)
               cGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTAccc                                                                                                                                       <  2:1168961/150‑1 (MQ=255)
                gCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTCCGGGGGCTTTACACg                                                                                                                                      >  2:712626/1‑150 (MQ=255)
                gCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCCCCGTAGGGGCGGGTTTAc                                                                                                                                         <  2:521086/147‑1 (MQ=39)
                gCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCCCCGTAGGGGCGGGTTTAc                                                                                                                                         >  1:521086/1‑147 (MQ=39)
                      cTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAaga                                                                                                                                >  2:510186/1‑150 (MQ=255)
                        gtgAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCGCCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAgg                                                                                                                              <  2:69488/150‑1 (MQ=255)
                            gggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTgg                                                                                                                          <  1:712626/150‑1 (MQ=255)
                                 tGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAg                                                                                                                     <  2:1502/150‑1 (MQ=255)
                                    ttattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCa                                                                                                                  >  1:853596/1‑150 (MQ=255)
                                       ttGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACgggggg                                                                                                                                              <  1:534582/119‑1 (MQ=255)
                                       ttGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACgggggg                                                                                                                                              >  2:534582/1‑119 (MQ=255)
                                                  tcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAgg                                                                                                    >  1:469661/1‑150 (MQ=255)
                                                    gCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGt                                                                                                  <  1:1125929/150‑1 (MQ=255)
                                                    gCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGt                                                                                                  >  1:597809/1‑150 (MQ=255)
                                                     cACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTc                                                                                                 >  1:63715/1‑150 (MQ=255)
                                                     cACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTc                                                                                                 >  2:351237/1‑150 (MQ=255)
                                                            aCTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGcc                                                                                          <  1:146190/150‑1 (MQ=255)
                                                              tCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGATGCTGGCACAGGCAACTACCGAACCAgg                                                                                                    <  1:195180/138‑1 (MQ=255)
                                                              tCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGGAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGATGCTGGCACAGGCAACTACCGAACCAgg                                                                                                    >  2:195180/1‑138 (MQ=255)
                                                                                      ttACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAccc                                                                >  1:842248/1‑150 (MQ=255)
                                                                                       tACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCa                                                               <  2:469661/150‑1 (MQ=255)
                                                                                              ccTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCtt                                                        <  1:510186/150‑1 (MQ=255)
                                                                                               cTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCtt                                                        <  2:1131643/149‑1 (MQ=255)
                                                                                               cTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCtt                                                        >  1:1131643/1‑149 (MQ=255)
                                                                                                      tGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCg                                                >  2:797880/1‑150 (MQ=255)
                                                                                                        gAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGaa                                              <  2:1271347/150‑1 (MQ=255)
                                                                                                          aaTATTTTTGCTGTTGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATc                                            <  2:415679/150‑1 (MQ=255)
                                                                                                          aaTATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCACACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATc                                            >  1:163671/1‑150 (MQ=255)
                                                                                                          aaTATTTTTGCTGTGACGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGCTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATc                                            <  2:894737/150‑1 (MQ=255)
                                                                                                              tttttGCTGTGGCGGCCTCTTCGCGTGTAGACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACAAGGGTCAAACGCCCCAACACGGGCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAgc                                        <  2:669646/150‑1 (MQ=255)
                                                                                                                  tGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCca                                    <  2:63715/150‑1 (MQ=255)
                                                                                                                    cTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCcagc                                  >  2:684357/1‑150 (MQ=255)
                                                                                                                                ctTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGcc                      <  2:597809/150‑1 (MQ=255)
                                                                                                                                  tcacgGGTAAACCCCCCGTAGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCgg                             >  1:621008/4‑139 (MQ=39)
                                                                                                                                  tcacgGGTAAACCCCCCGTAGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCgg                             <  2:621008/136‑1 (MQ=255)
                                                                                                                                      gggTAAACCCGCCCCTACGGGGGGTTTACCCGGGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCggg                >  1:594464/1‑150 (MQ=255)
                                                                                                                                            aCCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGTCAACTCCCGAACCAGGGTCAAGCGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTTCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTg          >  2:922640/1‑150 (MQ=255)
                                                                                                                                             cccGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGc         <  2:842248/150‑1 (MQ=255)
                                                                                                                                              ccGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTa                                                                                                              >  2:708703/1‑48 (MQ=255)
                                                                                                                                              ccGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTa                                                                                                              <  1:708703/48‑1 (MQ=255)
                                                                                                                                              ccGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGcc        >  2:1253683/1‑150 (MQ=255)
                                                                                                                                                    cTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAAGGCCCCAACACGGTCTGGCTTTAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGAc  >  1:1119640/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GGTGAACTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGAC  >  NC_002947/5229095‑5229392

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: