New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 2037696103 (0.650)3 (0.020) 3/134 17.6 2.9% coding (435/1176 nt) wzy O antigen polymerase (Wzy)
?REL606 = 2037796 110 (0.780)coding (335/1176 nt) wzy O antigen polymerase (Wzy)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TAAATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/2037614‑2037696
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑taataatgAGCGCGGAGGCGGAGTAATTTCAGCGTTATTTGTATTACTGAATGTTGATGTACTTGTTTTAATATTTCTAACAC  <  REL606/2037796‑2037722
                                                                                                                                                              
TAAATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGAT                                                                                    >  1:2461377/1‑76
TAAATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGAT                                                                                    >  2:129402/1‑76
 AAATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATA                                                                                   <  2:657153/76‑1
 AAATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATA                                                                                   >  1:5077568/1‑76
  AATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAA                                                                                  <  1:4649521/76‑1
  AATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAA                                                                                  <  1:4687662/76‑1
   ATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAAT                                                                                 >  2:3492934/1‑76
   ATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAAT                                                                                 <  1:3725911/76‑1
     GAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAA                                                                               <  1:2848775/76‑1
     GAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAA                                                                               <  1:5082225/76‑1
     GAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAA                                                                               >  2:811479/1‑76
      AGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAAT                                                                              <  1:3328195/76‑1
       GATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATG                                                                             >  1:4451776/1‑76
                                           ACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGAGCGCGGAGGCGGAGTAATTTCAGCGTTATTTGTAT                                         >  2:3146624/1‑76
                                                 AAGAAAAGAAATAATCCACAACATGATAATAATGAGCGCGGAGGCGGAGTAATTTCAGCGTTATTTGTATTACTGA                                   >  2:3594026/1‑76
                                                                        TGATAATAATGAGCGCGGAGGCGGAGTAATTTCAGCGTTATTTGTATTACTGAATGTTGATGTACTTGTTTTAATA            >  2:776815/1‑76
                                                                           TAATAATGAGCGCGGAGGCGGAGTAATTTCAGCGTTATTTGTATTACTGAATGTTGATGTACTTGTTTTAATATTT         <  2:4802995/76‑1
                                                                            AATAATGAGCGCGGAGGCGGAGTAATTTCAGCGTTATTTGTATTACTGAATGTTGATGTACTTGTTTTAATATTTC        >  2:34974/1‑76
                                                                              TAATGAGCGCGGAGGCGGAGTAATTTCAGCGTTATTTGTATTACTGAATGTTGATGTACTTGTTTTAATATTTCTA      <  2:634074/76‑1
                                                                              TAATGAGCGCGGAGGCGGAGTAATTTCAGCGTTATTTGTATTACTGAATGTTGATGTACTTGTTTTAATATTTCTA      <  2:1656448/76‑1
                                                                              TAATGAGCGCGGAGGCGGAGTAATTTCAGCGTTATTTGTATTACTGAATGTTGATGTACTTGTTTTAATATTTCTA      <  2:1431012/76‑1
                                                                              TAATGAGCGCGGAGGCGGAGTAATTTCAGCGTTATTTGTATTACTGAATGTTGATGTACTTGTTTTAATATTTCTA      <  1:3741079/76‑1
                                                                              TAATGAGCGCGGAGGCGGAGTAATTTCAGCGTTATTTGTATTACTGAATGTTGATGTACTTGTTTTAATATTTCTA      <  1:3711942/76‑1
                                                                               AATGAGCGCGGAGGCGGAGTAATTTCAGCGTTATTTGTATTACTGAATGTTGATGTACTTGTTTTAATATTTCTAA     >  2:4672381/1‑76
                                                                               AATGAGCGCGGAGGCGGAGTAATTTCAGCGTTATTTGTATTACTGAATGTTGATGTACTTGTTTTAATATTTCTAA     <  2:4057132/76‑1
                                                                               AATGAGCGCGGAGGCGGAGTAATTTCAGCGTTATTTGTATTACTGAATGTTGATGTACTTGTTTTAATATTTCTAA     <  2:1305519/76‑1
                                                                                  GAGCGCGGAGGCGGAGTAATTTCAGCGTTATTTGTATTACTGAATGTTGATGTACTTGTTTTAATATTTCTAACAC  >  2:567003/1‑76
                                                                                  GAGCGCGGAGGCGGAGTAATTTCAGCGTTATTTGTATTACTGAATGTTGATGTACTTGTTTTAATATTTCTAACAC  >  2:976718/1‑76
                                                                                                                                                              
TAAATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/2037614‑2037696
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑taataatgAGCGCGGAGGCGGAGTAATTTCAGCGTTATTTGTATTACTGAATGTTGATGTACTTGTTTTAATATTTCTAACAC  <  REL606/2037796‑2037722

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