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.

TGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCAGGCTCGATTTATGGCCTATACCGCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAG  >  NZ_CP009273/3516199‑3516376                                                                                         |                                                                                          ttgagcaatcagctcgccatcggcgctgtaaatctgcatcggaatttgcaggcgaacatcttttaatgtcgccacatccggcagttgtgg                                                                                          <  2:286705‑M1/2‑1 (MQ=255) tGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCAgg                                                                                          >  2:494420/1‑90 (MQ=255) tGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCAgg                                                                                          >  1:57876/1‑90 (MQ=255)  gCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCAGGc                                                                                         >  2:270142/1‑90 (MQ=255)  gCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCAGGc                                                                                         >  2:267040/1‑90 (MQ=255)        atcagctcgccatcggcgctgtaaatctgcatcggaatttgcaggcgaacatcttttaatgtcgccacatccggcagttgtGGCTc                                                                                       <  2:99299‑M1/5‑1 (MQ=255)        atcagctcgccatcggcgctgtaaatctgcatcggaatttgcaggcgaacatcttttaatgtcgccacatccggcagttgtGGCTc                                                                                       >  1:99299‑M1/82‑86 (MQ=255)          gAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCAGGCTCGATTTa                                                                                 <  1:509184/90‑1 (MQ=255)            aaCTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCAGGCTCGATTTATg                                                                               <  2:401328/90‑1 (MQ=255)                  aTGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCAGGCTCGATTTATGGCCtat                                                                         <  1:494420/90‑1 (MQ=255)                             ccAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCAGGCTCGATTTATGGCCTATACCGCTACATc                                                              <  2:44312/90‑1 (MQ=255)                             ccAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCAGGCTCGATTTATGGCCTATACCGCTACATc                                                              <  2:275686/90‑1 (MQ=255)                                      ttCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCAGGCTCGATTTATGGCCTATACCGCTACATCGAGCCACaa                                                     <  1:312921/90‑1 (MQ=255)                                                        acatcttttaatgtcgccacatccggcagttgtGGCTCGATTTATGGCCTATACCGCTACACCGGACCCCCACTTCCGGGTGTGGGGACa                                   <  1:202317‑M1/57‑1 (MQ=255)                                                          ccTTGCAGTCTGTTGCATTCTGCTGGGAGCAGGCTCGATTTATGGCCTATACCGCTACATCGAGCCACAACTGCCGGATGTGGCGACAtt                                 >  1:399526/1‑90 (MQ=255)                                                             tGCAGTCTGTTGCATTCTGCTGGGAGCAGGCTCGATTTATGGcc                                                                            >  2:241211/1‑44 (MQ=255)                                                             tGCAGTCTGTTGCATTCTGCTGGGAGCAGGCTCGATTTATGGcc                                                                            <  1:241211/44‑1 (MQ=255)                                                                   cTGTTGCATTCTGCTGGGAGCAGGCTCGATTTATGGCCTATACCGCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGt                        >  2:261601/1‑90 (MQ=255)                                                                     tcgccacatccggcagttgtGGCTCGATTTATGGCCtat                                                                         <  1:44948‑M1/19‑1 (MQ=255)                                                                     tcgccacatccggcagttgtGGCTCGATTTATGGCCtat                                                                         >  2:44948‑M1/21‑39 (MQ=255)                                                                        gCATTCTGCTGGGAGCAGGCTCGATTTATGGCCTATACCGCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGcc                   <  1:395112/90‑1 (MQ=255)                                                                          aTTCTGCTGGGAGCAGGCTCGATTTATGGCCTATACCGCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTg                 >  1:193621/1‑90 (MQ=255)                                                                             ctgctgGGAGCAGGCTCGATTTATGGCCTATACCGCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCaa              >  2:341763/1‑90 (MQ=255)                                                                             ctgctgGGAGCAGGCTCGATTTATGGCCTATACCGCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCaa              >  2:561186/1‑90 (MQ=255)                                                                                         GGCTCGATTTATGGCCTATACCGCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAg  <  2:124702‑M1/90‑1 (MQ=255)                                                                                         |                                                                                          TGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCAGGCTCGATTTATGGCCTATACCGCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAG  >  NZ_CP009273/3516199‑3516376

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