New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1979271 =	0 (0.000)	15 (1.480)	13/298	0.0	100%	intergenic (‑56/‑482)	insA5/uspC	IS1 family protein InsA/universal stress protein C
	?	NC_000913	= 2102947	NA (NA)					pseudogene (443/446 nt)	wbbL	interrupted rhamnosyltransferase WbbL

GCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGAT  >  NC_000913/2102830‑2103113                                                                                                                      |                                                                                                                                                                       gcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATAttt                                   >  1:44999/1‑251 (MQ=255)                                  gTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatgat  <  1:116325/251‑1 (MQ=255)                                    tgatgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCATAAAGCACAGGAAATTTTc                                                                                                                        <  1:73605/131‑1 (MQ=255)                                                                                            tCAGGNAGTTGNTNNAACCTTCCTTAATCCCTAANATAAATGACACAATANAATCNGAAAGNACNg                                                                                                                                 >  2:99681/1‑66 (MQ=255)                                                                                            tCAGGGACTTGTTCGCACCTTCCGTAATCCCTAACATAAAGGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGCAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAgg                                                                            <  1:99681/119‑1 (MQ=255)                                                                                                                      |                                                                                                                                                                       GCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGAT  >  NC_000913/2102830‑2103113

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