New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2041054 =50 (0.770)6 (0.100) 4/136 7.6 10.2% coding (330/1113 nt) vioA VioA, involved in dTDP‑N‑acetylviosamine synthesis
?REL606 2041468 = 61 (1.030)coding (1353/1428 nt) wzx O antigen translocase (Wzx)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTTGAGTAGCAATGGACAGGTAAAATTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATAGCATGTGAAGTTGCAAC  >  REL606/2040980‑2041128
                                                                          |                                                                          
gtTGAGTAGCAATGGACAGGTAAAATTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAatat                                                                           >  1:1115154/1‑76 (MQ=255)
  tGAGTAGCAATGGACAGGTAAAATTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAAtatta                                                                         <  1:630705/76‑1 (MQ=255)
  tGAGTAGCAATGGACAGGTAAAATTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAAtatta                                                                         >  2:444735/1‑76 (MQ=255)
      tAGCAATGGACAGGTAAAATTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATc                                                                     >  1:205770/1‑76 (MQ=255)
      tAGCAATGGACAGGTAAAATTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATc                                                                     >  1:1496687/1‑76 (MQ=255)
       aGCAATGGACAGGTAAAATTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATcc                                                                    <  1:13298/76‑1 (MQ=255)
          aaTGGACAGGTAAAATTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCcatc                                                                 >  2:130841/1‑76 (MQ=255)
            tGGACAGGTAAAATTGCTGAAGTTTTTGGTGTAATAGCTTGTTCTATTTTTCTGTAATCAATATTATATCcatcat                                                               >  2:2030953/1‑76 (MQ=255)
             ggACAGGTAAAATTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCAtt                                                              <  2:1081647/76‑1 (MQ=255)
             ggACAGGTAAAATTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCAtt                                                              <  2:1081635/76‑1 (MQ=255)
                    tAAAATTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAata                                                       >  1:1355723/1‑76 (MQ=255)
                     aaaaTTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAatat                                                      >  2:628732/1‑76 (MQ=255)
                     aaaaTTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAatat                                                      >  2:162339/1‑76 (MQ=255)
                        aTTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATcaa                                                   <  2:952514/76‑1 (MQ=255)
                          tGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATcaaca                                                 <  2:913885/76‑1 (MQ=255)
                              gAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGa                                             <  1:1332611/76‑1 (MQ=255)
                               aaGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGAc                                            <  1:1471822/76‑1 (MQ=255)
                               aaGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGAc                                            <  1:1471821/76‑1 (MQ=255)
                                aGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACa                                           >  1:1716503/1‑76 (MQ=255)
                                aGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACa                                           >  1:201428/1‑76 (MQ=255)
                                aGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACa                                           <  1:939568/76‑1 (MQ=255)
                                 gTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACan                                          <  1:1837565/76‑2 (MQ=255)
                                 gTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAg                                          <  1:540028/76‑1 (MQ=255)
                                   ttttGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGc                                        >  1:861747/1‑76 (MQ=255)
                                    tttGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCg                                       >  1:1213877/1‑76 (MQ=255)
                                    tttGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCg                                       >  1:1213899/1‑76 (MQ=255)
                                     ttGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGt                                      <  2:1310669/76‑1 (MQ=255)
                                          gtAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGtc                                 >  1:1486651/1‑76 (MQ=255)
                                           tAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGtcc                                <  1:120755/76‑1 (MQ=255)
                                             aTAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGtccat                              >  2:245264/1‑76 (MQ=255)
                                                  cTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACATAGACAGGCGTAAGtccattccat                         >  2:1320103/1‑76 (MQ=255)
                                                  cTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGtccattccat                         >  2:1373318/1‑76 (MQ=255)
                                                  cTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGtccattccat                         >  2:1373345/1‑76 (MQ=255)
                                                  cTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGtccattccat                         >  1:1654906/1‑76 (MQ=255)
                                                   tGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATa                        >  1:279705/1‑76 (MQ=255)
                                                    gTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATaa                       <  2:677794/76‑1 (MQ=255)
                                                     ttCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATaaa                      >  2:1059792/1‑76 (MQ=255)
                                                       cTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAAt                    <  1:1735183/76‑1 (MQ=255)
                                                        tATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATa                   >  2:1542444/1‑76 (MQ=255)
                                                        tATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATa                   >  2:456922/1‑76 (MQ=255)
                                                          tttttCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATAGc                 >  1:118227/1‑76 (MQ=255)
                                                           ttttCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATAGCa                >  2:95770/1‑76 (MQ=255)
                                                               cTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATAGCAtgtg            >  2:1178288/1‑76 (MQ=255)
                                                                tGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATAGCATGTGa           <  1:506310/76‑1 (MQ=255)
                                                                  tAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATAGCATGTGAAg         >  2:973493/1‑76 (MQ=255)
                                                                   aaTCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATAGCATGTGAAGt        >  1:741336/1‑76 (MQ=255)
                                                                      cAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATAGCATGTGAAGTTGc     <  2:332854/76‑1 (MQ=255)
                                                                      cAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATAGCATGTGAAGTTGc     >  1:1945280/1‑76 (MQ=255)
                                                                      cAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATAGCATGTGAAGTTGc     <  1:1716264/76‑1 (MQ=255)
                                                                         tattatATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATAGCATGTGAAGTTGcaac  >  1:1703141/1‑76 (MQ=255)
                                                                          |                                                                          
GTTGAGTAGCAATGGACAGGTAAAATTGCTGAAGTTTTTGGTGTAATAGCCTGTTCTATTTTTCTGTAATCAATATTATATCCATCATTTTCAATATCAACAAAGACAGGCGTAAGTCCATTCCATAAAATAGCATGTGAAGTTGCAAC  >  REL606/2040980‑2041128

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