Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,229,249:1	+C	15.4%	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,249	1	.	C	15.4%	77.6 / 11.0	26	intergenic (‑143/+33)	PP_4606/PP_4607	ferric siderophore receptor/transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base . (10/12);  new base C (2/2);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-01

TTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCTTGCG  >  NC_002947/5229102‑5229398                                                                                                                                                     |                                                                                                                                                      ttCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggg‑gg                                                                                                                                                     >  2:307839/1‑150 (MQ=255)          gCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCGACTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCC‑CCCGTAGGGGCGGGTTTACCCg                                                                                                                                            <  1:31521/150‑1 (MQ=38)           cAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGt                                                                                                                                           <  1:307839/150‑1 (MQ=255)           cAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGt                                                                                                                                           <  2:89380/150‑1 (MQ=255)           cAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGt                                                                                                                                           <  2:366406/150‑1 (MQ=255)            aGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGCGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTg                                                                                                                                          <  1:436865/150‑1 (MQ=255)                   tgAGGGGGTGATTATTGTGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACTCGTGGAGAGGc                                                                                                                                   >  2:255313/1‑150 (MQ=255)                   tgAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGc                                                                                                                                   >  2:301719/1‑150 (MQ=255)                       ggggTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGc                                                                                                                               >  1:215101/1‑150 (MQ=255)                                       atcaatcGCACTTTCACTCATTACACGGNTGCCACGTGCTGTTANCNCNCCTGTGCCTGGCAATATTTTTGCTGTGGGGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGGAGAGGCTGGCACGGGCAACTACCGaa                                                                                                               >  1:256294/1‑150 (MQ=255)                                          aatcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCAGGTGAAGCCGCCCCTACGCG‑AGGATTACCCGAGAAGAGGCTGGCACAGGCAACTACCGAACCa                                                                                                            >  1:228562/1‑150 (MQ=255)                                           atcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTgg                                                                                                                                >  2:435551/1‑129 (MQ=255)                                           atcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTgg                                                                                                                                <  1:435551/129‑1 (MQ=255)                                            ccGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAgg                                                                                                          <  2:256294/149‑1 (MQ=255)                                             cGCACTTGCACTTATGACACGGAGGCGACGTTCTGTCACCTTACCTTTTGCTGTAAATAGTTTTGCTGTGGCGGCCTCTTCGCGGGTAATCCCGCCCTTACGGT‑GGGTTTACCGGTGAAGAGGCGGGCACATGCAACTTCCGATCCAgtg                                                                                                         <  2:92338/150‑3 (MQ=255)                                                        tCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCC‑CCCGTAGGGGCGGGTTTACCCGCGAAGAGGc                                                                                                                                   <  2:163142/113‑1 (MQ=21)                                                        tCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCC‑CCCGTAGGGGCGGGTTTACCCGCGAAGAGGc                                                                                                                                   >  1:163142/1‑113 (MQ=25)                                                                cGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGt                                                                                      >  1:432921/1‑150 (MQ=255)                                                                  gATGCCACGTGCTGTTACCTAACCTTTCCCGTGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCt                                                                                    <  2:215101/150‑1 (MQ=255)                                                                                               cTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGc                                                       >  2:402335/1‑150 (MQ=255)                                                                                                    aaTATTTTTGCTGTGGCGGCCTCTTCGCGGGGAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATc                                                  >  1:436225/1‑150 (MQ=255)                                                                                                         ttttGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAgcg                                             >  1:338149/1‑150 (MQ=255)                                                                                                                            tCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCAc                          <  1:426528/150‑1 (MQ=255)                                                                                                                               cgGGTAAACCC‑CCCGTAGGGGCGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCgg                       >  2:52501/1‑150 (MQ=39)                                                                                                                                         cgtgccTACGGG‑TGGTTTCCACGTGAAGAGGCTGGCACATGCACGTACTGAACCAGGGTCAAACGCCCCGACTCGGTCTGCCGTGACAAGGTCAAGACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCg             <  2:109524/146‑1 (MQ=255)                                                                                                                                             ccTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGa         <  1:402335/150‑1 (MQ=255)                                                                                                                                                 cGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCtt     <  2:432921/150‑1 (MQ=255)                                                                                                                                                    g‑gggTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCTTgcg  <  1:48612/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      TTCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCTTGCG  >  NC_002947/5229102‑5229398

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