Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,586,035	+G	intergenic (+118/+97)	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,033	1	.	G	100.0%	71.5 / NA	23	intergenic (+116/+99)	PP_4061/PP_4063	hypothetical protein/long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (8/15);  total (8/15)

GCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCAGGGA‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCGG  >  NC_002947/4585887‑4586118                                                                                                                                                    |                                                                                          gccgAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGC‑gg                                                                                         <  5:271874/149‑1 (MQ=255)           gAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGTCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGc                                                                                <  1:234318/149‑1 (MQ=255)                aCCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTa                                                                           >  8:269907/1‑149 (MQ=255)                    gTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAg                                                                       <  7:122683/149‑1 (MQ=255)                       aaGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCCGGGa                                                                    <  6:237350/149‑1 (MQ=39)                       aaGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGa                                                                    <  6:364164/149‑1 (MQ=255)                            ttAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAggggcagcg                                                               >  6:15444/1‑143 (MQ=38)                              aaGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGA‑ccgtg                                                              <  1:359487/147‑5 (MQ=39)                               aGGGTCGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGCCCGGGTGCTTGGCAAGGTGAAGCTCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCg                                                           <  2:359328/149‑1 (MQ=38)                                gggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCgg                                                          >  3:67891/1‑149 (MQ=39)                                 ggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGc                                                         <  8:315113/149‑1 (MQ=39)                                        cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCtt                                                      <  8:33056/145‑1 (MQ=38)                                        cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCtt                                                      >  7:33056/1‑145 (MQ=38)                                                                  tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTTCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTATACCCGCTCATACCAGGGACC‑GTGCCGGCATTGAATG‑CGGCGCTGCCCCTGGTAGGAGC‑ggg                      <  6:261640/149‑1 (MQ=255)                                                                    cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGG‑CAGCGCC‑GCATTCAAGGCCGGCACGGTCCCTGGTAGGAGc                          <  5:15444/144‑1 (MQ=14)                                                                    cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGC‑ggg                      <  4:67891/147‑1 (MQ=38)                                                                    cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGC‑ggg                      <  4:65280/147‑1 (MQ=38)                                                                      ctgctCCATTGAAAGGTGGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGG‑CAGCGCC‑GCATTCAAGGCCGGCACGGTCCCTGGTAGGAGC‑GGGTTTa                  <  8:139236/149‑1 (MQ=14)                                                                            cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTTAACCCCCTACAAACAAGGACC‑CTGGCGGGCGTGTAAT‑CGGCGCGGCCCCGGGGAGGAGCGGGGTTTAcccc              >  1:163791/1‑147 (MQ=255)                                                                            cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCcg              <  6:355116/147‑1 (MQ=38)                                                                            cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCcg              >  5:355116/1‑147 (MQ=38)                                                                                     gTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCGAGAATG‑CGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCCGCG‑AGAggccgg  >  8:347236/1‑149 (MQ=37)                                                                                                ccGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGa                            <  1:70700/114‑1 (MQ=37)                                                                                                ccGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGa                            >  2:70700/1‑114 (MQ=37)                                                                                                                                                    |                                                                                          GCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCAGGGA‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCGG  >  NC_002947/4585887‑4586118

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