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–1300695	1300697	3–1199	23 [0]	[0] 23	insH21	insH21

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1299498	0 (0.000)	23 (0.540)	22/478	0.9	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

TGAAAAACGATAACCAACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGCCAGACGATTCTTTATTGGTAAT  >  NC_000913/1300456‑1300942                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     agaaacggagctgccatgctccgtttacttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgctattattaattccactgcctttggtagaggaaagagctaaATAATAATc                                                                                                                                                                                                                                              >  2:257848‑M2/243‑251 (MQ=255)                  gctccgtttacttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATtgt                                                                                                                                                                                                                             >  2:92382‑M2/226‑251 (MQ=255)                                                 gttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGaaa                                                                                                                                                                                               <  1:217634‑M2/56‑1 (MQ=255)                                                    cccgtgtagttatcaaaagcaaaaggagtaggttaaaatattcttctcaaattacagtgagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacagcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGTAACTGTAAGCAGAAAAGAt                                                                                                                                                                                           <  2:324252‑M2/60‑1 (MQ=255)                                                                        aaaaggaataggtaaaaatattcttctcaaattacggttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAAtctc                                                                                                                                                                       >  1:311707‑M2/172‑251 (MQ=255)                                                                                                 ctcaaattacagttagttataaggatttccttaactgattctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtgttt                                                                                                                                              >  1:332305‑M2/147‑251 (MQ=255)                                                                                                            gttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATg                                                                                                                                   >  1:374182‑M2/136‑251 (MQ=255)                                                                                                              tagttataaggatttccttaactgcttctcctcaccatcatgttagtttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCGCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATgg                                                                                                                                  <  1:47488‑M2/117‑1 (MQ=255)                                                                                                                         atttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagactaattgaaatgatattataaattccactgcatttggtagaggacagtgctaaATAATAATAAATTGTTAAATTATTGTGCATTTCACTACTGAAACTGTGATCAGAAAAGATAGACATGCTTAGCAAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGTCATGGATGCAAGCtat                                                                                                                       >  2:418087‑M2/123‑250 (MQ=255)                                                                                                                                tacctgcttccccgcaccagcatgttattttcgacacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAAGCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAAATAGAAAAGATAGACATGCTTAGCCAATCTCTAGTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTc                                                                                                                 <  2:275383‑M2/134‑1 (MQ=255)                                                                                                                                       ttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggtaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGt                                                                                                        >  1:352061‑M2/109‑251 (MQ=255)                                                                                                                                             tcaccatcatgttattttcgccacatcatcatcatgggcttgctgaaaaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCACTTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAAACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGTTGTTTGTTACGGCCTGGATGCAAACTATAGATTCTGATACGGTCAATa                                                                                                   >  2:457037‑M2/103‑250 (MQ=255)                                                                                                                                                    catgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAgaga                                                                                            >  1:57433‑M2/96‑250 (MQ=255)                                                                                                                                                        ttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTAATGTAACTGTAATCAGAAAGGATAGACATGCTTAGCCACTCTCTATTTGATTGAATGGAAAGATGTTTGATAAGGCATGAATGCATGCTATAGATTCTGATACGTACAATAAAAGAGACTTGc                                                                                       >  2:17190‑M2/92‑251 (MQ=255)                                                                                                                                                               cgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACa                                                                                 <  1:374340‑M2/166‑1 (MQ=255)                                                                                                                                                                        ataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCaa                                                                        >  2:335627‑M2/76‑250 (MQ=255)                                                                                                                                                                                                     aaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTgg                                          <  1:126795‑M2/205‑1 (MQ=255)                                                                                                                                                                                                         gatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGgtaat                                      >  2:465956‑M2/43‑249 (MQ=255)                                                                                                                                                                                                          atattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTACGGCATGGATGCAAGCTATAGATTCTGCTACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGAACtt                                     >  2:193492‑M2/42‑251 (MQ=255)                                                                                                                                                                                                                ttaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCggg                               >  1:216867‑M2/36‑251 (MQ=255)                                                                                                                                                                                                                  aattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCa                             <  1:465956‑M2/218‑1 (MQ=255)                                                                                                                                                                                                                                      aggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTCTTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTATGCAAAATGTATTGGCGAGTAAAAAACGCATTTGGTACTTTCCGGGCAACCACCCGACGATTCTTTAt         >  2:142858‑M2/14‑251 (MQ=255)                                                                                                                                                                                                                                              gctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGCCAGACGATTCTTTATTGGTAAt  <  2:137535‑M2/245‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     TGAAAAACGATAACCAACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGCCAGACGATTCTTTATTGGTAAT  >  NC_000913/1300456‑1300942

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