Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_002947	3,827,335	Δ51 bp	100%	coding (626‑676/1020 nt)	ptxS ←	2‑ketogluconate utilization repressor

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_002947	3827335	3827385	51	33 [0]	[0] 33	ptxS	2‑ketogluconate utilization repressor

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3827334	0 (0.000)	31 (1.830)	21/230	NT	100%	coding (677/1020 nt)	ptxS	2‑ketogluconate utilization repressor
	?	NC_002947	3827386 =	0 (0.000)					coding (625/1020 nt)	ptxS	2‑ketogluconate utilization repressor

GCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGATCGACGCGCCGAATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTG  >  NC_002947/3827207‑3827467                                                                                                                                |                                                                                                                                      gcgcCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGGAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgc                                                                                                                               <  2:538847‑M1/136‑9 (MQ=255) gcgcCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgc                                                                                                                               <  2:389340‑M1/136‑9 (MQ=255) gcgcCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgc                                                                                                                               >  2:411189‑M1/1‑128 (MQ=255)      cGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCGCTGCAGCCCGGCGCCGAatgcctgcacccg                                                                                                                          <  1:342134‑M1/136‑14 (MQ=255)           tCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctcca                                                                                                                     >  1:215865‑M1/1‑118 (MQ=255)                                gCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgg                                                                                                >  1:91614‑M1/1‑97 (MQ=255)                                gCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcga                                                                                                >  1:141274‑M1/1‑97 (MQ=255)                                  cAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgagg                                                                                              <  2:215865‑M1/136‑42 (MQ=255)                                   aaGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggg                                                                                             >  1:229905‑M1/1‑94 (MQ=255)                                                  ttGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgcca                                                                              <  2:39938‑M1/136‑58 (MQ=255)                                                  ttGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgcca                                                                              <  2:330176‑M1/136‑58 (MQ=255)                                                  ttGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgcca                                                                              <  2:325053‑M1/136‑58 (MQ=255)                                                   tGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccag                                                                             >  1:429666‑M1/1‑78 (MQ=255)                                                   tGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccag                                                                             >  1:84229‑M1/1‑78 (MQ=255)                                                     aaGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccagga                                                                           <  1:240395‑M1/136‑61 (MQ=255)                                                     aaGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccagga                                                                           >  2:27311‑M1/1‑76 (MQ=255)                                                          gAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcg                                                                      <  1:475066‑M1/136‑66 (MQ=255)                                                          gAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgagggggtcgctgaccgccaggatgtcg                                                                      <  2:565186‑M1/136‑66 (MQ=255)                                                          caagATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcg                                                                      <  1:531642‑M1/135‑66 (MQ=255)                                                                ggCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtag                                                                <  1:411189‑M1/136‑72 (MQ=255)                                                                 gCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagc                                                               <  2:403635‑M1/136‑73 (MQ=255)                                                                            aTGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgc                                                    <  1:96442‑M1/136‑84 (MQ=255)                                                                            aTGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgc                                                    <  1:538189‑M1/136‑84 (MQ=255)                                                                                ccACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCTAatgcctgcacccgctccaggcgtgaactggtgacatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgcaggt                                                >  2:50893‑M1/1‑49 (MQ=255)                                                                                           aaGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgcaggtggtccagggcc                                     >  1:276514‑M1/1‑38 (MQ=255)                                                                                                 cGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgcaggtggtccagggcctgctcg                               >  2:338073‑M1/1‑32 (MQ=255)                                                                                                 cGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgcaggtggtccagggcctgctcg                               >  2:341441‑M1/1‑32 (MQ=255)                                                                                                 cGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgcaggtggtccagggcctgctcg                               >  2:528518‑M1/1‑32 (MQ=255)                                                                                                          ttGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgcaggtggtccagggcctgctcgactgcgtca                      <  2:190172‑M1/136‑114 (MQ=255)                                                                                                            ggCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttagctgacctccagcatttcgccgtaaccctggccttgcaggtggtccaggtcctcctcgactccgacacc                    >  2:225166‑M1/1‑21 (MQ=255)                                                                                                               ccTGCAGCCCGGCGCCTAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgctgtagccctgggcttgcaggtggtccagggcctgctcgactgcgtcagcatt                 <  1:50893‑M1/136‑119 (MQ=255)                                                                                                                          gcgcCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgcaggtggtccagggcctgctcgactgcgtcagcattgtccaggccga      <  2:137589‑M1/136‑130 (MQ=255)                                                                                                                              cGActgccggcccccgctccaggcgtgaactggtgccatcgcggggttcgctgaccgccaggatgtcgcggtagccctgggcttgcaggtggtccagggcctgctcgactgcgtcagcattgtccaggccgacaag  <  1:421851‑M1/136‑134 (MQ=255)                                                                                                                                |                                                                                                                                      GCGCCCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGATCGACGCGCCGAATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTG  >  NC_002947/3827207‑3827467

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