New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	2685639 =	22 (1.010)	5 (0.260)	5/152	1.2	19.7%	coding (3264/3888 nt)	purL	phosphoribosylformylglycinamidine synthase
	?	NZ_CP009273	2685762 =	21 (1.080)					coding (3141/3888 nt)	purL	phosphoribosylformylglycinamidine synthase
Rejected: Frequency below/above cutoff threshold.

TCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGC  >  NZ_CP009273/2685550‑2685716                                                                                          |                                                                              atgtaaaactgtcgttcgatatcaacgaagatgtggcagcaccgtatattgccactggcgcacgtccgaaagttgccgtactgcgtgagc                                                                               >  1:464546‑M1/90‑90 (MQ=255)   gtaaaactgtcgttcgatatcaacgaagatgtggcagcaccgtatattgccactggcgcacgtccgaaagttgccgtactgcgtgagCAg                                                                             <  2:66993‑M1/3‑1 (MQ=255)            gTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGaa                                                                    >  2:129783/1‑90 (MQ=255)            gTCATTGAACAGGATTGACTTCGCCCAACCTGCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGaa                                                                    >  2:302692/1‑90 (MQ=255)             tCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGaaa                                                                   <  1:52934/90‑1 (MQ=255)              cATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAAt                                                                  <  1:431164/90‑1 (MQ=255)                      aGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAgg                                                          >  2:3873/1‑90 (MQ=255)                            gACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCa                                                                              >  2:574979/1‑64 (MQ=255)                            gACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCa                                                                              <  1:574979/64‑1 (MQ=255)                             agatgtggcagcaccgtatattgccactggcgcacgtccgaaagttgccgtactgcgtgagCAGGGCGGGGAAATCTTCCAGGCCCGTGc                                                   <  1:241646‑M1/29‑1 (MQ=255)                                 cGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTcc                                               >  2:421610/1‑90 (MQ=255)                                 cGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTcc                                               >  2:288966/1‑90 (MQ=255)                                  gCCCAACCTTCACCGGCACCCAGCACATCCCCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCg                                              <  1:16544/90‑1 (MQ=255)                                           tCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAgg                                     >  2:352342/1‑90 (MQ=255)                                             aCCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTc                                   <  1:237961/90‑1 (MQ=255)                                             aCCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTc                                   <  1:561115/90‑1 (MQ=255)                                               cGGCACCCAGCACATCACCGTAGGAGAAACCCCCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCAc                                 <  1:352342/90‑1 (MQ=255)                                                 ttgccactggcgcacgtccgaaagttgccgtactgcgtgagCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTc                               <  1:189276‑M1/49‑1 (MQ=255)                                                 gCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGaaa                                                                   >  2:266689/1‑54 (MQ=255)                                                 gCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGaaa                                                                   <  1:266689/54‑1 (MQ=255)                                                 gCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTc                               >  2:555181/1‑90 (MQ=255)                                                 gCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTc                               <  2:66430/90‑1 (MQ=255)                                                      cAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATg                          >  2:447977/1‑90 (MQ=255)                                                      actggcgcacgtccgaaagttgccgtactgcgtgagCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATg                          <  2:420019‑M1/54‑1 (MQ=255)                                                             tCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCg                   >  2:240497/1‑90 (MQ=255)                                                              acgtccgaaagttgccgtactgcgtgagCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGa                  <  1:279814‑M1/62‑1 (MQ=255)                                                                 cGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAg               <  2:590017/90‑1 (MQ=255)                                                                         ttgccgtactgcgtgagCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAg       <  1:73278‑M1/73‑1 (MQ=255)                                                                              accGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGc  <  1:447977/90‑1 (MQ=255)                                                                                          |                                                                              TCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGC  >  NZ_CP009273/2685550‑2685716

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