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%	76.7 / NA	25	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 (10/15);  total (10/15)

GCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCA‑GGGA‑C‑GTGCC‑G‑GCCTTGAA‑TG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCG  >  NC_002947/4585887‑4586139                                                                                                                                                    |                                                                                                                  gccgAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGC‑gg                                                                                                                 <  8:29881/149‑1 (MQ=255) gccgAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGC‑gg                                                                                                                 <  8:184467/149‑1 (MQ=255)                  cGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGATAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAcc                                                                                                 <  2:64282/149‑1 (MQ=255)                            ttAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑ggggcagcg                                                                                      >  3:434254/1‑143 (MQ=38)                                 ggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTcc                                                                                                     <  1:406076/130‑1 (MQ=255)                                 ggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTcc                                                                                                     >  2:406076/1‑130 (MQ=255)                                     gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTTGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCC‑G‑GCCtt                                                                           <  8:182651/148‑1 (MQ=38)                                        cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCC‑G‑GCCTTGAA‑t                                                                      >  6:146117/1‑149 (MQ=39)                                                        aCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCC‑G‑GCCTTGAA‑TG‑CGGCGCTGCCCCTgg                                                     <  7:132655/149‑1 (MQ=39)                                                               aTCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCC‑G‑GCCTTGAA‑TG‑CGGCGCTGCCCCTGGTAGGa                                                <  2:301549/147‑1 (MQ=39)                                                                  tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCATGGCGGGTAAAACCCCAACAAACA‑GGGACC‑CTGCC‑C‑GCCAAGGA‑AC‑CGGCGCAGCCCCAGGGAGGCGCggg                                           >  6:239176/1‑149 (MQ=255)                                                                   gcgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCC‑G‑GCCTTGAA‑TG‑CGGCGCTGCCCCTGGTAGGAGCggg                                           >  8:157958/1‑148 (MQ=38)                                                                   gcgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCC‑G‑GCCTTGAA‑TG‑CGGCGCTGCCCCTGGTAGGAGCggg                                           <  7:157958/148‑1 (MQ=38)                                                                    cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGG‑CAGCGCC‑‑‑GCATTCAA‑GGCCGGCACGGTCCCTGGTAGGAGCGGGt                                          <  4:434254/148‑1 (MQ=14)                                                                         ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCC‑G‑GCCTTGAACTG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTAc                                      <  1:52811/148‑1 (MQ=37)                                                                            cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCC‑G‑GCCTTGAA‑TG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCcgcg                                 <  5:146117/149‑1 (MQ=38)                                                                              ttGAAAGGTTGAGGATGACCGGTTGCTTGGCAAGGTGAAGATCTTGCATGGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCC‑G‑GCCTTGAA‑TG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGaa                               <  8:302248/149‑1 (MQ=37)                                                                                gAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAc                                                                                                  >  6:159663/1‑86 (MQ=37)                                                                                gAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAc                                                                                                  <  5:159663/86‑1 (MQ=37)                                                                                 aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCC‑G‑GCCTTGAA‑TG‑CGGCGCGGCCCCTGGTAGGAGCGGGTTTACCCGCGAaga                             <  5:305045/148‑1 (MQ=37)                                                                                           gATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCC‑G‑GCCTTGAA‑TG‑CGGCGCTGCCCCTGGTAGGAGCggg                                           >  6:362630/1‑124 (MQ=37)                                                                                           gATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCC‑G‑GCCTTGAA‑TG‑CGGCGCTGCCCCTGGTAGGAGCggg                                           <  5:362630/124‑1 (MQ=37)                                                                                                cggggTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑ggg                                                                                            <  6:108448/73‑1 (MQ=21)                                                                                                ccGGGTCCTTGGCAAGGTGAAGATCTTGCATTCATCGCGCCGGCCTCTACGCGGGTAAACCCGCTCCTACCA‑ggg                                                                                            >  5:108448/1‑75 (MQ=18)                                                                                                                  gAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTTAAACCCCGCCGACCACGGGC‑CGGGGCCTGTGCCGTGAA‑TG‑CGGCGCTGCCCCGGGTAGGAGCGGGTTTACCCGCGAAGCGGCCGGCCGCGCCAGCACATCAACCCg  >  1:350018/1‑145 (MQ=16)                                                                                                                                 aTCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCC‑G‑GCCTTGAA‑TG‑CGGCGCTGc                                                           >  2:429073/1‑70 (MQ=17)                                                                                                                                 aTCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCC‑G‑GCCTTGAA‑TG‑CGGCGCTGc                                                           <  1:429073/70‑1 (MQ=17)                                                                                                                                                    |                                                                                                                  GCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCA‑GGGA‑C‑GTGCC‑G‑GCCTTGAA‑TG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCG  >  NC_002947/4585887‑4586139

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