New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	5822892 =	69 (1.010)	4 (0.060)	4/314	NT	6.2%	coding (421/597 nt)	rdgB	dITP/XTP pyrophosphatase
	?	AE015451	5823034 =	57 (0.890)					coding (563/597 nt)	rdgB	dITP/XTP pyrophosphatase

TGTACAGCGCCAGCGGCGGCAGGCTGCTAAAGCCTGCCGCGCCGGGGGGCGACAGCATTTCGATCATGCCAGGCCCAGACGTTTGCGCAGCAGGGCCATGGCACGGGCGCGGTGGCTGAGCTGGTTCTTGTCTGCTGGGGCTAGATCGGCGCTGGAGCAGTTGCGCTCCGGTACCCAGAACAGCGGATCGTAGCCAAAGCCCTGGTCGCCGCTGGCCTCGAACAGGATGCGCCCGTGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/5823131‑5822892 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgcccgtgccaTGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTA  >  AE015451/5823034‑5823216                                                                                                                                                                                                                                                                                                                                                                                                                                         TGTACAGCGCCAGCGGCGGCAGGCTGCTAAAGCCTGCCGCGCCGGGGGGCGACAGCATTTCGATCATGCCAGGCCCAGACGTTTGCGCAGCAGGGCCATGGCACGGGCGCGGTGGCTGAGCTGGTTCTTGTCTGCTGGGGCTAGATCGGCGCTGGAGCAGTTGCGCTCCGGTACCCAGAACAGCGGATCGTAGCCAAAGCCCTGGTCACCGCTGGCCTCGAACAGGATGCGCC                                                                                                                                                                                                >  2:676000/1‑233                                cCCTGCCGCGCCGGGGGGCGACAGCATTTCGATCATGCCAGGCCCAGACGTTTGCGCAGCAGGGCCATGGCACGGGCGCGGTGGCTGAGCTGGTTCTTGTCTGCTGGGGCTAGATCGGCGCTGGAGCAGTTGCGCTCCGGTACCCAGAACAGCGGATCGTAGCCAAAGCCCTGGTCGCCGCTGGCCTCGAACAGGATGCGC                                                                                                                                                                                                 <  1:115803/200‑1                                                                         GCCCAGACGTTTGCGCAGCAGGGCCATGGCACGGGCGCGGTGGCTGAGCTGGTTCTTGTCTGCTGGGGCTAGATCGGCGCTGGAGCAGTTGCGCTCCGGTACCCAGAACAGCGGATCGTAGCCAAAGCCCTGGTCGCCGCTGGCCTCGAACAGGATGCGCCCGTG                                                                                                                                                                                            <  2:302874/165‑1                                                                                                    GGCACGGGCGCGGTGGCTGAGCTGGTTCTTGTCTGCTGGGGCTAGATCGGCGCTGGAGCAGTTGCGCTCCGGTACCCAGAACAGCGGATCGTAGCCAAAGCCCTGGTCGCCGCTGGCCTCGAACAGGATGCG                                                                                                                                                                                                  <  1:873141/132‑1                                                                                                                                             CTAGATCCGCGCTGGAGCAGTTGCGCTCCGGTACCCAGAACAGCGGATCGTAGCCAAAGCCCTGGTCGCCGCTGGCCTCGAACAGGACGCGCCC                                                                                                                                                                                               >  2:587902/1‑94                                                                                                                                                                         CGGTACCCAGAACAGCGGATCGTAGCCAAAGCCCTGGTCGCCGCTGGCCTCGAACAGGATGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTG                                                                                                                                                                  <  2:816232/95‑1                                                                                                                                                                         CGGTACCCAGAACAGCGGATCGTAGCCAAAGCCCTGGTCGCCGCTGGCCTCGAACAGGATG                                                                                                                                                                                                    <  2:664316/61‑1                                                                                                                                                                            TACCCAGAACAGCGGATCGTAGCCAAAGCCCTGGTCGCCGCTGGCCTCGAACAGGATGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCC                                                                                                                                                              <  2:1184315/96‑1                                                                                                                                                                                 AGAACAGCGGATCGTAGCCAAAGCCCTGGTCGCCGCTGGCCTCGAACAGGATGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTT                                            <  1:599540/205‑1                                                                                                                                                                                 AGAACAGCGGATCGTAGCCAAAGCCCTGGTCGCCGCTGGCCTCGAACAGGATGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTT                                            >  2:599537/1‑205                                                                                                                                                                                                                                       GCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCCGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTTGTG                                                                       >  2:316820/1‑124                                                                                                                                                                                                                                       GCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCACATGCCCTTATTGCGACTTCAACTCCAACGCAGCCGGGCCTGAACTACCGGAAGCTGCCTA  >  1:49748/1‑193                                                                                                                                                                                                                                           GTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGT                                                                    <  2:1228277/123‑1                                                                                                                                                                                                                                                                                                                                                                                                                                         TGTACAGCGCCAGCGGCGGCAGGCTGCTAAAGCCTGCCGCGCCGGGGGGCGACAGCATTTCGATCATGCCAGGCCCAGACGTTTGCGCAGCAGGGCCATGGCACGGGCGCGGTGGCTGAGCTGGTTCTTGTCTGCTGGGGCTAGATCGGCGCTGGAGCAGTTGCGCTCCGGTACCCAGAACAGCGGATCGTAGCCAAAGCCCTGGTCGCCGCTGGCCTCGAACAGGATGCGCCCGTGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/5823131‑5822892 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgcccgtgccaTGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTA  >  AE015451/5823034‑5823216

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