New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 2707188	20 (0.920)	3 (0.160)	3/148	1.6	13.2%	coding (934/1335 nt)	srmB	ATP‑dependent RNA helicase SrmB
	?	NZ_CP009273	= 2707267	22 (1.160)					coding (1013/1335 nt)	srmB	ATP‑dependent RNA helicase SrmB
Rejected: Frequency below/above cutoff threshold.

TACAGGGCAAGCGTAACGAAGCGATCAAGCGTTTGACCGAAGGTCGCGTAAACGTACTGGTCGCAACCGATGTTGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGG  >  NZ_CP009273/2707102‑2707264                                                                                       |                                                                             tACAGGGCAAGCGTAACGAAGCGATCAAGCGTTTGACCGAAGGTCGCGTAAACGTACTGGTCGCAACCGATGTTGCCGCGCGCGGTATCg                                                                           >  1:496501/1‑90 (MQ=255) tACAGGGCAAGCGTAACGAAGCGATCAAGCGTTTGACCGAAGGTCGCGTAAACGTACTGGTCGCAACCGATGTTGCCGCGCGCGGTATCg                                                                           >  2:376815/1‑90 (MQ=255)     gggCAAGCGTAACGAAGCGATCAAGCGTTTGACCGAAGGTCGCGTAAACGTACTGGTCGCAACCGATGTTGCCGCGCGCGGTATCGACAt                                                                       >  1:239575/1‑90 (MQ=255)         aaGCGTAACGAAGCGATCAAGCGTTTGACCGAAGGTCGCGTAAACGTACTGGTCGCAACCGATGTTGCCGCGCGCGGTAcgtccgatacg                                                                   >  2:356570‑M1/1‑79 (MQ=255)         aaGCGTAACGAAGCGATCAAGCGTTTGACCGAAGGTCGCGTAAACGTACTGGTCGCAACCGATGTTGCCGCGCGCGGTATCGACATTCCt                                                                   >  1:573129/1‑90 (MQ=255)          aGCGTAACGAAGCGATCAAGCGTTTGACCGAAGGTCGCGTAAACGTACTGGTCGCAACCGATGTTGCCGCGCGCGGTATCGACATTCCTg                                                                  >  2:208546/1‑90 (MQ=255)                    aGCGATCAAGCGTTTGACCGAATGGCGCGTAAACGTACTGGTCGCAACCGAGGTTGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCa                                                        <  1:131407/90‑1 (MQ=255)                    aGCGATCAAGCGTTTGACCGAAGGTCGCGTAAACGTACTGGTCGCAACCGATGTTGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCa                                                        <  2:189006/90‑1 (MQ=255)                            aGCGTTTGACCGAAGGTCGCGTAAACGTACTGGTCGCAACCGATGTTGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCACGTCTTTa                                                <  2:573129/90‑1 (MQ=255)                                        aaGGTCGCGTAAACGTACTGGTCGCAACCGATGTTGCCGCGCGCGGTAcgtccgatacggtgcaaataagtatcgccactgcgcggcata                                    >  1:83687‑M1/1‑48 (MQ=255)                                        aaGGTCGCGTAAACGTACTGGTCGCAACCGATGTTGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGc                                    <  2:592229/90‑1 (MQ=255)                                            tCGCGTAAACGTACTGGTCGCAACCGATGTTGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCcgcg                                >  2:616076/1‑90 (MQ=255)                                              gcgTAAACGTACTGGTCGCAACCGATGTTGCCGCGCGCGGTAcgtccgatacggtgcaaataagtatcgccactgcgcggcatatcgaag                              >  1:47693‑M1/1‑42 (MQ=255)                                              gcgTAAACGTACTGGTCGCAACCGATGTTGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCa                              <  2:507935/90‑1 (MQ=255)                                               cgTAAACGTACTGGTCGCAACCGATGTTGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAg                             <  1:51884/90‑1 (MQ=255)                                                       tACTGGTCGCAACCGATGTTGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCa                              >  2:122269/1‑81 (MQ=255)                                                       tACTGGTCGCAACCGATGTTGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCa                              <  1:122269/81‑1 (MQ=255)                                                        aCTGGTCGCAACCGATGTTGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATAc                    <  1:5180/90‑1 (MQ=255)                                                           ggTCGCAACCGATGTTGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTa                 >  1:123621/1‑90 (MQ=255)                                                                  aCCGATGTTGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCAc          <  2:168684/90‑1 (MQ=255)                                                                   ccGATGTTGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCAcc         <  1:65353/90‑1 (MQ=255)                                                                     gATGTTGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGt       >  2:42053/1‑90 (MQ=255)                                                                          tGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCgg  >  2:269126/1‑90 (MQ=255)                                                                                       |                                                                             TACAGGGCAAGCGTAACGAAGCGATCAAGCGTTTGACCGAAGGTCGCGTAAACGTACTGGTCGCAACCGATGTTGCCGCGCGCGGTATCGACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGG  >  NZ_CP009273/2707102‑2707264

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