New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2038724 =	51 (0.820)	12 (0.220)	6/134	5.8	20.3%	coding (225/822 nt)	wbbB	putative glycosyltransferase WbbB
	?	REL606	2039214 =	49 (0.890)					coding (481/747 nt)	wbbA	putative glycosyltransferase WbbA
Rejected: Coverage evenness skew score above cutoff.

ATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTTATAT  >  REL606/2039142‑2039288                                                                         |                                                                           aTGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGa                                                                         <  1:1391620/76‑1 (MQ=255)  ngTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGaa                                                                        >  1:1765466/2‑76 (MQ=255)       ttGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAg                                                                   <  1:1922509/76‑1 (MQ=255)       ttGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAg                                                                   <  1:1922490/76‑1 (MQ=255)        tGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGa                                                                  <  1:478378/76‑1 (MQ=255)        tGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGa                                                                  <  1:750723/76‑1 (MQ=255)         tgTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGAt                                                                 <  1:433021/75‑1 (MQ=255)         ggTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGAt                                                                 <  1:433004/76‑1 (MQ=255)         ggTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGAt                                                                 <  1:433010/76‑1 (MQ=255)          gTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATg                                                                <  1:341154/76‑1 (MQ=255)           tGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGt                                                               <  1:42390/76‑1 (MQ=255)               aTTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTcac                                                           <  1:1080478/76‑1 (MQ=255)                  ttCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAat                                                        <  1:1124059/76‑1 (MQ=255)                  ttCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAat                                                        <  1:179520/76‑1 (MQ=255)                  ttCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAat                                                        <  1:1124051/76‑1 (MQ=255)                  ttCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAat                                                        <  1:1124032/76‑1 (MQ=255)                            tACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACtata                                              <  1:1246964/76‑1 (MQ=255)                            tACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACtata                                              <  1:1246956/76‑1 (MQ=255)                            tACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTACACAATATTCACtata                                              <  1:1134637/76‑1 (MQ=255)                               tCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATAttt                                           >  1:818365/1‑76 (MQ=255)                                        aaaaTACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGgatg                                  >  1:1329851/1‑76 (MQ=255)                                           ttACAATAGTTTTTTTGTTTTTTTTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATc                               >  2:1716966/2‑76 (MQ=255)                                           aTACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATc                               >  1:1716966/1‑76 (MQ=255)                                            tACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCa                              <  1:1273855/76‑1 (MQ=255)                                            tACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCa                              <  1:899701/76‑1 (MQ=255)                                             aCAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCAt                             >  1:854803/1‑76 (MQ=255)                                             aCAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCAt                             >  2:836387/1‑76 (MQ=255)                                             aCAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCAt                             >  1:1700502/1‑76 (MQ=255)                                             aCAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCAt                             >  1:1700486/1‑76 (MQ=255)                                               aaTAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATaa                           <  2:543851/76‑1 (MQ=255)                                               aaTAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATaa                           <  1:1448605/76‑1 (MQ=255)                                                aTAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATaan                          <  1:719825/76‑2 (MQ=255)                                                aTAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATaaa                          >  1:260610/1‑76 (MQ=255)                                                aTAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATaaa                          >  2:759381/1‑76 (MQ=255)                                                aTAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATaaa                          >  2:759393/1‑76 (MQ=255)                                                  aGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGc                        >  1:729322/1‑76 (MQ=255)                                                  aGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGc                        <  1:1768429/76‑1 (MQ=255)                                                  aGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGc                        >  1:121876/1‑76 (MQ=255)                                                  aGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGc                        >  1:841846/1‑76 (MQ=255)                                                    tttttttGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGt                      <  1:660112/76‑1 (MQ=255)                                                    tttttttGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGt                      <  1:660124/76‑1 (MQ=255)                                                     ttttttGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGtt                     >  1:1679816/1‑76 (MQ=255)                                                      tttttGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTg                    <  2:67926/76‑1 (MQ=255)                                                       ttttGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGc                   >  1:1953037/1‑76 (MQ=255)                                                        tttGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATACACTATATTTATCTGGATGATCATAAAGCGTTGcc                  >  1:236591/1‑76 (MQ=255)                                                          tGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCat                >  1:784171/1‑76 (MQ=255)                                                          tGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCat                >  1:784180/1‑76 (MQ=255)                                                           gTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCata               <  2:1192309/76‑1 (MQ=255)                                                                tATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCa          <  1:1111234/76‑1 (MQ=255)                                                                  ttttttAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATc        >  1:1402020/1‑76 (MQ=255)                                                                    ttttAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCtt      <  1:792602/76‑1 (MQ=255)                                                                     tttAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTta     <  1:1205100/76‑1 (MQ=255)                                                                      ttAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTtat    >  2:536818/1‑76 (MQ=255)                                                                      ttAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTtat    >  2:588086/1‑76 (MQ=255)                                                                      ttAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTtat    <  2:626206/76‑1 (MQ=255)                                                                      ttAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTtat    >  2:70620/1‑76 (MQ=255)                                                                        aaTGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTTatat  >  1:1159677/1‑76 (MQ=255)                                                                        aaTGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTTatat  <  1:1192579/76‑1 (MQ=255)                                                                         |                                                                           ATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTTTGTTTTTATTTTTTAATGAATCAAGATGTTCACAATATTCACTATATTTATCTGGATGATCATAAAGCGTTGCCATATCCCATCTTATAT  >  REL606/2039142‑2039288

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