New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 1684417	11 (0.540)	4 (0.220) +CAGTGACAAC	3/146	1.5	21.6%	coding (309/1509 nt)	ydgA	YdgA family protein
	?	NZ_CP009273	= 3404280	22 (1.080)					coding (642/966 nt)	dusB	tRNA dihydrouridine synthase DusB
Rejected: Coverage evenness skew score above cutoff.

GAAAATCCGTGGATTAAAAGCGGTCAGAGCGTCATCTTCAACGAATCGGTTGATCATGGTCCCTTCCCGCTTGCCCAGCTTAAAAAACTGAACCTGATCCCGTCGATGGCATCAATTCAAACCACGCTGGTTAATAACGAAGTAAGCAAACCACTGTTTG  >  NZ_CP009273/1684346‑1684505                                                                        |                                                                                         gaaaATCCGTGGATTAAACGCGGTCAGAGCGGCATCTTCAACGAATCGGTTGATAATGGTCCCTTCCCGCCTGCCCAGCTGAAAAACCTg                                                                        <  1:597775/90‑1 (MQ=255) gaaaATCCGTGGATTAAAAGCGGTCAGAGCGTCATCTTCAACGAATCGGTTGATCATGGTCCCTTCCCGCTTGCCCAGCTTAAAAAACTg                                                                        <  1:351169/90‑1 (MQ=255)                                      tCAACGAATCGGTTGATCATGGTCCCTTCCCGCTTGCCCAGCTTAAAAAACTGAACCTGATCCCGTCGATGGCATCAATTCAAACCACGc                                   <  2:374604/90‑1 (MQ=255)                                           gAATCGGTTGATCATGGTCCCTTCCCGCTTGCCCAGCTTAAAAAACTGAACCTGATCCCGTCGATGGCATCAATTCAAACCACGCTGGtt                              <  1:492069/90‑1 (MQ=255)                                            aaTCGGTTGATCATGGTCCCTTCCCGCTTGCCCAGCTTAAAAAACTGAACCTGATCCCGTCGATGGCATCAATTCAAACCACGCTGGTta                             >  2:545171/1‑90 (MQ=255)                                                               cTTCCCGCTTGCCCAGCTTAAAAAACTGAACCTGATCCCGTCGATGGCATCAATTCAAAc                                        <  1:325672/60‑1 (MQ=255)                                                               cTTCCCGCTTGCCCAGCTTAAAAAACTGAACCTGATCCCGTCGATGGCATCAATTCAAAc                                        >  2:325672/1‑60 (MQ=255)                                                               cTTCCCGCTTGCCCAGCTTAAAAAACTGAACCTGATCCCGTCGATGGCATCAATTCAAACCACGCTGGTTAATAACGAAGTAAGCAAAcc          <  1:411355/90‑1 (MQ=255)                                                               cTTCCCGCTTGCCCAGCTTAAAAAACTGAACCTGATCCCGTCGATGGCATCAATTCAAACCACGCTGGTTAATAACGAAGTAAGCAAAcc          <  1:456929/90‑1 (MQ=255)                                                               cTTCCCGCTTGCCCAGCTTAAAAAACTGAACCTGATCCCGTCGATGGCATCAATTCAAACCACGCTGGTTAATAACGAAGTAAGCAAAcc          >  2:160883/1‑90 (MQ=255)                                                                       ttGCCCAGCTTAAAAAACTGAACCTGATCCCGTCGATGGCATCAATTCAAACCACGCTGGTTAATAACGAAGTAAGCAAACCACTGTTTg  >  2:494153/1‑90 (MQ=255)                                                                        |                                                                                         GAAAATCCGTGGATTAAAAGCGGTCAGAGCGTCATCTTCAACGAATCGGTTGATCATGGTCCCTTCCCGCTTGCCCAGCTTAAAAAACTGAACCTGATCCCGTCGATGGCATCAATTCAAACCACGCTGGTTAATAACGAAGTAAGCAAACCACTGTTTG  >  NZ_CP009273/1684346‑1684505

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