Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_000913	1,299,499	Δ1,199 bp	100%		insH21	insH21

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1299499–1300693	1300697	5–1199	28 [0]	[2] 30	insH21	insH21

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1299498	0 (0.000)	24 (0.940)	22/256	NT	100%	intergenic (+253/‑33)	ychE/insH21	putative inner membrane protein/IS5 transposase and trans‑activator
	?	NC_000913	1300698 =	0 (0.000)					intergenic (+151/‑484)	insH21/oppA	IS5 transposase and trans‑activator/oligopeptide ABC transporter periplasmic binding protein

ACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTT  >  NC_000913/1299361‑1299637                                                                                                                                          |                                                                                                                                            aCAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataat                                                                                                                                        >  1:252524‑M1/1‑138 (MQ=255)         ttaTAAGGATTTCCTTAACTGCTTCTCCTCACCATTATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaatt                                                                                                                                <  1:280171‑M1/143‑14 (MQ=255)               ggATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTAGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaa                                                                                                                          <  2:435464‑M1/143‑20 (MQ=255)                gATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaat                                                                                                                         <  1:288445‑M1/143‑21 (MQ=255)                 aTTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGCAAGTGCTAAataataatcaattgttaaatt                                                                                                                        >  1:111686‑M1/1‑122 (MQ=255)                            tGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcattt                                                                                                             >  1:503263‑M1/1‑111 (MQ=255)                              cTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCCGGGCTTGCTGAAGAATCATTGAAATGATATTAATAATCTGACTTCTTTTGGTAGAGGAACGTGCTCAataataatcaattattaacttattgtgcatttca                                                                                                           >  2:5181‑M1/1‑109 (MQ=255)                                              tGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaat                                                                                           <  2:435886‑M1/143‑51 (MQ=255)                                              tGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGAAGAGGAAAGTGATCAaaaataatcaattgttcaattatagtgcatttcactactagaactgtaat                                                                                           >  2:203785‑M1/1‑93 (MQ=255)                                              tGTTATTTTCGCCACATCATAAGCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaat                                                                                           <  2:99108‑M1/143‑51 (MQ=255)                                                ttATTTTCGACACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatca                                                                                         <  2:442563‑M1/143‑53 (MQ=255)                                                     ttCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaa                                                                                    <  1:487831‑M1/143‑58 (MQ=255)                                                        gCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagat                                                                                 <  2:457004‑M1/143‑61 (MQ=255)                                                              tcatAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacat                                                                           <  1:64857‑M1/143‑67 (MQ=255)                                                                atAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCTACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgttcatttaactactggtactttcatccgacaagatagacattt                                                                         >  2:20833‑M1/1‑75 (MQ=255)                                                                atAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgc                                                                         >  1:450667‑M1/1‑75 (MQ=255)                                                                                          ttGAAATGATATTATTAATTCCACTGCCGTGTTTAGACGAAAGTTCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaatt                                               <  1:456834‑M1/143‑95 (MQ=255)                                                                                             aaaTGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaa                                            <  1:414379‑M1/143‑98 (MQ=255)                                                                                                           aTTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaag                              <  1:224004‑M1/143‑112 (MQ=255)                                                                                                            ttCCACTGACTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatttccatttgattgaattgaaagatgtttgttaagg                             <  1:241523‑M1/143‑113 (MQ=255)                                                                                                                 cTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatgg                        >  1:507281‑M1/1‑26 (MQ=255)                                                                                                                    ccTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatg                     >  2:269300‑M1/1‑23 (MQ=255)                                                                                                                     cTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgc                    <  1:86785‑M1/143‑122 (MQ=255)                                                                                                                         ggTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagc                <  1:463353‑M1/143‑126 (MQ=255)                                                                                                                            agagGAAAGTGCTAAataataatcaattgttaaattgttgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctat             <  1:503119‑M1/143‑129 (MQ=255)                                                                                                                                  aaGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattc       <  2:441058‑M1/143‑135 (MQ=255)                                                                                                                                     tGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgattgttaaggcatggatgcaagctatagattctga    >  2:190447‑M1/1‑6 (MQ=255)                                                                                                                                       cTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgata  <  1:395198‑M1/143‑140 (MQ=255)                                                                                                                                          |                                                                                                                                            ACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTT  >  NC_000913/1299361‑1299637

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