New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2038074 =55 (0.920)8 (0.150) 7/134 5.6 12.7% coding (57/1176 nt) wzy O antigen polymerase (Wzy)
?REL606 2038286 = 61 (1.140)coding (663/822 nt) wbbB putative glycosyltransferase WbbB
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AATCAAGAGGAAGATTATAATTAAATGAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTAATATT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/2038154‑2038074
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttaatattAATTCGTTGTTTATCATTATGTAACAATTTAGCTTGATACCAAACCGTATCTGTGTCTATATCAATTAAATG  >  REL606/2038286‑2038357
                                                                                                                                                         
AATCAAGAGGAAGATTATAATTAAATGAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTA                                                                               <  1:1370734/76‑1
AATCAAGAGGAAGATTATAATTAAATGAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTA                                                                               <  1:812067/76‑1
 ATCAAGAGGAAGATTATAATTAAATGAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTAA                                                                              <  1:1425359/76‑1
 ATCAAGAGGAAGATTATAATTAAATGAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTAA                                                                              <  1:1425370/76‑1
 ATCAAGAGGAAGATTATAATTAAATGAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTAA                                                                              >  1:73179/1‑76
   CAAGAGGAAGATTATAATTAAATGAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTAATA                                                                            >  1:1396547/1‑76
   CAAGAGGAAGATTATAATTAAATGAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTAATA                                                                            <  1:81121/76‑1
    AAGAGGAAGATTATAATTAAATGAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTAATAT                                                                           <  1:713255/76‑1
     AGAGGAAGATTATAATTAAATGAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTAATATT                                                                          >  1:1207019/1‑76
     AGAGGAAGATTATAATTAAATGAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTAATATT                                                                          >  2:397209/1‑76
     AGAGGAAGATTATAATTAAATGAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTAATATT                                                                          >  1:1206993/1‑76
         GAAGATTATAATTAAATGAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTAATATTAATT                                                                      <  2:1734659/76‑1
         GAAGATTATAATTAAATGAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTAATATTAATT                                                                      <  2:75577/76‑1
          AAGATTATAATTAAATGAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTAATATTAATTC                                                                     <  2:1657069/76‑1
           AGATTATAATTAAATGAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTAATATTAATTCG                                                                    <  2:1228693/76‑1
                      AAATGAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTAATATTAATTCGTTGTTTATCAT                                                         <  2:1836861/76‑1
                       AATGAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTAATATTAATTCGTTGTTTATCATT                                                        <  2:1084354/76‑1
                          GAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTAATATTAATTCGTTGTTTATCATTATG                                                     <  2:565964/76‑1
                           AGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTAATATTAATTCGTTGTTTATCATTATGT                                                    <  2:466109/76‑1
                                                                          TAATATTAATTCGTTGTTTATCATTATGTAACAATTTAGCTTGATACCAAACCGTATCTGTGTCTATATCAATTAA     >  1:845962/1‑76
                                                                            ATATTAATTCGTTGTTTATCATTATGTAACAATTTAGCTTGATACCAAACCGTATCTGTGTCTATATCAATTAAAT   <  1:1694082/76‑1
                                                                            ATATTAATTCGTTGTTTATCATTATGTAACAATTTAGCTTGATACCAAACCGTATCTGTGTCTATATCAATTAAAT   <  1:1443068/76‑1
                                                                             TATTAATTCGTTGTTTATCATTATGTAACAATTTAGCTTGATACCAAACCGTATCTGTGTCTATATCAATTAAATG  >  1:324162/1‑76
                                                                                                                                                         
AATCAAGAGGAAGATTATAATTAAATGAGAAAAAATTTTTACACGACACCAATCACATTACTAATTCTTACGCTTAATATT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/2038154‑2038074
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttaatattAATTCGTTGTTTATCATTATGTAACAATTTAGCTTGATACCAAACCGTATCTGTGTCTATATCAATTAAATG  >  REL606/2038286‑2038357

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