New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 718897	19 (0.480)	4 (0.120)	3/132	5.1	15.5%	intergenic (+43/‑88)	yhhI/ybfD	predicted transposase; b3484_1/hypothetical protein
	?	REL606	= 719050	27 (0.780)					coding (66/762 nt)	ybfD	hypothetical protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTT  >  REL606/718824‑718934                                                                          |                                      gTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTaa                                     <  1:333203/76‑1 (MQ=255)    aGTCCTTGCGGGGAGCGGGCTTTCGTAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTaaaaa                                  <  1:339593/76‑1 (MQ=255)     gTCCTTGCGGGGAGCGGGCTTTCGTAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAAt                                 <  2:584858/76‑1 (MQ=255)                  gCGGGCTTTCGTAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCtt                    <  2:425326/76‑1 (MQ=255)                    gggCTTTCGTAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCttat                  <  1:394407/76‑1 (MQ=255)                    gggCTTTCGTAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCttat                  <  1:394412/76‑1 (MQ=255)                        tttCGTAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTAcc              <  1:1048297/76‑1 (MQ=255)                             tAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAg         <  2:1005358/76‑1 (MQ=255)                             tAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAg         >  1:293911/1‑76 (MQ=255)                             tAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAg         <  2:449268/76‑1 (MQ=255)                              aaTCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGc        >  1:800944/1‑76 (MQ=255)                              aaTCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGc        >  2:1001224/1‑76 (MQ=255)                              aaTCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGc        >  2:1032369/1‑76 (MQ=255)                              aaTCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGc        >  2:1032382/1‑76 (MQ=255)                                 cTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGt     >  1:1124023/1‑76 (MQ=255)                                    tCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGtttt  >  2:832294/1‑76 (MQ=39)                                    tCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGtttt  >  2:832267/1‑76 (MQ=39)                                    tCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGtttt  <  1:929848/76‑1 (MQ=39)                                    tCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGtttt  <  1:929833/76‑1 (MQ=39)                                                                          |                                      GTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATCTTTCCCTGCTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCTGAAGCTGTTTT  >  REL606/718824‑718934

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