Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,586,030	2 bp→TC	intergenic (+113/+101)	PP_4061 → / ← PP_4063	hypothetical protein/long‑chain fatty acid‑‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,586,030	0	C	T	100.0%	49.2 / NA	16	intergenic (+113/+102)	PP_4061/PP_4063	hypothetical protein/long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (5/11);  total (5/11)
*	NC_002947	4,586,031	0	T	C	100.0%	48.8 / NA	16	intergenic (+114/+101)	PP_4061/PP_4063	hypothetical protein/long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (5/11);  total (5/11)

CACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCAGGGA‑CGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGC  >  NC_002947/4585926‑4586107                                                                                                         ||                                                                                cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAA‑t                                     <  4:178084/149‑1 (MQ=39)    ggCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAA‑Tgcg                                  >  8:114541/1‑149 (MQ=39)                aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTg                      <  1:362081/149‑1 (MQ=39)                aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTg                      <  5:1821/149‑1 (MQ=39)                             cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCggg           <  2:365210/147‑1 (MQ=38)                             cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGt          <  7:14794/148‑1 (MQ=38)                                  ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAcc     <  7:114541/148‑1 (MQ=38)                                   tCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGt                    >  7:253881/1‑132 (MQ=38)                                   tCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGt                    <  8:253881/132‑1 (MQ=38)                                          aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAA‑TGCGGCGCTGc                           >  1:27156/1‑118 (MQ=37)                                          aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAA‑TGCGGCGCTGc                           <  2:27156/118‑1 (MQ=37)                                          aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAc      <  5:103520/139‑1 (MQ=38)                                          aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAc      >  6:103520/1‑139 (MQ=38)                                              gTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCgg            <  1:26232/129‑1 (MQ=38)                                              gTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCgg            >  2:26232/1‑129 (MQ=38)                                                       gACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATTGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCcgc  <  1:252142/131‑1 (MQ=37)                                                                                                         ||                                                                                CACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCAGGGA‑CGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGC  >  NC_002947/4585926‑4586107

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