New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1759424 =30 (0.940)6 (0.200) 5/242 2.9 17.6% coding (1100/1221 nt) sufS cysteine desulfurase, stimulated by SufE; selenocysteine lyase, PLP‑dependent
?NC_000913 1759462 = 28 (0.930)coding (1062/1221 nt) sufS cysteine desulfurase, stimulated by SufE; selenocysteine lyase, PLP‑dependent
Rejected: Frequency below/above cutoff threshold.

TGGAATCTGTACCGGATCTCACTCTCTATGGCCCACAAAACAGGCTTGGCGTTATTGCTTTTAATCTCGGTAAACACCACGCCTATGATGTTGGCAGTTTTCTCGATAATTACGGCATTGCTGTGCGTACCGGACATCACTGCGCAATGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1759574‑1759424
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcaatgccGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGT  >  NC_000913/1759462‑1759606
                                                                                                                                                                                                                                                                                                        
TGGAATCTGTACCGGATCTCACTCTCTATGGCCCACAAAACAGGCTTGGCGTTATTGCTTTTAATCTCGGTAAACACCACGCCTATGATGTTGGCAGTTTTCTCGATAATTACGGCATTGCTGTGCGTACCGGACATCACTGCGCAATG                                                                                                                                                     >  4:84256/1‑149
  GAATCTGTACCGGATCTCACTCTCTATGGCCCACAAAACAGGCTTGGCGTTATTGCTTTTAATCTCGGTAAACACCACGCCTATGATGTTGGCAGTTTTCTCGATAATTACGGCATTGCTGTGCGTACCGGACATCACTGCGCAATGC                                                                                                                                                    <  6:96352/148‑1
                                                                                                        GATAATTACGGCATTGCTGTGCGTACCGGACATCACTGCGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGG                                             >  4:73841/1‑149
                                                                                                                CGGCATTGCTGTGAGTACCGGACATCACTGCGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTCCCGAGATTAAAAGCAATAACGCAAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATT                                     >  6:106470/1‑149
                                                                                                                 GGCATTGCTGTGCGTACCGGACATCACTGCGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGA                                       <  1:133888/146‑1
                                                                                                                 GGCATTGCTGTGCGTACCGGACATCACTGCGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGA                                       >  2:133888/1‑146
                                                                                                                                          ACTGCGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAATCAATAACTCCAAGCCTGTTTTGTGGGACATAGAGAGTGATATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATC           <  5:106470/149‑1
                                                                                                                                          ACTGCGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATC           <  3:73841/149‑1
                                                                                                                                                  ATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTG   >  1:22401/1‑149
                                                                                                                                                   TGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGT  <  4:78787/149‑1
                                                                                                                                                   TGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGT  <  6:100202/149‑1
                                                                                                                                                                                                                                                                                                        
TGGAATCTGTACCGGATCTCACTCTCTATGGCCCACAAAACAGGCTTGGCGTTATTGCTTTTAATCTCGGTAAACACCACGCCTATGATGTTGGCAGTTTTCTCGATAATTACGGCATTGCTGTGCGTACCGGACATCACTGCGCAATGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1759574‑1759424
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcaatgccGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGT  >  NC_000913/1759462‑1759606

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.