Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,618,992	Δ6,688 bp		rhsB–[yhhJ]	rhsB, yhhH, yhhI, [yhhJ]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3618992	3625679	6688	21 [0]	[0] 21	rhsB–[yhhJ]	rhsB,yhhH,yhhI,[yhhJ]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3618991	0 (0.000)	21 (0.580)	18/94	0.4	100%	intergenic (+2/‑201)	nikR/rhsB	DNA‑binding transcriptional repressor NikR/rhs element protein RhsB
	?	NC_000913	3625680 =	NA (NA)					coding (1124/1125 nt)	yhhJ	ABC transporter family protein YhhJ

GGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGC  >  NC_000913/3618945‑3619029                                               |                                       ggcgTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTta                                      >  1:2132895‑M1/1‑47 (MQ=255)    gTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgc                                   >  1:2144604‑M1/1‑44 (MQ=255)    gTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgc                                   >  1:2243579‑M1/1‑44 (MQ=255)       cGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccat                                >  1:42273‑M1/1‑41 (MQ=255)          cACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgt                             >  1:12398‑M1/1‑38 (MQ=255)             gggCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcc                           >  1:732102‑M1/1‑35 (MQ=255)                 aTTTGCAGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatc                      >  1:1489238‑M1/1‑31 (MQ=255)                       aGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttc                <  1:1580362‑M1/49‑25 (MQ=255)                       aGTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttc                >  1:330215‑M1/1‑25 (MQ=255)                        gTGCTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcc               <  1:1476976‑M1/49‑26 (MQ=255)                           cTTGCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctg             <  1:2549799‑M1/48‑28 (MQ=255)                              gCCGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatc         >  1:2352124‑M1/1‑18 (MQ=255)                               ccGAAGGAAGATTGAGTtacgccattgtcccaatcgtcttcctgaatcg        <  1:2530448‑M1/49‑33 (MQ=255)                                 gaaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgca      <  1:2130431‑M1/49‑35 (MQ=255)                                 gaaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcg        >  1:2421733‑M1/1‑15 (MQ=255)                                 gaaggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcg        >  1:2712301‑M1/1‑15 (MQ=255)                                   aggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcag     >  1:2201239‑M1/1‑13 (MQ=255)                                    ggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagca   <  1:2539136‑M1/49‑38 (MQ=255)                                    ggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagca   <  1:1884926‑M1/49‑38 (MQ=255)                                    ggaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcag     >  1:913222‑M1/1‑12 (MQ=255)                                     gaagATTGAGTtacgccattgtcccaatcgtcttcctgaatcgcagcag  >  1:2089977‑M1/1‑11 (MQ=255)                                               |                                       GGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGC  >  NC_000913/3618945‑3619029

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.