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	33 [0]	[0] 33	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)	28 (1.230)	22/80	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

TCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGC  >  NC_000913/257866‑257953                                          |                                               tCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaa                                         <  1:484733‑M1/49‑8 (MQ=255) tCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgca                                          >  1:53179‑M1/1‑42 (MQ=255)   aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgca                                          >  1:66212‑M1/1‑40 (MQ=255)   aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgca                                          >  1:306480‑M1/1‑40 (MQ=255)    aaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgca                                          >  1:2122125‑M1/1‑39 (MQ=255)       aaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataa                                   <  1:1318462‑M1/49‑14 (MQ=255)       aaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataa                                   <  1:1481853‑M1/49‑14 (MQ=255)        aTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataat                                  <  1:73733‑M1/49‑15 (MQ=255)          ttACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatc                                 >  1:1462307‑M1/1‑33 (MQ=255)              cGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattc                            >  1:77599‑M1/1‑29 (MQ=255)               gCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattct                           >  1:1447510‑M1/1‑28 (MQ=255)                 aCTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattcttt                         >  1:1564315‑M1/1‑26 (MQ=255)                  cTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttt                        <  1:496774‑M1/49‑25 (MQ=255)                     ggCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcg                     <  1:226800‑M1/49‑28 (MQ=255)                     ggCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttc                      <  1:1135159‑M1/48‑27 (MQ=255)                     ggCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttc                      <  1:2115672‑M1/48‑27 (MQ=255)                     ggCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttc                      <  1:1643114‑M1/48‑27 (MQ=255)                      gCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcga                    >  1:557816‑M1/1‑21 (MQ=255)                      gCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcga                    <  1:1240712‑M1/49‑29 (MQ=255)                       cccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgat                   <  1:409651‑M1/49‑30 (MQ=255)                        ccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgatt                  >  1:1233613‑M1/1‑19 (MQ=255)                        ccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgatt                  >  1:295085‑M1/1‑19 (MQ=255)                        ccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgatt                  >  1:1761350‑M1/1‑19 (MQ=255)                          gTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgt                >  1:605899‑M1/1‑17 (MQ=255)                             tatTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctg             >  1:546837‑M1/1‑14 (MQ=255)                             tatTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtct              <  1:1060504‑M1/48‑35 (MQ=255)                                  gTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgt        <  1:1837174‑M1/49‑41 (MQ=255)                                  gTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctg         <  1:1872422‑M1/48‑40 (MQ=255)                                   tGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgta       <  1:1944911‑M1/49‑42 (MQ=255)                                   tGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgta       <  1:1138196‑M1/49‑42 (MQ=255)                                     aaGGTAagtgcaaagataatcgattctttttcgattgtctggctgtat      <  1:1965829‑M1/48‑43 (MQ=255)                                      aGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgc    >  1:1978068‑M1/1‑5 (MQ=255)                                         tAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgt  >  1:32431‑M1/1‑2 (MQ=255)                                          |                                               TCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGC  >  NC_000913/257866‑257953

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ 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.