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–1300696	1300697	2–1199	33 [0]	[0] 33	insH21	insH21

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1299498	0 (0.000)	33 (0.700)	29/488	0.4	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

ACGATAACCAACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGCCAGACGATTCTTTATTGG  >  NC_000913/1300462‑1300938                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 ggagctgccatgctccgtttacttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGt                                                                                                                                                                                                                                    >  1:151321‑M2/237‑251 (MQ=255)            gctccgtttacttcgtcattattttttctttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATtgt                                                                                                                                                                                                                         <  2:294769‑M2/26‑1 (MQ=255)                     acttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCac                                                                                                                                                                                                                >  1:101364‑M2/217‑251 (MQ=255)                                                cgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAg                                                                                                                                                                                     <  1:60787‑M2/62‑1 (MQ=255)                                                      gttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGc                                                                                                                                                                               <  2:101364‑M2/68‑1 (MQ=255)                                                           caaaagcaaaaggcataggtaaaaatattcttctcaaattacagttagttctaaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGc                                                                                                                                                                          >  2:109686‑M2/179‑251 (MQ=255)                                                                   aaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATctct                                                                                                                                                                  <  2:310805‑M2/81‑1 (MQ=255)                                                                                 aaatattcttctcaaatcacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGattgaattg                                                                                                                                                    >  1:70228‑M2/157‑251 (MQ=255)                                                                                                   acagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCAATTTGATTGAATTGAAAGATGTTTGTTAAGGc                                                                                                                                  <  2:299120‑M2/113‑1 (MQ=255)                                                                                                        tagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATgg                                                                                                                              <  1:151109‑M2/117‑1 (MQ=255)                                                                                                                           acctgcttctcctctccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatagtattaattccactgcctttggtagaggaaagtgcgaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGTATGCAAGCTATATATTCTg                                                                                                           <  2:355733‑M2/136‑1 (MQ=255)                                                                                                                                     cctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTc                                                                                                   <  1:280494‑M2/144‑1 (MQ=255)                                                                                                                                         accatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATaaa                                                                                             >  2:377721‑M2/101‑250 (MQ=255)                                                                                                                                            atcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACt                                                                                                                                                                                                      <  2:310594‑M2/45‑1 (MQ=255)                                                                                                                                            atcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACt                                                                                                                                                                                                      >  1:310594‑M2/98‑142 (MQ=255)                                                                                                                                                        tcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAAc                                                                              <  1:274091‑M2/165‑1 (MQ=255)                                                                                                                                                                        ctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTAt                                                             >  1:297494‑M2/70‑251 (MQ=255)                                                                                                                                                                          gggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTg                                                           >  1:353683‑M2/68‑251 (MQ=255)                                                                                                                                                                                   gaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGTACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTa                                                    >  2:323735‑M2/59‑249 (MQ=255)                                                                                                                                                                                              gaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTg                                       >  1:208375‑M2/48‑251 (MQ=255)                                                                                                                                                                                              gaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTg                                       >  1:173281‑M2/48‑251 (MQ=255)                                                                                                                                                                                                 atgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGt                                     >  1:58167‑M2/45‑250 (MQ=255)                                                                                                                                                                                                   gatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACt                                  >  1:190899‑M2/43‑251 (MQ=255)                                                                                                                                                                                                   gatattattaattccactgcctttggtagaggaaagagcttaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACt                                  >  1:242857‑M2/43‑251 (MQ=255)                                                                                                                                                                                                                 cactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGc                    <  1:126804‑M2/223‑1 (MQ=255)                                                                                                                                                                                                                       cttttgtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCATCCGCCAGACg              >  2:78314‑M2/23‑251 (MQ=255)                                                                                                                                                                                                                         ttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGCCAGACGa             <  1:167232‑M2/230‑1 (MQ=255)                                                                                                                                                                                                                         ttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGCCAGACGAt            >  2:369473‑M2/21‑251 (MQ=255)                                                                                                                                                                                                                              agaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTTTTAAGGCATGGATGCAAGCTATAGATGCTGATACGGTCAATAAAAGAGAATTGCTTAACAGTTTTTCAAAATGTATTGGCGAGGAAGAACCGCATTTGGAACTTTCCGGGCAACCGCCAGACGATTCttt       <  2:109067‑M2/236‑1 (MQ=255)                                                                                                                                                                                                                              agaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATTTATTTGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGCCAGACGATTCttt       <  1:109686‑M2/236‑1 (MQ=255)                                                                                                                                                                                                                              agaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGCCAGACGATTCttt       <  2:190899‑M2/236‑1 (MQ=255)                                                                                                                                                                                                                                  gaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGCCAGACGATTCTTTATTg   <  2:382043‑M2/240‑1 (MQ=255)                                                                                                                                                                                                                                   aaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGCCAGACGATTCTTTATTgg  <  2:73483‑M2/241‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 ACGATAACCAACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGCCAGACGATTCTTTATTGG  >  NC_000913/1300462‑1300938

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