Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,229,201:1	+C	12.8%	intergenic (‑143/+33)	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_002947	5,229,201	1	.	C	12.8%	159.6 / 20.2	47	intergenic (‑143/+33)	PP_4606/PP_4607	ferric siderophore receptor/transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base . (22/19);  new base C (3/3);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

CTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCTTG  >  NC_002947/5229053‑5229348                                                                                                                                                       |                                                                                                                                                    cTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggg                                                                                                                                                      >  2:158319/1‑149 (MQ=255)     aTGGACGCAGACCTGCGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTCTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggg‑gg                                                                                                                                                   <  1:58827/148‑1 (MQ=255)     aTGGACGCAGACCTGCGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggg‑gg                                                                                                                                                   >  2:58827/1‑147 (MQ=255)         aCGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggg‑ggg                                                                                                                                                  >  4:124031/1‑144 (MQ=255)         aCGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggg‑ggg                                                                                                                                                  <  3:124031/144‑1 (MQ=255)          cGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTAcc                                                                                                                                            <  7:94006/149‑1 (MQ=255)               aCCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGa                                                                                                                                       <  3:236595/149‑1 (MQ=255)                    tgAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCC‑CCCGTAGGGGCGGGTTTACCCGCGAagag                                                                                                                                   <  2:174234/148‑1 (MQ=38)                              attattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTAc                                                                                                                                             >  4:249632/1‑128 (MQ=255)                              attattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTAc                                                                                                                                             <  3:249632/128‑1 (MQ=255)                              agtcttGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGGGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAg                                                                                                                         <  5:83809/144‑1 (MQ=255)                               ttattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAgg                                                                                                                        >  7:111144/1‑148 (MQ=255)                                     cGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTAc                                                                                                                 <  8:279573/149‑1 (MQ=255)                                     cGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCGGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGGGAAGAGGCTGGCACAGGCAACTAc                                                                                                                 <  5:265976/149‑1 (MQ=255)                                            atcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAAc                                                                                                            >  7:223474/1‑147 (MQ=255)                                                   tttCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTgg                                                                                                                              <  2:302594/122‑1 (MQ=255)                                                   tttCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACCGG‑GGGTTTACCCGTGAAGAGGCTgg                                                                                                                              >  1:302594/1‑122 (MQ=255)                                                         tCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGccc                                                                                             >  8:173655/1‑149 (MQ=255)                                                           aTTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTAAAGAGTCTGGCAAAGGCAACTACCGAACCAGGGTCAAACGcccc                                                                                            >  6:307839/1‑148 (MQ=255)                                                                  ggATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGt                                                                                    >  8:64393/1‑149 (MQ=255)                                                                    tctccACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTc                                                                                   <  3:67585/145‑1 (MQ=255)                                                                        cACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAg                                                                                                                         >  5:124282/1‑106 (MQ=255)                                                                        cACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAg                                                                                                                         <  6:124282/106‑1 (MQ=255)                                                                        cACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTc                                                                                   >  4:67585/1‑144 (MQ=255)                                                                        cACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCg                                                                              >  4:202632/1‑149 (MQ=255)                                                                         aCGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGt                                                                             <  1:158319/149‑1 (MQ=255)                                                                               tGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGcac                                                                                                                           <  8:203291/97‑1 (MQ=255)                                                                               tGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGcac                                                                                                                           >  7:203291/1‑97 (MQ=255)                                                                                 ttACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAcc                                                                     <  3:202632/149‑1 (MQ=255)                                                                                            gtgCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGc                                                          >  6:292323/1‑149 (MQ=255)                                                                                               ccTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGcc                                                         >  5:160213/1‑147 (MQ=255)                                                                                                cTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAg                                                      >  2:272044/1‑149 (MQ=255)                                                                                                  ggAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGc                                                     >  5:38954/1‑148 (MQ=255)                                                                                                   gAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGc                                                     <  5:307839/147‑1 (MQ=255)                                                                                                   gAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCg                                                    >  7:61876/1‑148 (MQ=255)                                                                                                     aaTATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAt                                                 >  8:41513/1‑149 (MQ=255)                                                                                                            t‑tGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAgcg                                           <  6:160213/148‑1 (MQ=255)                                                                                                                cTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCcag                                       <  6:172580/149‑1 (MQ=255)                                                                                                                        gcCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCt                               >  8:205441/1‑149 (MQ=255)                                                                                                                              tcacgGGTAAACCC‑CCCGTAGGGGCGGGTTTACCCGCGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCa                         >  1:194240/4‑149 (MQ=37)                                                                                                                                gcgGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCAcc                       <  7:41513/149‑1 (MQ=255)                                                                                                                                acgGGTAAACCC‑CCCGTAGGGGCGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCAcc                       >  5:189684/2‑149 (MQ=38)                                                                                                                                 cgGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCg                      <  5:292323/149‑1 (MQ=255)                                                                                                                                 cgGGTAAACCC‑CCCGTAGGGGCGGGTTTACCCGCGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCg                      <  2:194240/149‑1 (MQ=38)                                                                                                                                   ggTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCgg                     >  3:108692/1‑148 (MQ=255)                                                                                                                                    gTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCggg                    >  4:171979/1‑148 (MQ=255)                                                                                                                                    gTAAACCC‑CCCGTAGGGGCGGGTTTACCCGCGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCggg                    >  2:38712/1‑148 (MQ=38)                                                                                                                                     tAAACCC‑CCCGTAGGGGCGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGaa                  <  6:189684/149‑1 (MQ=39)                                                                                                                                        aCCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCt               >  1:3881/1‑149 (MQ=255)                                                                                                                                           cGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCCCCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGcc            >  6:33517/1‑149 (MQ=255)                                                                                                                                             ccccTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGc          <  8:61876/149‑1 (MQ=255)                                                                                                                                               ccTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTg        <  8:111144/149‑1 (MQ=255)                                                                                                                                                     gg‑gggTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCTTg  <  7:64393/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                    CTTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTT‑TGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCTTG  >  NC_002947/5229053‑5229348

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