New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2038432 =23 (1.040)9 (0.460) 4/134 2.2 24.2% coding (517/822 nt) wbbB putative glycosyltransferase WbbB
?REL606 2038941 = 36 (1.820)coding (8/822 nt) wbbB putative glycosyltransferase WbbB

ATTATATAATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCGG  >  REL606/2038868‑2039012
                                                                         |                                                                       
tgtatgctaataaatgaaaaatctgcaaaaataattgttagttgctacaatgagcttccaaaatatttttttaCaa                                                                       <  2:345220‑M2/3‑1 (MQ=255)
tgtatgctaataaatgaaaaatctgcaaaaataattgttagttgctacaatgagcttccaaaatatttttttaCaa                                                                       <  2:345221‑M2/3‑1 (MQ=255)
atTATATAATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACaa                                                                       <  1:345221/76‑1 (MQ=255)
atTATATAATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACaa                                                                       <  1:345220/76‑1 (MQ=255)
     ataATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCAt                                                                  >  1:151640/1‑76 (MQ=255)
       aaTTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATta                                                                <  1:575403/76‑1 (MQ=255)
            attaattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAat                                                           <  2:110030/76‑1 (MQ=255)
             ttaattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatc                                                          <  1:71357/76‑1 (MQ=255)
             aaaattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatc                                                          >  2:361375/3‑76 (MQ=255)
              taattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatca                                                         >  2:90087/1‑76 (MQ=255)
              taattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTACAAGCCATTATAAatcat                                                        <  1:382259/76‑1 (MQ=255)
               gaaaaatctgcaaaaataattgttagttgctacaatgatcttccaaaatatttttttaCAATCCATTATAAatcat                                                        <  2:382259‑M2/18‑1 (MQ=255)
               gaaaaatctgcaaaaataattgttagttgctacaatgagcttccaaaatatttttttaCAATCCATTATAAatcat                                                        <  2:507917‑M2/18‑1 (MQ=255)
               gaaaaatctgcaaaaataattgttagttgctacaatgagcttccaaaatatttttttaCAATCCATTATAAatcat                                                        <  2:507930‑M2/18‑1 (MQ=255)
               gaaaaatctgcaaaaagaattgttagttgctacaatgagcttccaaaatatttttttaCAATCCATTATAAatcat                                                        <  2:507923‑M2/18‑1 (MQ=255)
               aattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatcat                                                        <  1:507923/76‑1 (MQ=255)
               aattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatcat                                                        <  1:507930/76‑1 (MQ=255)
               aattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatcat                                                        <  1:507917/76‑1 (MQ=255)
                 ttaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATcc                                                      <  1:517681/76‑1 (MQ=255)
                   aTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCa                                                    >  1:391629/1‑76 (MQ=255)
                    tGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCaa                                                   >  2:175358/1‑76 (MQ=255)
                     tctgcaaaaataattgttagttgctacaatgagcttccaaaatatttttttaCAATCCATTATAAATCATCCCAat                                                  <  2:421857‑M2/24‑1 (MQ=255)
                     tctgcaaaaataatggttagttgctacaatgagcttccaaaatatttttttaCAATCCATTATAAATCATCCCAat                                                  <  2:421853‑M2/24‑1 (MQ=255)
                     gAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAat                                                  <  1:421853/76‑1 (MQ=255)
                     gAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCAAAATTTTTTTTTACAATCCATTATAAATCATCCCAat                                                  <  1:421857/76‑1 (MQ=255)
                      aaTTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTCAATCCATTATAAATTATCCCCattt                                                >  1:257773/1‑75 (MQ=255)
                         tCTGGCTTTAAAGTGTAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAAGCATCCCAattgat                                              <  1:213200/76‑1 (MQ=255)
                         caaaaataattgttagttgctacaatgagcttccaaaatatttttttaCAATCCATTATAAATCATCCCAattgat                                              <  2:213200‑M2/28‑1 (MQ=255)
                              cTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAg                                         >  1:366837/1‑76 (MQ=255)
                              cTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAg                                         <  1:40840/76‑1 (MQ=255)
                                 tAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTg                                      >  2:555289/1‑76 (MQ=255)
                                 tAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTg                                      >  2:555293/1‑76 (MQ=255)
                                        taAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACaa                               <  1:619508/76‑1 (MQ=255)
                                        taAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACaa                               <  1:715641/76‑1 (MQ=255)
                                        taAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACaa                               <  1:619514/76‑1 (MQ=255)
                                        taAGAACAGGCCACAATGCATAATTTTTTTTTACAAGCCATTATAAATCATCCCAATTGATTGAAGGTGATAACaa                               <  1:619529/76‑1 (MQ=255)
                                        taAGAACAGCCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACaa                               >  2:349487/1‑76 (MQ=255)
                                              cAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGc                         >  2:448851/1‑76 (MQ=255)
                                                  ccACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATa                     >  1:471149/1‑76 (MQ=255)
                                                    acaATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACa                   <  1:343496/76‑1 (MQ=255)
                                                    acaATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACa                   >  1:243679/1‑76 (MQ=255)
                                                     caATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACaa                  >  1:204924/1‑76 (MQ=255)
                                                     caATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACaa                  >  1:204919/1‑76 (MQ=255)
                                                          naTAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGt             >  1:265161/2‑76 (MQ=255)
                                                                ttttttttACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAag       <  2:185894/76‑1 (MQ=255)
                                                                     tttaCAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCgg  >  1:278100‑M2/5‑76 (MQ=255)
                                                                         |                                                                       
ATTATATAATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCGG  >  REL606/2038868‑2039012

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