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.

ATTATGTAACAATTTAGCTTGATACCAAACCGTATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAG  >  REL606/2038301‑2038450
                                                                          |                                                                           
aTTATGTAACAATTTAGCTTGATACCAAACCGTATCTGTGTCTATATCAATTAAATGACTACGAATAAATtctctc                                                                            >  1:1201175/1‑76 (MQ=255)
 ttATGTAACAATTTAGCTTGATACCAAACCGTATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCa                                                                           >  1:1318652/1‑76 (MQ=255)
  tATGTAACAATTTAGCTTGATACCAAACCGTATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAg                                                                          <  2:1495852/76‑1 (MQ=255)
   aTGTAACAATTTAGCTTGATACCAAACCGTATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGt                                                                         <  2:268806/76‑1 (MQ=255)
         cAATTTAGCTTGATACCAAACCGTATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGAt                                                                   <  1:250358/76‑1 (MQ=255)
              tAGCTTGATACCAAACCGTATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAAt                                                              >  1:763086/1‑76 (MQ=255)
                 cTTGATACCAAACCGTATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTg                                                           >  1:1882934/1‑76 (MQ=255)
                 cTTGATACCAAACCGTATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTg                                                           <  1:1979562/76‑1 (MQ=255)
                  ttGATACCAAACCGTATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGc                                                          >  1:960257/1‑76 (MQ=255)
                          aaaCCGTATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCa                                                  >  1:853391/1‑76 (MQ=255)
                            aCCGTATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCAtt                                                <  2:840637/76‑1 (MQ=255)
                            aCCGTATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCAtt                                                >  1:947773/1‑76 (MQ=255)
                            aCCGTATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCAtt                                                >  1:947764/1‑76 (MQ=255)
                             ccGTATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCAttt                                               >  1:443011/1‑76 (MQ=255)
                             ccGTATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCAttt                                               >  1:1092051/1‑76 (MQ=255)
                              cGTATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTg                                              <  2:1357752/76‑1 (MQ=255)
                              cGTATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTg                                              <  2:1087913/76‑1 (MQ=255)
                                tATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTggg                                            <  1:1284480/76‑1 (MQ=255)
                                tATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTggg                                            >  1:1329793/1‑76 (MQ=255)
                                tATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTggg                                            >  1:1329802/1‑76 (MQ=255)
                                 aTCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTgggg                                           <  1:1556349/76‑1 (MQ=255)
                                      tgtCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCaat                                      >  1:1575602/1‑76 (MQ=255)
                                        tCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCaataa                                    <  1:210828/76‑1 (MQ=255)
                                        tCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCaataa                                    <  2:2055484/76‑1 (MQ=255)
                                        tCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCaataa                                    <  1:501461/76‑1 (MQ=255)
                                           atatCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATaatat                                 >  1:1122225/1‑76 (MQ=255)
                                            tatCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATaatatc                                >  2:1133957/1‑76 (MQ=255)
                                              tCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATaatatcaa                              >  1:648698/1‑76 (MQ=255)
                                               cAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATaatatcaat                             >  1:1895887/1‑76 (MQ=255)
                                                aaTTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATaatatcaata                            <  1:1335578/76‑1 (MQ=255)
                                                 aTTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATaatatcaatat                           >  1:374627/1‑76 (MQ=255)
                                                 aTTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATaatatcaatat                           >  1:23779/1‑76 (MQ=255)
                                                   tAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAatatcaatatca                         >  2:1700346/1‑76 (MQ=255)
                                                      aTGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTa                      >  1:1720410/1‑76 (MQ=255)
                                                           tACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTaaa                 <  2:1061067/76‑1 (MQ=255)
                                                               aaTAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTaaaaaaa             <  1:2018292/76‑1 (MQ=255)
                                                                     ttCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATAtttt       >  1:1848414/1‑76 (MQ=255)
                                                                      tctctcAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTg      <  2:1422261/76‑1 (MQ=255)
                                                                       ctctcAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTgg     >  1:820073/1‑76 (MQ=255)
                                                                       ctctcAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTgg     >  1:99070/1‑76 (MQ=255)
                                                                          tcAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAg  >  1:719194/1‑76 (MQ=255)
                                                                          |                                                                           
ATTATGTAACAATTTAGCTTGATACCAAACCGTATCTGTGTCTATATCAATTAAATGACTACGAATAAATTCTCTCAGTCTTGATTCAATCTGCATTCCCCATTTGGGGTCAATAATATCAATATCATTACTAAAAAAATATTTTGGAAG  >  REL606/2038301‑2038450

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