Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,490,116	Δ11 bp	intergenic (‑53/+15)	yjgR ← / ← idnR	DUF853 family protein with NTPase fold/transcriptional repressor, 5‑gluconate‑binding

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4490116	4490126	11	42 [0]	[0] 42	yjgR/idnR	DUF853 family protein with NTPase fold/transcriptional repressor, 5‑gluconate‑binding

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4490115	0 (0.000)	42 (0.950)	35/216	0.2	100%	intergenic (‑52/+26)	yjgR/idnR	DUF853 family protein with NTPase fold/transcriptional repressor, 5‑gluconate‑binding
	?	NC_000913	4490127 =	0 (0.000)					intergenic (‑64/+14)	yjgR/idnR	DUF853 family protein with NTPase fold/transcriptional repressor, 5‑gluconate‑binding

AGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACGGGG  >  NC_000913/4490050‑4490260                                                                              |                                                                                                                                      gggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGaaa                                                                                                                               <  1:408171‑M2/9‑1 (MQ=255) gggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGaaa                                                                                                                               >  2:408171‑M2/78‑86 (MQ=255)                gttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagagaaga                                                                                                                                   >  2:412269‑M2/63‑67 (MQ=255)                gttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagagaaga                                                                                                                                   <  1:412269‑M2/5‑1 (MQ=255)                 ttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGAt                                                                                                      <  1:32910‑M2/34‑1 (MQ=255)                 ttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGAt                                                                                                      >  2:32910‑M2/62‑95 (MQ=255)                 ttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTa                                                          >  1:167407‑M2/62‑139 (MQ=255)                   cactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTgttgtggttg                                                                   >  2:671469‑M2/60‑128 (MQ=255)                   cactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTgttgtggttg                                                                   <  1:671469‑M2/69‑1 (MQ=255)                     ctcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGGGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTAtttt                                                      <  2:880923‑M2/82‑1 (MQ=255)                      tcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATAccc                                                                                         <  2:519223‑M2/47‑1 (MQ=255)                      tcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATAccc                                                                                         >  1:519223‑M2/57‑103 (MQ=255)                                tcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTaaa                                                                                                                     >  2:451639‑M2/47‑65 (MQ=255)                                tcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTaaa                                                                                                                     <  1:451639‑M2/19‑1 (MQ=255)                                       gtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAgt                                                                             >  2:943717‑M2/40‑98 (MQ=255)                                       gtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAgt                                                                             <  1:943717‑M2/59‑1 (MQ=255)                                       gtttcctacctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATa                                                                                                     <  1:396396‑M2/35‑1 (MQ=255)                                       gtttcctacctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATa                                                                                                     >  2:396396‑M2/40‑74 (MQ=255)                                        tttcctgcctggagtatagtcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAg                                                                              <  1:721377‑M2/58‑1 (MQ=255)                                        tttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAg                                                                              >  2:721377‑M2/39‑96 (MQ=255)                                          tcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACAttt                                 >  2:556163‑M2/37‑139 (MQ=255)                                           cctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTg                                >  1:698441‑M2/36‑139 (MQ=255)                                               cctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGggcg                            <  1:556163‑M2/108‑1 (MQ=255)                                                 tggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAg                                         >  2:392791‑M2/30‑124 (MQ=255)                                                 tggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAg                                         <  1:392791‑M2/95‑1 (MQ=255)                                                    agtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCAt                       >  2:717571‑M2/27‑139 (MQ=255)                                                        taggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGa                   >  2:101769‑M2/23‑139 (MQ=255)                                                             aacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAat              >  2:266243‑M2/18‑139 (MQ=255)                                                              acgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAg                                                                              >  2:748648‑M2/17‑74 (MQ=255)                                                              acgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAg                                                                              <  1:748648‑M2/58‑1 (MQ=255)                                                               cgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGtt                                                           <  2:343329‑M2/77‑1 (MQ=255)                                                               cgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGtt                                                           >  1:343329‑M2/16‑92 (MQ=255)                                                               cgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGAcc                 >  1:171859‑M2/16‑134 (MQ=255)                                                               cgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGAcc                 <  2:171859‑M2/119‑1 (MQ=255)                                                               cgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAatat            >  1:52628‑M2/16‑139 (MQ=255)                                                               cgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAatat            >  2:304928‑M2/16‑139 (MQ=255)                                                                gtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATc           >  2:720192‑M2/15‑139 (MQ=255)                                                                gtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATc           >  1:110611‑M2/15‑139 (MQ=255)                                                                     gcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGtgttg                                                                         <  1:383673‑M2/63‑1 (MQ=255)                                                                     gcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGtgttg                                                                         >  2:383673‑M2/10‑72 (MQ=255)                                                                     gcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAAc      <  2:110611‑M2/130‑1 (MQ=255)                                                                         agagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACgggg  <  1:717571‑M2/134‑1 (MQ=255)                                                                              |                                                                                                                                      AGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACGGGG  >  NC_000913/4490050‑4490260

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