Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,299,499	Δ1,199 bp		insH21	insH21

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1299499–1300697	1300697	1–1199	16 [0]	[0] 16	insH21	insH21

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1299498	0 (0.000)	16 (0.520)	16/478	0.7	100%	intergenic (+253/‑69)	ychE/insH21	MarC family putative inner membrane protein YchE/IS5 family transposase and trans‑activator
	?	NC_000913	1300698 =	0 (0.000)					intergenic (+151/‑484)	insH21/oppA	IS5 family transposase and trans‑activator/oligopeptide ABC transporter periplasmic binding protein

AAAAACGATAACCAACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTA  >  NC_000913/1300458‑1300860                                                                                                                                                                                                                                                 |                                                                                                                                                                   aaacggagctgccatgctccgtttacttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCa                                                                                                                                                           >  1:233478‑M2/241‑250 (MQ=255)  aacggagctgccatgctccgtttacttcgtcattatttttactttgttcccgcgccgttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgtaatgatattattacttccactgcctttggtagaggaaagtgctaaCTAATACTCAAt                                                                                                                                                         >  2:117434‑M2/240‑251 (MQ=255)     ggagctgccatgctccgtttacttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGt                                                                                                                                                      <  1:30788‑M2/15‑1 (MQ=255)                            tcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCact                                                                                                                                 <  1:240210‑M2/36‑1 (MQ=255)                                    atttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACt                                                                                                                        <  2:313681‑M2/45‑1 (MQ=255)                                                                                    aaaatactcttctcaaattacagttagttataaggatttctttaactgcttctcctcaccaacattttattttcgccacatcataatcctgggcttgatgaagaataattgaaatgatattattaattccactgactttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCCTTGCCAATCTCTATTTGattgaatt                                                                       <  2:58139‑M2/94‑1 (MQ=255)                                                                                                 caaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCAGTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtgtt                                                             <  2:6786‑M2/104‑1 (MQ=255)                                                                                                  aaattacagttagttataaggatttccttaactgcttctcctcactatcatgtcattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatcttattaattccactgcctttggcagaggaaagtgctacATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGAAATCAGAAAAGATAGACATGCTTAGCCACTCTCTACTTGCTTGAATTGAAAGAtgtttgtt                                                         >  1:51251‑M2/144‑251 (MQ=255)                                                                                                                    aggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCt                                       <  1:148484‑M2/126‑1 (MQ=255)                                                                                                                       atttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCtat                                     >  1:290621‑M2/123‑250 (MQ=255)                                                                                                                       atttcctgatcggcttctcttctccatcatgttattttcgccacatcatattcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTGAAAGTATTGTGCATTTCACTACGGGACCTGTAATCAGAAAGGATGGACATGCTTAGACAATCTCTATTTGATTGAATTGAAAGATGGTTGTTAAGGCATGGATGCAAGCtata                                    <  1:168645‑M2/129‑1 (MQ=255)                                                                                                                                      tctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTgg                                                                                                                            <  1:3707‑M2/41‑1 (MQ=255)                                                                                                                                      tctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTgg                                                                                                                            >  2:3707‑M2/108‑148 (MQ=255)                                                                                                                                               catcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGAATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATaaa               >  2:197645‑M2/99‑248 (MQ=255)                                                                                                                                                      ttattttcgccacatcatactcctgggcttgctgaagaataattgacttgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATAAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGATAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGAGTGGATGCAAGCTATAGATTATTATACGGTCCATAAAAGAgaatggc     >  2:42648‑M2/92‑248 (MQ=255)                                                                                                                                                          ttgcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaacgtgctaaATAATAATCAATTGTTAAATTATTGTGCATTGCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTa  <  1:117434‑M2/163‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                   AAAAACGATAACCAACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTA  >  NC_000913/1300458‑1300860

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