New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 24614666 (1.850)4 (0.120) 4/224 3.2 6.5% coding (94/279 nt) yafQ mRNA interferase toxin of toxin‑antitoxin pair YafQ/DinJ
?NC_000913 = 246167 56 (1.740)coding (73/279 nt) yafQ mRNA interferase toxin of toxin‑antitoxin pair YafQ/DinJ
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AATCAGGATCCAGTCCGGTTCGACATGAGCATCGCGATAACCTTTCCATGAACCTTGCAGCGGGTGGTCTTTATAAACAGCTGGAAGCGGTAAAGTATTATTGATAAGAAGCGTCATAAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/246027‑246146
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aagcgtcataagGATATGAATAAATTGAAATATCTTATGACGCTTCTTATCAATAATACTTTACCGCTTCCAGCTGTTTATAAAGACCACCCGCTGCAAGGTTCATGGAAAGGTTATCGCGATGCTCATGTCGAACCGGACTGGATCCTGATTTA  <  NC_000913/246167‑246025
                                                                                                                                                                                                                                                                       
AATCAGGATCCAGTCCGGTTCGACATGAGCATCGCGATAACCTTTCCATGAACCTTGCAGCGGGTGGTCTTTATAAACAGCTGGAAGCGGTAAAGTATTATTGATAAGAAGCGTCATAAGGATATGAATAAATTGAAATATCTTATGAC                                                                                                                    >  6:176785/1‑149
      GATCCAGTCCGGTTCGACATGAGCATCGCGATAACCTTTCCATGAACCTTGCAGCGGGTGGTCTTTATAAACAGCTGGAAGCGGTAAAGTATTATTGATAAGAAGCGTCATAAGGATATGAATAAATTGAAATATCTTATGACGCTTCT                                                                                                              <  5:176785/149‑1
                                                  AACCTTGCAGCGGGTGGTCTTTATAAACAGCTGGAAGCGGTAAAGTATTATTGATAAGAAGCGTCATAAGGATATGAATAAATTGAAATATCTTATGACGC                                                                                                                  <  5:2601/101‑1
                                                  AACCTTGCAGCGGGTGGTCTTTATAAACAGCTGGAAGCGGTAAAGTATTATTGATAAGAAGCGTCATAAGGATATGAATAAATTGAAATATCTTATGACGC                                                                                                                  >  6:2601/1‑101
                                                                                                                  CATAAGGATATGAATAAATTGAAATATCTTATGACGCTTCTTATCAATAATACTTTACCGCTTCCAGCTGTTTATAAAGACCACCCGCTGCAAGGTTCATGGAAAGGTTATCGCGATGCTCATGTCGAACCGGACTGGATCCTGATTTA  >  3:125155/1‑149
                                                                                                                   ATAAGGATATGAATAAATTGAAATATCTTATGACGCTTCTTATCAATAATACTTTACCGCTTCCAGCTGTTTATAAAGACCACCCGCTGCAAGGTTCATGGAAAGGTTATCGCGATGCTCATGTCGAACCGGACTGGATCCTGATTTA  >  2:199818/1‑148
                                                                                                                                                                                                                                                                       
AATCAGGATCCAGTCCGGTTCGACATGAGCATCGCGATAACCTTTCCATGAACCTTGCAGCGGGTGGTCTTTATAAACAGCTGGAAGCGGTAAAGTATTATTGATAAGAAGCGTCATAAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/246027‑246146
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aagcgtcataagGATATGAATAAATTGAAATATCTTATGACGCTTCTTATCAATAATACTTTACCGCTTCCAGCTGTTTATAAAGACCACCCGCTGCAAGGTTCATGGAAAGGTTATCGCGATGCTCATGTCGAACCGGACTGGATCCTGATTTA  <  NC_000913/246167‑246025

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.