New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 233137033 (1.120)4 (0.140) 4/246 2.9 10.9% pseudogene (3033/4605 nt) yfaS pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
?NC_000913 = 2331375 34 (1.220)pseudogene (3028/4605 nt) yfaS pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Rejected: Frequency below/above cutoff threshold.

CACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2331318‑2331370
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cacaaaaCTGTCTTTTGTGGATAACAGCTGGCCCGTTGAACAGCAGAAAAATGTCATGCTCGGTGGTGGCGATAACGCGCTGATGTTGCCCGAGCAGGCGAGCAATATCCGGCTACAAAGTAGTGAAAC  <  NC_000913/2331375‑2331254
                                                                                                                                                                               
CACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAACTGTCTTTTGTGGATAAC                                                                                                          >  7:116965/1‑71
CACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAACTGTCTTTTGTGGATAAC                                                                                                          <  8:116965/71‑1
           ACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGGATAACAGCTGGCCCGTTGAACAGCAGAAAAATGTCATGCTCGGTGGTGGCGATAACGCGCTGATGTTGCCCGAGCAGGCGAGCAATATCCGGCT                 >  1:85441/1‑149
                           CAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGGATAACAGCTGGCCCGTTGAACAGCAGAAAAATGTCATGCTCGGTGGTGGCGATAACGCGCTGATGTTGCCCGAGCAGGCGAGCAATATCCGGCTACAAAGTAGTGAAAC  <  2:85441/148‑1
                                              CACAAAACTGTCTTTTGTGGATAACAGCTGGCCCGTTGAACAGCAGAAAAATGTCATGCTCGGTGGTGGCGATAACGCGCTGATGTTGCCCGAGCAGGCGAGCAATATC                      >  3:28313/1‑109
                                              CACAAAACTGTCTTTTGTGGATAACAGCTGGCCCGTTGAACAGCAGAAAAATGTCATGCTCGGTGGTGGCGATAACGCGCTGATGTTGCCCGAGCAGGCGAGCAATATC                      <  4:28313/109‑1
                                               ACAAAACTGTCTTTTGTGGATAACAGCTGGCCCGTTGAACAGCAGAAAAATGTCATGCTCGGTGGTGGCGATAACGCGCTGA                                                <  3:12188/82‑1
                                               ACAAAACTGTCTTTTGTGGATAACAGCTGGCCCGTTGAACAGCAGAAAAATGTCATGCTCGGTGGTGGCGATAACGCGCTGA                                                >  4:12188/1‑82
                                                CAAAACTGTCTTTTGTGGATAACAGCTGGCCCGTTGAACAGCAGAAAAATGTCATGCTCGGTGGTGGCGATAACGCGCTGATGTTGCCCGAGCAGGCGAGCAATATCCGGCTACAAAGTAGTGAAA   >  3:124095/1‑126
                                                CAAAACTGTCTTTTGTGGATAACAGCTGGCCCGTTGAACAGCAGAAAAATGTCATGCTCGGTGGTGGCGATAACGCGCTGATGTTGCCCGAGCAGGCGAGCAATATCCGGCTACAAAGTAGTGAAA   <  4:124095/126‑1
                                                                                                                                                                               
CACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2331318‑2331370
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cacaaaaCTGTCTTTTGTGGATAACAGCTGGCCCGTTGAACAGCAGAAAAATGTCATGCTCGGTGGTGGCGATAACGCGCTGATGTTGCCCGAGCAGGCGAGCAATATCCGGCTACAAAGTAGTGAAAC  <  NC_000913/2331375‑2331254

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