New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AM260479	3854038 =	11 (0.430)	5 (0.210)	4/200	1.9	33.0%	intergenic (+52/+36)	h16_A3587/h16_A3588	Hypothetical protein/putative peptidase, M20A subfamily
	?	AM260479	3854056 =	10 (0.420)					intergenic (+70/+18)	h16_A3587/h16_A3588	Hypothetical protein/putative peptidase, M20A subfamily

GAGTCCGCCGCGTGGTGTGCCGCGATATTGCGGAAAGCTTGAGAACAAGAAAGCAAAGCGCCACGGGGTCTCCGTGGCGCTTTTGCGTTGGCCGACCGCCTGTTGGCGGTCCGGCTCAATCGAAATGACTTACTTCATCAGGTCCGACAGCGCCAGCGCCACCACCTTGGTCGCGCGGTTCAGGTCGTTCAGGCGCAGGTTCTCGTCCGAGTTATGGCCGCGGG  >  AM260479/3853945‑3854168                                                                                                                |                                                                                                                 gagaacctgcgcctgaacgacctgaaccgcgcgaccaaggtggtggcgctggcgctgtcggacctgatgaagtaagtcatttcgattgagccggaccgccaacaggcggtcCGGCTCAATCGAAATGacttact                                                                                            <  2:662618‑M2/23‑1 (MQ=255) gagTCCGCCGCGTGGTGTGCCGCGATATTGCGGAAAGCTTGAGAACAAGAAAGCAAAGCGCCACGGGGTCTCCGTGGCGCTTTTGCGTTGGCCGACCGCCTGTTGGCGGTCCGGCTCAATCGAAATGacttact                                                                                            >  2:289410/1‑134 (MQ=255) gagTCCGCCGCGTGGTGTGCCGCGATATTGCGGAAAGCTTGAGAACAAGAAAGCAAAGCGCCACGGGGTCTCCGTGGCGCTTTTGCGTTGGCCGACCGCCTGTTGGCGGTCCGGCTCAATCGAAATGacttact                                                                                            >  2:917217/1‑134 (MQ=255)  agaacctgcgcctgaacgacctgaaccgcgcgaccaaggtggtggcgctggcgctgtcggacctgatgaagtaagtcatttcgattgagccggaccgccaacaggcggtcCGGCTCAATCGAAATGacttactt                                                                                           <  1:74377‑M2/24‑1 (MQ=255)  agaacctgcgcctgaacgacctgaaccgcgcgaccaaggtggtggcgctggcgctgtcggacctgatgaagtaagtcatttcgattgagccggaccgccaacaggcggtcCGGCTCAATCGAAATGacttactt                                                                                           >  2:74377‑M2/111‑134 (MQ=255)               gtgtGCCGCGATATTGCGGAAAGCTTGAGAACAAGAAAGCAAAGCGCCACGGGGTCTCCGTGGCGCTTTTGCGTTGGCCGACCGCCTGTTGGCGGTCCGGCTCAATCGAAATGACTTACTTCATCAGGTCCGa                                                                               >  2:868802/1‑133 (MQ=255)                           atTGCGGAAAGCTTGAGAACAAGAAAGCAAAGCGCCACGGGGTCTCCGTGGCGCTTTTGCGTTGGCCGACCGCCTGTTGGCGGTCCGGCTCAATCGAAATGACTTACTTCATCAGGTCCGAcagcgccagcgc                                                                   >  2:718193/1‑133 (MQ=255)                              gCGGAAAGCTTGAGAACAAGAAAGCAAAGCGCCACGGGGTCTCCGTGGCGCTTTTGCGTTGGCCGACCGCCTGTTGGCGGTCCGGCTCAATCGAAATGACTTACTTCATCAGGTCCGACAGCGCCAGCGccac                                                                >  1:836854/1‑133 (MQ=255)                                  accaaggtggtggcgctggcgctgtcggacctgatgaagtaagtcatttcgattgagccggaccgccaacaggcggtcCGGCTCAATCGAAATGacttactt                                                                                           <  1:205923‑M2/24‑1 (MQ=255)                                  accaaggtggtggcgctggcgctgtcggacctgatgaagtaagtcatttcgattgagccggaccgccaacaggcggtcCGGCTCAATCGAAATGacttactt                                                                                           >  2:205923‑M2/79‑102 (MQ=255)                                        tGAGAACAAGAAAGCAAAGCGCCACGGGGTCTCCGTGGCGCTTTTGCGTTGGCCGACCGCCTGTTGGCGGTCCGGCTCAATCGAAATGACTTACTTCATCAGGTCCGAcagc                                                                           <  2:325976/112‑1 (MQ=255)                                        tGAGAACAAGAAAGCAAAGCGCCACGGGGTCTCCGTGGCGCTTTTGCGTTGGCCGACCGCCTGTTGGCGGTCCGGCTCAATCGAAATGACTTACTTCATCAGGTCCGAcagc                                                                           >  1:325976/1‑112 (MQ=255)                                          agaACAAGAAAGCAAAGCGCCACGGGGTCTCCGTGGCGCTTTTGCGTTGGCCGACCGCCTGTTGGCGGTCCGGCTCAATCGAAATGACTTACTTCATCAGGTCCGACAGCGCCAGCGCCACCACCTTGGTcgcg                                                   >  1:652427/1‑134 (MQ=255)                                                                      aagtaagtcatttcgattgagccgGACCGCCTGTTGGCGGTCCGGCTCAATCGAAATGACTTACTTCATCAGGTCCGACAGCGCCAGCGCCACCACCTTGGt                                                       <  1:205923‑M1/78‑1 (MQ=255)                                                                      aagtaagtcatttcgattgagccgGACCGCCTGTTGGCGGTCCGGCTCAATCGAAATGACTTACTTCATCAGGTCCGACAGCGCCAGCGCCACCACCTTGGt                                                       >  2:205923‑M1/25‑102 (MQ=255)                                                                      aagtaagtcatttcgattgagccgGACCGCCTGTTGGCGGTCCGGCTCAATCGAAATGACTTACTTCATCAGGTCCGACAGCGCCAGCGCCACCACCTTGGTCGCGCGGTTCAGGTCGTTCAGGCGCAGGTtct                       >  2:74377‑M1/25‑134 (MQ=255)                                                                      aagtaagtcatttcgattgagccgGACCGCCTGTTGGCGGTCCGGCTCAATCGAAATGACTTACTTCATCAGGTCCGACAGCGCCAGCGCCACCACCTTGGTCGCGCGGTTCAGGTCGTTCAGGCGCAGGTtct                       <  1:74377‑M1/110‑1 (MQ=255)                                                                       agtaagtcatttcgattgagccgGACCGCCTGTTGGCGGTCCGGCTCAATCGAAATGACTTACTTCATCAGGTCCGACAGCGCCAGCGCCACCACCTTGGTCGCGCGGTTCAGGTCGTTCAGGCGCAGGTtctc                      <  2:662618‑M1/111‑1 (MQ=255)                                                                         cGTGGCGCTTTTGCGTTGGCCGACCGCCTGTTGGCGGTCCGGCTCAATCGAAATGACTTACTTCATCAGGTCCGACAGCGCCAGCGCCACCACCTTGGTCGCGCGGTTCAGGTCGTTCAGGCGCAGGTtct                       <  2:836854/131‑1 (MQ=255)                                                                                           gCCGACCGCCTGTTGGCGGTCCGGCTCAATCGAAATGACTTACTTCATCAGGTCCGACAGCGCCAGCGCCACCACCTTGGTCGCGCGGTTCAGGTCGTTCAGGCGCAGGTTCTCGTCCGAGTTATGGCCGCggg  >  1:735142/1‑134 (MQ=255)                                                                                                                |                                                                                                                 GAGTCCGCCGCGTGGTGTGCCGCGATATTGCGGAAAGCTTGAGAACAAGAAAGCAAAGCGCCACGGGGTCTCCGTGGCGCTTTTGCGTTGGCCGACCGCCTGTTGGCGGTCCGGCTCAATCGAAATGACTTACTTCATCAGGTCCGACAGCGCCAGCGCCACCACCTTGGTCGCGCGGTTCAGGTCGTTCAGGCGCAGGTTCTCGTCCGAGTTATGGCCGCGGG  >  AM260479/3853945‑3854168

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