New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =14 (0.420)4 (0.150) 3/128 2.7 25.7% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 12 (0.460)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

AATACCACCAATCATAACCGGTATAGCTATGAAACCGGTGCGTCCTAATGTAGCCGTTGCTGCATTGAGTTTTGCTCCTTCTCCCAAACCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1207881‑1207790
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTG  >  NC_000913/1209619‑1209694
                                                                                                                                                                        
AATACCACCAATCATAACCGGTATAGCTATGAAACCGGTGCGTCCTAATGTAGCCGTTGCTGCATTGAGTTTTGCTCCT                                                                                           <  1:239683/79‑1
  TACCACCAATCATAACCGGTATAGCTATGAAACCGGTGCGTCCTAATGTAGCCGTTGCTGCATTGAGTTTTGCTCCTTCTC                                                                                       <  1:239678/81‑1
      ACCAATCATAACCGGTATAGCTATGAAACCGGTGCGTCCTAATGTAGCCGTTGCTGCATTGAGTTTTGCTCCTTCTCCCA                                                                                    >  1:239662/1‑80
                      ATAGCTATGAAACCGGTGCGTCCTAATGTAGCCGTTGCTGCATTGAGTTTTGCTCCTTCTCCCAAACCAACGTTTATGAA                                                                    >  2:239889/1‑80
                                     gcGCGTCCTAATGTAGCCGTTGCTGCATTGAGTTTTGCTCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAAC                                                      >  2:239881/3‑79
                                                       GTTGCTGCATTGAGTTTTGCTCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCT                                  >  1:239891/1‑81
                                                           CTGCATTGAGTTTTGCTCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTT                              >  2:239896/1‑81
                                                                             CTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAC            <  1:239871/81‑1
                                                                                  CCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCT       >  2:239894/1‑81
                                                                                     AAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAT    >  2:239917/1‑81
                                                                                       ACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTG  <  1:239872/81‑1
                                                                                                                                                                        
AATACCACCAATCATAACCGGTATAGCTATGAAACCGGTGCGTCCTAATGTAGCCGTTGCTGCATTGAGTTTTGCTCCTTCTCCCAAACCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1207881‑1207790
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTG  >  NC_000913/1209619‑1209694

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