New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_017635 2314099 =13 (0.320)17 (0.610) 13/68 0.6 54.0% coding (38/414 nt) ECW_RS28245 hypothetical protein
?NC_017635 2315940 = 20 (0.720)pseudogene (1/459 nt) ECW_RS28250 phage tail protein

AGAAGAACGGAACGCCAATAGCTGGAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCTTTTCCACTCAAATGAGTCAGCGTC  >  NC_017635/2315895‑2315989
                                             |                                                 
agaagaACGGAACGCCAATAGCTGGAGCACCTTCTCCCAAACCAAGGTAt                                               <  2:751001/50‑1 (MQ=255)
agaagaACGGAACGCCAATAGCTGGAGCACCTTCTCCCAAACCAAGGTAt                                               <  1:283615/50‑1 (MQ=255)
agaagaACGGAACGCCAATAGCTGGAGCACCTTCTCCCAAACCAAGGTAt                                               >  1:212750/1‑50 (MQ=255)
    gaacggaacgCCAATAGCTGGAGCACCTTCTCCCAAACCAAGGTATGCga                                           <  2:1995906/50‑1 (MQ=255)
    gaacggaacgCCAATAGCTGGAGCACCTTCTCCCAAACCAAGGTATGCga                                           <  2:1750583/50‑1 (MQ=255)
      cgttacccacatcccgttttgccgcttctcccaaaccaaGGTATGCgag                                          >  1:798802‑M2/40‑49 (MQ=255)
       gttacccacatcccgttttgccgcttctcccaaaccaaGGTATGCgag                                          >  1:745072‑M2/39‑48 (MQ=255)
        ttacccacatcccgttttgccgcttctcccaaaccaaGGTATGCGagaag                                       <  2:667953‑M2/13‑1 (MQ=255)
         tacccacatcccgttttgccgcttctcccaaaccaaGGTATGCGagaaga                                      >  2:1208621‑M2/37‑50 (MQ=255)
           cccacatcccgttttgccgcttctcccaaaccaaGGTATGCGAGAAGAc                                     >  2:67557‑M2/35‑49 (MQ=255)
           acgCCAATAGCTGGAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGAcc                                    <  1:1411849/50‑1 (MQ=255)
            cgCCAATAGCTGGAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCa                                   <  1:1421064/50‑1 (MQ=255)
            ccacatcccgttttgccgcttctcccaaaccaaGGTATGCGAGAAGACCa                                   <  2:1394518‑M2/17‑1 (MQ=255)
            ccacatcccgttttgccgcttctcccaaaccaaGGTATGCGAGAAGACCa                                   <  1:1935009‑M2/17‑1 (MQ=255)
            ccacatcccgttttgccgcttctcccaaaccaaGGTATGCGAGAAGACCa                                   <  2:1401892‑M2/17‑1 (MQ=255)
             cacatcccgttttgccgcttctcccaaaccaaGGTATGCGAGAAGACCAg                                  >  1:1060392‑M2/33‑50 (MQ=255)
             cacatcccgttttgccgcttctcccaaaccaaGGTATGCGAGAAGACCAg                                  >  1:633821‑M2/33‑50 (MQ=255)
              ccAATAGCTGGAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGc                                 <  1:1047356/50‑1 (MQ=255)
               catcccgttttgccgcttctcccaaaccaaGGTATGCGAGAAGACCAGcc                                <  1:89740‑M2/20‑1 (MQ=255)
               catcccgttttgccgcttctcccaaaccaaGGTATGCGAGAAGACCAGc                                 >  2:704051‑M2/31‑49 (MQ=255)
                 aTAGCTGGAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCcac                              <  1:734039/50‑1 (MQ=255)
                    gCTGGAGCACCTTCTCCAAAACAAAGGTATGCGAGAAGACCAGCCACATc                           <  2:269553/50‑1 (MQ=18)
                     cTGGAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCt                          <  2:277974/50‑1 (MQ=255)
                     cTGGAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCt                          <  1:1462674/50‑1 (MQ=255)
                      ttttgccgcttctcccaaaccaaGGTATGCGAGAAGACCAGCCACATCtt                         <  1:1030646‑M2/27‑1 (MQ=255)
                      ttttgccgcttctcccaaaccaaGGTATGCGAGAAGACCAGCCACATCtt                         >  2:1352630‑M2/24‑50 (MQ=255)
                      ttttgccgcttctcccaaaccaaGGTATGCGAGAAGACCAGCCACATCtt                         >  1:2032255‑M2/24‑50 (MQ=255)
                       ggAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCttt                        >  2:962680/1‑50 (MQ=255)
                       ggAGCACCTTCTCCCAAACCAAGGAATGCGAGAAGACCAGCCACATCttt                        >  2:642542/1‑50 (MQ=255)
                        ttgccgcttctcccaaaccaaGGTATGCGAGAAGACCAGCCACATCtttt                       >  2:448213‑M2/22‑50 (MQ=255)
                        gAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCtttt                       >  2:371237/1‑50 (MQ=255)
                        gAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCtttt                       >  2:1707333/1‑50 (MQ=255)
                          gCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCTTTTcc                     <  1:1659267/50‑1 (MQ=255)
                          gCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCTTTTcc                     <  1:101414/50‑1 (MQ=255)
                           ccgcttcccccaaaccaaGGTATGCGAGAAGACCAGCCACATCTTTTcc                     <  2:1794654‑M2/31‑1 (MQ=255)
                           ccccTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCTTTTCCa                    <  2:1826922/48‑1 (MQ=255)
                           ccccTTCCCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCTTTTCCa                    <  2:1569146/48‑1 (MQ=17)
                                 ctcccaaaccaaGGTATGCGAGAAGACCAGCCACATCTTTTCCACTCaa               <  1:128619‑M2/37‑1 (MQ=255)
                                   cccaaaccaaGGTATGCGAGAAGACCAGCCACATCTTTTCCACTCAAAt             >  1:1499303‑M2/11‑49 (MQ=255)
                                    ccaacccaaGGTATGCGAGAAGACCAGCCACATCTTTTCCACTCAAATGa           <  2:992787‑M2/41‑1 (MQ=255)
                                    ccaaaccaaGTTATGCGAGAAGACCAGCCACATCTTTTCCACTCAAATGa           <  2:947669‑M2/41‑1 (MQ=255)
                                    ccaaaccaaGGTATGCGAGAAGACCAGCTACATCTTTTCCACTCAAATg            <  1:725154‑M2/40‑1 (MQ=255)
                                    ccaaaccaaGGTATGCGAGAAGACCAGCCACATCTTTTCCACTCAAATg            <  1:348967‑M2/40‑1 (MQ=255)
                                      aaaccaaGGTATGCGAGAAGACCAGCCACATCTTTTCCACTCAAATGAGt         <  1:326756‑M2/43‑1 (MQ=255)
                                       aaccaaGTTATGCGAGAAGACCAGCCACATCTTTTCCACTCAAATGAGTc        <  2:1775547‑M2/44‑1 (MQ=255)
                                       aaccaaGGTATGCGAGAAGACCAGCCACATCTTTTCCACTCAAATGAGTc        <  2:1214304‑M2/44‑1 (MQ=255)
                                       aaccaaGGTATGCGAGAAGACCAGCCACATCTTTTCCACTCAAATGAGTc        <  1:925239‑M2/44‑1 (MQ=255)
                                          caaGGTATGCGAGAAGACCAGCCACATCTTTTCCACTCAAATGAGTCAGc     <  2:1745604‑M2/47‑1 (MQ=255)
                                          caaGGTATGCGAGAAGACCAGCCACATCTTTTCCACTCAAATGAGTCAGc     <  1:538800‑M2/47‑1 (MQ=255)
                                           aaGGTATGCGAGAAGACCAGCCACATCTTTTCCACTCAAATGAGTCAGCg    >  2:1809634‑M2/3‑50 (MQ=255)
                                             ggTATGCGAGAAGACCAGCCACATCTTTTCCACTCAAATGAGTCAGCGtc  >  1:2051788/1‑50 (MQ=255)
                                             |                                                 
AGAAGAACGGAACGCCAATAGCTGGAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCTTTTCCACTCAAATGAGTCAGCGTC  >  NC_017635/2315895‑2315989

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