New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2038934 =	44 (0.950)	8 (0.190)	4/136	4.6	15.9%	coding (15/822 nt)	wbbB	putative glycosyltransferase WbbB
	?	REL606	2039357 =	45 (1.070)					coding (338/747 nt)	wbbA	putative glycosyltransferase WbbA
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTCCAACTCGACGATCTTCATCTACTAGTTTTACGCATTCACGCCAT  >  REL606/2039282‑2039428                                                                            |                                                                        cTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAat                                                                         <  1:550382/76‑1 (MQ=255) cTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAat                                                                         <  1:487540/76‑1 (MQ=255) cTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAat                                                                         >  1:243908/1‑76 (MQ=255) cTTATATATGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAat                                                                         >  1:963062/1‑76 (MQ=255)  ttATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAata                                                                        <  2:1458580/76‑1 (MQ=255)    atatTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAatata                                                                      <  2:105192/76‑1 (MQ=255)    atatTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAatata                                                                      >  1:263936/1‑76 (MQ=255)     tatTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAatatat                                                                     <  1:589319/76‑1 (MQ=255)     tatTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAatatat                                                                     <  1:589313/76‑1 (MQ=255)       tttgatttgaCCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTa                                                                   <  1:1052230/76‑1 (MQ=255)        ttgatttgaCCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTaa                                                                  <  2:911513/76‑1 (MQ=255)         tgatttgaCCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTaaa                                                                 <  1:102520/76‑1 (MQ=255)         tgatttgaCCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTaaa                                                                 <  1:1402288/76‑1 (MQ=255)          gatttgaCCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAAt                                                                <  1:1461865/76‑1 (MQ=255)            tttgaCCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCa                                                              <  2:648681/76‑1 (MQ=255)              tgaCCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCAtt                                                            >  1:637014/1‑76 (MQ=255)                  cAAGTTAAATATCTCATCAATCCATTTATACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTa                                                        <  1:1203597/76‑1 (MQ=255)                  cAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTa                                                        <  1:986256/76‑1 (MQ=255)                  cAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTa                                                        <  1:986241/76‑1 (MQ=255)                  cAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTa                                                        <  1:1203577/76‑1 (MQ=255)                   aaGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTAt                                                       <  1:1455914/76‑1 (MQ=255)                    aGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATc                                                      <  2:5953/76‑1 (MQ=255)                     gTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCa                                                     <  2:997420/76‑1 (MQ=255)                      ttAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTANATCATTTTTATCaa                                                    <  1:813604/76‑1 (MQ=255)                      ttAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCaa                                                    >  1:1331538/1‑76 (MQ=255)                      ttAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCaa                                                    >  1:1300991/1‑76 (MQ=255)                      accaaTATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCaa                                                    >  2:280646/4‑76 (MQ=255)                      accaaTATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCaa                                                    >  2:280631/4‑76 (MQ=255)                        aaaTATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAAtt                                                  <  1:1372213/76‑1 (MQ=255)                            atCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTc                                              >  1:1115939/1‑76 (MQ=255)                              ctcATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTCCa                                            >  1:968665/1‑76 (MQ=255)                                 atcaATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTCCAACt                                         <  1:1024239/76‑1 (MQ=255)                                     aTCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTCCAACTcgac                                     >  1:936315/1‑76 (MQ=255)                                       ccATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTCCAACTcgacga                                   >  2:484957/1‑76 (MQ=255)                                            tAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTCCAACTCGACGATCTTc                              >  1:591390/1‑76 (MQ=255)                                               aCTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTCCAACTCGACGATCTTCATc                           >  1:1166281/1‑76 (MQ=255)                                               aCTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTCCAACGCGACGATCTTCATc                           <  1:357081/76‑1 (MQ=255)                                                 tCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTCCAACTCGACGATCTTCATcta                         >  1:994106/1‑76 (MQ=255)                                                    tGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTCCAACTCGACGATCTTCATctacta                      <  2:746464/76‑1 (MQ=255)                                                    tGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTCCAACTCGACGATCTTCATctacta                      <  2:746463/76‑1 (MQ=255)                                                    tGCAAATAATCATTTTCTAAAAAATATATTAAACCATTTTTATCAATTTTTCCAACTCGACGATCTTCATctacta                      <  2:746461/76‑1 (MQ=255)                                                              cATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTCCAACTCGACGATCTTCATCTACTAGTTTTACGCa            >  2:138420/1‑76 (MQ=255)                                                               aTTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTCCAACTCGACGATCTTCATCTACTAGTTTTACGCAt           >  1:1350698/1‑76 (MQ=255)                                                                    cTAAAAAATATATTAAATCATTTTTATCAATTTTTCCAACTCGACGATCTTCATCTACTAGTTTTACGCATTCACg      >  1:1397995/1‑76 (MQ=255)                                                                    cTAAAAAATATATTAAATCATTTTTATCAATTTTTCCAACTCGACGATCTTCATCTACTAGTTTTACGCATTCACg      >  1:1398016/1‑76 (MQ=255)                                                                     tAAAAAATATATTAAATCATTTTTATCAATTTTTCCAACTCGACGATCTTCATCTACTAGTTTTACGCATTCACGc     <  1:821477/76‑1 (MQ=255)                                                                        aaaaTATATTAAATCATTTTTATCAATTTTTCCAACTCGACGATCTTCATCTACTAGTTTTACGCATTCACGCCAt  >  1:181043/1‑76 (MQ=255)                                                                            |                                                                        CTTATATTTGATTTGACCAAGTTAAATATCTCATCAATCCATTTAGACTCATGCAAATAATCATTTTCTAAAAAATATATTAAATCATTTTTATCAATTTTTCCAACTCGACGATCTTCATCTACTAGTTTTACGCATTCACGCCAT  >  REL606/2039282‑2039428

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