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

CCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGATCGACGCGCCGAATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTG  >  NC_002947/3827211‑3827467                                                                                                                            |                                                                                                                                      ccGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcaccc                                                                                                                           >  1:556565‑M1/1‑124 (MQ=255)  cGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccg                                                                                                                          <  2:577576‑M1/136‑14 (MQ=255)                  ccGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactg                                                                                                          >  2:561442‑M1/1‑107 (MQ=255)                  ccGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactg                                                                                                          >  2:329523‑M1/1‑107 (MQ=255)                  ccGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactg                                                                                                          >  2:653641‑M1/1‑107 (MQ=255)                       tGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgcc                                                                                                     >  2:573712‑M1/1‑102 (MQ=255)                         aCCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccat                                                                                                   <  2:43307‑M1/136‑37 (MQ=255)                            gCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcga                                                                                                >  1:659854‑M1/1‑97 (MQ=255)                                     gCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgc                                                                                       <  2:973768‑M1/136‑49 (MQ=255)                                     gCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgc                                                                                       <  2:808851‑M1/136‑49 (MQ=255)                                       aaCGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctg                                                                                     >  2:225660‑M1/1‑86 (MQ=255)                                        aCGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctga                                                                                    >  2:667036‑M1/1‑85 (MQ=255)                                               tGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccag                                                                             >  1:272462‑M1/1‑78 (MQ=255)                                                      gAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcg                                                                      <  1:989842‑M1/136‑66 (MQ=255)                                                            ggCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtag                                                                <  1:361881‑M1/136‑72 (MQ=255)                                                            ggCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtag                                                                <  1:329523‑M1/136‑72 (MQ=255)                                                                 gCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctg                                                           >  1:468966‑M1/1‑60 (MQ=255)                                                                 gCGGACCATGACCACTGTGGGCAAGAAACGAAGACAGTTGGCCCTGCAGCCCGGCGCCGCatgcctgaacccgctccaggcgtgaactggtgccatagaggggttcgcttaccgccaggatgtcgcggtagcactg                                                           >  2:907511‑M1/1‑60 (MQ=255)                                                                   ggACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctggg                                                         <  2:999862‑M1/136‑79 (MQ=255)                                                                       cATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttg                                                     <  2:556565‑M1/136‑83 (MQ=255)                                                                        aTGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgc                                                    <  1:436672‑M1/136‑84 (MQ=255)                                                                         tGACCACTGTGGGCAAGAAACGAAGCGAGGTGGCCCTGCAGCCCGGCTCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgca                                                   <  2:484173‑M1/136‑85 (MQ=255)                                                                          gACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagc                                                               <  1:975026‑M1/123‑73 (MQ=255)                                                                          gACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagc                                                               >  2:975026‑M1/1‑51 (MQ=255)                                                                                          aaaCGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgcaggtggtccagggcctgc                                  >  1:243503‑M1/1‑35 (MQ=255)                                                                                            aCGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgcaggtggtccagggcctgctc                                <  1:568432‑M1/136‑104 (MQ=255)                                                                                             cGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgcaggtggtccagggcctgctcg                               >  2:981761‑M1/1‑32 (MQ=255)                                                                                                aGCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgcaggtggtccagggcctgctcgact                            <  1:364045‑M1/136‑108 (MQ=255)                                                                                                 ggcagTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgcaggtggtccagggcctgctcgactg                           <  1:863443‑M1/134‑109 (MQ=255)                                                                                                 gCCAGTTGGCCCTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgcaggtggtccagggcctgctcgactg                           <  1:152964‑M1/136‑109 (MQ=255)                                                                                                          cccTGCAGCCCGGCGCCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgcaggtggtccagggcctgctcgactgcgtcagcat                  <  2:879670‑M1/136‑118 (MQ=255)                                                                                                                       cgcCGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgcaggtggtccagggcctgctcgactgcgtcagcattgtccaggccgac     >  2:954111‑M1/1‑6 (MQ=255)                                                                                                                          cGAatgcctgcacccgctccaggcgtgaactggtgccatcgaggggttcgctgaccgccaggatgtcgcggtagccctgggcttgcaggtggtccagggcctgctcgactgcgtcagcattgtccaggccgacaag  <  1:667036‑M1/136‑134 (MQ=255)                                                                                                                            |                                                                                                                                      CCGGCCTCGAGCAGCGCCCGGGTGACCGCCAAGGTGGCAACGCCGTTGAAGGTGAAGATGGCTTGCGGACCATGACCACTGTGGGCAAGAAACGAAGCCAGTTGGCCCTGCAGCCCGGCGCCGATTTCCAGTACCTGCTGACGCATGCCTGGCCGTCGGCTGATCGACGCGCCGAATGCCTGCACCCGCTCCAGGCGTGAACTGGTGCCATCGAGGGGTTCGCTGACCGCCAGGATGTCGCGGTAGCCCTGGGCTTG  >  NC_002947/3827211‑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: ‑