New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1979271 =	0 (0.000)	9 (0.900)	8/322	0.1	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

GGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCT‑CGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCA  >  NC_000913/2102740‑2103147                                                                                                                                                                                                                |                                                                                                                                                                                                          ggTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACCAGAAAAAATCTTTTCTGGTCTCCCGGCGCTTACTGCTGAATTACCTGTCGGCGAAGGTAAGATGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGaccaagaaaggc                                                                                                                                                                 >  1:104779/1‑239 (MQ=14)          tcTCGATGAATTCCACCATGTTTTGCCATGATAGAATCTGCCACATGCGGGACAAGAACATCTCTTTTCTGGTCTGACGACGCTTACTGCTGAATTCACTGTCGGCGAAGGTACGTTGATGACTCATGATGAACCATGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTAGCAACATCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTAc                                                                                                                                                       >  1:185934/2‑251 (MQ=17)                         ccaTGGTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACAAGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTa                                                                                                                                         >  1:70420/1‑250 (MQ=255)                                            aaTCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCCCCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATGTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCT‑Cgg                                                                                                                    <  1:187486/251‑1 (MQ=255)                                                                                                            cTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGaa                                                                                                                                                   >  1:16041/1‑157 (MQ=255)                                                                                                              gtcgtcGAAGTTAATTTGATGACTCATGAGAAACCCTGTTCTATGGCTCCAGTTGACAAACATGATTTCATATCAGGGACTTGTTCGCACCTTCTTTAATCCCTAACATAAAGGACACAAGAAAGTTAGTAAGCACAGAAATTTTTCTTTCGGAATAGCCATGTAAAGTTTTCG‑TTGCGCTGCGGAACAGTCTTTACGGAC‑GAAAGCATATCTCTTACTTTCGACATATTTAATATATGTCAGCAAatcat                                                 <  1:108033/251‑1 (MQ=255)                                                                                                                                                                                ntcATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCT‑CGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatcatca  >  1:24248/2‑233 (MQ=255)                                                                                                                                                                                                                |                                                                                                                                                                                                          GGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCT‑CGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCA  >  NC_000913/2102740‑2103147

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