New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2768888 =86 (0.990)3 (0.040) 3/250 NT 3.9% coding (224/822 nt) yfjQ CP4‑57 prophage; uncharacterized protein
?NC_000913 2768933 = 65 (0.770)coding (269/822 nt) yfjQ CP4‑57 prophage; uncharacterized protein
Rejected: Frequency below/above cutoff threshold.

GCACACCGCCCTGAACATGCCCGGCAGCATCTGGTACGAACTGGTGCCATCGTGAGAGTTAAGCAGGATAATTTCCGGAACCTGTTTACCGGTAATTTGCCCCTCCCGACGTAGGCGCAGCATATGCTTAGTATGCTCACGACGACGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/2769036‑2768888
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGTTCAGGGCGGTGTGCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGGG  >  NC_000913/2768933‑2769079
                                                                                                                                                                                                                                                                                                        
GCACACCGCCCTGAACATGCCCGGCAGCATCTGGTACGAACTGGTGCCATCGTGAGAGTTAAGCAGGATAATTTCCGGAACCTGTTTACCGGTAATTTGCCCCTCCCGACGTAGGCGCAGCATATGCTTAGTATGCTCACGACGACGG                                                                                                                                                      <  6:97171/148‑1
                                                                                                     CCTCCCGACGTAGGCGCAGCATATGCTTAGTATGCTCACGACGACGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGTTCAGGGCGGTG                                                >  4:357509/1‑149
                                                                                                       TCCCGACGTAGGCGCAGCATATGCTTAGTATGCTCACGACGACGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCA                                                              <  3:29812/133‑1
                                                                                                       TCCCGACGTAGGCGCAGCATATGCTTAGTATGCTCACGACGACGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCA                                                              >  4:29812/1‑133
                                                                                                                                                  GGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGTTCAGGGCGGTGTGCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGG   <  7:121868/149‑1
                                                                                                                                                   GGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGTTCAGGGCGGTGTGCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGGG  <  3:217333/149‑1
                                                                                                                                                   GGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACAAGTTCGTACCAGATGCTGCCGGGCATGTTCAGGGCGGTGTGACAGAACTGTCTGGTCTGCGGCGAGTCGTTTGGCTAGGTGCGGG  <  8:235340/149‑1
                                                                                                                                                                                                                                                                                                        
GCACACCGCCCTGAACATGCCCGGCAGCATCTGGTACGAACTGGTGCCATCGTGAGAGTTAAGCAGGATAATTTCCGGAACCTGTTTACCGGTAATTTGCCCCTCCCGACGTAGGCGCAGCATATGCTTAGTATGCTCACGACGACGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/2769036‑2768888
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGTTCAGGGCGGTGTGCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGGG  >  NC_000913/2768933‑2769079

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