New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	548215 =	23 (1.160)	4 (0.250)	4/132	1.1	18.8%	coding (342/510 nt)	purE	5‑(carboxyamino)imidazole ribonucleotide mutase
	?	NZ_CP009273	548231 =	16 (1.010)					coding (326/510 nt)	purE	5‑(carboxyamino)imidazole ribonucleotide mutase
Rejected: Frequency below/above cutoff threshold.

AAAACGCTGGTGCCGGTGCTGGGCGTGCCAGTACAGAGCGCCGCACTGAGCGGTGTCGATAGCCTCTACTCCATCGTACAAATGCCGCGCGGCATT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/548310‑548215 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aatgccgcgcggcattTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCA  >  NZ_CP009273/548231‑548285                                                                                                                                                         AAAACGCTGGTGCCGGTGCTGGGCGTGCCAGTACAGAGCGCCGCACTGAGCGGTGTCGATAGCCTCTACTCCATCGTACAAATGCCGCGC                                                               <  2:23972/90‑1                          TGCCAGTACAGAGCGCCGCACTGAGCGGTGTCGATAGCCTCTACTCCATCGTACAAATGCCGCGCGGCATTTGTACGATGGAGTAGAGG                                       <  2:376131/89‑1                          TGCCAGTACAGAGCGCCGCACTGAGCGGTGTCGATAGCCTCTACTCCATCGTACAAATGCCGCGCGGCATTTGTACGAT                                                 <  1:307192/79‑1                          TGCCAGTACAGAGCGCCGCACTGAGCGGTGTCGATAGCCTCTACTCCATCGTACAAATGCCGCGCGGCATTTGTACGAT                                                 >  2:307192/1‑79                                    GAGCGCCGCACTGAGCGGTGTCGATAGCCTCTACTCCATCGTACAAATGCCGCGCGGCAT                                                          >  1:580635/1‑60                                    GAGCGCCGCACTGAGCGGTGTCGATAGCCTCTACTCCATCGTACAAATGCCGCGCGGCAT                                                          <  2:580635/60‑1                                                               CCTCTACTCCATCGTACAAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCA  <  1:376131/89‑1                                                                                       GCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCA  <  1:42998/65‑1                                                                                       GCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCA  >  2:42998/1‑65                                                                                                                                                         AAAACGCTGGTGCCGGTGCTGGGCGTGCCAGTACAGAGCGCCGCACTGAGCGGTGTCGATAGCCTCTACTCCATCGTACAAATGCCGCGCGGCATT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/548310‑548215 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aatgccgcgcggcattTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCA  >  NZ_CP009273/548231‑548285

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