New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2038432 =51 (1.030)11 (0.250) 6/134 4.6 22.1% coding (517/822 nt) wbbB putative glycosyltransferase WbbB
?REL606 2038941 = 32 (0.720)coding (8/822 nt) wbbB putative glycosyltransferase WbbB
Rejected: Coverage evenness skew score above cutoff.

CTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCATACAAA  >  REL606/2038359‑2038506
                                                                         |                                                                          
cTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTaa                                                                          >  1:1457133/1‑76 (MQ=255)
cTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTaa                                                                          <  1:69047/76‑1 (MQ=255)
cTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTaa                                                                          <  1:69062/76‑1 (MQ=255)
 tACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTaaa                                                                         <  1:1004623/76‑1 (MQ=255)
 tACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTaaa                                                                         <  2:1548304/76‑1 (MQ=255)
        aaaTTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAatat                                                                  <  1:562017/76‑1 (MQ=255)
        aaaTTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAatat                                                                  <  1:562033/76‑1 (MQ=255)
             ctctcAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTgg                                                             <  1:1348297/76‑1 (MQ=255)
                tcAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGaan                                                          <  1:11533/76‑2 (MQ=255)
                tcAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAg                                                          >  1:236934/1‑76 (MQ=255)
                 cAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGc                                                         <  1:833676/76‑1 (MQ=255)
                 cAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGc                                                         <  1:643407/76‑1 (MQ=255)
                 cAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGc                                                         <  1:1079275/76‑1 (MQ=255)
                 cAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGc                                                         <  1:1079269/76‑1 (MQ=255)
                 cAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGc                                                         >  1:538422/1‑76 (MQ=255)
                   gTCGTGATTCAATCTGCATTCCCCATGTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTc                                                       <  1:522984/76‑1 (MQ=255)
                          ttCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAg                                                <  2:1587916/76‑1 (MQ=255)
                                cTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTa                                          >  1:148664/1‑76 (MQ=255)
                                 tGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTaa                                         <  1:836224/76‑1 (MQ=255)
                                 tGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTaa                                         <  2:557897/76‑1 (MQ=255)
                                   cATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACa                                       <  1:181719/76‑1 (MQ=255)
                                   cATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACa                                       <  1:1100129/76‑1 (MQ=255)
                                   cATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACa                                       <  1:61822/76‑1 (MQ=255)
                                      tCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAatt                                    <  1:1493210/76‑1 (MQ=255)
                                          cATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTAttt                                <  1:114942/76‑1 (MQ=255)
                                           aTTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTAtttt                               <  2:146937/76‑1 (MQ=255)
                                           aTTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTAtttt                               >  1:6328/1‑76 (MQ=255)
                                               ggggTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCa                           <  2:1116518/76‑1 (MQ=255)
                                                 ggTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGa                         <  1:152724/76‑1 (MQ=255)
                                                 ggTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGa                         >  1:1112864/1‑76 (MQ=255)
                                                        aatatcaatatcATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCa                  >  1:359096/1‑76 (MQ=255)
                                                         atatcaatatcaTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCAt                 >  1:1483437/1‑76 (MQ=255)
                                                          tatcaatatcaTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCAtt                >  1:429770/1‑76 (MQ=255)
                                                          tatcaatatcaTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCAtt                >  1:138982/1‑76 (MQ=255)
                                                          tatcaatatcaTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCAtt                >  1:429718/1‑76 (MQ=255)
                                                          tatcaatatcaTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCAtt                >  1:429744/1‑76 (MQ=255)
                                                          tatcaatatcaTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCAtt                >  1:429752/1‑76 (MQ=255)
                                                             caatatcaTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATttat             >  1:194018/1‑76 (MQ=255)
                                                               atatcaTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATttatta           <  1:147741/76‑1 (MQ=255)
                                                                tatcaTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAg          >  1:863475/1‑76 (MQ=255)
                                                                tatcaTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAg          >  1:863488/1‑76 (MQ=255)
                                                                  tcaTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCa        >  1:1525078/1‑76 (MQ=255)
                                                                  tcaTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCa        >  1:1525059/1‑76 (MQ=255)
                                                                  tcaTTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCa        >  1:1525060/1‑76 (MQ=255)
                                                                      tACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTCTTCATTTATTAGCATACa    >  2:1512454/1‑76 (MQ=255)
                                                                        cTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCATACaaa  <  1:1333344/76‑1 (MQ=255)
                                                                        cTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCATACaaa  <  1:1289847/76‑1 (MQ=255)
                                                                        cTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCATACaaa  <  1:1318128/76‑1 (MQ=255)
                                                                        cTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCATACaaa  <  1:527746/76‑1 (MQ=255)
                                                                        cTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCATACaaa  <  1:527775/76‑1 (MQ=255)
                                                                        cTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCATACaaa  <  1:527763/76‑1 (MQ=255)
                                                                         |                                                                          
CTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAGCTCATTGTAGCAACTAACAATTATTTTTGCAGATTTTTCATTTATTAGCATACAAA  >  REL606/2038359‑2038506

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