New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1979271 =	0 (0.000)	19 (1.080)	15/292	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

GCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTT  >  NC_000913/2102723‑2103166                                                                                                                                                                                                                                 |                                                                                                                                                                                                                            nccATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATaaa                                                                                                                                                                                                    >  1:130357/2‑250 (MQ=255)          aGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAAtaa                                                                                                                                                                                                >  1:37097/1‑243 (MQ=255)                                                tttGCCGTGGCAGAATCTGCTCCATGNGGGACANNAANATATCTTTTCTGGGCTAACNGCNNNTNNNGNGGAATNCANNGTCGACGAANGTAAGTTGATGGCTCATGATGAACACTGTTCTATGGTTCGAGATGACCGACATNATCTCATATNNGGGCCTTGTTCGCACCTTCCTNAAtnnanntncccnca                                                                                                                                                                                                               >  2:103115/1‑179 (MQ=255)                                                  tGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTc                                                                                                                                                   <  1:37395/250‑1 (MQ=255)                                                   gCCATGGCAGCATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATACATGACACAATAAAATCAGACAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGc                                                                                                                                                 >  1:105138/1‑251 (MQ=255)                                                                             naCAAGCAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGca                                                                                                                                                                                         >  1:359578/2‑185 (MQ=255)                                                                                                                 ctgctgAATTCACTGACGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCAACGGCTCCAGATGACAAACATGATCTCATATAAGGGACTTGTTAGAACCTTCCTTACTCCCTACCATCAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCg                                                                                                                                                  >  1:300058/1‑188 (MQ=21)                                                                                                                   nctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGaataat                                                                                                                                                                 >  1:44685/2‑171 (MQ=255)                                                                                                                                                     cTCATGATGAACCCTGTTCTATGGCTCCAGATGACAACCATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATc                                                                                                                                                                                                 >  1:288328/1‑105 (MQ=255)                                                                                                                                                               aCCCTGNNCTATGGCGCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCa                                                                                                                                                                                                <  2:37097/96‑1 (MQ=255)                                                                                                                                                                  ctggtCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATTACACAATAAAATCAGAAAGCACAGGAAATGTTCTTACGGAATAATCATGTAAAGATTTCGCTTCAtctag                                                                                                                                        <  1:430182/145‑3 (MQ=255)                                                                                                                                                                           ggCTCCAGATGACAAACATGATCTCATATCCGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACACTAAAATCCGAAAGCACAGGAACTTTTCTTACGGACTAATCATg                                                                                                                                                             >  1:393372/1‑119 (MQ=25)                                                                                                                                                                                                   cATATCAGGGACTTGTTCACACCTTCCATAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGTCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCATCAAATCATCAGGCTGCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATAttttn  <  1:91840/250‑2 (MQ=255)                                                                                                                                                                                                            gACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAatata                                                                                  <  1:224082/161‑1 (MQ=255)                                                                                                                                                                                                             aCTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAAGTATGTAAAGATTTCGCTTCATCTCGGGTCAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAg                                                                            <  1:425848/166‑1 (MQ=255)                                                                                                                                                                                                                         nttccttAATCCCTAACATAAATGACACAATAAAATCAGNAAGCNCNNGANANTTNCTTAc                                                                                                                                                                         <  2:204448/60‑1 (MQ=255)                                                                                                                                                                                                                         cttccttAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTAc                                                                                                                                                                         >  1:204448/1‑61 (MQ=255)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                            GCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTT  >  NC_000913/2102723‑2103166

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