Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,229,213	T→C	8.3%	intergenic (‑155/+21)	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,213	0	T	C	8.3%	132.1 / 3.0	49	intergenic (‑155/+21)	PP_4606/PP_4607	ferric siderophore receptor/transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base T (23/21);  new base C (2/2);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

ACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCTTGCGC  >  NC_002947/5229067‑5229351                                                                                                                                                    |                                                                                                                                           aCCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGa                                                                                                                                          <  3:236595/149‑1 (MQ=255)      tgAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCC‑CCCGTAGGGGCGGGTTTACCCGCGAagag                                                                                                                                      <  2:174234/148‑1 (MQ=38)                agtcttGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGGGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAg                                                                                                                            <  5:83809/144‑1 (MQ=255)                 ttattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAgg                                                                                                                           >  7:111144/1‑148 (MQ=255)                       cGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTAc                                                                                                                    <  8:279573/149‑1 (MQ=255)                       cGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCGGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGGGAAGAGGCTGGCACAGGCAACTAc                                                                                                                    <  5:265976/149‑1 (MQ=255)                              atcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAAc                                                                                                               >  7:223474/1‑147 (MQ=255)                                     tttCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTgg                                                                                                                                 <  2:302594/122‑1 (MQ=255)                                     tttCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACCGG‑GGGTTTACCCGTGAAGAGGCTgg                                                                                                                                 >  1:302594/1‑122 (MQ=255)                                           tCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGccc                                                                                                >  8:173655/1‑149 (MQ=255)                                             aTTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTAAAGAGTCTGGCAAAGGCAACTACCGAACCAGGGTCAAACGcccc                                                                                               >  6:307839/1‑148 (MQ=255)                                                    ggATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGt                                                                                       >  8:64393/1‑149 (MQ=255)                                                      tctccACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTc                                                                                      <  3:67585/145‑1 (MQ=255)                                                          cACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAg                                                                                                                            >  5:124282/1‑106 (MQ=255)                                                          cACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAg                                                                                                                            <  6:124282/106‑1 (MQ=255)                                                          cACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTc                                                                                      >  4:67585/1‑144 (MQ=255)                                                          cACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCg                                                                                 >  4:202632/1‑149 (MQ=255)                                                           aCGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGt                                                                                <  1:158319/149‑1 (MQ=255)                                                                 tGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGcac                                                                                                                              >  7:203291/1‑97 (MQ=255)                                                                 tGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGcac                                                                                                                              <  8:203291/97‑1 (MQ=255)                                                                   ttACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAcc                                                                        <  3:202632/149‑1 (MQ=255)                                                                              gtgCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGc                                                             >  6:292323/1‑149 (MQ=255)                                                                                 ccTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGcc                                                            >  5:160213/1‑147 (MQ=255)                                                                                  cTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAg                                                         >  2:272044/1‑149 (MQ=255)                                                                                    ggAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGc                                                        >  5:38954/1‑148 (MQ=255)                                                                                     gAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGc                                                        <  5:307839/147‑1 (MQ=255)                                                                                     gAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCg                                                       >  7:61876/1‑148 (MQ=255)                                                                                       aaTATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAAt                                                    >  8:41513/1‑149 (MQ=255)                                                                                              ttGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑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)                                                                                                                                        ggtgTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCTTgc    <  8:223474/146‑1 (MQ=255)                                                                                                                                          ggTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCTTgcgc  <  6:38954/149‑1 (MQ=255)                                                                                                                                                   gTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCgg                        >  4:118497/1‑118 (MQ=255)                                                                                                                                                   gTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCgg                        <  3:118497/118‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                           ACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCTTGCGC  >  NC_002947/5229067‑5229351

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