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

CCGACAGGATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCTCGGGCCCATCGACTGGGCACTGCCCCCGCTGTGAGAAACCGCGGTGGGGTACGCAATCGCTCAGGGCGTCATGCTCATGCGCGGCTCGCTGCGCTGGCCAATGGCAATCCCCCGCTCATGGTTCTGGCGCACGTACTCCACCAAGGCTTCAAGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCT  >  AE015451/1612652‑1613099                                                                                                                                                                                                                                        |                                                                                                                                                                                                                         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)                                                    tGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCTCGGGCCCATCGACTGGGCACTgccc                                                                                                                                                                                                                      >  2:315983‑M2/1‑181 (MQ=255)                                                           gggggCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCTCGGGCCCATCGACTGGGCACTGCCCCCGCTGTGAGAAACCGCGGTGGGGTACGCAATCGCTCAGGGCGTCATGCTCATGcgcg                                                                                                                                                            >  1:475064/1‑236 (MQ=255)                                                                                                     cAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCTCGGGCCCATCGACTGGGCACTGCCCCCGCTGTGAGAAACCGCGGTGGGGTACGCAATCGCTCAGGGCGtcatgctcat                                                                                                                                                                 >  2:402651/1‑189 (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)                                                                                                                                                                                           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)                                                                                                                                                                                                               gggCTCGGGCCCATCGACTGGGCACTGCCCCCGCTGTGAGAAACCGCGGTGGGGTACGCAATCGCGCAGGGCGTCATGCTCATGCGCGGCTCGCTGCGCTGGCCAATGGCAATCCCCCGCTCATGGTTCTGGCGCACGTACTccac                                                                                                  >  1:323366/1‑146 (MQ=255)                                                                                                                                                                                                                     gggCCCATCGACTGGGCACTgcccccgctgctcgaggaagctgttggcacggctgaaatgcccgcaaccgaggaacccgcgatacgccgacagcggcgaaggattgaccgacttcaaacccagatgcttggtaccgtcgatcagcttctgcttgctctgtgcatgcgccccccatagcaggaacaccacgttcgg                                           >  1:158141‑M2/1‑20 (MQ=255)                                                                                                                                                                                                                         ccATCGACTGGGCACTGCCCCGGCTGGGAGAAACCGCGGTGGGGTACGCAATCGCTCAGGGCGTCATGCTCATGCGCGGCTCGCTGCGCTGGCCAATGGCAATCCCCCGCTCATGGTTCTGGCGCACGTACTCCACCAAGGCTTCAAGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCt  <  1:218192/232‑1 (MQ=255)                                                                                                                                                                                                                               aCTGGGCACTGCCCCCGCTGTGAGAAACCGCCGGGGGGTACGCAATCGCTCCGGGCGccatgctca                                                                                                                                                                  >  2:61468/1‑66 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                         CCGACAGGATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCGGCGTATCGCGGGTTCATCGGTTGCGGGCATTTCAGCCGTGCCAACAGCTTCCTCGAGCAGCGGGGGCTCGGGCCCATCGACTGGGCACTGCCCCCGCTGTGAGAAACCGCGGTGGGGTACGCAATCGCTCAGGGCGTCATGCTCATGCGCGGCTCGCTGCGCTGGCCAATGGCAATCCCCCGCTCATGGTTCTGGCGCACGTACTCCACCAAGGCTTCAAGGTCGGTGATGCCGGTGTCCTGGCGTCGCCTTCCTTCTTTGAAAGTGGCCAAGCCACCGCCACCATCGGCCAGGAAAAGGTTGCT  >  AE015451/1612652‑1613099

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