Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	257,908	Δ776 bp		insB9–[crl]	insB9, insA9, [crl]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	257908	258683	776	35 [0]	[0] 35	insB9–[crl]	insB9,insA9,[crl]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 257907	0 (0.000)	32 (1.400)	26/82	0.0	100%	intergenic (+8/+16)	crl/insB9	RNA polymerase holoenzyme assembly factor Crl/IS1 transposase B
	?	NC_000913	258684 =	0 (0.000)					pseudogene (9/331 nt)	crl	RNA polymerase holoenzyme assembly factor Crl

ACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCG  >  NC_000913/258636‑258728                                                 |                                             gattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAg                                             >  1:1572555‑M2/49‑50 (MQ=255)  attgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGt                                            <  1:1648579‑M2/3‑1 (MQ=255)    tgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTg                                           >  1:1024930‑M2/46‑49 (MQ=255)     gatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGc                                          >  1:1286792‑M2/45‑49 (MQ=255)       tcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCaaa                                       <  1:1460494‑M2/8‑1 (MQ=255)       tcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCaaa                                       <  1:2029348‑M2/8‑1 (MQ=255)         aaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAg                                      <  1:737517‑M2/9‑1 (MQ=255)         aaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAg                                      <  1:606016‑M2/9‑1 (MQ=255)           aaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATa                                   <  1:1195488‑M2/12‑1 (MQ=255)            aaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATaa                                  >  1:1437‑M2/38‑50 (MQ=255)            aaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATaa                                  <  1:1947631‑M2/13‑1 (MQ=255)            aaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATaa                                  <  1:1868635‑M2/13‑1 (MQ=255)              atttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATc                                <  1:2236674‑M2/15‑1 (MQ=255)               tttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATc                                >  1:2165970‑M2/35‑49 (MQ=255)                   ccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGAtt                            <  1:1969164‑M2/19‑1 (MQ=255)                   ccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGAtt                            <  1:1229659‑M2/19‑1 (MQ=255)                   ccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTc                           <  1:372467‑M2/20‑1 (MQ=255)                      cactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCttt                        >  1:710653‑M2/28‑50 (MQ=255)                       actaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCtttt                       <  1:2292684‑M2/24‑1 (MQ=255)                           ggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTc                     <  1:651688‑M2/26‑1 (MQ=255)                           ggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCg                    <  1:652639‑M2/27‑1 (MQ=255)                           ggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGa                   >  1:1557233‑M2/23‑50 (MQ=255)                            gcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGAt                  <  1:498604‑M2/29‑1 (MQ=255)                             cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGAtt                 <  1:145758‑M2/30‑1 (MQ=255)                                 tatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCt             >  1:663097‑M2/17‑50 (MQ=255)                                  atattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGACTGTCTg            >  1:773295‑M2/16‑50 (MQ=255)                                   tattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTgg           <  1:1611898‑M2/36‑1 (MQ=255)                                   tattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTgg           <  1:2219519‑M2/36‑1 (MQ=255)                                   tattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTg            <  1:2250959‑M2/35‑1 (MQ=255)                                     ttcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCt         <  1:57362‑M2/38‑1 (MQ=255)                                     ttcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCt         <  1:1238250‑M2/38‑1 (MQ=255)                                        gtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTa      <  1:63987‑M2/41‑1 (MQ=255)                                         tgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTAt     >  1:1703679‑M2/9‑50 (MQ=255)                                          gaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATg    >  1:840746‑M2/8‑50 (MQ=255)                                            atgtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCg  >  1:717122‑M2/6‑50 (MQ=255)                                                 |                                             ACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCG  >  NC_000913/258636‑258728

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.