New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2038432 =	23 (1.040)	9 (0.460)	4/134	2.2	24.2%	coding (517/822 nt)	wbbB	putative glycosyltransferase WbbB
	?	REL606	2038941 =	36 (1.820)					coding (8/822 nt)	wbbB	putative glycosyltransferase WbbB

GACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCATACAAA  >  REL606/2038357‑2038506                                                                            |                                                                           gACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACt                                                                            >  2:237900/1‑76 (MQ=255)   cTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTaa                                                                          >  1:148505/1‑76 (MQ=255)    ccgctctttcaactcattgtatggcttctttgttatcaccttcaatcaattgggatgatttataatggattgTaaa                                                                         >  1:278100‑M1/73‑76 (MQ=255)     aCGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTaaaa                                                                        <  1:560627/76‑1 (MQ=255)      cGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTaaaaa                                                                       >  1:401972/1‑76 (MQ=255)       gAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTaaaaaa                                                                      <  1:209961/76‑1 (MQ=255)                        cTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCAt                                                     >  1:239333/1‑76 (MQ=255)                            cttCAATCTGCATTCCCCATTTGGGGTCAATAATATCTATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTa                                                 >  2:40840/2‑76 (MQ=255)                                aaTCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCaa                                             <  2:283307/76‑1 (MQ=255)                                         tCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACAAACAatt                                    <  1:460848/76‑1 (MQ=255)                                               tttGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTAttttt                              <  1:164526/76‑1 (MQ=255)                                                atcaattgggatgatttataatggattgTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTg                             <  2:213200‑M1/48‑1 (MQ=255)                                                 tGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGc                            <  1:628834/76‑1 (MQ=255)                                                    attgggatgatttataatggattgTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACCATTATTTTTGCAGa                         <  2:421853‑M1/52‑1 (MQ=255)                                                    attgggatgatttataatggattgTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGa                         <  2:421857‑M1/52‑1 (MQ=255)                                                       cAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGAttt                      >  1:192888/1‑76 (MQ=255)                                                         ataataTCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGAttttt                    >  1:295572/1‑76 (MQ=255)                                                         ataataTCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGAttttt                    >  1:295573/1‑76 (MQ=255)                                                          atgatttataatggattgTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTCTTTTTGCAGATTTTTc                   <  2:507923‑M1/58‑1 (MQ=255)                                                          atgatttataatggattgTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTc                   <  2:507917‑M1/58‑1 (MQ=255)                                                          atgatttataatggattgTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTc                   <  2:507930‑M1/58‑1 (MQ=255)                                                          atgatttataatggattgTAAAAAAATATTTTGGAAGATCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTc                   <  2:382259‑M1/58‑1 (MQ=255)                                                                 aatatcaTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATttatt            >  2:47372/1‑76 (MQ=255)                                                                  atatcaTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATttatta           >  1:182040/1‑76 (MQ=255)                                                                     tcaTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCa        >  1:23536/1‑76 (MQ=255)                                                                     tcaTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCa        >  1:23531/1‑76 (MQ=255)                                                                     tcaTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCa        <  1:115862/76‑1 (MQ=255)                                                                      caTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCAt       <  1:27403/76‑1 (MQ=255)                                                                      caTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTNCATTTATTAGCAt       <  1:485174/76‑1 (MQ=255)                                                                         ttgTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCATACa    <  2:345220‑M1/73‑1 (MQ=255)                                                                         ttgTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCATACa    <  2:345221‑M1/73‑1 (MQ=255)                                                                           cTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCATACaaa  <  2:337887/76‑1 (MQ=255)                                                                           cTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCATACaaa  <  1:279527/76‑1 (MQ=255)                                                                            |                                                                           GACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCATACAAA  >  REL606/2038357‑2038506

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