New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4341778 =	16 (0.800)	5 (0.350)	3/96	1.4	39.1%	intergenic (‑150/‑133)	melR/melA	melibiose operon transcriptional regulator; autoregulator/alpha‑galactosidase, NAD(P)‑binding
	?	NC_000913	4341834 =	4 (0.280)					intergenic (‑206/‑77)	melR/melA	melibiose operon transcriptional regulator; autoregulator/alpha‑galactosidase, NAD(P)‑binding

CCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGATTTACTGCTGCTTCACGCAGGATCTGAGTTTATGGGAATGCTCAACCTGGA  >  NC_000913/4341714‑4341845                                                                 |                                                                    cccAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTc                                                                  >  2:163851/1‑68 (MQ=255) cccAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAAATTATTc                                                                  >  1:298039/1‑68 (MQ=255)    gcttgaataacttcatcatggcaggcgaatctgcagaaaacctccggcttccaggttgagcATTCCCa                                                               >  2:504268‑M1/62‑68 (MQ=255)    gcttgaataacttcatcatggcaggcgaatctgcagaaaacctccggcttccaggttgagcATTCCCa                                                               >  1:193468‑M1/62‑68 (MQ=255)       ggggCGAGGGGACTCTCCGAGTATCATGAGGCCGAAGACGCTGCTTTTCAGGTAATTTATTCCCATaa                                                            >  1:216094/1‑68 (MQ=255)        gggCGAGGTGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATaaa                                                           <  2:497729/68‑1 (MQ=255)        gggCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATaaa                                                           <  1:179444/68‑1 (MQ=255)              ggggACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGa                                                     <  2:338212/68‑1 (MQ=255)              ggggACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGa                                                     <  2:327934/68‑1 (MQ=255)               ttcatcatggcaggcgaatctgcagaaaacctccggcttccaggttgagcATTCCCATAAACTCAGAt                                                    <  2:216191‑M1/18‑1 (MQ=255)               gggACTCTCCGAGTATCATGAGGCCGAATACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGAt                                                    >  2:200860/1‑68 (MQ=255)                 gACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGAttt                                                  >  1:260287/1‑68 (MQ=255)                    tctcCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGATTTAct                                               <  2:236126/68‑1 (MQ=255)                            tATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGATTTACTGCTGCTTc                                       >  2:435142/1‑68 (MQ=255)                               aatctgctgaatacctccggcttccatgttgagcATTCCCATAATCTCAGATTTACTGCTGcttctcg                                    >  2:135564‑M1/35‑65 (MQ=255)                               aatctgcagaaaacctccggcttccaggttgagcATTCCCATAAACTCAGATTTACTGCTGCTTCACg                                    >  1:670205‑M1/35‑68 (MQ=255)                                atctgcagaaaacctccggcttccaggttgagcATTCCCATAAACTCAGATTTACTGCTGCTTCACg                                    <  2:293381‑M1/34‑1 (MQ=255)                                atctgcagaaaacctccggcttccaggttgagcATTCCCATAAACTCAGATTTACTGCTGCTTCACg                                    >  1:293381‑M1/34‑67 (MQ=255)                                 tGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGATTTACTGCTGCTTCACGCa                                  <  2:310518/68‑1 (MQ=255)                                  ctgcagaaaacctccggcttccaggttgagcATTCCCATAAACTCAGATTTACTGCTGCTTCACGCAg                                 <  1:135564‑M1/37‑1 (MQ=255)                                   aGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGATTTACTGCTGCTTCACGCAgg                                <  1:107526/68‑1 (MQ=255)                                    ggCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGATTTACTGCTGCTTCACGCAGGa                               <  1:347681/68‑1 (MQ=255)                                           aaCTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGATTTACTGCTGCTTCACGCAGGATCTGAGt                        <  1:96838/68‑1 (MQ=255)                                                               ttATTCCCATAAACTCAGATTTACTGCTGCTTCACGCAGGATCTGAGTTTATGGGAATGCTCAACCTg    >  1:544178/1‑68 (MQ=255)                                                                 ATTCCCATAAACTCAGATTTACTGCTGCTTCACGCAGGATCTGAGTTTATGGGAATGATCAACCTGGa  <  2:39648‑M1/68‑1 (MQ=255)                                                                 |                                                                    CCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGATTTACTGCTGCTTCACGCAGGATCTGAGTTTATGGGAATGCTCAACCTGGA  >  NC_000913/4341714‑4341845

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