New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 381473	5 (0.030)	6 (0.040)	4/446	12.8	55.7%	noncoding (214/1331 nt)	IS2	repeat region
	?	NC_000913	= 381839	NA (NA)					noncoding (580/1331 nt)	IS2	repeat region
Rejected: Coverage evenness skew score above cutoff.

CCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCA  >  NC_000913/381235‑381505                                                                                                                                                                                                                                               |                                 ccAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGccagccag                      <  2:1387280/251‑1 (MQ=39)          aTAAACGCATAATTACTAGACTGGCCACCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGAAATTGTTCAGCAGAGATTTGAACCGGGGATGACGGTCTCCCGCGTTGCCCGGCAACATGGGGTAGCAGCCAGCCAGTTATTtctc             >  1:863272/1‑251 (MQ=255)               aagacagttACTTGACTGGCCCCCTGAATCTCCAGACAACCAATATAACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCg        <  1:851835/244‑1 (MQ=11)                   tAATTACTGTACTGGCCCCCTCAATCTCCAAATAACGAATATTGCTTAAAAAAGTGATAGTCTTAGTGCTAGTTATGAGACTAGTCAGTGGAAAACAGATGATGGACGTCTTAGGGCCGGAGAAACGCAGGCGGCGTACCACACAGGAAGAGAGCGCAAGTGTTCAGCAGAGCTTTGAACCGGGGAGGACGGGCTCCCTCGGTGACCGGCAACATGGTGGAGGAGCCAGCCAGTTATTTCTCTGGCGTAAg    <  1:784145/251‑1 (MQ=11)                     aTTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCa  <  2:863070/251‑1 (MQ=31)                                                                                                                                                                                                                                               |                                 CCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCA  >  NC_000913/381235‑381505

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