Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,586,057	+C	intergenic (+140/+75)	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,056	1	.	C	100.0%	95.5 / NA	29	intergenic (+139/+76)	PP_4061/PP_4063	hypothetical protein/long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (14/15);  total (14/15)

CTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCCTACCAGGGA‑CGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTG  >  NC_002947/4585913‑4586198                                                                                                                                                    |                                                                                                                                                cTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGA‑ccgt                                                                                                                                             >  2:117322/1‑146 (MQ=39)     aGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCg                                                                                                                                          <  3:329877/149‑1 (MQ=39)            aCCACGGCACCGGCGTGAACGGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTg                                                                                                                                   >  3:297874/1‑148 (MQ=38)              cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAAt                                                                                                                                <  3:267345/149‑1 (MQ=39)                        gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGc                                                                                                                      <  1:117322/149‑1 (MQ=39)                        gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGGGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGc                                                                                                                      <  8:64901/149‑1 (MQ=38)                          gTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGccc                                                                                                                    <  4:124133/149‑1 (MQ=39)                             aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTg                                                                                                                 <  2:112474/149‑1 (MQ=39)                                    gATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGAG‑GGTAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGa                                                                                                                                  <  4:105436/125‑1 (MQ=38)                                    gATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGAAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGAG‑GGTAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGa                                                                                                                                  >  3:105436/1‑125 (MQ=38)                                               ctCCATTGATAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAg                                                                                                             <  2:229073/135‑1 (MQ=38)                                               ctCCATTGAAAGGTTGAGGATGACCGGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAg                                                                                                             >  1:229073/1‑136 (MQ=38)                                                  cATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCcgcg                                                                                            <  4:297874/149‑1 (MQ=38)                                                       aacggTTGCGCATCACC‑GGGGTCTTGGCACGGGTAATATTTTGCATTGATGGGGCCGGCCTCTTCGGG‑GGTAACCCCGCTCCTCCCTGGCACCGTGCCGG‑CCTTTAATGCGGCGCTGCCCCTGGTTGTAGCGGGTTTACCCGCGAaga                                                                                        <  4:332242/145‑1 (MQ=255)                                                       aaaGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAaga                                                                                        <  6:206151/148‑1 (MQ=38)                                                                        ggtgagCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAggc                                                                                     <  6:183771/130‑1 (MQ=37)                                                                                  gCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCCTTGAATGCGGCGCTGCCCCTGGTAGGa                                                                                                           <  1:18278/104‑1 (MQ=25)                                                                                  gCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGa                                                                                                           >  2:18278/1‑103 (MQ=25)                                                                                                ttGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCt                                                                                                                  >  7:137434/1‑82 (MQ=21)                                                                                                ttGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCt                                                                                                                  <  8:137434/82‑1 (MQ=21)                                                                                                                          gcgggTAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCa                      >  1:220536/4‑149 (MQ=255)                                                                                                                            gggTAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATg                    >  8:337372/2‑149 (MQ=255)                                                                                                                                  aCCCGCTCCTACCAGCGACCGTGCCGG‑CCTTGAATGCGGCGCTCCCCCTGGAAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAAc              >  6:211864/1‑149 (MQ=39)                                                                                                                                   cccGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAggccggc                                                                                 >  5:388572/1‑81 (MQ=255)                                                                                                                                   cccGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAggccggc                                                                                 <  6:388572/81‑1 (MQ=255)                                                                                                                                          cTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCAc      >  3:283665/1‑149 (MQ=255)                                                                                                                                            aCCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCt    >  1:182270/1‑149 (MQ=255)                                                                                                                                             ccAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCtt   >  3:281044/1‑149 (MQ=255)                                                                                                                                             ccAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCtt   >  3:153610/1‑149 (MQ=255)                                                                                                                                              cAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTg  >  3:390544/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                CTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCCTACCAGGGA‑CGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTG  >  NC_002947/4585913‑4586198

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