New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 2038833	30 (0.550)	7 (0.140)	6/134	5.3	21.6%	coding (116/822 nt)	wbbB	putative glycosyltransferase WbbB
	?	REL606	= 2039157	24 (0.490)					coding (538/747 nt)	wbbA	putative glycosyltransferase WbbA
Rejected: Coverage evenness skew score above cutoff.

CAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAATTTATTAATTATGAATTCTGGCTTTAAAGTA  >  REL606/2038759‑2038907                                                                           |                                                                           cAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGAAAt                                                                           <  1:1671545/76‑1 (MQ=255)  aaaaaGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGAAATc                                                                          <  1:832520/76‑1 (MQ=255)   aaaaGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGAAATca                                                                         <  1:976243/76‑1 (MQ=255)       gagaATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGAAATCacaac                                                                     <  2:496085/76‑1 (MQ=255)           aTTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGAAATCACAACACCa                                                                 >  1:1335624/1‑76 (MQ=255)                 cccGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATCGGAAATCACAACACCATCTAGt                                                           >  2:1105220/1‑76 (MQ=255)                           atatGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGAAATCACAACACCATCTAGTTCATAATCTg                                                 <  1:1111585/76‑1 (MQ=255)                               gAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTaa                                             <  2:1583524/76‑1 (MQ=255)                                          aTTGTAATCTGTTCTACTATATGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATtatata                                  >  1:74368/1‑76 (MQ=255)                                          aTTGTAATCTGTTCTACTATATGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATtatata                                  >  1:74377/1‑76 (MQ=255)                                           ttGTAATCTGTTCTACTATATGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATaa                                 <  1:106465/76‑1 (MQ=255)                                            tGTAATCTGTTCTACTATATGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAAt                                >  1:1291395/1‑76 (MQ=255)                                            tGTAATCTGTTCTACTATATGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAAt                                >  1:1478939/1‑76 (MQ=255)                                                   tGTTCTACTATATGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAATTTattaa                         >  1:147757/1‑76 (MQ=255)                                                   tGTTCTACTATATGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAATTTattaa                         >  1:993219/1‑76 (MQ=255)                                                   tGTTCTACTATATGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAATTTattaa                         >  1:993217/1‑76 (MQ=255)                                                     ttCTACTATATGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAATTTattaatt                       >  1:162009/1‑76 (MQ=255)                                                     ttCTACTATATGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAATTTattaatt                       >  1:851960/1‑76 (MQ=255)                                                      tCTACTATATGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAATTTattaatta                      >  1:1076733/1‑76 (MQ=255)                                                       ctactaTATGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAATTTATTAATTAt                     >  1:867569/1‑76 (MQ=255)                                                       ctactaTATGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAATTTATTAATTAt                     >  1:902923/1‑76 (MQ=255)                                                        tactaTATGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAATTTATTAATTATg                    >  1:1565828/1‑76 (MQ=255)                                                         actaTATGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAATTTATTAATTATGa                   >  1:1725203/1‑76 (MQ=255)                                                           tataTGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAATTTATTAATTATGAAt                 >  1:207295/1‑76 (MQ=255)                                                            atatGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAATTTATTAATTATGAAtt                >  1:822238/1‑76 (MQ=255)                                                              atGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAATTTATTAATTATGAATTCt              >  1:471871/1‑76 (MQ=255)                                                                     ttGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAATTTATTAATTATGAATTCTGGCTTTa       >  1:1501927/1‑76 (MQ=255)                                                                     ttGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAATTTATTAATTATGAATTCTGGCTTTa       >  1:1501938/1‑76 (MQ=255)                                                                       ggAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAATTTATTAATTATGAATTTTGGCTTTaaa     >  2:478911/1‑76 (MQ=255)                                                                          aaTCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAATTTATTAATTATGAATTCTGGCTTTAAAGta  <  1:751394/76‑1 (MQ=255)                                                                           |                                                                           CAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGAAATCACAACACCATCTAGTTCATAATCTGGTAAATAATTATATAATTTATTAATTATGAATTCTGGCTTTAAAGTA  >  REL606/2038759‑2038907

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