New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AM260480	= 1197838	9 (0.410)	4 (0.200)	4/186	1.6	30.5%	coding (970/1437 nt)	h16_B1061	Nucleotidyltransferase/DNA polymerase involved in DNA repair
	?	AM260480	= 1197865	10 (0.500)					coding (997/1437 nt)	h16_B1061	Nucleotidyltransferase/DNA polymerase involved in DNA repair

CGCGCACCTGTCGCGGCTGCTGCGCGAGAAGCTCGACAAGCTCAATTTCCACGCCCCGGTGACCGGCCTGGCGCTGCGCGTTGAGGCGATGGCCGCATGCGTGCCGCACAGCGATGCGCTGTTCCCCGAGCCCGGCGGCACGCCGGAAGACCTGGGCCGCCTGCTCGATACGCTGATGGCGCGGCTGGGCCGTGACAACGTGCTGCAGCCGCGGCCGCTGGCCGACCACCGGCCCGAGCGCGCCAA  >  AM260480/1197732‑1197977                                                                                                                                      |                                                                                                                 cgcgCACCTGTCGCGGCTGCTGCGCGAGAAGCTCGACAAGCTCAATTTCCACGCCCCGGTGACCGGCCTGGCGCTGCGCGTTGAGGCGATGGCCGCATGCGTGCCGCACAGCGATGCGCTGTTCCCCGAGCCcg                                                                                                                  >  2:431622/1‑134 (MQ=255)                         cgAGAAGCTCGACAAGCTCAATTTCCACGCCCCGGTGACCGGCCTGGCGCTGCGCGTTGAGGCGATGGCCGCATGCGTGCCGCACAGCGATGCGCTGTTCCCCGAGCCCGGCGGCACGCCGGAAGACCTGGgcc                                                                                          <  2:913691/134‑1 (MQ=255)                                            aaTTTCCACGCCCCGGTGACCGGCCTGGCGCTGCGCGTTGAGGCGATGGCCGCATGCGTGCCGCACAGCGATGCGCTGTTCCCCGAGCCCGGCGGCACGCCGGAAGACCTGGGCCGCCTGCTCGATACGCTg                                                                         <  1:1013848/132‑1 (MQ=255)                                                            tGACCGGCCTGGCGCTGCGCGTTGAGGCGATGGCCGCATGCGTGCCGCACAGCGATGCGCTGTTCCCCGAGCCCGGCGGCACGCCGGAAGACCTGGGCCGCCTGCTCGATACGCTGATGGCGCGGCTGGGCCGt                                                       <  1:431622/134‑1 (MQ=255)                                                                       gcgcTGCGCGTTGAGGCGATGGCCGCATGCGTGCCGCACAGCGATGCGCTGTTCCCCGAGCCCGGCGGCACGCCGGAAGACCTGGGCCGCCTGCTCGATACGCTGATGGCGCGGCTGGGCCGTGACAACGtgc                                             <  2:210255/133‑1 (MQ=255)                                                                        cgcTGCGCGTTGAGGCGATGGCCGCATGCGTGCCGCACAGCGATGCGCTGTTCCCCGAGCCCGGCGGCACGCCGGAAGACCTGGGCCGCCTGCTCGATACGCTGATGGCGCGGCTGGGCCGTGACAACGtgctg                                           <  1:1057289/134‑1 (MQ=255)                                                                                                         cGCACAGCGATGCGCTGTTCCCCGAGCCCGGCGGCACGCCGGAAGACCTGGGCCGCCTGCTCGATACGCTGATGGCGCGGCTGGGCCGTGACAACGTGCTGCAGCCGCGGCCGCTGGCCGACCACCGGCCCGAg          >  2:861095/1‑134 (MQ=255)                                                                                                          gCACAGCGATGCGCTGTTCCCCGAGCCCGgcggcacgcatgcggccatcgcctcaa                                                                                       <  1:438851‑M2/56‑28 (MQ=255)                                                                                                          gCACAGCGATGCGCTGTTCCCCGAGCCCGgcggcacgcatgcggccatcgcctcaa                                                                                       >  2:438851‑M2/1‑29 (MQ=255)                                                                                                            acaGCGATGCGCTGTTCCCCGAGCCCGGCGGCACGCCGGAAGACCTggg                                                                                            >  1:145460/1‑49 (MQ=255)                                                                                                            acaGCGATGCGCTGTTCCCCGAGCCCGGCGGCACGCCGGAAGACCTggg                                                                                            <  2:145460/49‑1 (MQ=255)                                                                                                                       cTGTTCCCCGAGCCCGgcggcacgcatgcggccatcgcatcaacgcgcagcgccaggccggtcaccggggcgtggaaattgagcttgtcgagcttctcgc                              <  1:274705‑M2/100‑85 (MQ=255)                                                                                                                       cTGTTCCCCGAGCCCGgcggcacgcatgcggccatcgcatcaacgcgcagcgccaggccggtcaccggggcgtggaaattgagcttgtcgagcttctcgc                              >  2:274705‑M2/1‑16 (MQ=255)                                                                                                                                gAGCCCGGCGGCACGCCGGAAGACCTGGGCCGCCTGCTCGATACGCTGATGGCGCGGCTGGGCCGTGACAACGTGCTGCAGCCGCGGCCGCTGGCCGACCACCGGCCCGAGCGCGCCaa  <  1:139194/119‑1 (MQ=255)                                                                                                                                gAGCCCGGCGGCACGCCGGAAGACCTGGGCCGCCTGCTCGATACGCTGATGGCGCGGCTGGGCCGTGACAACGTGCTGCAGCCGCGGCCGCTGGCCGACCACCGGCCCGAGCGCGCCaa  >  2:139194/1‑119 (MQ=255)                                                                                                                                      |                                                                                                                 CGCGCACCTGTCGCGGCTGCTGCGCGAGAAGCTCGACAAGCTCAATTTCCACGCCCCGGTGACCGGCCTGGCGCTGCGCGTTGAGGCGATGGCCGCATGCGTGCCGCACAGCGATGCGCTGTTCCCCGAGCCCGGCGGCACGCCGGAAGACCTGGGCCGCCTGCTCGATACGCTGATGGCGCGGCTGGGCCGTGACAACGTGCTGCAGCCGCGGCCGCTGGCCGACCACCGGCCCGAGCGCGCCAA  >  AM260480/1197732‑1197977

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