New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 203837541 (0.630)3 (0.050) 3/132 8.0 6.8% coding (574/822 nt) wbbB putative glycosyltransferase WbbB
?REL606 = 2038763 46 (0.800)coding (186/822 nt) wbbB putative glycosyltransferase WbbB
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGAAATCAC  >  REL606/2038696‑2038837
                                                                   |                                                                          
cTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAAtt                                                                    <  1:155994/76‑1 (MQ=255)
 tCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAAttt                                                                   <  2:2050199/76‑1 (MQ=255)
 tCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAAttt                                                                   <  1:1600376/76‑1 (MQ=255)
    tGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTc                                                                <  1:1303917/76‑1 (MQ=255)
      aTGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTccc                                                              <  2:22793/76‑1 (MQ=255)
         gTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGtt                                                           >  1:495397/1‑76 (MQ=255)
         gTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGtt                                                           <  2:2020258/76‑1 (MQ=255)
             taatGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCg                                                       <  2:178069/76‑1 (MQ=255)
               atGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGat                                                     >  1:1975523/1‑76 (MQ=255)
               atGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGat                                                     >  2:1135850/1‑76 (MQ=255)
                tGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGata                                                    >  2:503229/1‑76 (MQ=255)
                  aGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATg                                                  <  1:2013385/76‑1 (MQ=255)
                    cTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGaa                                                >  2:1745199/1‑76 (MQ=255)
                      aTCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAAt                                              <  1:1458128/76‑1 (MQ=255)
                      aTCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAAt                                              >  2:157791/1‑76 (MQ=255)
                      aTCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAAt                                              >  2:1599168/1‑76 (MQ=255)
                      aTCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAAt                                              >  2:1699695/1‑76 (MQ=255)
                      aTCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAAt                                              >  1:1395127/1‑76 (MQ=255)
                       tcttctAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATc                                             <  1:1280224/76‑1 (MQ=255)
                        cttctAAAAAAGCTATTCCTTCCTTATAGGCAAAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATcc                                            <  1:427820/76‑1 (MQ=255)
                         ttctAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATccc                                           <  2:772407/76‑1 (MQ=255)
                           ctAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCtt                                         <  2:122299/76‑1 (MQ=255)
                            tAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTc                                        <  1:808968/76‑1 (MQ=255)
                                 aaTCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGt                                   <  1:78679/76‑1 (MQ=255)
                                  aTCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTa                                  <  1:1994808/76‑1 (MQ=255)
                                  aTCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTa                                  >  1:1845120/1‑76 (MQ=255)
                                    cTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAAt                                >  2:98009/1‑76 (MQ=255)
                                    cTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAAt                                >  1:1910866/1‑76 (MQ=255)
                                     tATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATc                               >  1:1088900/1‑76 (MQ=255)
                                       ttccttccTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTg                             >  1:734782/1‑76 (MQ=255)
                                         ccttccttATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGtt                           >  1:1353840/1‑76 (MQ=255)
                                            tccttATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTcta                        >  2:653508/1‑76 (MQ=255)
                                            tccttATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTcta                        <  2:645426/76‑1 (MQ=255)
                                               ttATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTctacta                     <  2:1798794/76‑1 (MQ=255)
                                                tataGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACtat                    >  2:117354/1‑76 (MQ=255)
                                                tataGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACtat                    >  1:783051/1‑76 (MQ=255)
                                                  taGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACtatan                  <  1:1236347/76‑2 (MQ=255)
                                                   aGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATg                 <  2:1891009/76‑1 (MQ=255)
                                                   aGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATg                 <  1:531020/76‑1 (MQ=255)
                                                             cTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGaa       >  1:1930520/1‑76 (MQ=255)
                                                             cTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGaa       >  1:1807921/1‑76 (MQ=255)
                                                             cTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGaa       <  1:518090/76‑1 (MQ=255)
                                                              tCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGaaa      <  1:205104/76‑1 (MQ=255)
                                                               cAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGAAAt     <  1:918042/76‑1 (MQ=255)
                                                                aaaaaGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGAAATc    <  2:2030077/76‑1 (MQ=255)
                                                                aaaaaGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGAAATc    <  2:58036/76‑1 (MQ=255)
                                                                  aaaGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTTTACTATATGATTATTGGAAATtac  >  1:896427/1‑76 (MQ=255)
                                                                  aaaGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGAAATcac  >  1:452076/1‑76 (MQ=255)
                                                                  aaaGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGAAATcac  >  2:1579945/1‑76 (MQ=255)
                                                                   |                                                                          
CTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTTCATTGTAATCTGTTCTACTATATGATTATTGGAAATCAC  >  REL606/2038696‑2038837

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