New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	3717977 =	7 (0.630)	3 (0.310)	3/144	0.9	33.3%	coding (702/912 nt)	glyQ	glycine‑‑tRNA ligase subunit alpha
	?	NZ_CP009273	3717989 =	6 (0.630)					coding (690/912 nt)	glyQ	glycine‑‑tRNA ligase subunit alpha

GCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTT  >  NZ_CP009273/3717900‑3718051                                                                                          |                                                               gcATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCtgctgc                                                                >  1:39543/1‑90 (MQ=255)                     gaatacgcggatgtggacttcctgttcacctgcttcgagcagtacgagaaagaagcgcagcagctgctgCGCTTCTTTCTCGTACTGCTc                                            <  1:222305‑M2/21‑1 (MQ=255)                     gaatacgcggatgtggactccctgttcccctgcttcgagcagtacgagaaagaagcgCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTc                                            >  1:314331‑M1/58‑90 (MQ=255)                              gatgtggacttcctgttcacctgcttcgagcagtacgagaaagaagcgcagcagctgctgCGCTTCTTTCTCGTACTGCTCGAAGCAGGt                                   <  2:181548‑M2/30‑1 (MQ=255)                                      tACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGa                           >  2:89847/1‑90 (MQ=255)                                                acctgcttcgagcagtacgagaaagaagcgCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATc                 <  2:181548‑M1/60‑1 (MQ=255)                                                    ggcagcggcagcGGATTTTCCAGCGCCAGCAGCTGCTGCGcttc                                                           <  1:304454/44‑1 (MQ=255)                                                    ggcagcggcagcGGATTTTCCAGCGCCAGCAGCTGCTGCGcttc                                                           >  2:304454/1‑44 (MQ=255)                                                    ggcagcggcagcGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCcgcg             >  2:175683/1‑90 (MQ=255)                                                         gagcagtacgagaaagaagcgcagcagctgctgCGCTTCTTTCTCGTACTGCTCGAAGCAGGGGAACAGGGAGTCCACATCCGCGTATTc        >  1:314331‑M2/34‑90 (MQ=255)                                                         gagcagtacgagaaagaagcgCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTc        <  1:222305‑M1/69‑1 (MQ=255)                                                               cggATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGtt  <  1:281814/90‑1 (MQ=255)                                                                                          |                                                               GCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTT  >  NZ_CP009273/3717900‑3718051

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