New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 203691756 (1.030)9 (0.180) 7/140 5.1 18.3% coding (42/1113 nt) wbbC putative glycosyltransferase WbbC
?REL606 = 2037192 28 (0.550)coding (939/1176 nt) wzy O antigen polymerase (Wzy)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTCCTAAGTG  >  REL606/2036844‑2036986
                                                                         |                                                                     
caatcCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTaa                                                                     <  2:463264/76‑1 (MQ=255)
 aatcCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAAt                                                                    >  2:1749033/1‑76 (MQ=255)
  atcCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAAtt                                                                   <  2:619159/76‑1 (MQ=255)
    cccATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCt                                                                 <  1:1298063/76‑1 (MQ=255)
     ccATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTg                                                                <  2:1733329/76‑1 (MQ=255)
      cATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGa                                                               >  1:352043/1‑76 (MQ=255)
       aTCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGaa                                                              >  1:1115027/1‑76 (MQ=255)
        tCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGaaa                                                             <  1:381596/76‑1 (MQ=255)
         cGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGaaaa                                                            >  1:312796/1‑76 (MQ=255)
         cGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGaaaa                                                            <  1:766877/76‑1 (MQ=255)
                ttCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGt                                                     >  1:138928/1‑76 (MQ=255)
                   aTTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTaaa                                                  <  1:1445979/76‑1 (MQ=255)
                   aTTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTaaa                                                  <  1:375713/76‑1 (MQ=255)
                   aTTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTaaa                                                  <  1:1306491/76‑1 (MQ=255)
                    ttAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAAt                                                 <  1:563632/76‑1 (MQ=255)
                      aaTATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTg                                               <  1:999698/76‑1 (MQ=255)
                      aaTATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTg                                               <  1:926075/76‑1 (MQ=255)
                        tatTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGaa                                             <  1:236766/76‑1 (MQ=255)
                        tatTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGaa                                             <  1:236773/76‑1 (MQ=255)
                        tatTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGaa                                             <  1:236793/76‑1 (MQ=255)
                         atTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGaaa                                            <  2:1499419/76‑1 (MQ=255)
                          ttttttGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGaaaa                                           >  1:1660291/1‑76 (MQ=255)
                            ttttGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACt                                         >  1:461059/1‑76 (MQ=255)
                            ttttGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACt                                         >  1:32719/1‑76 (MQ=255)
                             tttGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTa                                        >  1:117273/1‑76 (MQ=255)
                             tttGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTa                                        <  2:138928/76‑1 (MQ=255)
                                    ttAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGt                                 <  1:1742409/76‑1 (MQ=255)
                                      aaaTCTCTAAAGTATACTCCAACCCCATAATATTTTTATTTTGAAAAAGAAAGTAAATTGAAAACTATTTTCGTCa                               >  1:868758/1‑76 (MQ=255)
                                      aaaTCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCa                               <  1:839313/76‑1 (MQ=255)
                                      aaaTCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCa                               >  1:761453/1‑76 (MQ=255)
                                       aaTCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCAt                              >  1:606492/1‑76 (MQ=255)
                                        aTCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATa                             >  1:498381/1‑76 (MQ=255)
                                        aTCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATa                             >  1:506719/1‑76 (MQ=255)
                                        aTCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATa                             >  1:710025/1‑76 (MQ=255)
                                           tctAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAg                          <  1:227005/76‑1 (MQ=255)
                                           tctAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAg                          >  2:201327/1‑76 (MQ=255)
                                           tctAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAg                          >  1:627451/1‑76 (MQ=255)
                                              aaaGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGtt                       <  2:1476126/76‑1 (MQ=255)
                                               aaGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGttt                      >  1:1611176/1‑76 (MQ=255)
                                                tggataCTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGtttt                     >  2:1555958/4‑76 (MQ=255)
                                                       tCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCttt              >  1:575761/1‑76 (MQ=255)
                                                        ccAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCtttt             >  1:1414689/1‑76 (MQ=255)
                                                         cAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTa            >  2:636509/1‑76 (MQ=255)
                                                           aCCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGt          >  1:729180/1‑76 (MQ=255)
                                                             cccATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTcc        <  1:1103219/76‑1 (MQ=255)
                                                              ccATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTCCt       >  1:575837/1‑76 (MQ=255)
                                                              ccATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTCCt       <  1:1170599/76‑1 (MQ=255)
                                                              ccATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTCCt       <  2:424503/76‑1 (MQ=255)
                                                              ccATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTCCt       <  2:424526/76‑1 (MQ=255)
                                                              ccATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTCCt       <  2:582662/76‑1 (MQ=255)
                                                               cATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTCCTa      >  1:481357/1‑76 (MQ=255)
                                                               cATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTCCTa      >  1:1495406/1‑76 (MQ=255)
                                                                   atatTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTCCTAAgtg  <  1:1129959/76‑1 (MQ=255)
                                                                   atatTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTCCTAAgtg  <  1:1129956/76‑1 (MQ=255)
                                                                   atatTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTCCTAAgtg  <  1:1129944/76‑1 (MQ=255)
                                                                   atatTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTCCTAAgtg  >  1:1077241/1‑76 (MQ=255)
                                                                         |                                                                     
CAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTCCTAAGTG  >  REL606/2036844‑2036986

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