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%	110.8 / NA	33	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 (18/15);  total (18/15)

GCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCCTACCAGGGA‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTT  >  NC_002947/4585912‑4586197                                                                                                                                                    |                                                                                                                                               gCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGA‑ccg                                                                                                                                             <  2:176774/149‑3 (MQ=255) gCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGA‑ccg                                                                                                                                             >  7:129060/1‑147 (MQ=39)   ttAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGA‑ccgtg                                                                                                                                           >  1:342332/1‑145 (MQ=39)          tAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCt                                                                                                                                     <  7:337957/149‑1 (MQ=39)             aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTg                                                                                                                                   >  7:238151/1‑148 (MQ=39)                   gCACCGGCTTGAACAGCTGATCTGCTCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcgg                                                                                                                            <  6:32279/149‑1 (MQ=38)                            tGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGcccc                                                                                                                   <  8:238151/149‑1 (MQ=39)                                   cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTTAAACCCCAACAAACAGGGACCCGGGCGGCCGTGTATTCGGCGCTGCCCCggggttg                                                                                                            >  5:46812/1‑145 (MQ=255)                                   cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAgg                                                                                                            <  6:401680/149‑1 (MQ=39)                                           cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑ggg                                                                                                     <  2:342332/147‑1 (MQ=38)                                           cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑gg                                                                                                      <  6:303184/146‑1 (MQ=38)                                            gctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑gg                                                                                                      >  1:400971/1‑145 (MQ=38)                                            gctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑gg                                                                                                      <  2:400971/145‑1 (MQ=38)                                                ctCCATTGACAGGTGGAGGATGACCGGGAGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑GGGTTTAccc                                                                                              <  7:191487/149‑1 (MQ=37)                                                ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCCCTCCTAACACGGACCCTGGCGGCCGTGAATGCGGCGCGGCCCCGGGTTGGAGCGGGGTTTAcc                                                                                               >  1:42785/1‑149 (MQ=255)                                                   cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCcg                                                                                             <  7:17271/147‑1 (MQ=38)                                                        aatgtttGAGTATGCCGGGGTGTTTGGCATGGTCAAGATCTTGCATTGATCGGGCCGGCCTCTTCGGG‑GTTAAACCCGCTCCAACCCGGCACCGTGCCGGCTTTTATTGCGGCGGTGCCCCGTTTAGGTGG‑GGTTTTACCCGCGAaga                                                                                       <  6:46812/143‑1 (MQ=255)                                                           ggTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑GGGtt                                                                                                   <  5:100887/133‑1 (MQ=38)                                                           ggTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑GGGtt                                                                                                   >  6:100887/1‑133 (MQ=38)                                                             tttcggATGACGGGGGGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAACCCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑GGTTTTACCCGCGAAGAggccgg                                                                                 <  2:42785/145‑1 (MQ=37)                                                             ttGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTAGCGGGTAAACCCGCTACTAACAGGGACCCTGCCGGCCGTGAATGCGGCGCTGGCCCTGGTAGGAGC‑GGGTTTAACCGCGAAGAggccgg                                                                                 >  3:347667/1‑149 (MQ=255)                                                                                                                         cgcgggTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGGTGCCCCTGGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCCCATAAACCCCCAGATGCAGCCGCAAACACCGCCACACTGATATCAcc                       >  6:378693/5‑147 (MQ=37)                                                                                                                         cgcgggTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCAGGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACg                       >  8:134263/5‑148 (MQ=255)                                                                                                                           cgggTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCGGCCGCGCCa                                                                         >  1:250458/3‑96 (MQ=255)                                                                                                                             ggTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATg                   >  1:389016/1‑148 (MQ=255)                                                                                                                             ggTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGc                  >  3:15692/1‑149 (MQ=255)                                                                                                                              gTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCg                 >  7:256285/1‑149 (MQ=255)                                                                                                                               tAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg                <  7:134263/149‑1 (MQ=255)                                                                                                                               tAAACCCCCGCCTACCAGGGACCGTGCCGGCCTTGAAGGCGGCGCTGCCCCTGTTAGGAGC‑GGGTTAACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg                <  5:378693/149‑1 (MQ=255)                                                                                                                                   cccGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACt            >  6:287311/1‑149 (MQ=255)                                                                                                                                    ccGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACtt           >  6:437899/1‑149 (MQ=255)                                                                                                                                    ccGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACt            >  6:9185/1‑148 (MQ=255)                                                                                                                                    ccGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACt            >  6:282803/1‑148 (MQ=255)                                                                                                                                     cGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTc          <  5:9185/149‑1 (MQ=255)                                                                                                                                          cTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCAc     >  1:240095/1‑149 (MQ=255)                                                                                                                                             ccAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCtt  >  3:355916/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                               GCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCCTACCAGGGA‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTT  >  NC_002947/4585912‑4586197

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