New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 4195934	21 (0.660)	3 (0.110)	3/146	2.4	12.2%	coding (1527/1590 nt)	purH	bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase
	?	NZ_CP009273	4196004 =	25 (0.900)					coding (1457/1590 nt)	purH	bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase
Rejected: Frequency below/above cutoff threshold.

TCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/4195844‑4195934 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcggcggcaTCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGCCGCTTTAATACCGG  >  NZ_CP009273/4196004‑4196091                                                                                                                                                                                     TCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGC                                                                                           >  1:478461/1‑90                                                          GAGCATCGCAATACCGTGCTCGGCGGCGGCGGCATCAATACCGTCGCGGAAAGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCAC                                  >  2:487673/1‑90                                                            GCATCGCAATACCGTGCTCGTCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAAC                                                               >  1:307790/1‑59                                                            GCATCGCAATACCGTGCTCGTCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAAC                                                               <  2:307790/59‑1                                                                                  GGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCAC                                  <  1:423142/66‑1                                                                                  GGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCAC                                  >  2:423142/1‑66                                                                                          CATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGCCGCTTTAATACCGG  >  1:187010/1‑90                                                                                          CATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGCCGCTTTAATACCGG  >  2:381586/1‑90                                                                                                                                                                                     TCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/4195844‑4195934 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcggcggcaTCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGCCGCTTTAATACCGG  >  NZ_CP009273/4196004‑4196091

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