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	20 [0]	[2] 22	insH21	insH21

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

AACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGCCAGACGATTCTT  >  NC_000913/1300471‑1300932                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           atgctccgtttacttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTCAattatt                                                                                                                                                                                                                     >  1:211973‑M2/228‑251 (MQ=255)                   cattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgcgaaATAATAATCAATTGTTAAATTATTGCGCATTTCACTACTgg                                                                                                                                                                                                    <  2:63102‑M2/41‑1 (MQ=255)                       atttttactttgatcccgcgcatttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataagtatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataatttacatgatattattaattcctctgcctttggtagtcgatagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACt                                                                                                                                                                                                <  1:74522‑M2/45‑1 (MQ=255)                             actttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataagtatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataagtgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAATTTATTGTGCATTTCACTACTGGAACTGTAATCa                                                                                                                                                                                         <  1:372944‑M2/52‑1 (MQ=255)                                ttgtgcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcatgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATATTAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGa                                                                                                                                                                                       >  1:355235‑M2/197‑250 (MQ=255)                                    tcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAg                                                                                                                                                                                   <  2:293834‑M2/58‑1 (MQ=255)                                      tcgcgccgttatcaatagcaaaaggaataggtcaaactattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatactgggcttgctgaagaataattgaaatgatattattaattacactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTCTGCATTTCACTACGGGAAATGTAATcagaaaaact                                                                                                                                                                                 >  2:257193‑M2/191‑247 (MQ=255)                                                                 taggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATttgat                                                                                                                                                     >  1:36766‑M2/164‑251 (MQ=255)                                                                                         tacagttagttataaggatttccttaactgcttctcctcaccataatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAgg                                                                                                                             >  1:130785‑M2/140‑251 (MQ=255)                                                                                                   tataaggctttcctaaactgcttctcctcaccatcatgttattttccccacatcatcatcctgggcttgctgaagaataattgaaatgatattattaattccactgactttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAACTAGACATGCTTAGCCAATCTCTCTTTGATTGAATTGAAAGATGTTTGTTAAGGAATGGATGc                                                                                                                    >  2:249586‑M2/130‑250 (MQ=255)                                                                                                            ttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCtat                                                                                                             >  1:374893‑M2/121‑248 (MQ=255)                                                                                                                                        gttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAAtt                                                                               <  1:67120‑M2/158‑1 (MQ=255)                                                                                                                                        gttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaagtccactgccttttgtagaggaaagtgctaaACAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTTATTGAATTGAAAGATGTTTGTAAAGGCATGGATGCAAGCTATAGATTCTTATACGGTCAATAAAAGAGAatg                                                                               <  2:221596‑M2/158‑2 (MQ=255)                                                                                                                                           attttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagcgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAAATGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGTTACGGTCAATAAAAGAGAATTGCtt                                                                           >  1:252052‑M2/90‑251 (MQ=255)                                                                                                                                                cgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTCAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACa                                                                       <  1:307722‑M2/166‑1 (MQ=255)                                                                                                                                                            atcctgtgcttgctgaagaataattgaaatgatattaataattccactgcctgtggtagatgaaagtgctaaAAAATAATCAATTGTTAAATTATTTTGCAGTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGATTAGCCAATCTCTATTTGATTGAATGGAAAGATTTTTGTTAAGGCATGGATGCAAGCTATAGATACTGATACGGTCAATATAAGATAATTGCTTAACAAATGTACAAAATg                                                          <  2:55224‑M2/179‑1 (MQ=255)                                                                                                                                                                   gcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATGCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGc                                                   <  1:33936‑M2/186‑1 (MQ=255)                                                                                                                                                                                                         actgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAACTTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTATATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTACCAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGcc             >  1:365495‑M2/28‑251 (MQ=255)                                                                                                                                                                                                           tgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAAACGcc             >  2:112983‑M2/26‑249 (MQ=255)                                                                                                                                                                                                                    tagaggaaagcgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGCCAGACGATTCtt  >  1:387288‑M2/17‑251 (MQ=255)                                                                                                                                                                                                                                   aataataatCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTg                                                                              <  2:124409/160‑1 (MQ=255)                                                                                                                                                                                                                                   aataataatCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTg                                                                              >  1:124409/1‑160 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           AACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGCCAGACGATTCTT  >  NC_000913/1300471‑1300932

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: ‑