New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3000957 =51 (0.930)5 (0.110) 5/134 2.9 11.7% coding (89/1197 nt) dpaL diaminopropionate ammonia‑lyase
?NZ_CP009273 3000973 = 34 (0.760)coding (105/1197 nt) dpaL diaminopropionate ammonia‑lyase
Rejected: Frequency below/above cutoff threshold.

TTCTTCACACCAAAAAGGTTTGCGAGATCGTCCAGCGCACAAAGCGGTGTTGGGCGATAACCAGCTATTTTCTGGTGGAACTGGCGCGCCAGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3001050‑3000957
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aactggcgcgccagttCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCAAACCTTTTTGGTGTGAAGAAAATTCTC  >  NZ_CP009273/3000973‑3001057
                                                                                                                                                                                   
TTCTTCACACCAAAAAGGTTTGCGAGATCGTCCAGCGCACAAAGCGGTGTTGGGCGATAACCAGCTATTTTCTGGTGGAACTGGCGCGCC                                                                                           >  1:1617244/1‑90
TTCTTCACACCAAAAAGGTTTGCGAGATCGTCCAGCGCACAAAGCGGTGTTGGGCGATAACCAGCTATTTTCTGGTGGAACTGGCGCGCC                                                                                           >  2:1426152/1‑90
    TCACACCAAAAAGGTTTGCGAGATCGTCCAGCGCACAAAGCGGTGTTGGGCGATAACCAGCTATTTTCTGGTGGAACTGGCGCGCCAGTT                                                                                       <  1:869691/90‑1
                               CCAGCGCACAAAGCGGTGTTGGGCGATAACCAGCTATTTTCTGGTGGAACTGGCGCGCCAGTTCCCCCAGAAAATAGCTGGTTATCGCCC                                                            <  2:326495/90‑1
                                         aaaCGGTGTTGGGCGATAACCATCTTTTTTCTGGTGGAACTGGCGCTCCCGTCACACCAGAAAATAGCTGGTTATCGCCCAACACCGCTT                                                  <  2:45041/87‑1
                                                        ATAACCAGCTATTTTCTGGTGTAACGGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGTGCTGGACGAT                                   >  2:62123/1‑90
                                                                    TTTCTGGTGGAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCAA                            >  1:1169045/1‑85
                                                                    TTTCTGGTGGAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCAA                            <  2:1169045/85‑1
                                                                                   GCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCAAACCTTTTTGGTGTGAAGAAA        >  1:19397/1‑90
                                                                                        CCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCAAACCTTTTTGGTGTGAAGAAAATTCT   <  1:1175550/90‑1
                                                                                         CAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCAAACCTTTTTGGTGTGAAGAAAATTCTC  <  2:199427/90‑1
                                                                                                                                                                                   
TTCTTCACACCAAAAAGGTTTGCGAGATCGTCCAGCGCACAAAGCGGTGTTGGGCGATAACCAGCTATTTTCTGGTGGAACTGGCGCGCCAGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3001050‑3000957
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aactggcgcgccagttCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCAAACCTTTTTGGTGTGAAGAAAATTCTC  >  NZ_CP009273/3000973‑3001057

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