New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 203856072 (1.230)13 (0.260) 5/126 5.9 21.1% coding (389/822 nt) wbbB putative glycosyltransferase WbbB
?REL606 = 2039125 37 (0.750)coding (570/747 nt) wbbA putative glycosyltransferase WbbA
Rejected: Coverage evenness skew score above cutoff.

AAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTT  >  REL606/2039050‑2039198
                                                                           |                                                                         
aagtaTCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTc                                                                           <  2:316582/76‑1 (MQ=255)
  aaaTCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCat                                                                         >  2:1192240/3‑76 (MQ=255)
   taTCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCaac                                                                        <  2:1321554/76‑3 (MQ=255)
    aTCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCataa                                                                       >  2:1841436/1‑76 (MQ=255)
           ataatTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGt                                                                <  2:1810089/76‑1 (MQ=255)
            taatTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGtt                                                               <  1:356773/76‑1 (MQ=255)
            taatTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGtt                                                               <  1:356751/76‑1 (MQ=255)
             aatTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTg                                                              <  1:1722820/76‑1 (MQ=255)
               tttGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGca                                                            >  2:265379/1‑76 (MQ=255)
                     aTCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTa                                                      >  1:237993/1‑76 (MQ=255)
                      tCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTAc                                                     <  1:1241615/76‑1 (MQ=255)
                         tAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTg                                                  >  2:441180/1‑76 (MQ=255)
                               gCTAAGTGTGGAATGATGTGTGTCCTAGCAAAAGCCCGGGCTTTCATAATAGATGTTGCACATGTACTTGGTGCTa                                            <  2:177357/76‑1 (MQ=255)
                                 tAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTAtt                                          >  1:681353/1‑76 (MQ=255)
                                 tAATGTTAGAATTATGTGTGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTAtt                                          <  2:1214190/76‑1 (MQ=255)
                                       tAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCaa                                    >  1:793929/1‑76 (MQ=255)
                                        aGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAAt                                   >  1:715902/1‑76 (MQ=255)
                                        aGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAAt                                   <  2:575843/76‑1 (MQ=255)
                                            ttattTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCt                               >  2:1169756/1‑76 (MQ=255)
                                               tttttGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACt                            <  1:456761/76‑1 (MQ=255)
                                                  ttGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTccc                         >  1:1371702/1‑76 (MQ=255)
                                                   tGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCa                        >  1:164695/1‑76 (MQ=255)
                                                    gTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAg                       <  1:625004/76‑1 (MQ=255)
                                                    gTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAg                       <  1:1163342/76‑1 (MQ=255)
                                                    gTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAg                       <  2:1534691/76‑1 (MQ=255)
                                                       cTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAAt                    <  1:539289/76‑1 (MQ=255)
                                                        tATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATa                   <  2:699110/76‑1 (MQ=255)
                                                          tCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATaaa                 <  1:329229/76‑1 (MQ=255)
                                                           cAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATaaaa                <  1:859404/76‑1 (MQ=255)
                                                           cAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATaaaa                <  1:859376/76‑1 (MQ=255)
                                                            aaaaaCCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAAt               >  1:392928/1‑76 (MQ=255)
                                                            aaaaaCCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAAt               <  2:309652/76‑1 (MQ=255)
                                                            aaaaaCCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAAt               <  1:1403523/76‑1 (MQ=255)
                                                            aaaaaCCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAAt               >  2:993255/1‑76 (MQ=255)
                                                              aaaCCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATAc             <  1:1246868/76‑1 (MQ=255)
                                                               aaCCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACa            <  2:1833215/76‑1 (MQ=255)
                                                               aaCCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACa            <  1:272887/76‑1 (MQ=255)
                                                                aCCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACaa           <  2:386666/76‑1 (MQ=255)
                                                                aCCCTCGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCACCATAAAAAACaa           >  2:1731503/1‑76 (MQ=255)
                                                                 cccTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAAt          <  1:912056/76‑1 (MQ=255)
                                                                  ccTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATa         <  1:1677275/76‑1 (MQ=255)
                                                                        tttCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGttttt   <  2:726793/76‑1 (MQ=255)
                                                                         ttCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCCCCATAACATACACTCAtttttt  >  1:1849918/1‑76 (MQ=255)
                                                                         ttCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGattttt  >  2:1849918/1‑76 (MQ=255)
                                                                           |                                                                         
AAGTATCAAAAATAATTTGTAATCGTAAATAGCTAATTTTAGAATTATTTTTGTCCTATCAAAAACCCTGGCTTTCATAATATATGTTGCACATGTACTTGGTGCTATTTTCCAATGCCTACTCCCAGAATAAAATACAATAGTTTTTT  >  REL606/2039050‑2039198

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.