New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 2035140NA (NA)7 (0.350) 4/134 2.2 20.0% coding (717/768 nt) wbbD putative glycosyltransferase WbbD
?REL606 = 2035556 28 (1.420)coding (301/768 nt) wbbD putative glycosyltransferase WbbD

ACTAAATCAGTAAGTTGGCAGCATCACCCAAATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGAAATGGCATATTTTATTCCACCTCGTCTTTTTGC  >  REL606/2035066‑2035214
                                                                          |                                                                          
acTAAATCAGTAAGTTGGCAGCATCACCCAAATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGatatat                                                                           >  1:541403/1‑76 (MQ=255)
acTAAATCAGTAAGTTGGCAGCATCACCCAAATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGatatat                                                                           >  2:551079/1‑76 (MQ=255)
   aaaTCAGTAAGTTGGCAGCATCACCCAAATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATATAAt                                                                        >  1:559743/1‑76 (MQ=255)
   aaaTCAGTAAGTTGGCAGCATCACCCAAATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATATAAt                                                                        <  1:26721/76‑1 (MQ=255)
      tCAGTAAGTTGGCAGCATCACCCAAATATCCGTGTGATAATATTTATTATAAAGTCTTGATATGATATATAATTTg                                                                     <  2:205482/76‑1 (MQ=255)
        aGTAAGTTGGCAGCATCACCCAAATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATATAATTTGCa                                                                   <  2:853/76‑1 (MQ=255)
          tAAGTTGGCAGCATCACCCAAATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATATAATTTGCACg                                                                 >  1:654402/1‑76 (MQ=255)
             gTTGGCAGCATCACCCAAATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATATAATTTGCACGATa                                                              >  2:19261/1‑76 (MQ=255)
               tGGCAGCATCACCCAAATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATATAATTTGCACGATaaa                                                            >  1:454399/1‑76 (MQ=255)
                     caTCACCCAAATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGAt                                                      >  1:528609/1‑76 (MQ=255)
                      aTCACCCAAATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATa                                                     <  1:415335/76‑1 (MQ=255)
                       tCACCCAAATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATAc                                                    >  1:174304/1‑76 (MQ=255)
                        cACCCAAATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACt                                                   >  2:341794/1‑76 (MQ=255)
                        cACCCAAATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACt                                                   >  1:290135/1‑76 (MQ=255)
                            cAAATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGc                                               >  1:302515/1‑76 (MQ=255)
                              aaTATCCTTGTGATAATATTTATTATTAAGTCTTGATATGATATAaatagaacaaatcgattcgaatgtatattta                                             >  2:672473‑M1/1‑45 (MQ=255)
                              aaTATCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATAaatagaacaaatcgattcgaatgtatacttg                                             >  2:672460‑M1/1‑45 (MQ=255)
                              aaTATCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCtt                                             >  1:672460/1‑76 (MQ=255)
                                tatCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTa                                           <  1:55378/76‑1 (MQ=255)
                                  tCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTAAt                                         >  2:111860/1‑76 (MQ=255)
                                  tCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTAAt                                         <  1:83699/76‑1 (MQ=255)
                                     ttGTGATAATATTTATTATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTc                                      >  1:127585/1‑76 (MQ=255)
                                          ataataTTTATTATAAAGTCTTGATATGATATAaatagaacaaatcgattcgaatgtatacttgattttatgacta                                 >  2:422554‑M1/1‑33 (MQ=255)
                                          ataataTTTATTATAAAGTCTTGATATGATATAaatagaacaaatcgattcgaatgtatacttgattttatgacta                                 >  2:422553‑M1/1‑33 (MQ=255)
                                          ataataTTTATTATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGaaa                                 >  1:422554/1‑76 (MQ=255)
                                          ataataTTTATTATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGaaa                                 >  1:422553/1‑76 (MQ=255)
                                          ataataTTTATTATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGaaa                                 <  2:193033/76‑1 (MQ=255)
                                            aataTTTATTATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGAAATg                               >  1:629436/1‑76 (MQ=255)
                                             atatTTATTATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGAAATgg                              >  2:77916/1‑76 (MQ=255)
                                               ctttattATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGAAATGGCa                            >  1:396449/2‑76 (MQ=255)
                                               atttattATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGAAATGGCa                            >  2:603402/1‑76 (MQ=255)
                                                tttattATAAAGTCTTGATATGATATAaatagaacaaatcgattcgaatgtatacttgattttatgactaagaatg                           >  2:95311‑M1/1‑27 (MQ=255)
                                                tttattATAAAGTCTTGATATGATATATAATTTGCACGGTAAAAATGATACTGAGCTTTAATTTCTGAAATGGCat                           >  1:244526/1‑76 (MQ=255)
                                                tttattATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGAAATGGCat                           >  1:95311/1‑76 (MQ=255)
                                                tttattATAAAGTCTTGATATGACATAaatagaaccaatcgattggaatgtatatttggttttatgactaagaatg                           >  2:244526‑M1/1‑27 (MQ=255)
                                                 ttattaTAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGAAATGGCata                          >  1:163317/1‑76 (MQ=255)
                                                     tataAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGAAATGGCATAtttt                      >  1:603332/1‑76 (MQ=255)
                                                       taAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGAAATGGCATATTTTAt                    >  2:293279/1‑76 (MQ=255)
                                                           gTCTTGATATGATATAaatagaacaaatcgattcgaatgtatacttgattttatgactaagaatggtgatgttcat                >  2:664065‑M1/1‑16 (MQ=255)
                                                           gTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGAAATGGCATATTTTATTCCa                >  1:664065/1‑76 (MQ=255)
                                                             cTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGAAATGGCATATTTTATTCCAcc              >  2:16664/1‑76 (MQ=255)
                                                             cTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGAAATGGCATATTTTATTCCAcc              >  1:443266/1‑76 (MQ=255)
                                                             cTTGACATGATATATAATTTGCACGATAAAAATGACACTGAGCTTTAATTTCTGAAATGGCATATTTTATTCCAcc              >  2:643307/1‑76 (MQ=255)
                                                                gatatgatatATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGAAATGGCATATTTTATTCCACCTCg           >  1:103242/1‑76 (MQ=255)
                                                                  tatgatatATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGAAATGGCATATTTTATTCCACCTCGTc         >  2:605669/1‑76 (MQ=255)
                                                                   atgatatATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGAAATGGCATATTTTATTCCACCTCGTCt        >  1:158208/1‑76 (MQ=255)
                                                                       tataaatagaacaaatcgattcgaatgtatacttgattttatgactaagaatggtgatgttcatattcttagctct    >  2:642574‑M1/1‑4 (MQ=255)
                                                                        atataATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGAAATGGCATATTTTATTCCACCTCGTCTTTTTg   <  2:336984/76‑1 (MQ=255)
                                                                        ataaatagaacaaatcgattcgaatgtatacttgattttatgactaagaatggtgatgttcatattcttagctctt   >  2:357080‑M1/1‑3 (MQ=255)
                                                                         tataATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGAAATGGCATATTTTATTCCACCTCGTCTTTTTGc  >  1:237871/1‑76 (MQ=255)
                                                                         tataATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGAAATGGCATATTTTATTCCACCTCGTCTTTTTGc  >  1:17237/1‑76 (MQ=255)
                                                                          |                                                                          
ACTAAATCAGTAAGTTGGCAGCATCACCCAAATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATATAATTTGCACGATAAAAATGATACTGAGCTTTAATTTCTGAAATGGCATATTTTATTCCACCTCGTCTTTTTGC  >  REL606/2035066‑2035214

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