New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2037554 =68 (1.090)3 (0.050) 3/136 6.6 4.9% coding (577/1176 nt) wzy O antigen polymerase (Wzy)
?REL606 2037718 = 54 (0.950)coding (413/1176 nt) wzy O antigen polymerase (Wzy)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACAAATAACGCTGAAATTACTCCGCCTCCGC  >  REL606/2037647‑2037793
                                                                       |                                                                           
atatCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACAtg                                                                         >  1:244502/1‑76 (MQ=255)
    ccAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGtt                                                                     >  2:1908880/1‑76 (MQ=255)
    ccAACCACGCTGAAGAAAAGAAATAATCCACAACAGGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGtt                                                                     >  2:1802400/1‑76 (MQ=255)
     cAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTa                                                                    <  1:1582164/76‑1 (MQ=255)
      aaCCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAg                                                                   <  1:1199996/76‑1 (MQ=255)
       aCCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGa                                                                  <  2:1307598/76‑1 (MQ=255)
          aCGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAat                                                               >  2:919848/1‑76 (MQ=255)
          aCGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAat                                                               <  1:1821698/76‑1 (MQ=255)
          aCGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAat                                                               >  2:1603540/1‑76 (MQ=255)
           cGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAata                                                              >  1:1315034/1‑76 (MQ=255)
            gCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAatat                                                             <  1:677500/76‑1 (MQ=255)
            gCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAatat                                                             <  2:1407522/76‑1 (MQ=255)
               gaagaaAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTaa                                                          <  1:1776616/76‑1 (MQ=255)
               gaagaaAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTaa                                                          >  2:1019713/1‑76 (MQ=255)
                aagaaaagaaATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTaaa                                                         <  1:159342/76‑1 (MQ=255)
                aagaaaagaaATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTaaa                                                         >  1:154293/1‑76 (MQ=255)
                 agaaaagaaaTAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTaaaa                                                        <  1:20467/76‑1 (MQ=255)
                   aaaagaaaTAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACa                                                      >  2:919914/1‑76 (MQ=255)
                   aaaagaaaTAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACa                                                      <  1:1866974/76‑1 (MQ=255)
                       gaaaTAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTa                                                  <  2:414584/76‑1 (MQ=255)
                         aataatCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACa                                                >  1:1790823/1‑76 (MQ=255)
                            aatCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCa                                             <  2:108262/76‑1 (MQ=255)
                            aatCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCa                                             >  1:1928687/1‑76 (MQ=255)
                             atCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCaa                                            <  1:167352/76‑1 (MQ=255)
                             atCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCaa                                            <  1:199140/76‑1 (MQ=255)
                             atCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCaa                                            <  1:199132/76‑1 (MQ=255)
                             atCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCaa                                            >  1:1708247/1‑76 (MQ=255)
                                 acaacaTGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACAtt                                        >  2:1554744/1‑76 (MQ=255)
                                 acaacaTGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACAtt                                        >  1:965813/1‑76 (MQ=255)
                                 acaacaTGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAACATTAAAACAAGAACATAAACAtc                                        >  2:978810/1‑75 (MQ=255)
                                 acaacaTAATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACAtt                                        >  1:1282174/1‑76 (MQ=255)
                                  caacaTGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTc                                       <  1:534238/76‑1 (MQ=255)
                                   aacaTGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCa                                      >  2:756257/1‑76 (MQ=255)
                                   aacaTGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCa                                      >  2:646163/1‑76 (MQ=255)
                                   aacaTGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCa                                      <  2:1375072/76‑1 (MQ=255)
                                    acaTGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAg                                     >  1:30537/1‑76 (MQ=255)
                                     caTGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGt                                    >  1:427765/1‑76 (MQ=255)
                                      aTGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTa                                   >  1:1576996/1‑76 (MQ=255)
                                       tGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTaa                                  <  1:1636445/76‑1 (MQ=255)
                                         ataataATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATa                                <  2:1332630/76‑1 (MQ=255)
                                         ataataATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATa                                <  2:1502549/76‑1 (MQ=255)
                                          taataaTGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATAc                               <  1:842595/76‑1 (MQ=255)
                                             taatGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACaaa                            >  1:1740307/1‑76 (MQ=255)
                                             taatGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACaaa                            <  1:1329262/76‑1 (MQ=255)
                                               atGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACAAATa                          <  1:1143016/76‑1 (MQ=255)
                                               atGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACAAATa                          >  2:1008968/1‑76 (MQ=255)
                                                    cTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACAAATAACGCt                     <  1:1802400/76‑1 (MQ=255)
                                                    cTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACAAATAACGCt                     <  2:807753/76‑1 (MQ=255)
                                                       gTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACAAATAACGCTGaa                  >  1:1916963/1‑76 (MQ=255)
                                                       gTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACAAATAACGCTGaa                  <  1:1308087/76‑1 (MQ=255)
                                                           ttacttGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACAAATAACGCTGAAATTa              >  2:1699906/1‑76 (MQ=255)
                                                            tacttGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACAAATAACGCTGAAATTAc             >  2:644148/1‑76 (MQ=255)
                                                            tacttGTATTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACAAATAACGCTGAAATTAc             >  1:1288233/1‑76 (MQ=255)
                                                               ttGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACAAATAACGCTGAAATTActcc          <  2:1751961/76‑1 (MQ=255)
                                                                   tttAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACAAATAACGCTGAAATTActccgcct      >  2:1395294/1‑76 (MQ=255)
                                                                   tttAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACAAATAACGCTGAAATTActccgcct      >  1:1011844/1‑76 (MQ=255)
                                                                   tttAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACAAATAACGCTGAAATTActccgcct      >  1:736391/1‑76 (MQ=255)
                                                                   tttAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACAAATAACGCTGAAATTActccgcct      >  1:1901678/1‑76 (MQ=255)
                                                                     tAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACAAATAACGCTGAAATTActccgcctcc    >  2:347156/1‑76 (MQ=255)
                                                                      aaCATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACAAATAACGCTGAAATTActccgcctccg   >  2:1671989/1‑76 (MQ=255)
                                                                       aCATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACAAATAACGCTGAAATTActccgcctccgc  >  2:1293093/1‑76 (MQ=255)
                                                                       |                                                                           
ATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGTAATACAAATAACGCTGAAATTACTCCGCCTCCGC  >  REL606/2037647‑2037793

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