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.

ATCGTGGCTGGCAAGCGCAATTTTTCGCGCGCGGCACAGGGCGGCGATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAGGGAGTATTTGCCCTGATAATATTCGCGATACTTCTCGCG  >  NZ_CP009273/4306346‑4306513                                                                                   |                                                                                      aTCGTGGCTGGCAAGCGCAATTTTTCGCGCGCGGCACAGGGCGGCGATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAg                                                                                >  2:533672/1‑90 (MQ=255) aTCGTGGCTGGCAAGCGCAATTTTTCGCGCGCGGCACAGGGCGGCGATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAg                                                                                >  2:571482/1‑90 (MQ=255)       gctggcAAGCGCAATTTTTCGCGCGCGGCACAGGGCGGCGATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGctc                                                                          <  1:475286/90‑1 (MQ=255)       gctggcAAGCGCAATTTTTCGCGCGCGGCACAGGGCGGCGATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGctc                                                                          <  2:554729/90‑1 (MQ=255)       gctggcAAGCGCAATTTTTCGCGCGCGGCACAGGGCGGCGATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGctc                                                                          <  2:124605/90‑1 (MQ=255)       acgctgccgctgtttgaaaaactggtgcagcgcgagccggtgacgctggtgtcgctgatggaccactcgccgggccAGCGCCAGTTGctc                                                                          <  2:430645‑M1/14‑1 (MQ=255)                  cAATTTTTCGCGCGCGGCACAGGGCGGCGATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTAct                                                               <  2:339204/90‑1 (MQ=255)                        ttCGCGCGCGGCACAGGGCGGCGATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCa                                                         <  2:447743/90‑1 (MQ=255)                          aactggtgcagcgcgagccggtgacgctggtgtcgctgatggaccactcgccgggccAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATc                                                       >  2:223138‑M1/58‑90 (MQ=255)                               gtgcagcgcgagccggtgacgctggtgtcgctgatggaccactcgccgggccAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTgcgc                                                  >  1:430645‑M1/53‑90 (MQ=255)                                cgGCACAGGGCGGCGATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCa                                                 <  2:209744/90‑1 (MQ=255)                                cgGCACAGGGCGGCGATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCa                                                 <  1:382731/90‑1 (MQ=255)                                  gCACAGGGCGGCGATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATc                                               <  1:563411/90‑1 (MQ=255)                                       gggCGGCGATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGa                                          <  1:372237/90‑1 (MQ=255)                                           ggcgATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAGGGa                                      <  1:545818/90‑1 (MQ=255)                                           ggcgATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAGGGa                                      <  1:134945/90‑1 (MQ=255)                                            gcgATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAGGGAg                                     >  2:452033/1‑90 (MQ=255)                                                 ggATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAGGGAGTAttt                                >  2:380624/1‑90 (MQ=255)                                                     tCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAGGGAGTATTTGccc                            <  1:533672/90‑1 (MQ=255)                                                          aTTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAGGGAGTATTTGCCCTGata                       <  1:571482/90‑1 (MQ=255)                                                           ttCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAGGGAGTATTTGCCCTGataa                      <  1:121723/90‑1 (MQ=255)                                                           ttCGGCTGCGACCAGCGTCCGGCGAGCGCCAGTTCCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAGGGAGTATTTGCCCTGataa                      <  2:420700/90‑1 (MQ=255)                                                             cGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAGGGAGTATTTGCCCTGataata                    >  2:524540/1‑90 (MQ=255)                                                                    gACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAGGGa                                      <  2:76220/65‑1 (MQ=255)                                                                    gACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAGGGa                                      >  1:76220/1‑65 (MQ=255)                                                                            tGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAGGGAGTATTTGCCCTGATAATATTCGCGATACTtctc     >  2:216674/1‑90 (MQ=255)                                                                               ggcgAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAGGGAGTATTTGCCCTGATAATATTCGCGATACTTCTcgcg  >  2:343033/1‑90 (MQ=255)                                                                                   |                                                                                      ATCGTGGCTGGCAAGCGCAATTTTTCGCGCGCGGCACAGGGCGGCGATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAGGGAGTATTTGCCCTGATAATATTCGCGATACTTCTCGCG  >  NZ_CP009273/4306346‑4306513

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