New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 1612861	11 (0.430)	3 (0.120)	3/298	NT	19.8%	coding (658/693 nt)	ung	Uracil‑DNA glycosylase
	?	AE015451	= 1612883	14 (0.580)					coding (680/693 nt)	ung	Uracil‑DNA glycosylase

CCGACAGGATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCTCGGGCCCATCGACTGGGCACTGCCCCCGCTGTGAGAAACCGCGGTGGGGTACGCAATCGCTCAGGGCGTCATGCTCATGCGCGGCTCGCTGCGCTGGCCAATGGCAATCCCCCGCTCATGGTTCTGGCGCACGTACTCCACCAAGGCTTCAAGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAG  >  AE015451/1612652‑1613093                                                                                                                                                                                                                  |                                                                                                                                                                                                                                         ccGACAGGATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCTCGGGCCCATCGACTGGGCACTgccc                                                                                                                                                                                                                <  1:16257‑M2/236‑5 (MQ=255)                   tggtCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCTCGGGCCCATCGACTGGGCACTGCCCCCGCTGTGAGAAACcg                                                                                                                                                                                                <  1:252450/234‑1 (MQ=255)                          tGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCTCGGGCCCATCGACTGGGCACTGCCCCCGCTGTGAGAAACCGCGGTgggg                                                                                                                                                                                        >  1:207291/1‑235 (MQ=255)                             gCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCTCGGGCCCATCGACTGGGCACTGCCCCCGCTGTGAGAAACCGCGGTGGGGTa                                                                                                                                                                                      <  1:428921/234‑1 (MQ=255)                                    ccGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGGTTTGAAGTCGGTCAATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCCTCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCagtgcccagtcgatgggccc                                                                                                                                                                                                                      >  1:158141‑M1/1‑175 (MQ=255)                                       aaCGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCTCgg                                                                                                                                                                                                                                      >  1:280051/1‑176 (MQ=255)                                                    tGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCTCGGGCCCATCGACTGGGCACTgccc                                                                                                                                                                                                                >  2:315983‑M2/1‑181 (MQ=255)                                                           gggggCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCTCGGGCCCATCGACTGGGCACTGCCCCCGCTGTGAGAAACCGCGGTGGGGTACGCAATCGCTCAGGGCGTCATGCTCATGcgcg                                                                                                                                                      >  1:475064/1‑236 (MQ=255)                                                                                                     cAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCTCGGGCCCATCGACTGGGCACTGCCCCCGCTGTGAGAAACCGCGGTGGGGTACGCAATCGCTCAGGGCGtcatgctcat                                                                                                                                                           >  2:402651/1‑189 (MQ=255)                                                                                                                                      gccgcTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCagtgcccagtcgatgggcccgagcccccgct                                                                                                                                                                                                           <  2:158137‑M1/108‑32 (MQ=255)                                                                                                                                         gcTGTCGGCGTAGCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCTCGGGCCCATCGACTGGGCACTGCCCCCGCTGTGAGAAACCGCGGTGGGGTACGCAATCGCTCAGGGCGTCATGCTCATGCGCGGCTcgctgcgctg                                                                                                                                         <  2:351626/171‑1 (MQ=255)                                                                                                                                                      cgcgGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCTCGGGCCCATCGACTGGGCACTGCCCCCGCTGTGAGAAACCGCGGTGGGGTACGCAATCGCTCAGGGCGTCATGCTCATGCGCGGCTCGCTGCGCTGGCCAATGGCAATCCCCCGCTCATGGTTc                                                                                                             <  2:409060/186‑1 (MQ=255)                                                                                                                                                                   gTTGCGGGCATCTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCTCGGGGCCATCGACTGGGCACTGCCCCCGCTGTGAGAAAc                                                                                                                                                                                                  <  2:50461/88‑1 (MQ=255)                                                                                                                                                                                        ccaccaGCTTCCTCGAGCAGCGGGGGCTCGGGCCCATCGACTGGGCACTGCCCCCGCTGTGAGAAACCGCGGTGGGGTACGCAATCGCTCAGGGCGTCATGCTCATGCGCGGCTCGCTGCGCTGGCCAATGGCAATCCCCCGCTCATGGTTCTGGCGCACGTACTCCACCAAGGCTTCAAGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTg                                  <  1:361054/227‑1 (MQ=255)                                                                                                                                                                                           aCAGCTTCCTCGAGCAGCGGGGGCagtgcccagtcgatgggcccgagcccccgctgctcga                                                                                                                                                                                                     >  2:202228‑M1/1‑24 (MQ=255)                                                                                                                                                                                           aCAGCTTCCTCGAGCAGCGGGGGCTCGGGCCCATCGACTGGGCACTGCCCCCGCTGTGAGAAACCGCGGTGGGGTACGCAATCGCTCAGGGCGTCATGCTCATGCGCGGCTCGCTGCGCTGGCCAATGGCAATCCCCCGCTCATGGTTCTGGCGCACGTACTCCACCAAGGCTTCAAGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAgccaccg                      >  1:126783/1‑236 (MQ=255)                                                                                                                                                                                                    tCGAGCAGCGGGGGCTCGGGCCCATCGACTGGGCACTgcccccgctgctcgaggaagctgt                                                                                                                                                                                            >  2:202228‑M2/1‑37 (MQ=255)                                                                                                                                                                                                          agcGGGGGCTCGGGCCCATCGACTGGGCACTgcccccgctgctcgaggaagctgttggcacggctgaaatgcccgcaaccgatgaacccgcgatacgccgacagcggc                                                                                                                                       <  2:158137‑M2/108‑78 (MQ=255)                                                                                                                                                                                                          agcGGGGGCTCGGGCCCATCGACTGGGCACTGCCCCCGCTGTGAGAAACcgc                                                                                                                                                                                               >  2:12379/1‑52 (MQ=255)                                                                                                                                                                                                               gggCagtgcccagtcgatgggcccgagcccccgctgctcgaggaagctgttggcacggctgaaatgcccgcaaccgatgaacccgcgatacgccgacagcggcgaaggatggaccgacttcaaaaccagatgcttggtaccgtcgatcagcttctgcttgctctgtgcatgcgccccccatagcaggaacaccacgttcgggcattgctcactgaccacctggataatcctgtcgg  <  1:16257‑M1/236‑233 (MQ=255)                                                                                                                                                                                                               gggCagtgcccagtcgatgggcccgagcccccgctgctcgaggaagctgttggcacggctgaaatgcccgcaaccgatgaacccgcgatacgccgacagcggcgaaggatggaccgacttcaaaaccagatgcttggtaccgtcgatcagcttctgcttgctctgtgcatgcgccccccatagca                                                     >  2:315983‑M1/1‑4 (MQ=255)                                                                                                                                                                                                               gggCTCGGGCCCATCGACTGGGCACTGCCCCCGCTGTGAGAAACCGCGGTGGGGTACGCAATCGCGCAGGGCGTCATGCTCATGCGCGGCTCGCTGCGCTGGCCAATGGCAATCCCCCGCTCATGGTTCTGGCGCACGTACTccac                                                                                            >  1:323366/1‑146 (MQ=255)                                                                                                                                                                                                                  |                                                                                                                                                                                                                                         CCGACAGGATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCTCGGGCCCATCGACTGGGCACTGCCCCCGCTGTGAGAAACCGCGGTGGGGTACGCAATCGCTCAGGGCGTCATGCTCATGCGCGGCTCGCTGCGCTGGCCAATGGCAATCCCCCGCTCATGGTTCTGGCGCACGTACTCCACCAAGGCTTCAAGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAG  >  AE015451/1612652‑1613093

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