New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	4306428 =	24 (1.100)	3 (0.160)	3/150	1.6	13.0%	coding (608/1137 nt)	phnM	alpha‑D‑ribose 1‑methylphosphonate 5‑triphosphate diphosphatase
	?	NZ_CP009273	4306531 =	19 (0.990)					coding (505/1137 nt)	phnM	alpha‑D‑ribose 1‑methylphosphonate 5‑triphosphate diphosphatase
Rejected: Frequency below/above cutoff threshold.

TTCGTACTGCTGCATCTGCGCATCAGTGAGGGAGTATTTGCCCTGATAATATTCGCGATACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTG  >  NZ_CP009273/4306445‑4306607                                                                                       |                                                                             ttcGTACTGCTGCATCTGCGCATCAGTGAGGGAGTATTTGCCCTGATAATATTCGCGATACTTCTCGCGGTTGGCGAACTGGCGCTGGcc                                                                           <  1:216674/90‑1 (MQ=255) ttcGTACTGCTGCATCTGCGCATCAGTGAGGGAGTATTTGCCCTGATAATATTCGCGATACTTCTCGCGGTTGGCGAACTGGCGCTGGcc                                                                           <  1:58868/90‑1 (MQ=255) ttcGTACTGCTGCATCTGCGAATCAGTGAGGGAGTATTTGCCCTGATAATATTCGCGATACTTCTCGCGGTTGGCGAACTGGCGCTGGcc                                                                           <  2:151132/90‑1 (MQ=255)       ctgctgCATCTGCGCATCAGTGAGGGAGTATTTGCCCTGATAATATTCGCGATACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGa                                                                     >  2:548847/1‑90 (MQ=255)         gctgcATCTGCGCATCAGTGAGGGAGTATTTGCCCTGATAATATTCGCGATACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAgg                                                                   >  1:601818/1‑89 (MQ=255)         gctgcATCTGCGCATCAGTGAGGGAGTATTTGCCCTGATAATATTCGCGATACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGt                                                                   >  1:142142/1‑90 (MQ=255)         gctgcATCTGCGCATCAGTGAGGGAGTATTTGCCCTGATAATATTCGCGATACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGt                                                                   <  1:548847/90‑1 (MQ=255)           tgcATCTGCGCATCAGTGAGGGAGTATTTGCCCTGATAATATTCGCGATACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGTgg                                                                 <  2:600103/90‑1 (MQ=255)                        cAGTGAGGGAGTATTTGCCCTGATAATATTCGCGATACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGTGGTCCATCAGCGaca                                                    >  1:367645/1‑90 (MQ=255)                                  gTATTTGCCCTGATAATATTCGCGATACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCAc                                          <  1:380624/90‑1 (MQ=255)                                      ttGCCCTGATAATATTCGCGATACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGc                                      >  1:120411/1‑90 (MQ=255)                                            tGATAATATTCGCGATACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGTGGTCAATCAGCGACACCAGCGTCACCGGCTCgcgc                                >  2:489207/1‑90 (MQ=255)                                             gATAATATTCGCGATACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCt                               >  2:254179/1‑90 (MQ=255)                                                 gcgcagatgcagcagtacgaagaagagcaactggcgctGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCAc                           >  1:430645‑M2/39‑90 (MQ=255)                                                 atatTCGCGATACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCAc                           <  2:121506/90‑1 (MQ=255)                                                      gatgcagcagtacgaagaagagcaactggcgctGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGtt                      >  2:223138‑M2/34‑90 (MQ=255)                                                      cgcgATACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAggt                      >  2:404881/1‑88 (MQ=255)                                                       gcgATACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCCCCGCACCCGttt                     >  1:542652/1‑90 (MQ=255)                                                       gcgATACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCCCCGCACCAGgtt                     >  1:408531/1‑90 (MQ=255)                                                            aCTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCaa                <  1:196130/90‑1 (MQ=255)                                                                         ggCGAACTGGCGCTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAG‑TTTTCAAACAGCGGCAggggg  >  2:405867/1‑86 (MQ=255)                                                                         gagcaactggcgctGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCgt   <  2:430645‑M2/76‑1 (MQ=255)                                                                                       |                                                                             TTCGTACTGCTGCATCTGCGCATCAGTGAGGGAGTATTTGCCCTGATAATATTCGCGATACTTCTCGCGGTTGGCGAACTGGCGCTGGCCCGGCGAGTGGTCCATCAGCGACACCAGCGTCACCGGCTCGCGCTGCACCAGTTTTTCAAACAGCGGCAGCGTG  >  NZ_CP009273/4306445‑4306607

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