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

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

CCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCC  >  NC_000913/1299248‑1299733                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                            ccTTTCAAATGAAACGGAGCGGCCAGGATCCGTTTACTTCGTCATTAGTTTTATTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCGCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCTCCATCATGTTATTGTCTCCACGTCATAATCCTGTGCTTGGTGAAGAATAATTGAAATGATAGGATTAATTCCACGTCCTTTTGTAGAGGAAAGTGCTaa                                                                                                                                                                                                                                             <  2:304332/251‑1 (MQ=255)           gAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatc                                                                                                                                                                                                                                    >  2:180999‑M1/1‑241 (MQ=255)            aaaCGGACGTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAAGTAAAAATATTCTTCTCAAATTACAGTTAGGTATAAGGATTTCCTTAACTGCTTCTCCTCACCAGCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGCGTAAAGTGCTAAataataatcaa                                                                                                                                                                                                                                  <  2:264910‑M1/251‑12 (MQ=255)                                     cTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcact                                                                                                                                                                                                         >  2:54626‑M1/1‑215 (MQ=255)                                          tCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactgg                                                                                                                                                                                                    >  2:43467‑M1/1‑210 (MQ=255)                                             ttattTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaac                                                                                                                                                                                                 >  1:12754‑M1/1‑207 (MQ=255)                                                 ttttACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTCAAAATATTCTTCTCAAATTCCAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTATTTCCACTGCCTTTGGTAGAGGCAAGTGCTAAatactaatcatttgttaaattattgtgcatttcactactggaactgt                                                                                                                                                                                              >  2:20025‑M1/1‑203 (MQ=255)                                                                                caaaagGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatct                                                                                                                                                              <  1:180999‑M1/251‑80 (MQ=255)                                                                                          aGGTAAAGATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattatggtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgatt                                                                                                                                                    <  1:150960‑M1/251‑90 (MQ=255)                                                                                              aaaaaTATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaat                                                                                                                                                >  1:29681‑M1/1‑158 (MQ=255)                                                                                                        ttctCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCATCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgacagatgt                                                                                                                                      >  2:243778‑M1/1‑148 (MQ=255)                                                                                                            cAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgt                                                                                                                                  <  1:86512‑M1/251‑108 (MQ=255)                                                                                                                                      ccTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagatt                                                                                                        <  2:29681‑M1/251‑134 (MQ=255)                                                                                                                                                                            acaTCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaatttt                                                                  <  2:131888‑M1/251‑172 (MQ=255)                                                                                                                                                                             catcatAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgtgaaggcatggatgcatgctatagattcggatacggtcaataaaagagaattgcttaacaattttg                                                                 <  1:292275‑M1/251‑173 (MQ=255)                                                                                                                                                                                    aTCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatg                                                          <  1:22006‑M1/251‑180 (MQ=255)                                                                                                                                                                                     tCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgt                                                         >  1:216563‑M1/1‑71 (MQ=255)                                                                                                                                                                                               gCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgtattggcgagt                                               <  2:216563‑M1/251‑191 (MQ=255)                                                                                                                                                                                                                              tCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagacaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatcaggtcaataaaagagaattgcttaacaattttgcaaaatgtattggcgagtaagaaccgcatttggtactttccgggcaacc                >  1:307073‑M1/1‑30 (MQ=255)                                                                                                                                                                                                                                      cTTTGGTAGAGGAAAGTGATAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgtattggcgagtaagaaccgcatttggtactttccgggcaaccgccagacg        >  2:211721‑M1/1‑22 (MQ=255)                                                                                                                                                                                                                                         tGGTAGAGGAAAGTGCTAAataatgatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgtattggcgagtaagaaccgcatttggtactttccgggcaaccgccagacgatt     >  2:233841‑M1/1‑19 (MQ=255)                                                                                                                                                                                                                                             agagGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgtattggcgagtaagaaccgcatttggtactttccgggcaaccgccagacgattctt  <  1:43467‑M1/250‑236 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                            CCTTTCAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCC  >  NC_000913/1299248‑1299733

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