New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 257907	2 (0.300)	5 (0.790)	4/238	0.3	84.2%	intergenic (+8/+16)	crl/insB9	RNA polymerase holoenzyme assembly factor Crl/IS1 transposase B
	?	NC_000913	258684 =	0 (0.000)					pseudogene (9/331 nt)	crl	RNA polymerase holoenzyme assembly factor Crl

AAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTG‑GACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGT  >  NC_000913/257781‑258007                                                                                                                                |                                                                                                     aaaTTTTTGCACTTGGGAAAAAAAGTGGC‑TCTAAACATGGGATAGCAATGACCTTACCG‑GTGTGAACCCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAaggtaagt                                                                                                   <  1:430536/129‑1 (MQ=255) aaaTTTTTCCATTTTGG‑AAAA‑AGTTTGCTCCCAAAAGGGGAAAAGAAAGACGGTTCCGGGTG‑GACACCCGAAGGGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAaggtaagtg                                                                                                  <  1:75603/129‑2 (MQ=255)                                      aGGAGATAGCAATGACGTTACCGAGTG‑GACACCCGAAGAGCAGATTGATCAAAAAATTGACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctg                                                               <  1:544507‑M1/129‑40 (MQ=255)                                               cAATGACGTTACCGAGTG‑GAACCCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtca                                                      <  1:550532‑M1/129‑49 (MQ=255)                                                                                         aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacttgaacccggcaccggaagtgcgtgacttctggggcttg             >  1:544545‑M1/1‑40 (MQ=255)                                                                                         aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggccccggaagggcgtgaattttggggctgg             >  1:320790‑M1/1‑40 (MQ=255)                                                                                                    cGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaatgcgtggaatcttgggccggggggtgggacct  >  1:329576‑M1/1‑29 (MQ=255)                                                                                                                                |                                                                                                     AAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTG‑GACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGT  >  NC_000913/257781‑258007

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