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

CATGGAACAGGGACACCAACGGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGCAGAGATAACGGCTAACTGGCATCATCCCCGGTT  >  NC_017635/987896‑987992
                                                 |                                               
cATGGAACAGGGACACCAACGGGTAATGCAGAGCCTTCTCCCAAACCAAc                                                 >  2:503193/1‑50 (MQ=255)
 aTGGAACAGGGACACCAACGGGTAATGCAGAGCCTTCTCCCAAACCAACg                                                >  1:9832/1‑50 (MQ=255)
 aTGGAACAGGGACACCAACGGGTAATGCAGAGCCTTCTCCCAAACCAACg                                                >  2:316096/1‑50 (MQ=255)
 aTGGAACAGGGACACCAACGGGTAATGCAGAGCCTTCTCCCAAACCAACg                                                >  1:1402719/1‑50 (MQ=255)
     aaCAGGGACACCAACGGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTa                                            >  1:504508/1‑50 (MQ=255)
      ttagctcctcgctgaatgacatatctggcttctcccaaaccaaCGTTTAt                                           <  2:859334‑M2/7‑1 (MQ=255)
      aCAGGGACACCAACGGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAt                                           >  1:1366826/1‑50 (MQ=255)
        agctcctcgctgaatgacatatctggcttctcccaaaccaaCGTTTATg                                          <  1:1246922‑M2/8‑1 (MQ=255)
           tcctcgctgaatgacatatctggcttctcccaaaccaaCGTTTATGaaaa                                      >  1:204239‑M2/39‑50 (MQ=255)
            cctcgctgaatgacatatctggcttctcccaaaccaaCGTTTATGAAAAt                                     >  1:1082693‑M2/38‑50 (MQ=255)
                cAACGGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                    >  1:1287504/1‑47 (MQ=255)
                 aaCGGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGCaga                                <  2:341514/50‑1 (MQ=255)
                 aaCGGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGCaga                                <  1:1384907/50‑1 (MQ=255)
                 aaCGGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGCaga                                <  1:850118/50‑1 (MQ=255)
                    gggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAGATGCAgaga                              >  2:1963785/1‑49 (MQ=255)
                    gggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGCAgaga                              >  2:1617841/1‑49 (MQ=255)
                    gggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGCAgaga                              >  2:1769898/1‑49 (MQ=255)
                    gggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGCAgaga                              >  1:569188/1‑49 (MQ=255)
                    aatgacatatctggcttctcccaaaccaaCGTTTATGAAAATGCAGAGAt                             <  1:1571814‑M2/21‑1 (MQ=255)
                     ggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGCAGAGAt                             >  1:2036196/1‑49 (MQ=255)
                      tgacatatctggcttctcccaaaccaaCGTTTATGAAAATGCAGAGATaa                           <  2:1630933‑M2/23‑1 (MQ=255)
                      tgacatatctggcttctcccaaaccaaCGTTTATGAAAATGCAGAGATaa                           <  1:214483‑M2/23‑1 (MQ=255)
                      gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGCAGAGATaa                           >  2:1542257/1‑50 (MQ=255)
                        acatatctggcttctcccaaaccaaCGTTTATGAAAATGCAGAGATAAc                          <  1:1661173‑M2/24‑1 (MQ=255)
                        aaTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGCAGAGATAACg                         >  1:686551/1‑50 (MQ=255)
                        aaTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGCAGAGATAACg                         >  1:87505/1‑50 (MQ=255)
                        aaTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGCAGAGATAACg                         >  2:1651323/1‑50 (MQ=255)
                        aaTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGCAGAGATAACg                         >  2:1356101/1‑50 (MQ=255)
                         aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGCAGAGATAACg                         <  2:2042182/49‑1 (MQ=255)
                             agagCCTTCTCCCAAACCAACGTTTATGAAAATGCAGAGATAACGGCTaa                    <  1:1243137/50‑1 (MQ=255)
                             agagCCTTCTCCCAAACCAACGTTTATGAAAATGCAGAGATAACGGCTaa                    >  2:761000/1‑50 (MQ=255)
                             agagCCTTCTCCCAAACCAACGTTTATGAAAATGCAGAGATAACGGCTaa                    >  1:1632167/1‑50 (MQ=255)
                               agCCTTCTCCCAAACCAACGTTTATGAAAATGCAGAGATAACGGCTAACt                  >  2:356654/1‑50 (MQ=255)
                               agCCTTCTCCCAAACCAACGTTTATGAAAATGCAGAGATAACGGCTAACt                  <  2:494329/50‑1 (MQ=255)
                                 gcttctcccaaaccaaCGTTTATGAAAATGCAGAGATAACGGCTAACTgg                <  2:305983‑M2/34‑1 (MQ=255)
                                  cttctcccaaaccaaCGTTTATGAAAATGCAGAGATAACGGCTAACTgg                >  1:689019‑M2/16‑49 (MQ=255)
                                       cccaaaccaaCGTTTATGAAAATGCAGAGATAACGGCTAACTGGcatcat          >  1:2060804‑M2/11‑50 (MQ=255)
                                       cccaaaccaaCGTTTATGAAAATGCAGAGATAACGGCTAACTGGcatcat          >  1:1646519‑M2/11‑50 (MQ=255)
                                         caaaccaaCGTTTATGAAAATGCAGAGATAACGGCTAACTGGCATCATcc        <  2:804070‑M2/42‑1 (MQ=255)
                                         caaaccaaCGTTTATGAAAATGCAGAGATAACGGCTAACTGGCATCATcc        <  1:1316614‑M2/42‑1 (MQ=255)
                                           aaccaaCGTTTATGAAAATGCAGAGATAACGGCTAACTGGCATCATcccc      <  2:1095065‑M2/44‑1 (MQ=255)
                                               aaCGTTTATGAAAATGCAGAGATAACGGCTAACTGGCATCATCCCCGGtt  <  1:872989‑M2/48‑1 (MQ=255)
                                               aaCGTTTATGAAAATGCAGAGATAACGGCTAACTGGCATCATCCCCGGtt  >  2:1907524‑M2/3‑50 (MQ=255)
                                               aaCGTTTATGAAAATGCAGAGATAACGGCTAACTGGCATCATCCCCGGtt  >  1:2049558‑M2/3‑50 (MQ=255)
                                               aaCGTTTATGAAAATGCAGAGATAACGGCTAACTGGCATCATCCCCGGtt  >  2:695092‑M2/3‑50 (MQ=255)
                                               aaCGTTTATGAAAATGCAGAGATAACGGCTAACTGGCATCATCCCCGGtt  >  1:155643‑M2/3‑50 (MQ=255)
                                                 |                                               
CATGGAACAGGGACACCAACGGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGCAGAGATAACGGCTAACTGGCATCATCCCCGGTT  >  NC_017635/987896‑987992

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