Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,229,246 C→G 15.4% intergenic (‑140/+36) 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,2460CG15.4% 59.2 / 5.2 26intergenic (‑140/+36)PP_4606/PP_4607ferric siderophore receptor/transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base C (11/11);  new base G (2/2);  total (13/13)
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.99e-01

TTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCTT  >  NC_002947/5229102‑5229395
                                                                                                                                                |                                                                                                                                                      
ttCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggggg                                                                                                                                                  >  2:307839/1‑150 (MQ=255)
         gCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCGACTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCCCCGTAGGGGCGGGTTTACCCg                                                                                                                                         <  1:31521/150‑1 (MQ=38)
          cAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGt                                                                                                                                        <  1:307839/150‑1 (MQ=255)
          cAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGt                                                                                                                                        <  2:89380/150‑1 (MQ=255)
          cAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGt                                                                                                                                        <  2:366406/150‑1 (MQ=255)
           aGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGCGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTg                                                                                                                                       <  1:436865/150‑1 (MQ=255)
                  tgAGGGGGTGATTATTGTGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACTCGTGGAGAGGc                                                                                                                                >  2:255313/1‑150 (MQ=255)
                  tgAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGc                                                                                                                                >  2:301719/1‑150 (MQ=255)
                      ggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGc                                                                                                                            >  1:215101/1‑150 (MQ=255)
                                      atcaatcGCACTTTCACTCATTACACGGNTGCCACGTGCTGTTANCNCNCCTGTGCCTGGCAATATTTTTGCTGTGGGGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGGAGAGGCTGGCACGGGCAACTACCGaa                                                                                                            >  1:256294/1‑150 (MQ=255)
                                         aatcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCAGGTGAAGCCGCCCCTACGCGAGGATTACCCGAGAAGAGGCTGGCACAGGCAACTACCGAACCa                                                                                                         >  1:228562/1‑150 (MQ=255)
                                          atcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTgg                                                                                                                             >  2:435551/1‑129 (MQ=255)
                                          atcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTgg                                                                                                                             <  1:435551/129‑1 (MQ=255)
                                           ccGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAgg                                                                                                       <  2:256294/149‑1 (MQ=255)
                                            cGCACTTGCACTTATGACACGGAGGCGACGTTCTGTCACCTTACCTTTTGCTGTAAATAGTTTTGCTGTGGCGGCCTCTTCGCGGGTAATCCCGCCCTTACGGTGGGTTTACCGGTGAAGAGGCGGGCACATGCAACTTCCGATCCAgtg                                                                                                      <  2:92338/150‑3 (MQ=255)
                                                       tCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCCCCGTAGGGGCGGGTTTACCCGCGAAGAGGc                                                                                                                                <  2:163142/113‑1 (MQ=21)
                                                       tCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCCCCGTAGGGGCGGGTTTACCCGCGAAGAGGc                                                                                                                                >  1:163142/1‑113 (MQ=25)
                                                               cGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGt                                                                                   >  1:432921/1‑150 (MQ=255)
                                                                 gATGCCACGTGCTGTTACCTAACCTTTCCCGTGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCt                                                                                 <  2:215101/150‑1 (MQ=255)
                                                                                              cTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGc                                                    >  2:402335/1‑150 (MQ=255)
                                                                                                   aaTATTTTTGCTGTGGCGGCCTCTTCGCGGGGAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATc                                               >  1:436225/1‑150 (MQ=255)
                                                                                                        ttttGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAgcg                                          >  1:338149/1‑150 (MQ=255)
                                                                                                                           tCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCAc                       <  1:426528/150‑1 (MQ=255)
                                                                                                                              cgGGTAAACCCCCCGTAGGGGCGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCgg                    >  2:52501/1‑150 (MQ=39)
                                                                                                                                        cgtgccTACGGGTGGTTTCCACGTGAAGAGGCTGGCACATGCACGTACTGAACCAGGGTCAAACGCCCCGACTCGGTCTGCCGTGACAAGGTCAAGACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCg          <  2:109524/146‑1 (MQ=255)
                                                                                                                                            ccTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGa      <  1:402335/150‑1 (MQ=255)
                                                                                                                                                cGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCtt  <  2:432921/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
TTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCTT  >  NC_002947/5229102‑5229395

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