New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2038688 =16 (0.310)5 (0.100) 3/138 6.8 25.8% coding (261/822 nt) wbbB putative glycosyltransferase WbbB
?REL606 2038966 = 14 (0.290)coding (729/747 nt) wbbA putative glycosyltransferase WbbA
Rejected: Coverage evenness skew score above cutoff.

TGAATAATAACTCACAATTTTCCTTAAAAAGAATTTTGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAA  >  REL606/2038613‑2038762
                                                                           |                                                                          
tgAATAATAACTCACAATTTTCCTTAAAAAGAATTTTGCATTATGTTTCAAGAACAATGTATCATTCAAAATgaga                                                                            >  1:1338177/1‑76 (MQ=255)
    taataaCTCACAATTTTCCTTAAAAAGAATTTTGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTg                                                                        <  1:233220/76‑1 (MQ=255)
      ataaCTCACAATTTTCCTTAAAAAGAATTTTGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGt                                                                      <  2:674237/76‑1 (MQ=255)
       taaCTCACAATTTTCCTTAAAAAGAATTTTGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTa                                                                     >  1:741939/1‑76 (MQ=255)
                           aaaGAATTTTGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAg                                                 <  2:269364/76‑1 (MQ=255)
                            aaGAATTTTGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGc                                                <  2:1339760/76‑1 (MQ=255)
                             aGAATTTTGTATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCt                                               >  1:927876/1‑76 (MQ=255)
                                      cATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTaa                                      <  1:1632647/76‑1 (MQ=255)
                                         tATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTaaaaa                                   >  1:942555/1‑76 (MQ=255)
                                         tATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTaaaaa                                   >  1:942575/1‑76 (MQ=255)
                                           tGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAAt                                 >  1:1437320/1‑76 (MQ=255)
                                                cAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCtgca                            <  2:130770/76‑4 (MQ=255)
                                                                       tGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTc     >  1:592188/1‑76 (MQ=255)
                                                                       tGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTc     >  1:592191/1‑76 (MQ=255)
                                                                         agaATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCaa   <  2:560152/76‑1 (MQ=255)
                                                                          gaATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCaaa  <  1:570664/76‑1 (MQ=255)
                                                                           |                                                                          
TGAATAATAACTCACAATTTTCCTTAAAAAGAATTTTGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAA  >  REL606/2038613‑2038762

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