New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2039283 =29 (0.530)10 (0.190) 7/142 5.2 23.1% coding (412/747 nt) wbbA putative glycosyltransferase WbbA
?REL606 2039547 = 39 (0.760)coding (148/747 nt) wbbA putative glycosyltransferase WbbA
Rejected: Coverage evenness skew score above cutoff.

TATATTTTTTAGAAAATGATTATTTGCATGAGTCTAAATGGATTGATGAGATATTTAACTTGGTCAAATCAAATATAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/2039360‑2039283
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ataaAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATA  >  REL606/2039547‑2039596
                                                                                                                                
TATATTTTTTAGAAAATGATTATTTGCATGAGTCTAAATGGATTGATGAGATATTTAACTTGGTCAAATCAAATAT                                                      >  1:822101/1‑76
TATATTTTTTAGAAAATGATTATTTGCATGAGTCTAAATGGATTGATGAGATATTTAACTTGGTCAAATCAAATAT                                                      >  2:1501926/1‑76
      TTTTAGAAAATGATTATTTGCATGAGTCTAAATGGATTGATGAGATATTTAACTTGGTCAAATCAAATATAAAAGT                                                <  2:104299/76‑1
           GAAAATGATTATTTGCATGAGTCTAAATGGATTGATGAGATATTTAACTTGGTCAAATCAAATATAAAAGTAAATG                                           <  2:1267132/76‑1
            AAAATGATTATTTGCATGAGTCTAAATGGATTGATGAGATATTTAACTTGGTCAAATCAAATATAAAAGTAAATGC                                          <  2:1464840/76‑1
              AATGATTATTTGCATGAGTCTAAATGGATTGATGAGATATTTAACTTGGTCAAATCAAATATAAAAGTAAATGCTA                                        <  2:1628628/76‑1
                           ATGAGTCTAAATGGATTGATGAGATATTTAACTTGGTCAAATCAAATATAAAAGTAAATGCTATTTGGTCTTTATA                           <  2:783653/76‑1
                            TGAGTCTAAATGGATTGATGAGATATTTAACTTGGTCAAATCAAATATAAAAGTAAATGCTATTTGGTCTTTATAT                          <  2:1239409/76‑1
                                                    ATTTAACTTGGTCAAATCAAATATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATA  <  2:497079/76‑1
                                                    ATTTAACTTGGTCAAATCAAATATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATA  <  2:497096/76‑1
                                                    ATTTAACTTGGTCAAATCAAATATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATA  <  2:877829/76‑1
                                                    ATTTAACTTGGTCAAATCAAAAATAAAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATA  <  2:877811/76‑1
                                                                                                                                
TATATTTTTTAGAAAATGATTATTTGCATGAGTCTAAATGGATTGATGAGATATTTAACTTGGTCAAATCAAATATAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/2039360‑2039283
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ataaAAGTAAATGCTATTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATA  >  REL606/2039547‑2039596

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