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.

GTTTCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2331254‑2331370
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cacaaaaCTGTCTTTTGTGGATAACAGCTGGCCCGTTGAACAGCAGAAAAATGTCATGCTCGGTGGTGGCGATAACGCGCTGATGTTGCCCGAGCAGGCGAGCAATATCCGGCTACAAAGTAGTGAAA  <  NC_000913/2331375‑2331255
                                                                                                                                                                                                                                              
GTTTCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAACTGTCTTTTGTGGATAACAGCTGGCCCGTTG                                                                                            >  2:85441/1‑148
               AGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAACTGTCTTTTGTGGATAACAGCTGGCCCGTTGAACAGCAGAAAAATGT                                                                            <  1:85441/149‑1
                                                                CACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAACTGTCTTTTGTGGATAAC                                                                                                         >  7:116965/1‑71
                                                                CACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAACTGTCTTTTGTGGATAAC                                                                                                         <  8:116965/71‑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
                                                                                                                                                                                                                                              
GTTTCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2331254‑2331370
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cacaaaaCTGTCTTTTGTGGATAACAGCTGGCCCGTTGAACAGCAGAAAAATGTCATGCTCGGTGGTGGCGATAACGCGCTGATGTTGCCCGAGCAGGCGAGCAATATCCGGCTACAAAGTAGTGAAA  <  NC_000913/2331375‑2331255

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.