Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,762,587	A→T	23.5%	Q631L (CAG→CTG)	PP_3325 →	ferric siderophore receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,762,587	0	A	T	23.5%	29.1 / 7.1	17	Q631L (CAG→CTG)	PP_3325	ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base A (6/7);  new base T (2/2);  total (8/9)
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.65e-01

GACGCAGCCCTGCCCAAAGGCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCT‑CCGGCCAGGCCGGCAGCGACTTCACCCTGGACGCCTACAGCACGGTCGACCTGCTGGCGCACTACAAGGCCAGCGAACAGGTCACCGTCGGCCTTAACCTCAACAACCTGTTCGACCGCAAGTACTA  >  NC_002947/3762442‑3762707                                                                                                                                                   |                                                                                                                         gACGCAGCCCTGCCCAAAGGCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTGCCGGCCTGGc                                                                                                                       >  2:73660/1‑150 (MQ=255) gACGCAGCCCTGCCCAAAGGCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCT‑CCGGCCAGGcc                                                                                                                      >  2:430534/1‑150 (MQ=255)          cTGCCCAAAGGCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCT‑CCGGCCAGGCCGGCAGCGAc                                                                                                             <  1:45263/150‑1 (MQ=255)          cTGCCCAAAGGCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCT‑CCGGCCAGGCCGGCAGCCAc                                                                                                             <  1:508575/150‑1 (MQ=255)           tGCCCAAAGGCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTGCCGGCCTGGCCggagcgttc                                                                                                             <  1:73660/150‑8 (MQ=255)              ccAAAGGCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCT‑CCGGCCAGGCCGGCAGCGACTTCa                                                                                                         >  1:229068/1‑150 (MQ=255)                 aaGGCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCT‑CCGGCCAGGCCGGCAGCGACTTCAccc                                                                                                      >  2:316799/1‑150 (MQ=255)                   ggCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCT‑CCGGCCAGGCCGGCAGCGACTTCACCCTg                                                                                                    <  1:430534/150‑1 (MQ=255)                                         cccGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCT‑CCGGCCAGGCCGGCAGCGACTTCACCCTGGACGCCTACAGCACGGTCGAcc                                                                              >  1:351113/1‑150 (MQ=255)                                                   cGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCT‑CCGGCCAGGCCGGCAGCGACTTCACCCTGGACGCCTACAGCACGGTCGACCTGCTGGCGCa                                                                    <  2:190905/150‑1 (MQ=255)                                                             aTGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTGCCGGCCTGGCCggagc                                                                                                                 <  2:306148/96‑4 (MQ=255)                                                             aTGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTGCCGGCCTGGCCggagc                                                                                                                 >  1:306148/1‑93 (MQ=255)                                                                    tgtgCGAAGTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGTGCAGCGGTGAACTACGTCGGCGAACGCT‑CCGGCCAGGCCGGCAGCGACTTCGCCCTGGACGCCTACAGCACGGTCGACCTGCTGGCGCACTACAAGGCCAGCGAAc                                                   <  1:455860/150‑1 (MQ=255)                                                                    tgtACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCT‑CCGGCCAGGCCGGCAGCGACTTCACCCTGGACGCCTACAGCACGGTCGACCTGCTGGCGCACTACa                                                               <  1:332475/138‑1 (MQ=255)                                                                    tgtACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCT‑CCGGCCAGGCCGGCAGCGACTTCACCCTGGACGCCTACAGCACGGTCGACCTGCTGGCGCACTACa                                                               >  2:332475/1‑138 (MQ=255)                                                                                          ccTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCT‑CCGGCCAGGCCGGCAGCGACTTCACCCTGGACGCCTACAGCACGGTCGACCTGCTGGCGCACTACAAGGCCAGCGAACAGGTCACCGTCGGCCTTAACCt                             <  2:212206/150‑1 (MQ=255)                                                                                                                     ggTGAACTACGTCGGCGAACGCT‑CCGGCCAGGCCGGCAGCGACTTCACCCTGGACGCCTACAGCACGGTCGACCTGCTGGCGCACTACAAGGCCAGCGAACAGGTCACCGTCGGCCTTAACCTCAACAACCTGTTCGACCGCAAGtacta  >  1:284036/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                         GACGCAGCCCTGCCCAAAGGCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCT‑CCGGCCAGGCCGGCAGCGACTTCACCCTGGACGCCTACAGCACGGTCGACCTGCTGGCGCACTACAAGGCCAGCGAACAGGTCACCGTCGGCCTTAACCTCAACAACCTGTTCGACCGCAAGTACTA  >  NC_002947/3762442‑3762707

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