New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 1432761NA (NA)7 (0.060) 5/170 NT NA coding (351/576 nt) tfaR Rac prophage; putative tail fiber assembly protein
?NC_000913 = 1432785 NA (NA)coding (375/576 nt) tfaR Rac prophage; putative tail fiber assembly protein

TATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGAAAAACTGTTCCGGATCCGGGAGGCGGAAGAAACAAAAAAC  >  NC_000913/1432709‑1432812
                                                                            |                           
tATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGAAAAACTGTTCCGCCTCCCGGATCCGgaacagttttt      <  2:2685149/100‑9 (MQ=1)
    gcaggAAGAGAAGTATAAGAAGTGGAACGGCACAGCCTGGGTGAAGTATACGGAAGCAGAAAAACTGTTCCGGATCCGGGATGaaaaaaaaaaaaaaaa   <  1:1664399/96‑1 (MQ=1)
      cgtaaGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGAAGAAAAACTGTTCCGGATCCGGGAGGCGGAAGAAACAAAAaac  <  1:1278752/95‑1 (MQ=0)
                                                                            |                           
TATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGCACAGCCTGGGTGAAGGATACGGAAGCAGAAAAACTGTTCCGGATCCGGGAGGCGGAAGAAACAAAAAAC  >  NC_000913/1432709‑1432812

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