New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2038432 =65 (1.110)16 (0.310) 6/134 5.4 21.6% coding (517/822 nt) wbbB putative glycosyltransferase WbbB
?REL606 2038941 = 58 (1.110)coding (8/822 nt) wbbB putative glycosyltransferase WbbB
Rejected: Coverage evenness skew score above cutoff.

TATAATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCGGGG  >  REL606/2038872‑2039014
                                                                     |                                                                         
tataATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCa                                                                     <  1:748024/76‑1 (MQ=255)
tataATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCa                                                                     <  1:748040/76‑1 (MQ=255)
tataACTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCa                                                                     <  1:748037/76‑1 (MQ=255)
 ataGTTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCAt                                                                    <  1:160380/76‑1 (MQ=255)
 ataATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCAt                                                                    <  1:160378/76‑1 (MQ=255)
 ataATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCAt                                                                    <  1:607580/76‑1 (MQ=255)
 ataATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCAt                                                                    <  1:607586/76‑1 (MQ=255)
 ataATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCAt                                                                    <  1:607599/76‑1 (MQ=255)
   aaTTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATta                                                                  <  1:864748/76‑1 (MQ=255)
   aaTTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATta                                                                  <  1:1441832/76‑1 (MQ=255)
   aaTTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATta                                                                  <  1:864759/76‑1 (MQ=255)
    aTTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATtat                                                                 <  1:1806656/76‑1 (MQ=255)
         ttaattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatc                                                            >  1:915657/1‑76 (MQ=255)
           aattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatcat                                                          <  1:569525/76‑1 (MQ=255)
           aattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatcat                                                          >  1:951566/1‑76 (MQ=255)
           aattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatcat                                                          <  1:569543/76‑1 (MQ=255)
            attaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatcatc                                                         <  1:1697913/76‑1 (MQ=255)
            attaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatcatc                                                         <  1:1697929/76‑1 (MQ=255)
                 gAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAat                                                    <  1:982572/76‑1 (MQ=255)
                  aaTTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAatt                                                   <  1:907281/76‑1 (MQ=255)
                    ttCTGGCTTTAAAGTATAAGAATAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAaatgg                                                 <  2:261963/76‑2 (MQ=255)
                    ttCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAattga                                                 <  1:1774172/76‑1 (MQ=255)
                    ttCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAattga                                                 <  1:1774162/76‑1 (MQ=255)
                    ttCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAaatgg                                                 <  2:261947/76‑2 (MQ=255)
                    ttCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAaatgg                                                 <  2:261962/76‑2 (MQ=255)
                     tCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAattgat                                                <  1:1098758/76‑1 (MQ=255)
                     tCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAattgat                                                <  1:1098765/76‑1 (MQ=255)
                      cTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAattgatt                                               >  1:1197873/1‑76 (MQ=255)
                      cTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAattgatt                                               >  1:1197866/1‑76 (MQ=255)
                       tGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAattgattg                                              <  1:31674/76‑1 (MQ=255)
                       tGGCTTTAAAGTATAAGAACACGCCACAATGCATAATTTTTGTTTACAATCCATTATAAATCATCCCAattgattg                                              <  1:31700/76‑1 (MQ=255)
                           tttAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAgg                                          <  1:1235861/76‑1 (MQ=255)
                           tttAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAgg                                          <  1:1455889/76‑1 (MQ=255)
                           tttAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAgg                                          <  1:1455866/76‑1 (MQ=255)
                             tAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTg                                        <  1:954148/76‑1 (MQ=255)
                             tAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTg                                        <  1:271589/76‑1 (MQ=255)
                             tAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTg                                        <  1:950927/76‑1 (MQ=255)
                             tAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTg                                        <  1:950913/76‑1 (MQ=255)
                             tAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTg                                        <  1:950902/76‑1 (MQ=255)
                               aaGTATAAGAACAGGCCCCAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGAt                                      <  1:1129995/76‑1 (MQ=255)
                               aaGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGAt                                      <  1:1129997/76‑1 (MQ=255)
                                aGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATa                                     <  1:701535/76‑1 (MQ=255)
                                  tataAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAAc                                   >  1:922528/1‑76 (MQ=255)
                                       gAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAaaga                              <  1:478802/76‑1 (MQ=255)
                                        aagaGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAaagaa                             <  1:92372/73‑1 (MQ=255)
                                          cAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGc                           <  1:647171/76‑1 (MQ=255)
                                          cAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGc                           <  1:1200836/76‑1 (MQ=255)
                                            ggCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCa                         <  1:1310232/76‑1 (MQ=255)
                                              ccACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATa                       >  1:1553165/1‑76 (MQ=255)
                                              ccACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATa                       >  1:1553144/1‑76 (MQ=255)
                                              ccACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATa                       >  1:1553135/1‑76 (MQ=255)
                                                acaATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATaat                     <  2:110360/76‑3 (MQ=255)
                                                acaATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACa                     <  1:943805/76‑1 (MQ=255)
                                                acaATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACa                     <  1:943807/76‑1 (MQ=255)
                                                 caATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACaa                    <  1:1195278/76‑1 (MQ=255)
                                                 caATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACaa                    <  1:147186/76‑1 (MQ=255)
                                                 caATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACaa                    <  1:1528428/76‑1 (MQ=255)
                                                           tttttttttACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGaaa          <  1:1027558/76‑1 (MQ=255)
                                                             tttttttACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAaga        <  1:1562446/76‑1 (MQ=255)
                                                              ttttttACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAagag       <  1:1153060/76‑1 (MQ=255)
                                                              ttttttACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAagag       <  1:1153071/76‑1 (MQ=255)
                                                               tttttACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGc      <  1:1452104/76‑1 (MQ=255)
                                                                  ttACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCggg   >  1:1800376/1‑76 (MQ=255)
                                                                   tACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCgggg  >  1:433449/1‑76 (MQ=255)
                                                                   tACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCgggg  >  1:433452/1‑76 (MQ=255)
                                                                   tACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCgggg  >  1:433465/1‑76 (MQ=255)
                                                                     |                                                                         
TATAATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCGGGG  >  REL606/2038872‑2039014

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