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.

GCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTGGCAATATACGGTGCTGCCACATCTTCGTTGATATCGAACGACAGTTTTACAT  >  NZ_CP009273/2685677‑2685850                                                                                      |                                                                                         gcacGTCACTCATATGAACGTCGATAGCATCAAACCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCAc                                                                                      <  1:215958/90‑1 (MQ=255)   agGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGc                                                                                    <  1:133420/90‑1 (MQ=255)             tatGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCaa                                                                          <  1:421610/90‑1 (MQ=255)             ctttgatgctatcgacgtgcatatgagtgacctgctgaccggacgcacgggcctggaagatttccacgccctgCTCACGCAGTACGGCaa                                                                          <  1:73278‑M2/17‑1 (MQ=255)                   cGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCAc                                                                                      <  2:473307/72‑1 (MQ=255)                   cGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCAc                                                                                      >  1:473307/1‑72 (MQ=255)                   cGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCg                                                                    >  1:122194/1‑90 (MQ=255)                   cGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCg                                                                    >  2:584502/1‑90 (MQ=255)                   cGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCg                                                                    >  2:583900/1‑90 (MQ=255)                   cGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCg                                                                    >  1:72852/1‑90 (MQ=255)                     tCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGa                                                                  >  2:581292/1‑90 (MQ=255)                       gATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACg                                                                >  2:231795/1‑90 (MQ=255)                        tcgacgtgcatatgagtgacctgctgaccggacgcacgggcctggaagatttccacgccctgCTCACGCAGTACGGCAACTTTCGGACGt                                                               <  1:279814‑M2/28‑1 (MQ=255)                           gCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTgcg                                                            >  1:9630/1‑90 (MQ=255)                             aTCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGcc                                                          >  2:120038/1‑90 (MQ=255)                              tCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCa                                                         <  1:288966/90‑1 (MQ=255)                               cAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAg                                                        >  1:528401/1‑90 (MQ=255)                                catatgagtgacctgctgaccggacgcacgggcctggaagatttccacgccctgCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGt                                                       <  2:420019‑M2/36‑1 (MQ=255)                                     gagtgacctgctgaccggacgcacgggcctggaagatttccacgccctgCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTGGCaa                                                  <  1:189276‑M2/41‑1 (MQ=255)                                                      gccgccATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTgg                                                     <  1:31816/70‑1 (MQ=255)                                                      gccgccATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTgg                                                     >  2:31816/1‑70 (MQ=255)                                                      gccgccATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTGGCAATATACGGTGCTGCcaca                                 >  1:90325/1‑90 (MQ=255)                                                      gccgccATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTGGCAATATACGGTGCTGCcaca                                 >  1:150170/1‑90 (MQ=255)                                                         gcacgggcctggaagatttccccgccctgCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTGGCAATATACGGTGCTGCCACATCt                              <  1:241646‑M2/61‑1 (MQ=255)                                                                      gCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTGGCAATATACGGTGCTGCCACATCTTCGTTGATATCGa                 <  2:377764/90‑1 (MQ=255)                                                                           ttCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTGGCAATATACGGTGCTGCCACATCTTCGTTGATATCGAACGAc            <  1:120038/90‑1 (MQ=255)                                                                                   ctgCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTGGCAATATACGGTGCTGCCACATCTTCGTTGATATCGAACGACAGTTTTAc    <  2:66993‑M2/87‑1 (MQ=255)                                                                                     gCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTGGCAATATACGGTGCTGCCACATCTTCGTTGATATCGAACGACAGTTTTACAt  >  1:464546‑M2/2‑90 (MQ=255)                                                                                      |                                                                                         GCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTGGCAATATACGGTGCTGCCACATCTTCGTTGATATCGAACGACAGTTTTACAT  >  NZ_CP009273/2685677‑2685850

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