New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 233362915 (0.490)4 (0.150) 4/112 2.1 18.7% pseudogene (774/4605 nt) yfaS pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
?NC_000913 = 2334464 22 (0.850)coding (496/624 nt) yfaT DUF1175 family protein, putative host defense protein
Rejected: Frequency below/above cutoff threshold.

GCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2333556‑2333629
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aagcgccacgAAAACTGACAACGGAATGCGCGCAGTCAGTCTGCAACAACTTATGACATGGAAGGACACCCGA  <  NC_000913/2334464‑2334402
                                                                                                                                         
GCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGC                                                                       <  2:356760/68‑1
     GCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCAC                                                                  >  1:201495/1‑68
      CTTCACCCTTTTTTTTACCCGCTGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCTCG                                                                 <  2:552555/68‑1
                                  CACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGAAAACTGCCAACGGCATGCGCGCAGTCA                                     >  1:549548/1‑68
                                    CCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGAAAACTGACAACGGAATGCGCGCA                                         >  1:637356/1‑62
                                    CCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGAAAACTGACAACGGAATGCGCGCA                                         <  2:637355/62‑1
                                         AGCTCTTTGCCTGACACTTTGCTAAGCGCCACGAAAACTGACAACGGAATGCGCGCAGTCAGTCTGCg                              >  2:392541/1‑67
                                                                AAGCGCCACGAAAACTGACAACGGAATGCGCGCAGTCAGTCTGCAACAACTTATGACATGGAAGGACA       <  1:658663/68‑1
                                                                 AGCGCCACGAAAACTGACAACGGAATGCGCGCAGTCAGTCTGCAACAACTTATGACATGGAAGGACAC      >  1:439592/1‑68
                                                                     CCACGAAAACTGACAACGGAATGCGCGCAGTCAGTCTGCAACAACTTATGACATGGAAGGACACCCGA  <  1:1152030/68‑1
                                                                     CCACGAAAACTGACAACGGAATGCGCGCAGTCAGTCTGCAACAACTTATGACATGGAAGGACACCCGA  <  1:1153544/68‑1
                                                                                                                                         
GCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2333556‑2333629
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aagcgccacgAAAACTGACAACGGAATGCGCGCAGTCAGTCTGCAACAACTTATGACATGGAAGGACACCCGA  <  NC_000913/2334464‑2334402

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