New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AM260479	2538231 =	15 (0.590)	4 (0.170)	4/192	1.8	23.4%	coding (552/900 nt)	h16_A2338	conserved hypothetical protein
	?	AM260479	2538284 =	12 (0.500)					coding (499/900 nt)	h16_A2338	conserved hypothetical protein

GTGGCGGCGCACCAGCGTCTTGCGCTCTTCGGCGCTGGCGTGGTCGGCGAGCGCGTCATACGTCATCACGCGGTAGGTCAGCCAGCCCCACAGCACCGGCGGGATCAGCGCAAAGAACGGCGGGATCAGCCACAGCGGCAGCGTCACCAGCACCAGCACCAGGAACACCAGCGTGCTGCCCAGCGACTGGAACACGCTGCCCAGGATGC  >  AM260479/2538198‑2538406                                                                                       |                                                                                                                           gTGGCGGCGCACCAGCGTCTTGCGCTCTTCGGCGCTGGCGTGGTCGGCGAGCGCGTCATACGTCATCACGCGGTAGGTCAGCCAGCCCCACAGCACCGGCGGGATCAGCGCAAAGAACGGCGGGATCAGCca                                                                               >  1:548305/1‑132 (MQ=255)           aCCAGCGTCTTGCGCTCTTCGGCGCTGGCGTGGTCGGCGAGCGCGTCATACGTCATCACGCGGTAGGTCAGCCAGCCCCACAGCACCGGCGGGATCAGCGCAAAGAACGGCGGGATCAGCCACAGCGGCAGCGt                                                                   >  1:109908/1‑134 (MQ=255)                      gcgcTCTTCGGCGCTGGCGTGGTCGGCGAGCGCGTCATACGTCATCACGCGGTAGGTCAGCCAGCCCCACAGCACCGGCggg                                                                                                            >  2:799586/1‑82 (MQ=255)                      gcgcTCTTCGGCGCTGGCGTGGTCGGCGAGCGCGTCATACGTCATCACGCGGTAGGTCAGCCAGCCCCACAGCACCGGCggg                                                                                                            <  1:799586/82‑1 (MQ=255)                                gcgcTGGCGTGGTCGGCGAGCGCGTCATACGTCATCACGCGGTAGGTCAGCCAGCCCCACAGCACCGGCGGGATCAGCGCAAAGAACGGCGGGAt                                                                                     <  1:956527/95‑1 (MQ=255)                                gcgcTGGCGTGGTCGGCGAGCGCGTCATACGTCATCACGCGGTAGGTCAGCCAGCCCCACAGCACCGGCGGGATCAGCGCAAAGAACGGCGGGAt                                                                                     >  2:956527/1‑95 (MQ=255)                                      GCGTGGTCGGCGAGCGCGTCATACGTCATCACGCGGTAGGTCAGCCAGCCCCACAGCACCGGCGGGATCAGCGCAAAGAACGGCGGGATCAGCCACAGCGGCAGCGTCACCAGCACCAGCACCAGGAACACCAg                                        >  2:109908‑M1/1‑134 (MQ=255)                                         gacctaccgcgtgatgacgtatgacgcgctcgccgaccacgccagcCCCACAGCACCGGCGGGATCAGCGCAAAGAACGGCGGGATCAGCCACAGCGGCAGCGTc                                                                  >  1:282699‑M2/47‑105 (MQ=255)                                         gacctaccgcgtgatgacgtatgacgcgctcgccgaccacgccagcCCCACAGCACCGGCGGGATCAGCGCAAAGAACGGCGGGATCAGCCACAGCGGCAGCGTc                                                                  <  2:282699‑M2/59‑1 (MQ=255)                                          ggTCGGCGAGCGCGTCATACGTCATCACGCGGTAGGTCAGCCAGCCCCACAGCACCGGCGGGATCAGCGCAAAGAAcggcg                                                                                         <  2:274218/81‑1 (MQ=255)                                          ggTCGGCGAGCGCGTCATACGTCATCACGCGGTAGGTCAGCCAGCCCCACAGCACCGGCGGGATCAGCGCAAAGAAcggcg                                                                                         >  1:274218/1‑81 (MQ=255)                                                          cgtatgacgcgctcgccgaccacgccagcCCCACAGCACCGGCGGGATCAGCGCAAAGAACGGCGGGATCAGCCACAGCGGc                                                                        <  2:123223‑M2/53‑1 (MQ=255)                                                          cgtatgacgcgctcgccgaccacgccagcCCCACAGCACCGGCGGGATCAGCGCAAAGAACGGCGGGATCAGCCACAGCGGc                                                                        >  1:123223‑M2/30‑82 (MQ=255)                                                           tACGTCATCACGCGGTAGGTCAGCCAGCCCCACAGCACCGGCGGGATCAGCGCAAAGAACGGCGGGATCAGCCACAGCGGCAGCGTCACCAGCACCAGCACCAGGAACACCAGCGTGCTGCCCAGCGACTGGaa                   <  2:548305/134‑1 (MQ=255)                                                                            ggTCAGCCAGCCCCACAGCACCGGCGGGATCAGCGCAAAGAACGGCGGGATCAGCCACAGCGGCAGCGTCACCAGCACCAGCACCAGGAACACCAGCGTGCTGCCCAGCGACTGGAACACGCTGCCCAGGAtgc  >  1:444897/1‑134 (MQ=255)                                                                              tCAGCCAGCCCCACAGCACCGGCGGGATCAGCGCAAAGAACGGCGGGATCAGCCACAGCGGCAGCGTcaccagcaccagcacca                                                  >  1:920118/1‑84 (MQ=255)                                                                              tCAGCCAGCCCCACAGCACCGGCGGGATCAGCGCAAAGAACGGCGGGATCAGCCACAGCGGCAGCGTcaccagcaccagcacca                                                  <  2:920118/84‑1 (MQ=255)                                                                                       |                                                                                                                           GTGGCGGCGCACCAGCGTCTTGCGCTCTTCGGCGCTGGCGTGGTCGGCGAGCGCGTCATACGTCATCACGCGGTAGGTCAGCCAGCCCCACAGCACCGGCGGGATCAGCGCAAAGAACGGCGGGATCAGCCACAGCGGCAGCGTCACCAGCACCAGCACCAGGAACACCAGCGTGCTGCCCAGCGACTGGAACACGCTGCCCAGGATGC  >  AM260479/2538198‑2538406

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