New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 1685981	78 (0.910)	4 (0.050)	4/342	NT	5.1%	coding (798/1425 nt)	patD	medium chain aldehyde dehydrogenase
	?	AE015451	= 1685988	75 (0.910)					coding (791/1425 nt)	patD	medium chain aldehyde dehydrogenase

GCCCGCTGCACCAGCCCGGCGACCTTGTCGCGGTGCTGGGCGCTGATCAGTGGCCCCAGTTCGGTTTCTTCAGACTCCTGCAGGCCCGGCTTGAGGCTGGCTACGGCTTTACCCAGGCGTTCGACAAAGCGCTCGTAGATCCCGTCCTGCACATACAGGCGGCAGGCTGCGGTGCAGTCCTGGCCGGCGTTGTAGAAGCCGAAGGTGCGGATGCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/1685766‑1685981 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gatgccgCCATCGACGGCATCCGCACCTTCGGCTTCTACAACGCCGGCCAGGACTGCACCGCAGCCTGCCGCCTGTATGTGCAGGACGGGATCTACGAGCGCTTTGTCGAACGCCTGGGTAAAGCCGTAGCCAGCCTCAAGCCGGGCCTGCAGGAGTCTGAAGAAACCGAACTGGGGCCACTGATCAGCGCCCAGCACCGCGACAAGGTCGCCGGGCTGGTGCAGCGGGCCATCGCCCAGCC  <  AE015451/1685988‑1685754                                                                                                                                                                                                                                                                                                                                                                                                                                                                     GCCCGCTGCACCAGCCCGGCGACCTTGTCGCGGTGCTGGGCGCTGCGCAGTGGCCCCAGTTCGGTTTCTTCAGACTCCTGCAGGCCCGGCTTGAGGCTGGCTACGGCTTTACCCAGGCGTTCGACAAAGCGCTCGTAGATCCCGTCCTGCACATACAGGCGGCAGGCTGCGGTGCAGTCCTGGCCGGCGTTGTAGAAGCCGAAGGTGCGG                                                                                                                                                                                                                                                   <  2:282647/210‑1     GCTGCACCAGCCCGGCGACCTTGTCGCGGTGCTGGGCGCTGATCAGTGGCCCCAGTTCGGTTTCTTCAGACTCCTGCAGGCCCGGCTTGAGGCTGGCTACGGCTTTACCCAGGCGTTCGACAAAGCGCTCGTAGATCCCGTCCTGCACATACAGGCGGCAGGCTGCGGTGCAGTCCTGGCCGGCGTTGTAGAAGCCGAAGGTGCGG                                                                                                                                                                                                                                                   <  1:1187424/206‑1                      cCCTTGTCGCGGTGCTGGGCGCTGATCAGTGGCCCCAGTTCGGTTTCTTCAGACTCCTGCAGGCCCGGCTTGAGGCTGGCTACGGCTTTACCCAGGCGTTCGACAAAGCGCTCGTAGATCCCGTCCTGCACATACAGGCGGCAGGCTGCGGTGCAGTCCTGGCCGGCGTTGTAGAAGCCGAAGGTGCGGATG                                                                                                                                                                                                                                                <  1:392529/191‑1                                               TCAGTGGCCCCAGTTCGGTTTCTTCAGACTCCTGCAGGCCCGGCTTGAGGCTGGCTACGGCTTTACCCAGGCGTTCGACAAAGCGCTCGTAGATCCCGTCCTGCACATACAGGCGGCAGGCTGCGGTGCAGTCCTGGCCGGCGTTGTAGAAGCCGAAGGTGCGGATGCCG                                                                                                                                                                                                                                             >  2:1136691/1‑170                                                                                  AGGCCCGGCTTGAGGCTGGCTACGGCTTTACCCAGGCGTTCGACAAAGCGCTCGTAGATCCCGTCCTGCACATACAGGCGGCAGGCTGCGGTGCAGTCCTGGCCGGCGTTGTAGAAGCCGAAGGTGCGGAc                                                                                                                                                                                                                                                 >  2:165803/1‑130                                                                                                                    GGCGTTCGACAAAGCGCGCGTAGATCCCGTCCTGCACATACAGGCGGCAGGCTGCGGTGCAGTCCTGGCCGGCGTTGTAGAAGCCGAAGGTGCGG                                                                                                                                                                                                                                                   <  2:313829/95‑1                                                                                                                                           tTCCCGTCCTGCACATCCAGGCGGCAGGCTGCGGTGCAGTCCTGGCCGGCGTTGTAGAAGCCGAAGGTGCGGATGCCGCCAGCGACGGCATCCGCACCTTCGGCTTCTACAACGCCGGCCAGGACTGCACCGCAGCCTGCCGCC                                                                                                                                                                           <  1:506641/143‑1                                                                                                                                                         ATACAGGCGGCAGGCTGCGGTGCAGTCCTGGCCGGCGTTGTAGAAGCCGAAGGTGCGGATGCCGCCATCGACGGCATCCGCACCTTCGGCTTCTACA                                                                                                                                                                                                            <  2:258226/97‑1                                                                                                                                                                            GTGCAGTCCTGGCCGGCGTTGTAGAAGCCGAAGGTGCGGATGCCGCCATCGACGGCATCCGCACCTTCGGCTTCTACAACGCC                                                                                                                                                                                                       <  1:1136694/83‑1                                                                                                                                                                                       GCCGGCGTTGTAGAAGCCGAAGGTGCGGATGCCGCCATCGACGGCATCCGCACCTTCGGCTTCTACAACGCCGGCCAGGACTGCACCGCAGCCTGCCGCCTGTATGTGCAGGACGGGATCTACGAGCGCTTTGTCGAACGCCTGGG                                                                                                                             >  2:955669/1‑146                                                                                                                                                                                                                      CCGCCATCGACGGCATCCGCACCTTCGGCTTCTACAACGCCGGCCAGGACTGCACCGCAGCCTGCCGCCTGTATGTGCAGGACGGGATCTACGAGCGCTTTGTCGAACGCCTGGGTAAAGCCGTAGCCAGCCTCAAGCCGGGCCTGCAGGAGTCTGAAGAAACCGAACTGGGGCCACTGATCAGCGCCCAGCACCGCGACAAGGTCGCCGGGCTGGgg                      >  1:100010/1‑216                                                                                                                                                                                                                        GCCATCGACGGCATCCGCACCTTCGGCTTCTACAACGCCGGCCAGGACTGCACCGCAGCCTGCCGCCTGTATGTGCAGGACGGGATCTACGAGCGCTTTGTCGAACGCCTGGGTAAAGCCGTAGCCAGCCTCAAGCCGGGCCTGCAGGAGTCTGAAGAAACCGAACTGGGGCCACTGATCAGCGCCCAGCACCGCGACAAGGTCGCCGGGCTGGTGCAGCGGGCCATCGCCCAGCC  <  1:572302/236‑1                                                                                                                                                                                                                                                                                                                                                                                                                                                                     GCCCGCTGCACCAGCCCGGCGACCTTGTCGCGGTGCTGGGCGCTGATCAGTGGCCCCAGTTCGGTTTCTTCAGACTCCTGCAGGCCCGGCTTGAGGCTGGCTACGGCTTTACCCAGGCGTTCGACAAAGCGCTCGTAGATCCCGTCCTGCACATACAGGCGGCAGGCTGCGGTGCAGTCCTGGCCGGCGTTGTAGAAGCCGAAGGTGCGGATGCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/1685766‑1685981 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gatgccgCCATCGACGGCATCCGCACCTTCGGCTTCTACAACGCCGGCCAGGACTGCACCGCAGCCTGCCGCCTGTATGTGCAGGACGGGATCTACGAGCGCTTTGTCGAACGCCTGGGTAAAGCCGTAGCCAGCCTCAAGCCGGGCCTGCAGGAGTCTGAAGAAACCGAACTGGGGCCACTGATCAGCGCCCAGCACCGCGACAAGGTCGCCGGGCTGGTGCAGCGGGCCATCGCCCAGCC  <  AE015451/1685988‑1685754

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