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

CATACACTAAATCAGTAAGTTGGCAGCATCACCCAAATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/2035061‑2035140
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atgatataAATAGAACAAATCGATTCGAATGTATACTTGATTTTATGACTAAGAATGGTGATGTTCATATTCTTAGCTCTT  <  REL606/2035556‑2035484
                                                                                                                                                         
CATACACTAAATCAGTAAGTTGGCAGCATCACCCAAATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGA                                                                               >  1:284901/1‑76
                                   AATATCCTTGTGATAATATTTATTATTAAGTCTTGATATGATATAAATAGAACAAATCGATTCGAATGTATATTTa                                            >  2:672473/1‑75
                                   AATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATAAATAGAACAAATCGATTCGAATGTATACTTG                                            >  2:672460/1‑76
                                               ATAATATTTATTATAAAGTCTTGATATGATATAAATAGAACAAATCGATTCGAATGTATACTTGATTTTATGACTA                                >  2:422553/1‑76
                                               ATAATATTTATTATAAAGTCTTGATATGATATAAATAGAACAAATCGATTCGAATGTATACTTGATTTTATGACTA                                >  2:422554/1‑76
                                                     TTTATTATAAAGTCTTGATATGATATAAATAGAACAAATCGATTCGAATGTATACTTGATTTTATGACTAAGAATG                          >  2:95311/1‑76
                                                     TTTATTATAAAGTCTTGATATGACATAAATAGAACCAATCGATTGGAATGTATATTTGGTTTTATGACTAAGAATG                          >  2:244526/1‑76
                                                                GTCTTGATATGATATAAATAGAACAAATCGATTCGAATGTATACTTGATTTTATGACTAAGAATGGTGATGTTCAT               >  2:664065/1‑76
                                                                            TATAAATAGAACAAATCGATTCGAATGTATACTTGATTTTATGACTAAGAATGGTGATGTTCATATTCTTAGCTCT   >  2:642574/1‑76
                                                                             ATAAATAGAACAAATCGATTCGAATGTATACTTGATTTTATGACTAAGAATGGTGATGTTCATATTCTTAGCTCTT  >  2:357080/1‑76
                                                                                                                                                         
CATACACTAAATCAGTAAGTTGGCAGCATCACCCAAATATCCTTGTGATAATATTTATTATAAAGTCTTGATATGATATA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/2035061‑2035140
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atgatataAATAGAACAAATCGATTCGAATGTATACTTGATTTTATGACTAAGAATGGTGATGTTCATATTCTTAGCTCTT  <  REL606/2035556‑2035484

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