New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4415329121 (0.990)5 (0.040) 5/244 NT 4.0% coding (1055/1626 nt) aidB DNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase
?NC_000913 = 4415338 121 (1.020)coding (1064/1626 nt) aidB DNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase
Rejected: Frequency below/above cutoff threshold.

GCCCTGGGTAGCCATGCCATGATGCGCCGTGCATTTTCGCTGGCGATTTATCATGCACATCAACGCCATGTTTTTGGTAATCCATTGATCCAACAGCCCCTTATGCGTCATGTCTTAAGTCGCATGGCACTTCAGCTTGAAGGGCAAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4415181‑4415329
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaacAGCAACGCCGTTTGCCCTTCAAGCTGAAGTGCCATGCGACTTAAGACATGACGCATAAGGGGCTGTTGGATCAATGGATTACCAAAAACATGGCGTTGATGTGCATGATAAATCGCCAGCGAAAATGCACGGCGCATCATGGCATGGC  <  NC_000913/4415338‑4415191
                                                                                                                                                                                                                                                                                                         
GCCCTGGGTAGCCATGCCATGATGCGCCGTGCATTTTCGCTGGCGATTTATCATGCACATCAACGCCATGTTTTTGGTAATCCATTGATCCAACAGCCCCTTATGCGTCATGTCTTAAGTCGCATGGCACTTCAGCTTGAAGGGCAAAC                                                                                                                                                      <  3:558029/149‑1
              TGCCATGATGCGCCGTGCATTTTCGCTGGCGATTTATCATGCACATCAACGCCATGTTTTTGGTAATCCATTGATCCAACAGCCCCTTATGCGTCATGTCTTAAGTCGCATGGCACTTCAGCTTGAAGGGCAAACAGCAACGCCGTTTG                                                                                                                                        <  4:1045851/149‑1
                                                                                                           TCATGTCTTAAGTCGCATGGCACTTCAGCTTGAAGGGCAAACAGCAACGCCGTTTGCC                                                                                                                                      >  3:175821/1‑58
                                                                                                           TCATGTCTTAAGTCGCATGGCACTTCAGCTTGAAGGGCAAACAGCAACGCCGTTTGCC                                                                                                                                      <  4:175821/58‑1
                                                                                                                                        TTGAAGGGCAAACGGCGTTGCTGTTTGCCCTTCAAGCTGAAGTGCCATGCGACTTAAGACATGACGCATAAGGGGCTGTTGGATCAAT                                                                           <  1:864699/88‑1
                                                                                                                                        TTGAAGGGCAAACGGCGTTGCTGTTTGCCCTTCAAGCTGAAGTGCCATGCGACTTAAGACATGACGCATAAGGGGCTGTTGGATCAAT                                                                           >  2:864699/1‑88
                                                                                                                                                  AACAGCAACGCCGTTTGCCCTTCAAGCTGAAGTGCCATGCGACTTAAGACATGACGCATAAGGGGCTGTTGGATCAATGGATTACCAAAAACATGGCGTTGATGTGCATGATAAATCGCCAGCGAAAATGCACGGCGCATCATGGCATG    <  3:430999/149‑1
                                                                                                                                                   ACAGCAACGCCGTTTGCCCTTCAAGCTGAAGTGCCATGCGACTTAAGACATGACGCATAAGGGGCTGTTGGATCAATGGATTACCAAAAACATGGCGTTGATGTGCATGATAAATCGCCAGCGAAAATGCACGGCGCATCATGGCATGG   >  4:581598/1‑149
                                                                                                                                                    CAGCAACGCCGTTTGCCCTTCAAGCTGAAGTGCCATGCGACTTAAGACATGACGCATAAGGGGCTGTTGGATCAATGGATTACCAAAAACATGGCGTTGATGTGCATGATAAATCGCCAGCGAAAATGCACGGCGCATCATGGCATGGC  >  3:339251/1‑149
                                                                                                                                                                                                                                                                                                         
GCCCTGGGTAGCCATGCCATGATGCGCCGTGCATTTTCGCTGGCGATTTATCATGCACATCAACGCCATGTTTTTGGTAATCCATTGATCCAACAGCCCCTTATGCGTCATGTCTTAAGTCGCATGGCACTTCAGCTTGAAGGGCAAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4415181‑4415329
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaacAGCAACGCCGTTTGCCCTTCAAGCTGAAGTGCCATGCGACTTAAGACATGACGCATAAGGGGCTGTTGGATCAATGGATTACCAAAAACATGGCGTTGATGTGCATGATAAATCGCCAGCGAAAATGCACGGCGCATCATGGCATGGC  <  NC_000913/4415338‑4415191

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