Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_002947	3,827,335	Δ51 bp	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 (0.840)	25/250	0.2	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

GTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGATCGACGCGCCGAATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGCTCGACTGCGTCAGCATTGTCCAGGCCGACAAG  >  NC_002947/3827262‑3827518                                                                                                                             |                                                                                                                                     ccggcctcgagcagcgcccgggtgaccgccaaggtggcaacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCAccc                                                                                                                           >  1:556565‑M2/125‑136 (MQ=255)  cggcctcgagcagcgcccgggtgaccgccaaggtggcaacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCg                                                                                                                          <  2:577576‑M2/13‑1 (MQ=255)                  ccgggtgaccgccaaggtggcaacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTg                                                                                                          >  2:561442‑M2/108‑136 (MQ=255)                  ccgggtgaccgccaaggtggcaacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTg                                                                                                          >  2:329523‑M2/108‑136 (MQ=255)                  ccgggtgaccgccaaggtggcaacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTg                                                                                                          >  2:653641‑M2/108‑136 (MQ=255)                       tgaccgccaaggtggcaacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGcc                                                                                                     >  2:573712‑M2/103‑136 (MQ=255)                         accgccaaggtggcaacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCAt                                                                                                   <  2:43307‑M2/36‑1 (MQ=255)                            gccaaggtggcaacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGa                                                                                                >  1:659854‑M2/98‑136 (MQ=255)                                     gcaacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGc                                                                                       <  2:973768‑M2/48‑1 (MQ=255)                                     gcaacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGc                                                                                       <  2:808851‑M2/48‑1 (MQ=255)                                       aacgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTg                                                                                     >  2:225660‑M2/87‑136 (MQ=255)                                        acgccgttgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGa                                                                                    >  2:667036‑M2/86‑136 (MQ=255)                                               tgaaggtgaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAg                                                                             >  1:272462‑M2/79‑136 (MQ=255)                                                      gaagatggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTcg                                                                      <  1:989842‑M2/65‑1 (MQ=255)                                                            ggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAg                                                                <  1:361881‑M2/71‑1 (MQ=255)                                                            ggcttgcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAg                                                                <  1:329523‑M2/71‑1 (MQ=255)                                                                 gcggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTg                                                           >  1:468966‑M2/61‑136 (MQ=255)                                                                 gcggaccatgaccactgtgggcaagaaacgaagacagttggccctgcagcccggcgccgcATGCCTGAACCCGCTCCAGGCGTGAACTGGTGCCATAGAGGGGTTCGCTTACCGCCAGGATGTCGCGGTAGCACTg                                                           >  2:907511‑M2/61‑136 (MQ=255)                                                                   ggaccatgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTggg                                                         <  2:999862‑M2/78‑1 (MQ=255)                                                                       catgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTg                                                     <  2:556565‑M2/82‑1 (MQ=255)                                                                        atgaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGc                                                    <  1:436672‑M2/83‑1 (MQ=255)                                                                         tgaccactgtgggcaagaaacgaagcgaggtggccctgcagcccggctccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCa                                                   <  2:484173‑M2/84‑1 (MQ=255)                                                                          gaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGc                                                               <  1:975026‑M2/72‑1 (MQ=255)                                                                          gaccactgtgggcaagaaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGc                                                               >  2:975026‑M2/52‑123 (MQ=255)                                                                                          aaacgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGc                                  >  1:243503‑M2/36‑136 (MQ=255)                                                                                            acgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGCTc                                <  1:568432‑M2/103‑1 (MQ=255)                                                                                             cgaagccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGCTCg                               >  2:981761‑M2/33‑136 (MQ=255)                                                                                                agccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGCTCGACt                            <  1:364045‑M2/107‑1 (MQ=255)                                                                                                 ggcagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGCTCGACTg                           <  1:863443‑M2/108‑1 (MQ=255)                                                                                                 gccagttggccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGCTCGACTg                           <  1:152964‑M2/108‑1 (MQ=255)                                                                                                          ccctgcagcccggcgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGCTCGACTGCGTCAGCAt                  <  2:879670‑M2/117‑1 (MQ=255)                                                                                                                       cgccgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGCTCGACTGCGTCAGCATTGTCCAGGCCGAc     >  2:954111‑M2/7‑136 (MQ=255)                                                                                                                          cgaATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGCTCGACTGCGTCAGCATTGTCCAGGCCGACAAg  <  1:667036‑M2/133‑1 (MQ=255)                                                                                                                             |                                                                                                                                     GTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGATCGACGCGCCGAATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTGCAGGTGGTCCAGGGCCTGCTCGACTGCGTCAGCATTGTCCAGGCCGACAAG  >  NC_002947/3827262‑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: ‑