New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AM260479 = 146542618 (0.680)4 (0.170) 4/194 1.9 17.3% coding (4779/4857 nt) gudB NAD‑specific glutamate dehydrogenase
?AM260479 = 1465451 22 (0.910)coding (4754/4857 nt) gudB NAD‑specific glutamate dehydrogenase
Rejected: Frequency below/above cutoff threshold.

GACGATGCGGACCTCGACACGGACGCGGAAGCAGACGAATAAGCCTGGGCTGGCCAGCGCTATGCGCGTTCCAGCATGCCGATCTCGCGCACCGCGACCGACAGCATCGACAGCCCCGCAGCCCCCGAGGCGCGCTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AM260479/1465289‑1465426
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgcgctggAACGCTTTACGCGCATGCTGGCCGACCAGCGCG  <  AM260479/1465451‑1465419
                                                                                                                                                                           
GACGATGCGGACCTCGACACGGACGCGGAAGCAGACGAATAAGCCTGGGCTGGCCAGCGCTATGCGCGTTCCAGCATGCCGATCTCGCGCACCGCGACCGACAGCATCGACAGCCCCGCAGCCCCCGAGGCGCG                                       <  2:971377/134‑1
                                                GCTGGCCAGCGCTATGCGCGTTCCAGCATGCCGATCTCGCGCACCGCGACCGACAGCATCGACAGCCCCGCAGCCCCCGAGGCGCGCTGGAACGCTTTACGCGCATGCTGGCCGACCAG      <  1:914042/119‑1
                                                GCTGGCCAGCGCTATGCGCGTTCCAGCATGCCGATCTCGCGCACCGCGACCGACAGCATCGACAGCCCCGCAGCCCCCGAGGCGCGCTGGAACGCTTTACGCGCATGCTGGCCGACCAG      >  2:914042/1‑119
                                                                                                    ACAGCATCGACAGCCCCGCAGCCCCCGAGGCGCGCTGGAACGCTTTACGCGCATGCTGGCCGACCAGCGCG  >  1:446566/1‑71
                                                                                                    ACAGCATCGACAGCCCCGCAGCCCCCGAGGCGCGCTGGAACGCTTTACGCGCATGCTGGCCGACCAGCGCG  <  2:446566/71‑1
                                                                                                                                                                           
GACGATGCGGACCTCGACACGGACGCGGAAGCAGACGAATAAGCCTGGGCTGGCCAGCGCTATGCGCGTTCCAGCATGCCGATCTCGCGCACCGCGACCGACAGCATCGACAGCCCCGCAGCCCCCGAGGCGCGCTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AM260479/1465289‑1465426
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgcgctggAACGCTTTACGCGCATGCTGGCCGACCAGCGCG  <  AM260479/1465451‑1465419

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