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

GAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGATCGACGCGCCGAATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGCTCGACTGCGTCAGCATTGTCCAGGCCGACAAG  >  NC_002947/3827258‑3827518                                                                                                                                 |                                                                                                                                     gcgcccggcctcgagcagcgcccgggtgaccgccaaggtggcaacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaaggaacgaagccagttggccctgcagcccggcgccgaATGCCTGc                                                                                                                               <  2:538847‑M2/8‑1 (MQ=255) gcgcccggcctcgagcagcgcccgggtgaccgccaaggtggcaacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGc                                                                                                                               <  2:389340‑M2/8‑1 (MQ=255) gcgcccggcctcgagcagcgcccgggtgaccgccaaggtggcaacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGc                                                                                                                               >  2:411189‑M2/129‑136 (MQ=255)      cggcctcgagcagcgcccgggtgaccgccaaggtggcaacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggcgctgcagcccggcgccgaATGCCTGCACCCg                                                                                                                          <  1:342134‑M2/13‑1 (MQ=255)           tcgagcagcgcccgggtgaccgccaaggtggcaacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCa                                                                                                                     >  1:215865‑M2/119‑136 (MQ=255)                                gccaaggtggcaacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATcgg                                                                                                >  1:91614‑M2/98‑135 (MQ=255)                                gccaaggtggcaacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGa                                                                                                >  1:141274‑M2/98‑136 (MQ=255)                                  caaggtggcaacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAgg                                                                                              <  2:215865‑M2/41‑1 (MQ=255)                                   aaggtggcaacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAggg                                                                                             >  1:229905‑M2/95‑136 (MQ=255)                                                  ttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCa                                                                              <  2:39938‑M2/57‑1 (MQ=255)                                                  ttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCa                                                                              <  2:330176‑M2/57‑1 (MQ=255)                                                  ttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCa                                                                              <  2:325053‑M2/57‑1 (MQ=255)                                                   tgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAg                                                                             >  1:429666‑M2/79‑136 (MQ=255)                                                   tgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAg                                                                             >  1:84229‑M2/79‑136 (MQ=255)                                                     aaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGa                                                                           <  1:240395‑M2/60‑1 (MQ=255)                                                     aaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGa                                                                           >  2:27311‑M2/77‑136 (MQ=255)                                                          gaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTcg                                                                      <  1:475066‑M2/65‑1 (MQ=255)                                                          gaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGGTCGCTGACCGCCAGGATGTcg                                                                      <  2:565186‑M2/65‑1 (MQ=255)                                                          caagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTcg                                                                      <  1:531642‑M2/65‑1 (MQ=255)                                                                ggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAg                                                                <  1:411189‑M2/71‑1 (MQ=255)                                                                 gcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGc                                                               <  2:403635‑M2/72‑1 (MQ=255)                                                                            atgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGc                                                    <  1:96442‑M2/83‑1 (MQ=255)                                                                            atgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGc                                                    <  1:538189‑M2/83‑1 (MQ=255)                                                                                ccactgtgggcaagaaacgaagccagttggccctgcagcccggcgcctaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGACATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAggt                                                >  2:50893‑M2/50‑136 (MQ=255)                                                                                           aagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGcc                                     >  1:276514‑M2/39‑136 (MQ=255)                                                                                                 cgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGCTCg                               >  2:338073‑M2/33‑136 (MQ=255)                                                                                                 cgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGCTCg                               >  2:341441‑M2/33‑136 (MQ=255)                                                                                                 cgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGCTCg                               >  2:528518‑M2/33‑136 (MQ=255)                                                                                                          ttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGCTCGACTGCGTCa                      <  2:190172‑M2/113‑1 (MQ=255)                                                                                                            ggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTAGCTGACCTCCAGCATTTCGCCGTAACCCTGGCCTTGCAGGTGGTCCAGGTCCTCCTCGACTCCGacacc                    >  2:225166‑M2/22‑134 (MQ=255)                                                                                                               cctgcagcccggcgcctaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCTGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGCTCGACTGCGTCAGCAtt                 <  1:50893‑M2/118‑1 (MQ=255)                                                                                                                          gcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGCTCGACTGCGTCAGCATTGTCCAGGCCGa      <  2:137589‑M2/129‑1 (MQ=255)                                                                                                                              cgaCTGCCGGCCCCCGCTCCAGGCGTGAACTGGTGCCATCGCGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGCTCGACTGCGTCAGCATTGTCCAGGCCGACAAg  <  1:421851‑M2/133‑1 (MQ=255)                                                                                                                                 |                                                                                                                                     GAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGATCGACGCGCCGAATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGCTCGACTGCGTCAGCATTGTCCAGGCCGACAAG  >  NC_002947/3827258‑3827518

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