New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1633763 =	NA (NA)	4 (0.170)	3/118	2.2	100%	coding (142/591 nt)	pinQ	Qin prophage; putative site‑specific recombinase
	?	NC_000913	1633785 =	0 (0.000)					coding (164/591 nt)	pinQ	Qin prophage; putative site‑specific recombinase

AGCAAATAATCGAAGAACACATTAGCGGCTCAGCAGCAACCAGTGAGCGTCCTGGTTTTAACCGGTTGCTTGCTCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAAT  >  NC_000913/1633722‑1633852                                                                |                                                                    cacaatcaattgatcaccacatttcaggcgagcaagcaaccggttaaaaccaggacgctcactGGttg                                                                 >  1:155178‑M2/64‑68 (MQ=255) cacaatcaattgatcaccacatttcaggcgagcaagcaaccggttaaaaccaggacgctcactGGttg                                                                 >  1:288614‑M2/64‑68 (MQ=255) cacaatcaattgatcaccacatttcaggcgagcaagcaaccggttaaaaccaggacgctcactGGttg                                                                 >  1:522861‑M2/64‑68 (MQ=255) cacaatcaattgatcaccacatttcaggcgagcaagcaaccggttaaaaccaggacgctcactGGttg                                                                 >  1:85297‑M2/64‑68 (MQ=255)   caatcaattgatcaccacatttcaggcgagcaagcaaccggttaaaaccaggacgctcactGGttgct                                                               <  2:154489‑M2/7‑1 (MQ=255)   caatcaattgatcaccacatttcaggcgagcaagcaaccggttaaaaccaggacgctcactGGttgct                                                               <  2:321456‑M2/7‑1 (MQ=255)                                                aaccaggacgctcactGGTTGCTTGCTCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGa                  >  1:438195‑M2/17‑68 (MQ=255)                                                 accaggacgctcactGGTTGCTTGCTCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGAt                 >  2:307538‑M2/16‑68 (MQ=255)                                                 accaggacgctcactGGTTGCTTGCTCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGAt                 >  2:415669‑M2/16‑68 (MQ=255)                                                  ccaggacgctcactGGTTGCTTGCTCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATc                >  1:536095‑M2/15‑68 (MQ=255)                                                                GGTTGCTTGCTCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTCGTAAt  <  2:711108‑M2/68‑1 (MQ=255)                                                                |                                                                    AGCAAATAATCGAAGAACACATTAGCGGCTCAGCAGCAACCAGTGAGCGTCCTGGTTTTAACCGGTTGCTTGCTCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAAT  >  NC_000913/1633722‑1633852

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