New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 1901548	87 (1.270)	3 (0.050)	3/320	NT	3.5%	coding (2099/4026 nt)	PP_1703	assimilatory nitrate reductase/sulfite reductase
	?	AE015451	= 1901589	84 (1.290)					coding (2140/4026 nt)	PP_1703	assimilatory nitrate reductase/sulfite reductase
Rejected: Frequency below/above cutoff threshold.

CGGCGCTGAACAAACTGGAGCCCGGCCCCTACCTTGAGATTCACCCCGAGGATGCCGCCCGGCTGGGTGTCGCTGAAAAGGACCAGGTGGCCGTTCGCTCGCGCCGGGGGCAGGCCGTGTTGCCGGCGCGCATCAGTGACCGGGTCATGCCGGGCAACTGCTTTGCGCCGTTTCACTGGAACGACCTGATCGGCGAGCAATTGGCGATCAATGCCGCGACCTGCGATGCGGTCGACCCGCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/1901308‑1901548 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgacccgctCCAGGGCTACGGCGCAGTGCTTGAAGGCAGGCTGCAGCGACAGCGGGTCGACCGCATCGCAGGTCGCGGCATTGATCGCCAATTGCTCGCCGATCAGGTCGTTCCAGTGAAACGGCGCAAAGCAGTTGCCCGGCATGACCCGGTCACTGATGCGCGCCGGCAACACGGCCTGCCCCCGGCGCGAGCGAACGGCCACCTGGTCCTTTTCAGCGACACCCAGCCGGGCGGCATCCTC  <  AE015451/1901589‑1901355                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              CGGCGCTGAACAAACTGGAGCCCGGCCCCTACCTTGAGATTCACCCCGAGGATGCCGCCCGGCTGGGTGTCGCTGAAAAGGACCAGGTGGCCGTTCGCTCGCGCCGGGGGCAGGCCGTGTTGCCGGCGCGCATCAGTGACCGGGTCATGCCGGGCAACTGCTTTGCGCCGTTTCACTGGAACGACCTGATCGGCGAGCAATTGGCGATCAATGCCGCGACCTGCGATGCGGTC                                                                                                                                                                                                                                                     >  2:410368/1‑233           CAAACTGGAGCCCGGCCCCTACCTTGAGATTCACCCCGAGGATGCCGCCCGGCTGGGTGTCGCTGAAAAGGACCAGGTGGCCGTTCGCTCGCGCCGGGGGCAGGCCGTGTTGCCGGCGCGCATCAGTGACCGGGTCATGCCGGGCAACTGCTTTGCGCCGTTTCACTGGAACGACCTGATCGGCGAGCAATTGGCGATCAATGCCGCGACCTGCGACGCGGTCGACCCGCT                                                                                                                                                                                                                                             >  1:118566/1‑231                                           ACCCCGAGGATGCCGCCCGGCTGGGTGTCGCTGAAAAGGACCAGGTGGCCGTTCGCTCGCGCCGGGGGCAGGCCGTGTTGCCGGCGCGCATCAGTGACCGGGTCATGCCGGGCAACTGCTTTGCGCCGTTTCACTGGAACGACCTGATCGGCGAGCAATTGGCGATCAATGCCGCGACCTGCGATGCGGTCGACCCG                                                                                                                                                                                                                                               <  1:1080622/197‑1                                                                               gGGACCAGGTGGCCGTTCGCTCGCGCCGGGGGCAGGCCGTGTTGCCGGCGCGCATCAGTGACCGGGTCATGCCGGGCAACTGCTTTGCGCCGTTTCACTGGAACGACCTGATCGGCGAGCAATTGGCGATCAATGCCGCGACCTGCGATGCGGTCGACCCGCTCCAGGGCTACGGCGCAGTGCTTGAAGGCAGGCTGCAGCGACAGCGGGTCGACCGCATCGCAGGTCGCGGCAT                                                                                                                                                                     <  1:548020/234‑1                                                                                              TTCGCTCGCGCCGGGGGCAGGCCGTGTTGCCGGCGCGCATCAGTGACCGGGTCATGCCGGGCAACTGCTTTGCGCCGTTTCACTGGAACGACCTGATCGGCGAGCAATTGGCGATCAATGCCGCGACCTGCGATGCGGTCGACCCGCTCCAGGGCTACGGCGCAGTGCTTGAAGGCAGGCTGCAGCGACAGCGGGTCGACCGCATCGCAGGTCGCGGCA                                                                                                                                                                      >  2:548017/1‑219                                                                                                            GGGCAGGCCGTGTTGCCGGCGCGCATCAGTGACCGGGTCATGCCGGGCAACTGCTTTGCGCCGTTTCACTGGAACGACCTGATCGGCGAGCAATTGGCGATCAATGCCGCGACCTGCGATGCGGTCG                                                                                                                                                                                                                                                    <  1:638128/127‑1                                                                                                            GGGCAGGCCGTGTTGCCGGCGCGCATCAGTGACCGGGTCATGCCGGGCAACTGCTTTGCGCCGTTTCACTGGAACGACCTGATCGGCGAGCAATTGGCGATCAATGCCGCGACCTGCGATGCGGTCG                                                                                                                                                                                                                                                    >  2:638125/1‑127                                                                                                                                                                                                     GCAATTGGCGATCAATGCCGCGACCTGCGATGCGGTCGACCCGCTCCAGGGCTACGGCGCAGTGCTTGAAGGCAGGCTGCAGCGACAGCGGGTCGACCGCATCGCAGGT                                                                                                                                                                             >  1:1046688/1‑109                                                                                                                                                                                                                                          GACCCGCTCCAGGGCTACGGCGCAGTGCTTGAAGGCAGGCTGCAGCGACAGCGGGTCGACCGCATCGCAGGTCGCGGCATTGATCGCCAATTGCTCGCCGATCAGGTCGTTCCAGTGAAACG                                                                                                                           <  2:1104368/122‑1                                                                                                                                                                                                                                              CGCTCCGGGGCTACGGCGCAGTGCTTGAAGGCAGGCTGCAGCGACAGCGGGTCGACCGCATCGCAGGTGGCGGCATTGATCGCCAATTGCTCGCCGATCAGGTCG                                                                                                                                        <  2:139351/105‑1                                                                                                                                                                                                                                                 TCCAGGGCTACGGCGCAGTGCTTGAAGGCAGGCTGCAGCGACAGCGGGTCGACCGCATCGCAGGTCGCGGCATTGATCGCCAATTGCTCGCCGATCAGGTCGTTCCAGTGAAACGGCGCAAAGCAGTTGCCCGGCATGACCCGGTCACTGATGCGCGCCGGCAACACGGCCTGCCCCCGGCGCGAGCGAACGGCCACCTGGTCCTTTTCAGCGACACCCAGCCGGGCGGCATCCTC  >  1:1059130/1‑236                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              CGGCGCTGAACAAACTGGAGCCCGGCCCCTACCTTGAGATTCACCCCGAGGATGCCGCCCGGCTGGGTGTCGCTGAAAAGGACCAGGTGGCCGTTCGCTCGCGCCGGGGGCAGGCCGTGTTGCCGGCGCGCATCAGTGACCGGGTCATGCCGGGCAACTGCTTTGCGCCGTTTCACTGGAACGACCTGATCGGCGAGCAATTGGCGATCAATGCCGCGACCTGCGATGCGGTCGACCCGCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/1901308‑1901548 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgacccgctCCAGGGCTACGGCGCAGTGCTTGAAGGCAGGCTGCAGCGACAGCGGGTCGACCGCATCGCAGGTCGCGGCATTGATCGCCAATTGCTCGCCGATCAGGTCGTTCCAGTGAAACGGCGCAAAGCAGTTGCCCGGCATGACCCGGTCACTGATGCGCGCCGGCAACACGGCCTGCCCCCGGCGCGAGCGAACGGCCACCTGGTCCTTTTCAGCGACACCCAGCCGGGCGGCATCCTC  <  AE015451/1901589‑1901355

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