New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3766956 =97 (0.940)4 (0.040) 4/248 NT 4.5% coding (620/843 nt) yibA putative immunity protein for polymorphic toxin RhsA; HEAT‑domain protein; lethality reduction protein
?NC_000913 3766985 = 77 (0.790)coding (649/843 nt) yibA putative immunity protein for polymorphic toxin RhsA; HEAT‑domain protein; lethality reduction protein
Rejected: Frequency below/above cutoff threshold.

CTAAAACAGGAAGTAGCGTTTTATCACCTAATTCACCCGCAGCTTCAATGATATCATCATAAACAGTATTTTTTTTTAACTCATCGCATAAAACAGATAAAACCCTTTTATCTTTTCTGTAGGAAAGTCCGATTATTGCTTCAATA‑‑CGGACTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3767108‑3766956
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggacttTCCTACAGAAAAGATAAAAGGGTTTTATCTGTTTTATGCGATGAGTTAAAAAAAAATACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGAT  >  NC_000913/3766985‑3767131
                                                                                                                                                                                                                                                                                                              
CTAAAACAGGAAGTAGCGTTTTATCACCTAATTCACCCGCAGCTTCAATGATATCATCATAAACAGTATTTTTTTTTAACTCATCGCATAAAACAGATAAAACCCTTTTATCTTTTCTGTAGGAAAGTCCGATTATTGCTTCAATA‑‑CGG                                                                                                                                                         >  4:143141/1‑149
CTAAAACAGGAAGTAGCGTTTTATCACCTAATTCACCCGCAGCTTCAATGATATCATCATAAACAGTATTTTTTTTTAACTCATCGCATAAAACAGATAAAACCCTTTTATCTTTTCTGTAGGAAAGTCCGATTATTGCTTCAATA‑‑CGG                                                                                                                                                         >  4:795132/1‑149
 TAAAACAGGAAGTAGCGTTTTATCACCTAATTCACCCGCAGCTTCAATGATATCATCATAAACAGTATTTTTTTTTAACTCATCGCATAAAACAGATAAAACCCTTTTATCTTTTCTGTAGGAAAGTCCGATTATTGCTTCAATA‑‑CGGA                                                                                                                                                        >  3:170629/1‑149
              AGCGTTTTATCACCTAATTCACCCGCAGCTTCAATGATATCATCATAAACAGTATTTTTTTTTAACTCATCGCATAAAACAGATAAAACCCTTTTATCTTTTCTGTAGGAAAGTCCG‑‑TATTGAAGCAATAATCGGACTTTCCTACAGAA                                                                                                                                           <  2:450600/149‑1
                                                                                                                            AAGTCCGATTATTGCTTCAATA‑‑CGGACTTTCCTACAGAAAAGATAAAAGGGTTTTATCTGTTTTATGCGATGAGTTAAAAAAAAATACTGTTTATGATGATATCAT                                                                        >  1:835639/1‑106
                                                                                                                            AAGTCCGATTATTGCTTCAATA‑‑CGGACTTTCCTACAGAAAAGATAAAAGGGTTTTATCTGTTTTATGCGATGAGTTAAAAAAAAATACTGTTTATGATGATATCAT                                                                        <  2:835639/106‑1
                                                                                                                             AGTCCGATTATTGCTTCAATA‑‑CGGACTTTCCTACAGAAAAGATAAAAGGGTTTTATCTGTTTTATGCGATGAGTTAAAAAAAAATACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTT                            <  1:450600/149‑1
                                                                                                                                                     GGACTTTCCTACAGAAAAGATAAAAGGGTTTTATCTGTTTTATGCGATGAGTTAAAAAAAAATACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTT      <  3:675430/149‑1
                                                                                                                                                     GGACTTTCCTACAGAAAAGATAAAAGGGTTTTATCTGTTTTATGCGATGAGTTAAAAAAAAATACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTT      >  4:228840/1‑149
                                                                                                                                                      GACTTTCCTACAGAAAAGATAAAAGGGTTTTATCTGTTTTATGCGATGAGTTAAAAAAAAATACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTT     >  2:134970/1‑149
                                                                                                                                                       ACTTTCTTACAGAAAAGATAAAAGGGTTTTATCTGTTTTATGCGATGAGTTAAAAAAAAATACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTG    <  4:305570/149‑1
                                                                                                                                                        CTTTCCTACAGAAAAGATAAAAGGGTTTTATCTGTTTTATGCGATGAGTTAAAAAAAAATACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGA   <  2:107420/149‑1
                                                                                                                                                        CTTTCCTACAGAAAAGATAAAAGGGTTTTATCTGTTTTATGCGATGAGTTAAAAAAAAATACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTCTACAAGTTTGA   <  1:489983/149‑1
                                                                                                                                                         TTTCCTACAGAAAAGATAAAAGGGTTTTATCTGTTTTATGCGATGAGTTAAAAAAAAATACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGAT  <  3:463143/149‑1
                                                                                                                                                         TTTCCTACAGAAAAGATAAAAGGGTTTTATCTGTTTTATGCGATGAGTTAAAAAAAAATACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGAT  >  3:817217/1‑149
                                                                                                                                                                                                                                                                                                              
CTAAAACAGGAAGTAGCGTTTTATCACCTAATTCACCCGCAGCTTCAATGATATCATCATAAACAGTATTTTTTTTTAACTCATCGCATAAAACAGATAAAACCCTTTTATCTTTTCTGTAGGAAAGTCCGATTATTGCTTCAATA‑‑CGGACTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3767108‑3766956
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggacttTCCTACAGAAAAGATAAAAGGGTTTTATCTGTTTTATGCGATGAGTTAAAAAAAAATACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGAT  >  NC_000913/3766985‑3767131

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