New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	274038 =	NA (NA)	3 (0.040)	3/166	NT	NA	noncoding (1112/1195 nt)	IS5	repeat region
	?	NC_000913	274083 =	NA (NA)					noncoding (1067/1195 nt)	IS5	repeat region

CTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAAT  >  NC_000913/273939‑274051                                                                                                    |              cTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGt               <  2:2233021/100‑1 (MQ=255)  tGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTa              <  2:872141/100‑1 (MQ=255)   gACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTaa             <  1:1921980/100‑1 (MQ=255)     cAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATGCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAAt           >  2:1201192/1‑100 (MQ=255)      aGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATg          <  2:1076837/100‑1 (MQ=38)       ncACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGc         >  2:83093/2‑100 (MQ=25)       gCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGc         >  2:19242/1‑100 (MQ=38)        cACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGcc        <  1:1965861/100‑1 (MQ=35)          cTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCtt      <  1:545162/100‑1 (MQ=32)           tGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTg     <  1:208374/100‑1 (MQ=32)           tGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTg     <  2:1589780/100‑1 (MQ=33)           tGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTg     <  2:1685762/100‑1 (MQ=33)           tGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTg     <  2:1776030/100‑1 (MQ=33)             ttATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTTATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGaa   <  2:528646/100‑1 (MQ=17)              tATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAAt  <  1:1396903/100‑1 (MQ=32)                                                                                                    |              CTGACAGCACTGTTATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAAT  >  NC_000913/273939‑274051

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