New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 2039480	32 (0.670)	6 (0.140)	6/134	4.3	20.5%	coding (215/747 nt)	wbbA	putative glycosyltransferase WbbA
	?	REL606	= 2039680	18 (0.420)					coding (15/747 nt)	wbbA	putative glycosyltransferase WbbA
Rejected: Coverage evenness skew score above cutoff.

AGTTTTACGCATTCACGCCATGCTTTGCGCTGATCTCCTCCATTGATTATGAATATCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATG  >  REL606/2039408‑2039555                                                                         |                                                                            aGTTTTACGCATTCACGCCATGCTTTGCGCTGATCTCCTCCATTGATTATGAATATCTTTTTATTATTCATATCGa                                                                          >  2:543146/1‑76 (MQ=255)             tCACGCCATGCTTTGCGCTGATCTCCTCCATTGATTATGAATATCTTTTTATTATTCATATCGATATTTTCAAAgt                                                              <  1:1271230/76‑1 (MQ=255)                      gCTTTGCGCTGATCTCCTCCATTGATTATGAATATCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGAtata                                                     >  2:38137/1‑76 (MQ=255)                       cTTTGCGCTGATCTCCTCCATTGATTATGAATATCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATATaa                                                    <  2:942998/76‑1 (MQ=255)                         ttGCGCTGATCTCCTCCATTGATTATGAATATCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATATAACg                                                  <  1:600255/76‑1 (MQ=255)                         ttGCGCTGATCTCCTCCATTGATTATGAATATCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATATAACg                                                  <  1:1381702/76‑1 (MQ=255)                         ttGCGCTGATCTCCTCCATTGATTATGAATATCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATATAACg                                                  <  1:1381691/76‑1 (MQ=255)                               tGATCTCCTCCATTGATTATGAATATCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTa                                            <  2:871368/76‑1 (MQ=255)                               tGATCTCCTCCATTGATTATGAATATCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTa                                            >  2:75687/1‑76 (MQ=255)                                 aTCTCCTCCATTGATTATGAATATCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAg                                          <  1:1262771/76‑1 (MQ=255)                                     cctccATTGATTATGAATATCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAg                                      <  1:734926/76‑1 (MQ=255)                                     cctccATTGATTATGAATATCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAg                                      <  1:1184493/76‑1 (MQ=255)                                        ccATTGATTATGAATATCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCaa                                   <  1:288138/76‑1 (MQ=255)                                          aTTGATTATGAATATCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACg                                 <  1:8503/76‑1 (MQ=255)                                          aTTGATTATGAATATCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATACAACGGATCTAAGGAACCAACg                                 <  1:8525/76‑1 (MQ=255)                                            tGATTATGAATATCTTTTTATTATACATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGtt                               >  1:616691/1‑76 (MQ=255)                                              aTTATGAATATCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTa                             >  1:1305723/1‑76 (MQ=255)                                               ttATGAATATCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTaa                            >  1:873295/1‑76 (MQ=255)                                                 aTGAATATCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATg                          >  1:324964/1‑76 (MQ=255)                                                  tGAATATCTTTTTATTACTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGa                         >  1:308449/1‑76 (MQ=255)                                                     atatCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATcc                      >  2:266597/1‑76 (MQ=255)                                                      tatCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCa                     >  1:704762/1‑76 (MQ=255)                                                       atCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCAt                    >  1:963314/1‑76 (MQ=255)                                                          tttttATTATTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCaa                 >  1:54456/1‑76 (MQ=255)                                                           ttttATTATTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCaaa                <  1:845636/76‑1 (MQ=255)                                                              tattatTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAAt             <  1:1037523/76‑1 (MQ=255)                                                                  attCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATaaaa         <  1:466474/76‑1 (MQ=255)                                                                    tCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGt       <  1:1136245/76‑1 (MQ=255)                                                                      atatCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTaa     <  1:187900/76‑1 (MQ=255)                                                                       tatCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTaaa    <  1:166697/76‑1 (MQ=255)                                                                        atCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAAt   <  1:861106/76‑1 (MQ=255)                                                                         tCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATg  <  1:1127290/76‑1 (MQ=255)                                                                         |                                                                            AGTTTTACGCATTCACGCCATGCTTTGCGCTGATCTCCTCCATTGATTATGAATATCTTTTTATTATTCATATCGATATTTTCAAAGTGTTGATATAACGGATCTAAGGAAGCAACGTTTAATGATCCATCAAAAATAAAAGTAAATG  >  REL606/2039408‑2039555

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