New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_017635 985116 =46 (1.140)9 (0.320) 8/68 1.1 23.8% coding (1049/1392 nt) ECW_RS04805 phage tail protein
?NC_017635 987945 = 26 (0.930)coding (37/504 nt) ECW_RS04830 phage tail protein

CGCTTACCAATTTGAGTGGAAAGGATGTAGCTGGTCTTCTCGCATACCTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCGAGGAGCTAATGCCAA  >  NC_017635/985068‑985164
                                                |                                                
cGCTTACCAATTTGAGTGGAAAGGATGTAGCTGGTCTTCTCGCATACCtt                                                 <  2:1235207/50‑1 (MQ=255)
aaccggggatgatgccagttagccgttatctctgcattttcataaacgtt                                                 <  1:872989‑M1/2‑1 (MQ=255)
aaccggggatgatgccagttagccgttatctctgcattttcataaacgtt                                                 >  2:1907524‑M1/49‑50 (MQ=255)
aaccggggatgatgccagttagccgttatctctgcattttcataaacgtt                                                 >  1:2049558‑M1/49‑50 (MQ=255)
aaccggggatgatgccagttagccgttatctctgcattttcataaacgtt                                                 >  2:695092‑M1/49‑50 (MQ=255)
aaccggggatgatgccagttagccgttatctctgcattttcataaacgtt                                                 >  1:155643‑M1/49‑50 (MQ=255)
 gCTTACCAATTTGAGTGGAAAGGATGTAGCTGGTCTTCTCGCATACCTTg                                                >  2:1996970/1‑50 (MQ=255)
  cTTACCAATTTGAGTGGAAAGGATGTAGCTGGTCTTCTCGCATACCTTgg                                               >  2:1977717/1‑50 (MQ=255)
  cTTACCAATTTGAGTGGAAAGGATGTAGCTGGTCTTCTCGCATACCTTgg                                               >  1:1597330/1‑50 (MQ=255)
   ttACCAATTTGAGTGGAAAGGATGTAGCTGGTCTTCTCGCATACCTTGGt                                              >  1:1432333/1‑50 (MQ=255)
    tACCAATTTGAGTGGAAAGGATGTAGCTGGTCTTCTCGCATACCTTGGtt                                             >  1:1968719/1‑50 (MQ=255)
    ggggatgatgccagttagccgttatctctgcattttcataaacgTTGGtt                                             <  2:1095065‑M1/6‑1 (MQ=255)
     aCCAATTTGAGTGGAAAGGATGTAGCTGGTCTTCTCGCATACCTTGGttt                                            >  2:521138/1‑50 (MQ=255)
      ggatgatgccagttagccgttatctctgcattttcataaacgTTGGTTTg                                           <  2:804070‑M1/8‑1 (MQ=255)
      ggatgatgccagttagccgttatctctgcattttcataaacgTTGGTTTg                                           <  1:1316614‑M1/8‑1 (MQ=255)
      ccAATTTGAGTGGAAAGGATGTAGCTGGTCTTCTCGCATACCTTGGttt                                            <  1:1916468/49‑1 (MQ=255)
        atgatgccagttagccgttatctctgcattttcataaacgTTGGTTTggg                                         >  1:1646519‑M1/41‑50 (MQ=255)
        atgatgccagttagccgttatctctgcattttcataaacgTTGGTTTggg                                         >  1:2060804‑M1/41‑50 (MQ=255)
        aaTTTGAGTGGAAAGGATGTAGCTGGTCTTCTCGCATACCTTGGTTTgg                                          >  1:1459634/1‑49 (MQ=255)
         aTTTGAGTGGAAAGGATGTAGCTGGTCTTCTCGCATACCTTGGTTTgg                                          >  2:1920885/1‑48 (MQ=255)
             gAGTGGAAAGGATGTAGCTGGTCTTCTCGCATACCTTGGTTTGGGAGAAg                                    >  2:1576962/1‑50 (MQ=255)
              ccagttagccgttatctctgcattttcataaacgTTGGTTTGGGAGAAg                                    >  1:689019‑M1/35‑49 (MQ=255)
              ccagttagccgttatctctgcattttcataaacgTTGGTTTGGGAGAAGc                                   <  2:305983‑M1/16‑1 (MQ=255)
              aGTGGAAAGGATGTAGCTGGTCTTCTCGCATACCTTGGTTTGGGAGAAg                                    <  2:1904043/49‑1 (MQ=255)
              aGTGGAAAGGATGTAGCTGGTCTTCTCGCATACCTTGGTTTGGGAGAAGc                                   >  2:527221/1‑50 (MQ=255)
              aGTGGAAAGGATGTAGCTGGTCTTCTCGCATACCTTGGTTTGGGAGAAGc                                   >  1:643258/1‑50 (MQ=255)
               gTGGAAAGGATGTAGCTGGTCTTCTCGCATACCTTGGTTTGGGAGAAggc                                  <  2:914426/50‑3 (MQ=255)
               gTGGAAAGGATGTAGCTGGTCTTCTCGCATACCTTGGTTTGGGAGAAggc                                  >  2:993120/1‑48 (MQ=255)
                 ggAAAGGATGTAGCTGGTCTTCTCGCATACCTTGGTTTGGGAGAAGc                                   >  2:1040798/1‑47 (MQ=255)
                   aaaGGATGTAGCTGGTCTTCTCGCATACCTTGGTTTGGGAGAAGCCAgtt                              >  2:1541446/1‑48 (MQ=255)
                   aaaGGATGTAGCTGGTCTTCTCGCATACCTTGGTTTGGGAGAAGCCAGat                              <  2:28387/50‑1 (MQ=255)
                    aaGGATGTAGCTGGTCTTCTCGCATACCTTGGTTTGGGAGAAGCCAGata                             <  2:347862/50‑1 (MQ=255)
                       gATGTAGCTGGTCTTCTCGCATACCTTGGTTTGGGAGAAGCCAGATATGt                          <  1:515916/50‑1 (MQ=255)
                       gATGTAGCTGGTCTTCTCGCATACCTTGGTTTGGGAGAAGCCAGATATGt                          <  1:618389/50‑1 (MQ=255)
                        gttatctctgcattttcataaacgTTGGTTTGGGAGAAGCCAGATATGt                          <  1:1661173‑M1/25‑1 (MQ=255)
                        aTGTAGCTGGTCTTCTCGCATACCTTGGTTTGGGAGAAGCCAGATATGt                          >  2:1511526/1‑49 (MQ=255)
                         ttatctctgcattttcataaacgTTGGTTTGGGAGAAGCCAGATATGTCa                        <  2:1630933‑M1/27‑1 (MQ=255)
                         ttatctctgcattttcataaacgTTGGTTTGGGAGAAGCCAGATATGTCa                        <  1:214483‑M1/27‑1 (MQ=255)
                          gTAGCTGGTCTTCTCGCATACCTTGGTTTGGGAGAAGCCAGATATGTCAt                       >  2:955006/1‑50 (MQ=255)
                          gTAGCTGGTCTTCTCGCATACCTTGGTTTGGGAGAAGCCAGATATGTCAt                       >  1:1507501/1‑50 (MQ=255)
                           tAGCTGGTCTTCTCGCATACCTTGGTTTGGGAGAAGCCAGATATGTCAtt                      >  1:1180880/1‑50 (MQ=255)
                           tAGCTGGTCTTCTCGCATACCTTGGTTTGGGAGAAGCCAGATATGTCAtt                      >  2:85635/1‑50 (MQ=255)
                           atctctgcattttcataaacgTTGGTTTGGGAGAAGCCAGATATGTCAtt                      <  1:1571814‑M1/29‑1 (MQ=255)
                              cTGGTCTTCTCGCATACCTTGGTTTGGGAGAAGCCAGATATGTCATTCAg                   >  2:774320/1‑50 (MQ=255)
                               tGGTCTTCTCGCATACCTTGGTTTGGGAGAAGCCAGATATGTCATTCAGc                  >  2:126735/1‑50 (MQ=255)
                                ggTCTTCTCGCATACCTTGGTTTGGGAGAAGCCAGATATGTCATTCAGc                  >  2:479103/1‑49 (MQ=255)
                                  tcttctCGCATACCTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCgag               >  2:425009/1‑50 (MQ=255)
                                  tcttctCGCATACCTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCgag               <  1:1260600/50‑1 (MQ=255)
                                  tcttctCGCATACCTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCgag               >  1:142348/1‑50 (MQ=255)
                                   cttctCGCATACCTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCgag               >  2:859166/1‑49 (MQ=255)
                                   attttcataaacgTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCgagg              >  1:1082693‑M1/14‑50 (MQ=255)
                                    ttttcataaacgTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCgagga             >  1:204239‑M1/13‑50 (MQ=255)
                                    ttctCGCATACCTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCgag               >  2:888148/1‑48 (MQ=255)
                                     tctcGCATACCTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCgaggag            >  2:1839972/1‑50 (MQ=255)
                                      ctcGCATACCTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCgaggag            <  2:1954424/49‑1 (MQ=255)
                                       tcGCATACCTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCGAGGAGc           >  1:1696471/1‑49 (MQ=255)
                                       tcGCATACCTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCGAGGAGc           >  2:1706397/1‑49 (MQ=255)
                                        cataaacgTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCGAGGAGCt          <  1:1246922‑M1/41‑1 (MQ=255)
                                         gCATACCTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCGAGGAGCTaa        >  2:1738676/1‑50 (MQ=255)
                                         ataaacgTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCGAGGAGCTaa        <  2:859334‑M1/43‑1 (MQ=255)
                                           aTACCTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCGAGGAGCTAATg      >  2:1796809/1‑50 (MQ=255)
                                           aTACCTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCGAGGAGCTAATg      >  2:926475/1‑50 (MQ=255)
                                            tACCTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCGAGGAGCTAATGc     >  1:1496224/1‑50 (MQ=255)
                                             aCCTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCGAGGAGCTAATGc     >  2:636978/1‑49 (MQ=255)
                                              ccTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCGAGGAGCTAATGCCa   <  1:400731/50‑1 (MQ=255)
                                               cTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCGAGGAGCTAATGCCaa  >  2:2066366/1‑50 (MQ=255)
                                                |                                                
CGCTTACCAATTTGAGTGGAAAGGATGTAGCTGGTCTTCTCGCATACCTTGGTTTGGGAGAAGCCAGATATGTCATTCAGCGAGGAGCTAATGCCAA  >  NC_017635/985068‑985164

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