New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1534318 =71 (0.690)12 (0.150) 4/38 NT 24.9% coding (550/1023 nt) yddR D‑ala‑D‑ala transporter subunit
?REL606 = 2036370 16 (0.200)coding (589/1113 nt) wbbC putative glycosyltransferase WbbC

GTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGAT  >  REL606/2036337‑2036404
                                 |                                  
gTTTATCTTCATGCAAAGCACCCACCGAGAAAAAa                                   <  1:17396908‑M2/35‑2 (MQ=255)
     tCTTCATGCAAAGCACCCACCGAGAAAAAAAAAgg                              <  1:5632973/35‑1 (MQ=255)
      cTTCATGCAAAGCACCCACCGAGAAAAAAAAAGGt                             >  1:9962516/1‑35 (MQ=255)
      cTTCATGCAAAGCACCCACCGAGAAAAAAAAAGGt                             >  1:20547159/1‑35 (MQ=255)
         cATGCAAAGCACCCACCGAGAAAAAAAAAGGTTca                          <  1:5092057/35‑1 (MQ=255)
           tGCAAAGCACCCACCGAGAAAAAAAAAGGTTcaca                        <  1:7260203/35‑1 (MQ=255)
           tGCAAAGCACCCACCGAGAAAAAAAAAGGTTcaca                        <  1:7205585/35‑1 (MQ=255)
            gCAAAGCACCCACCGAGAAAAAAAAAGGTTCACaa                       >  1:14101442/1‑35 (MQ=255)
              aaaGCACCCACCGAGAAAAAAAAAGGTTCACAATc                     <  1:18050836/35‑1 (MQ=255)
               aaGCACCCACCGAGAAAAAAAAAGGTTCACAATCa                    >  1:1326609/1‑35 (MQ=255)
               aaGCACCCACCGAGAAAAAAAAAGGTTCACAATCa                    >  1:10642294/1‑35 (MQ=255)
                  cccccaccGAGAAAAAAAAAGGTTCACAATCatat                 <  1:9965452/33‑1 (MQ=255)
                    cccaccGAGAAAAAAAAAGGTTCACAATCATATTg               >  1:15771617/1‑35 (MQ=255)
                     gcagcgAGAAAAAaaaagccggtaacgtgcgttgg              >  1:7319422‑M2/2‑13 (MQ=255)
                     gcaccaAGAAAAAaaaagccggtaacgtgcgttgg              >  1:14878415‑M2/2‑13 (MQ=255)
                     gcaacgAGAAAAAaaaagccggtaacgtgcgttgg              >  1:7796676‑M2/2‑13 (MQ=255)
                      caccGAGAAAAAAAAAGGTTCACAATCATATTGgt             <  1:20454440/35‑1 (MQ=255)
                       accGAAAAAAAaaaagccggtaacgtgcgttggtg            >  1:3454113‑M2/1‑11 (MQ=255)
                        gcgAAAAAAAaaaagccggtaacgtgcgttggtgg           >  1:16611306‑M2/2‑10 (MQ=255)
                        ccGAGAAAAAAAAAGGTTCACAATCATATTGgtgt           <  1:19079918/35‑1 (MQ=255)
                        ccGAGAAAAAAAAAGGTTCACAATCATATTGgtgt           >  1:13169272/1‑35 (MQ=255)
                         ggagaaAAAaaaagccggtaacgtgcgttggtgga          >  1:13390894‑M2/2‑9 (MQ=255)
                         cGAGAAAAAaaaagccggtaacgtgcgttggtgga          >  1:8514367‑M2/1‑9 (MQ=255)
                         cGAGAAAAAAAAAGGTTCACAATCATATTGGTGta          <  1:10177228/35‑1 (MQ=255)
                         cGAGAAAAAAAAAGGTTCACAATCATATTGGTGta          >  1:9210303/1‑35 (MQ=255)
                          gagaAAAAAAAATGTTCACAATCATATTGGTGtat         >  1:14472584/1‑35 (MQ=255)
                          aagaaAAAaaaagccggtaacgtgcgttggtggat         >  1:5416206‑M2/2‑8 (MQ=255)
                           agaAGAAaaaagccggtaacgtgcgttggtggatc        >  1:18516955‑M2/1‑7 (MQ=255)
                           agaAAAAaaaagccggtaacgtgcgttggtggatc        >  1:5664242‑M2/1‑7 (MQ=255)
                           agaAAAAaaaagccggtaacgtgcgttggtggatc        >  1:14162110‑M2/1‑7 (MQ=255)
                           agaAAAAAAAAGGTTCACAATAATATTGGTGtata        <  1:12059627/35‑1 (MQ=255)
                            gaAAAAaaaagccggtaacgtgcgttggtggatcc       >  1:18963262‑M2/1‑6 (MQ=255)
                            gaAAAAAAAAGGTTCACAATCATATTGGTGTATag       >  1:19824118/1‑35 (MQ=255)
                             aaaaaaaaagccggtaacgtgcgttggtggatcca      >  1:16322120‑M2/1‑5 (MQ=255)
                             aaaaaaaaaGGTTCACAATCATATTGGTGTATaga      <  1:13749163/35‑1 (MQ=255)
                             aaaaaaaaaGGTTCACAATCATATTGGTGTATaga      >  1:17650263/1‑35 (MQ=255)
                              aaaaaaaagccggtaacgtgcgttggtggatccag     >  1:11425654‑M2/1‑4 (MQ=255)
                              aaaaaaaagccggtaacgtgcgttggtggatccag     >  1:8961321‑M2/1‑4 (MQ=255)
                              aaaaaaaagccggtaacgtgcgttggtggatccag     >  1:9761052‑M2/1‑4 (MQ=255)
                              aaaaaaaaGGTTCACAATCATATTGGTGTATagag     >  1:20005452/1‑35 (MQ=255)
                              aaaaaaaaGGTTCACAATCATATTGGTGTATagag     <  1:6429631/35‑1 (MQ=255)
                               aaaaaaagccggtaacgtgcgttggtggatccagc    >  1:12339290‑M2/1‑3 (MQ=255)
                                aaaaaagccggtaacgtgcgttggtggatccagcc   >  1:8274339‑M2/1‑2 (MQ=255)
                                 aaaaagccggtaacgtgcgttggtggatccagcca  >  1:15980558‑M2/1‑1 (MQ=255)
                                 aaaaagccggtaacgtgcgttggtggatccagcca  >  1:11059650‑M2/1‑1 (MQ=255)
                                 aaaaaGGTTCACAATCATATTGGTGTATAGAGGAt  <  1:6321103/35‑1 (MQ=255)
                                 aaaaaGGTTCACAATCATATTGGTGTATAGAGGAt  <  1:19420843/35‑1 (MQ=255)
                                 |                                  
GTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGAT  >  REL606/2036337‑2036404

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