New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	548215 =	23 (1.160)	4 (0.250)	4/132	1.1	18.8%	coding (342/510 nt)	purE	5‑(carboxyamino)imidazole ribonucleotide mutase
	?	NZ_CP009273	548231 =	16 (1.010)					coding (326/510 nt)	purE	5‑(carboxyamino)imidazole ribonucleotide mutase
Rejected: Frequency below/above cutoff threshold.

CATGAGTCGCAAGAATTTGTGCTGCCAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCACGCCCAG  >  NZ_CP009273/548136‑548292                                                                                |                                                                              cATGAGTCGCAAGAATTTGTGCTGCCAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCcgcgcg                                                                     <  2:340858/90‑1 (MQ=255) cATGAGTCGCAAGAATTTGTGCTGCCAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCcgcgcg                                                                     <  1:503902/90‑1 (MQ=255)  aTGAGTCGCAAGAATTTGTGCTGCCAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCgg                                                                    >  2:370245/1‑90 (MQ=255)  aTGAGTCGCAAGAATTTGTGCTGCCAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCgg                                                                    >  2:212728/1‑90 (MQ=255)  aTGAGTCGCAAGAATTTGTGCTGCCAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCgg                                                                    >  1:595806/1‑90 (MQ=255)    gAGTCGCAAGAATTTGTGCTGCCAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCa                                                                  >  2:510059/1‑90 (MQ=255)    gAGTCGCAAGAATTTGTGCTGCCAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCa                                                                  <  2:176127/90‑1 (MQ=255)              aaTTTGTGCTGCCAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGAt                                                        <  2:18283/90‑1 (MQ=255)                tttGTGCTGCCAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGATgg                                                      <  2:140444/90‑1 (MQ=255)                      ctgcCAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGATGGAGTaga                                                <  1:59087/90‑1 (MQ=255)                      ctgcCAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGAATTTGTACGATGGAGTaga                                                <  1:212728/90‑1 (MQ=255)                         tgccagtacagagcgccgcactgagcggtgtcgatagcctctactccatcgtacaAATGCcgcgc                                                                      <  1:42998‑M1/10‑1 (MQ=255)                         tgccagtacagagcgccgcactgagcggtgtcgatagcctctactccatcgtacaAATGCcgcgc                                                                      >  2:42998‑M1/56‑65 (MQ=255)                         tgccagtacagagcgccgcactgagcggtgtcgatagcctctactccatcgtacaAATGCCGCGCGGCATTTGTACGATGGAGTAGAgg                                              <  1:376131‑M1/34‑1 (MQ=255)                         ccAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGc                                             <  2:2896/90‑1 (MQ=255)                            gTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTA‑CCACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATc                                         <  1:553458/90‑1 (MQ=255)                             tAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATc                                         <  2:595806/90‑1 (MQ=255)                                  cctcgTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGAt                                                        <  1:329645/67‑1 (MQ=255)                                  cctcgTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGAt                                                        >  2:329645/4‑70 (MQ=255)                                                   gCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGctc                   <  2:514313/90‑1 (MQ=255)                                                             tCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGa                                                         >  2:318792/1‑42 (MQ=255)                                                             tCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGa                                                         <  1:318792/42‑1 (MQ=255)                                                             tCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCa         <  1:370245/90‑1 (MQ=255)                                                              cctctactccatcgtacaAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCa         <  2:376131‑M1/71‑1 (MQ=255)                                                                    cGTACCCACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCACGCCCAg  <  2:538117/90‑1 (MQ=255)                                                                        atcgtacaAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCa         >  2:307192‑M1/9‑79 (MQ=255)                                                                        atcgtacaAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCa         <  1:307192‑M1/71‑1 (MQ=255)                                                                          cACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATc                                         <  1:143638/45‑1 (MQ=255)                                                                          cACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATc                                         >  2:143638/1‑45 (MQ=255)                                                                                |                                                                              CATGAGTCGCAAGAATTTGTGCTGCCAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCACGCCCAG  >  NZ_CP009273/548136‑548292

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