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	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.330)	22/482	1.5	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

GAAAAACGATAACCAACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGG  >  NC_000913/1300457‑1300912                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                       gaaacggagctgccatgctccgtttacttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATc                                                                                                                                                                                                                <  2:451387‑M2/9‑1 (MQ=255)   aacggagctgccatgctccgtttacttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttagaaggatttccttaactgcttctcctcaccatcatgttattgtcgtcacagcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAAt                                                                                                                                                                                                             <  2:211424‑M2/12‑1 (MQ=255)        agctgccatgctccgtttacttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTa                                                                                                                                                                                                        <  1:137737‑M2/17‑1 (MQ=255)                      gtttacttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCAtt                                                                                                                                                                                          >  1:44940‑M2/221‑251 (MQ=255)                      gtttacttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCAtt                                                                                                                                                                                          >  1:526659‑M2/221‑251 (MQ=255)                       tttacttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCAttt                                                                                                                                                                                         >  1:513306‑M2/220‑251 (MQ=255)                              cgtcattatttgtactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgcctcatcataatcctgggcttgctgaagaataattgaaatgatattattaattcctctgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCactac                                                                                                                                                                                   <  2:166417‑M2/38‑1 (MQ=255)                                     atttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACt                                                                                                                                                                            <  2:223054‑M2/45‑1 (MQ=255)                                                   cccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGa                                                                                                                                                              <  2:393371‑M2/59‑1 (MQ=255)                                                       cgcagttatcaaaagcaaaaggaataggtaaaactattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGAc                                                                                                                                                         >  1:582263‑M2/188‑251 (MQ=255)                                                                 aaaagcaaaaggaatagggaaaaatattattctctcgttacaactagttataagaatttacttaactgctactcatcaccatcatgttcttttcgccacatcataatcctggccttgctgacgaataattgaaatgatattattaattccactgcctttggtagaggcaagtgctaaATAATAATCACTTGTTAAATTATTGTGCATTTCACTACAGGAACTGTAATCAGAAAAGATAGACATGCTTAGcc                                                                                                                                               >  2:271405‑M2/178‑251 (MQ=255)                                                                      caaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAAtct                                                                                                                                          >  1:344530‑M2/173‑251 (MQ=255)                                                                                aggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATttgatt                                                                                                                                <  1:198638‑M2/89‑1 (MQ=255)                                                                                     aaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGattgaatt                                                                                                                           >  1:157553‑M2/158‑251 (MQ=255)                                                                                                  caaatgtcagttagttataaggatttccttaactgcttctcctcaccatcattttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATCATAATCAATTGTAAAATTATTGTGCATTTCACTACTGGAATTTTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtgtttg                                                                                                               <  2:612848‑M2/106‑1 (MQ=255)                                                                                                    aattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTa                                                                                                            <  2:542590‑M2/109‑1 (MQ=255)                                                                                                                          ttctttaactgcttctcctcgccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgccgttggtagaggaaagtgctaaATAATAATCAATTGTTAAAGTATTGGGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGGCAACCTCTATTTGATTGAATTGGAAGATGTTTGTTAAGGCATGGATTCAAGCTATAg                                                                                       <  2:332835‑M2/130‑1 (MQ=255)                                                                                                                                                      gttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAAtt                                                           <  2:44940‑M2/158‑1 (MQ=255)                                                                                                                                                                                      ctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGt                           >  2:512932‑M2/61‑250 (MQ=255)                                                                                                                                                                                             ataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCg                   >  2:105654‑M2/54‑251 (MQ=255)                                                                                                                                                                                                  tgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCAttt              >  2:243884‑M2/49‑251 (MQ=255)                                                                                                                                                                                                        gatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTAc         <  2:513306‑M2/208‑1 (MQ=255)                                                                                                                                                                                                              attaattccactgcctttgggagaggaaagtgataaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGACAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGCTTCTGATACGGTCACTAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCgg  >  2:337778‑M2/37‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                       GAAAAACGATAACCAACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGG  >  NC_000913/1300457‑1300912

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