New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2039076 =44 (0.750)11 (0.210) 7/136 4.5 16.9% coding (619/747 nt) wbbA putative glycosyltransferase WbbA
?REL606 2039562 = 68 (1.280)coding (133/747 nt) wbbA putative glycosyltransferase WbbA
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGAAAGAGCGGGGATAGGTGAAAGAAGTCTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTG  >  REL606/2039002‑2039150
                                                                          |                                                                          
tGAAAGAGCGGGGATAGGTGAAAGAAGTCTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTaa                                                                           <  2:1400548/76‑1 (MQ=255)
tGAAAGAGCGGGGATAGGTGAAAGAAGTCTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTaa                                                                           >  1:380636/1‑76 (MQ=255)
  aaaGAGCGGGGATAGGTGAAAGAAGTCTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAAt                                                                         <  1:133935/76‑1 (MQ=255)
  aaaGAGCGGGGATAGGTGAAAGAAGTCTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAAt                                                                         >  2:1438452/1‑76 (MQ=255)
    agagCGGGGATAGGTGAAAGAAGTCTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAg                                                                       >  2:848638/1‑76 (MQ=255)
               aGGTGAAAGAAGTCTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGa                                                            >  2:1082075/1‑76 (MQ=255)
                ggTGAAAGAAGTCTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGaa                                                           <  2:995487/76‑1 (MQ=255)
                  tGAAAGAAGTCTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAatt                                                         >  1:205659/1‑76 (MQ=255)
                   gAAAGAAGTCTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAatta                                                        >  2:367025/1‑76 (MQ=255)
                   gAAAGAAGTCTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAatta                                                        >  1:327444/1‑76 (MQ=255)
                   gAAAGAAGTCTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAatta                                                        >  2:367027/1‑76 (MQ=255)
                     aagaagTCTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAattatt                                                      <  2:325546/76‑1 (MQ=255)
                     aagaagTCTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAattatt                                                      >  2:133694/1‑76 (MQ=255)
                     aagaagTCTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAattatt                                                      >  2:635386/1‑76 (MQ=255)
                     aagaagTCTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAattatt                                                      >  2:1049337/1‑76 (MQ=255)
                       gaagTCTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTAtttt                                                    >  1:684629/1‑76 (MQ=255)
                           tCTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTc                                                >  1:1461808/1‑76 (MQ=255)
                            cTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTcc                                               >  2:1057071/1‑76 (MQ=255)
                            cTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTcc                                               >  2:1057091/1‑76 (MQ=255)
                               aTTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTAt                                            >  1:1844077/1‑76 (MQ=255)
                               aTTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTAt                                            <  1:693094/76‑1 (MQ=255)
                                ttCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATc                                           >  1:890199/1‑76 (MQ=255)
                                ttCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATc                                           <  1:402196/76‑1 (MQ=255)
                                   tAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCaaa                                        >  1:801361/1‑76 (MQ=255)
                                   tAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCaaa                                        >  1:1248439/1‑76 (MQ=255)
                                   tAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCaaa                                        >  1:124413/1‑76 (MQ=255)
                                    aaaTATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCaaaa                                       >  2:1566983/1‑76 (MQ=255)
                                     aaTATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCaaaaa                                      <  1:1380706/76‑1 (MQ=255)
                                             agtaagtaTCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCt                              >  2:1052261/1‑76 (MQ=255)
                                             agtaagtaTCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCt                              >  1:264530/1‑76 (MQ=255)
                                              gtaagtaTCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCtt                             <  2:1262143/76‑1 (MQ=255)
                                               taagtaTCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCttt                            <  2:11080/76‑1 (MQ=255)
                                               taagtaTCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCttt                            >  1:54671/1‑76 (MQ=255)
                                                aagtaTCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTc                           <  2:765248/76‑1 (MQ=255)
                                                 agtaTCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCa                          <  1:1366068/76‑1 (MQ=255)
                                                 agtaTCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCa                          <  2:762358/76‑1 (MQ=255)
                                                 agtaTCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCa                          <  2:762359/76‑1 (MQ=255)
                                                      cAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCataata                     <  1:1540943/76‑1 (MQ=255)
                                                      cAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCataata                     >  2:1084231/1‑76 (MQ=255)
                                                        aaaaTAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAatata                   <  2:1819415/76‑1 (MQ=255)
                                                         aaaTAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAatatat                  >  2:1831314/1‑76 (MQ=255)
                                                                ttGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGcac           >  1:1321893/1‑76 (MQ=255)
                                                                    aaTCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGt       <  2:202374/76‑1 (MQ=255)
                                                                         tAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTg  >  1:632529/1‑76 (MQ=255)
                                                                          |                                                                          
TGAAAGAGCGGGGATAGGTGAAAGAAGTCTAATTCTAAATATTTTAGTAAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTG  >  REL606/2039002‑2039150

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