New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	3516287 =	18 (0.870)	3 (0.160)	3/154	1.6	17.0%	coding (58/2553 nt)	mrcA	peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA
	?	NZ_CP009273	3516322 =	13 (0.700)					coding (93/2553 nt)	mrcA	peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA
Rejected: Frequency below/above cutoff threshold.

CTGCATCGGAATTTGCAGGCGAACATCTTTTAATGTCGCCACATCCGGCAGTTGTGGCTCGATGTAGCGGTATAGGCCATAAATCGAGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3516376‑3516287 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atcgagccACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAGATTTACAGCGCCGATGGCGAGCTGATTGCTCAA  >  NZ_CP009273/3516322‑3516409                                                                                                                                                                                    CTGCATCGGAATTTGCAGGCGAACATCTTTTAATGTCGCCACATCCGGCAGTTGTGGCTCGATGTAGCGGTATAGGCCATAAATCGAGCC                                                                                          <  2:124702/90‑1                                  TGTCCCCACACCCGGAAGTGGGGGTCCGGTGTAGCGGTATAGGCCATAAATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGT                                                         <  1:202317/90‑1                                                                        ATAGGCCATAAATCGAGCCACAACTGCCGGATGTGGCGA                                                                      <  1:44948/39‑1                                                                        ATAGGCCATAAATCGAGCCACAACTGCCGGATGTGGCGA                                                                      >  2:44948/1‑39                                                                                      GAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAGATTTACAGCGCCGATGGCGAGCTGAT         >  1:99299/1‑86                                                                                      GAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAGATTTACAGCGCCGATGGCGAGCTGAT         <  2:99299/86‑1                                                                                         CCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAGATTTACAGCGCCGATGGCGAGCTGATTGCTCAA  <  2:286705/90‑1                                                                                                                                                                                    CTGCATCGGAATTTGCAGGCGAACATCTTTTAATGTCGCCACATCCGGCAGTTGTGGCTCGATGTAGCGGTATAGGCCATAAATCGAGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3516376‑3516287 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atcgagccACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAGATTTACAGCGCCGATGGCGAGCTGATTGCTCAA  >  NZ_CP009273/3516322‑3516409

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