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

AAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCT  >  NC_000913/1299259‑1299661                                                                                                                                                                                                                                                |                                                                                                                                                                    aaaCGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatca                                                                                                                                                           >  1:233478‑M1/1‑240 (MQ=255)  aaCGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCCGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGTAATGATATTATTACTTCCACTGCCTTTGGTAGAGGAAAGTGCTAActaatactcaat                                                                                                                                                         >  2:117434‑M1/1‑239 (MQ=255)     ggAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgt                                                                                                                                                      <  1:30788‑M1/251‑16 (MQ=255)                            tCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcact                                                                                                                                 <  1:240210‑M1/249‑37 (MQ=255)                                    attTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaact                                                                                                                        <  2:313681‑M1/250‑46 (MQ=255)                                                                                    aaaaTACTCTTCTCAAATTACAGTTAGTTATAAGGATTTCTTTAACTGCTTCTCCTCACCAACATTTTATTTTCGCCACATCATAATCCTGGGCTTGATGAAGAATAATTGAAATGATATTATTAATTCCACTGACTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgccttgccaatctctatttgattgaatt                                                                       <  2:58139‑M1/251‑95 (MQ=255)                                                                                                 cAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcagttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtt                                                             <  2:6786‑M1/248‑105 (MQ=255)                                                                                                  aaaTTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACTATCATGTCATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATCTTATTAATTCCACTGCCTTTGGCAGAGGAAAGTGCTACataataatcaattgttaaattattgtgcatttcactactggaactgaaatcagaaaagatagacatgcttagccactctctacttgcttgaattgaaagatgtttgtt                                                         >  1:51251‑M1/1‑143 (MQ=255)                                                                                                                    aGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagct                                       <  1:148484‑M1/251‑127 (MQ=255)                                                                                                                       atttcctgatcggcttctcttctccATCATGTTATTTTCGCCACATCATATTCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgtgaaagtattgtgcatttcactacgggacctgtaatcagaaaggatggacatgcttagacaatctctatttgattgaattgaaagatggttgttaaggcatggatgcaagctata                                    <  1:168645‑M1/239‑130 (MQ=255)                                                                                                                       aTTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctat                                     >  1:290621‑M1/1‑122 (MQ=255)                                                                                                                                      tctcCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactgg                                                                                                                            <  1:3707‑M1/148‑42 (MQ=255)                                                                                                                                      tctcCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactgg                                                                                                                            >  2:3707‑M1/1‑107 (MQ=255)                                                                                                                                               catcatGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgaatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaa               >  2:197645‑M1/1‑98 (MQ=255)                                                                                                                                                      ttATTTTCGCCACATCATACTCCTGGGCTTGCTGAAGAATAATTGACTTGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataataaattgttaaattattgtgcatttcactactggaactgtaatcagataagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggagtggatgcaagctatagattattatacggtccataaaagagaatggc     >  2:42648‑M1/1‑91 (MQ=255)                                                                                                                                                          ttgcgccACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAACGTGCTAAataataatcaattgttaaattattgtgcattgcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgctta  <  1:117434‑M1/247‑164 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                    AAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCT  >  NC_000913/1299259‑1299661

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