New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 3001820	27 (0.870)	3 (0.110)	3/152	2.2	9.6%	coding (952/1197 nt)	dpaL	diaminopropionate ammonia‑lyase
	?	NZ_CP009273	= 3001851	32 (1.140)					coding (983/1197 nt)	dpaL	diaminopropionate ammonia‑lyase
Rejected: Frequency below/above cutoff threshold.

GGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/3001726‑3001820 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgttgccgTACGGATTACCCAGCACGCGCATACCTAATGCGGCAACGCTGTCCTGGCAGGAGATGAATTGGGTGGCGCAGTTACGTAGGAT  <  NZ_CP009273/3001851‑3001769                                                                                                                                                                                    GGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGT                                                                                          <  2:447221/90‑1   CAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTG                                                                                        <  2:617019/90‑1                           ACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGTACGGATTACCCAGCACGCGC                                                                >  1:761780/1‑90                               GCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGTACGGATTACCCAGCACGCGCATAC                                                            >  2:176456/1‑90                                                      ACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGC                                                                                       >  1:462870/1‑40                                                      ACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGC                                                                                       <  2:462870/40‑1                                                                                       GCGTTGCCGTACGGATTACCCAGCACGCGCATACCTAATGCGGCAACGCTGTCCTGGCAGGAGATGAATTGGGTGGCGCAGTTACGTAGG    >  2:419196/1‑90                                                                                         GTTGCCGTACGGATTACCCAGCACGCGCATACCTAATGCGGCAACGCTGTCCTGGCAGGAGATGAATTGGGTGGCGCAGTTACGTAGGAT  <  2:49275/90‑1                                                                                                                                                                                    GGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/3001726‑3001820 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgttgccgTACGGATTACCCAGCACGCGCATACCTAATGCGGCAACGCTGTCCTGGCAGGAGATGAATTGGGTGGCGCAGTTACGTAGGAT  <  NZ_CP009273/3001851‑3001769

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