New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2038943 =	57 (0.970)	30 (0.540)	12/142	4.0	30.3%	coding (6/822 nt)	wbbB	putative glycosyltransferase WbbB
	?	REL606	2039511 =	84 (1.520)					coding (184/747 nt)	wbbA	putative glycosyltransferase WbbA
Rejected: Coverage evenness skew score above cutoff.

ATTATATAATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCGGGGATA  >  REL606/2038868‑2039017                                                                            |                                                                           atTATATAATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACaa                                                                            <  1:962725/76‑1 (MQ=255) atTATATAATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACaa                                                                            <  1:671816/76‑1 (MQ=255) atTATATAATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACaa                                                                            <  1:671800/76‑1 (MQ=255) atTATATAATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACaa                                                                            <  1:19812/76‑1 (MQ=255)   tatataATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATc                                                                          >  2:1722820/1‑76 (MQ=255)   tatataATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATc                                                                          >  1:726793/1‑76 (MQ=255)    atataATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATcc                                                                         <  2:854865/76‑1 (MQ=255)    atataATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATcc                                                                         <  1:499106/76‑1 (MQ=255)     tataATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCa                                                                        <  2:549007/76‑1 (MQ=255)     tataATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCa                                                                        <  1:1335879/76‑1 (MQ=255)     tataATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCa                                                                        <  1:1339777/76‑1 (MQ=255)     tataATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCa                                                                        <  1:1339793/76‑1 (MQ=255)     tataATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCa                                                                        >  2:272887/1‑76 (MQ=255)       taATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCAtt                                                                      >  2:539289/1‑76 (MQ=255)         aTTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATtat                                                                    >  2:870556/1‑76 (MQ=255)            tattaATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATaaa                                                                 >  1:699110/1‑76 (MQ=255)             attaattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAat                                                                <  1:524785/76‑1 (MQ=255)                aattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatcat                                                             <  1:271113/76‑1 (MQ=255)                aattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatcat                                                             <  1:767379/76‑1 (MQ=255)                    aTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAGTCATCCCa                                                         <  1:1280207/76‑1 (MQ=255)                      gAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATTCATTATAAATCATCCCAat                                                       >  2:1186142/1‑76 (MQ=255)                      gAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAat                                                       >  1:1810089/1‑76 (MQ=255)                         ttCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAattct                                                    <  2:47828/76‑3 (MQ=255)                            tGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAattgattg                                                 <  1:295187/76‑1 (MQ=255)                             ggCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAAttgattga                                                <  1:126185/76‑1 (MQ=255)                             ggCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAAttgattga                                                >  1:1833215/1‑76 (MQ=255)                               cTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAg                                              >  2:329229/1‑76 (MQ=255)                               cTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAg                                              <  1:1566671/76‑1 (MQ=255)                               cTTTAAAGTATAAAAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAg                                              <  1:1566677/76‑1 (MQ=255)                                    aaGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGAt                                         <  1:644399/76‑1 (MQ=255)                                        ataAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACa                                     <  2:205220/76‑1 (MQ=255)                                         taAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACaa                                    <  1:1780967/76‑1 (MQ=255)                                          aaGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACaaa                                   >  2:1628774/1‑76 (MQ=255)                                          aaGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACaaa                                   >  2:1666140/1‑76 (MQ=255)                                              aCAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAaagaag                               <  1:989189/76‑1 (MQ=255)                                               cAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGc                              >  1:562214/1‑76 (MQ=255)                                               cAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGc                              >  1:309652/1‑76 (MQ=255)                                                aGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGcc                             >  1:1065347/1‑76 (MQ=255)                                                   ccACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATa                          <  2:1554434/76‑1 (MQ=255)                                                      caATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACaa                       >  2:1476659/1‑76 (MQ=255)                                                        aTGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATg                     <  1:1851380/76‑1 (MQ=255)                                                         tGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACACTGa                    >  1:983880/1‑76 (MQ=255)                                                         tGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGa                    >  1:983863/1‑76 (MQ=255)                                                          gCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAg                   >  2:423902/1‑76 (MQ=255)                                                              aaTTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGa               >  1:1356778/1‑76 (MQ=255)                                                               aTTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGaa              <  1:688888/76‑1 (MQ=255)                                                               aTTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGaa              <  1:344825/76‑1 (MQ=255)                                                                 ttttttttACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCGATACAATGAGTTGAAag            >  1:273812/1‑76 (MQ=255)                                                                 ttttttttACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAag            <  2:1803129/76‑1 (MQ=255)                                                                   ttttttACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAagag          <  1:576879/76‑1 (MQ=255)                                                                    tttttACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGc         <  1:409065/76‑1 (MQ=255)                                                                     ttttACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCg        >  2:112001/1‑76 (MQ=255)                                                                     ttttACAATCCATTAAAAATCATCCAAATTGATTGAAGGTGATAAAAAAAAAGCCATACAATGAGTTGAAAGAGCg        <  1:843712/76‑1 (MQ=255)                                                                      tttACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCgg       >  1:577425/1‑76 (MQ=255)                                                                      tttACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCgg       >  1:577423/1‑76 (MQ=255)                                                                       ttACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCggg      >  1:882473/1‑76 (MQ=255)                                                                           aaTCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCGGGGATa  >  1:1685896/1‑76 (MQ=255)                                                                            |                                                                           ATTATATAATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCGGGGATA  >  REL606/2038868‑2039017

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