New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2037389 =30 (0.760)3 (0.080) 3/142 4.7 7.9% coding (742/1176 nt) wzy O antigen polymerase (Wzy)
?REL606 2037953 = 42 (1.120)coding (178/1176 nt) wzy O antigen polymerase (Wzy)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CGGTCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAACCCCATTAAA  >  REL606/2037879‑2038026
                                                                          |                                                                         
cGGTCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGa                                                                          >  1:1185637/1‑76 (MQ=255)
 ggTCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGaa                                                                         <  1:1040287/76‑1 (MQ=255)
 ggTCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGaa                                                                         >  1:369181/1‑76 (MQ=255)
 ggTCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGaa                                                                         <  1:31687/76‑1 (MQ=255)
   tCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGTTAGAaat                                                                       <  1:454593/76‑1 (MQ=255)
   tCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAaat                                                                       >  1:78336/1‑76 (MQ=255)
    ctactaTTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAaata                                                                      >  1:928914/1‑76 (MQ=255)
      actaTTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAaataat                                                                    <  1:925758/76‑1 (MQ=255)
          ttAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGc                                                                <  1:765047/76‑1 (MQ=255)
             gAATCCGCTAATTTTGTTGTAAGTGGTTTTATTTTTGCCAAGCCATTATATAATACCGATAGAAATAATAAGCTaa                                                             <  1:635318/76‑1 (MQ=255)
             gAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTaa                                                             >  2:59176/1‑76 (MQ=255)
             gAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTaa                                                             <  1:341434/76‑1 (MQ=255)
               aTCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAaga                                                           <  1:297421/76‑1 (MQ=255)
                 cAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAgaga                                                         <  1:1071635/76‑1 (MQ=255)
                     tAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGatata                                                     <  1:518229/76‑1 (MQ=255)
                     tAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGatata                                                     <  1:733835/76‑1 (MQ=255)
                      aaTTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGatgtat                                                    >  1:717595/1‑76 (MQ=255)
                      aaTTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGatatat                                                    >  2:762761/1‑76 (MQ=255)
                       aTTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGAtatata                                                   <  1:236028/76‑1 (MQ=255)
                       aTTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGAtatata                                                   >  2:117846/1‑76 (MQ=255)
                        ttttGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATaa                                                  <  1:22093/76‑1 (MQ=255)
                          ttgttgTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAAt                                                <  1:862887/76‑1 (MQ=255)
                             ttgtAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTag                                             >  2:618339/1‑76 (MQ=255)
                             ttgtAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTag                                             >  2:533728/1‑76 (MQ=255)
                                   gTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTg                                       <  1:912500/76‑1 (MQ=255)
                                   gTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTg                                       >  1:444850/1‑76 (MQ=255)
                                      gTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAat                                    <  2:185989/76‑1 (MQ=255)
                                       ttttATTTTTGCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAatt                                   <  1:295811/76‑1 (MQ=255)
                                        tttATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAatta                                  <  1:243905/76‑1 (MQ=255)
                                           aTTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAAttatta                               <  1:13538/76‑1 (MQ=255)
                                           aTTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAAttatta                               >  1:1172150/1‑76 (MQ=255)
                                            ttttttCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTAc                              <  2:1078718/76‑1 (MQ=255)
                                             tttttCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACg                             >  1:313056/1‑76 (MQ=255)
                                              ttttCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGa                            >  2:92974/1‑76 (MQ=255)
                                                ttCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACa                          <  1:360978/76‑1 (MQ=255)
                                                  ccAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGt                        <  1:92974/76‑1 (MQ=255)
                                                  ccAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGt                        >  1:1089776/1‑76 (MQ=255)
                                                        cATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGa                  <  1:289937/76‑1 (MQ=255)
                                                           tatataATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGAttt               >  1:1011256/1‑76 (MQ=255)
                                                           tatataATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGAttt               <  1:1164252/76‑1 (MQ=255)
                                                              ataataCCGATAGCAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTaga            >  1:842763/1‑74 (MQ=255)
                                                                   aCCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAACCCCa       >  1:355966/1‑76 (MQ=255)
                                                                       aTAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAACCCCATTaa   <  1:167671/76‑1 (MQ=255)
                                                                        tAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAACCCCATTaaa  <  2:369181/76‑1 (MQ=255)
                                                                          |                                                                         
CGGTCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAACCCCATTAAA  >  REL606/2037879‑2038026

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