New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	3953277 =	62 (1.000)	4 (0.070)	4/330	NT	6.2%	coding (1246/1977 nt)	PP_5606	conserved exported protein of unknown function with YVTN beta‑propeller repeat
	?	AE015451	3953322 =	62 (1.040)					coding (1201/1977 nt)	PP_5606	conserved exported protein of unknown function with YVTN beta‑propeller repeat

TCCAGGCTGAACAGGGCCTGGGGCCGATGCGCTTCAGCCATGATGGCCGCTACGGCATCGTGCTCAATACCCTTGAGAACCAGGCACTGGTCATCGACGCCAGTACCGACAAGCTGATCCACCACATCCCGGTGGCGGCCGAACCGTACCAGCTGACGTTCACCAAAGGCTATGCCTATGTGCGTGGCCTGGCCTCGCCGAAAGTGAGCATGATCAATCTGGCCAGCCTGGGCGAGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/3953513‑3953277 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgaggCCAGGCCACGCACATAGGCATAGCCTTTGGTGAACGTCAGCTGGTACGGTTCGGCCGCCACCGGGATGTGGTGGATCAGCTTGTCGGTACTGGCGTCGA  >  AE015451/3953322‑3953420                                                                                                                                                                                                                                                                                                                                                  TCCAGGCTGAACAGGGCCTGGGGCCGATGCGCTTCAGCCATGATGGCCGCTACGGCATCGTGCTCAATACCCTTGAGAACCAGGCACTGGTCATCGACGCCAGTACCGACAAGCTGATCCACCACATCCCGGTGGCGGCCGAACCGTACCAGCTGACGTTCACCAAAGGCTATGCCTATGTGCGTGGCCTGGCCTCGCCGAAAGTGAGCATGATCAATCTGGCCAGCCTGGGCGAG                                                                                                      <  1:523267/236‑1      GCTGAACAGGGCCTGGGGCCGATGCGCTTCAGCCATGATGGCCGCTACGGCATCGTGCTCAATACCCTTGAGAACCAGGCACTGGTCATCGACGCCAGTACCGACAAGCTGATCCACCACATCCCGGTGGCGGCCGAACCGTACCAGCTGACGTTCACCAAAGGCTATGCCTATGTGCGTGGCCTGGCCTCGCCGAAAGTGAGCATGATCAATCTGGCCAGCCTGGGCGAGGCC                                                                                                   <  2:411868/234‑1       CTGAACAGGGCCTGGGGCCGATGCGCTTCAGCCATGATGGCCGCTACGGCATCGTGCTCAATACCCTTGAGAACCAGGCACTGGTCATCGACGCCAGTACCGACAAGCTGATCCACCACATCCCGGTGGCGGCCGAACCGTACCAGCTGACGTTCACCAAAGGCTATGCCTATGTGCGTGGCCTGGCCTCGCCGAAAGTGAGCATGATCAATCTGGCCAGCCTGGGCGAGGCCAGG                                                                                                <  1:1083532/236‑1                                    AGCCATGATGGCCGCTACGGCATCGTGCTCAATACCCTTGAGAACCAGGCACTGGTCATCGACGCCAGTACCGACAAGCTGATCCACCACATCCCGGTGGCGGCCGAACCGTACCAGCTGACGTTCACCAAAGGCTATGCCTATGTGCGTGGCCTGGCCTCGCCGAAAGTGAGCATGATCAATCTGGCCAGCCTGG                                                                                                           >  2:1057577/1‑196                                                                                                                                 CCGGTGGCGGCCGAACCGTACCAGCTGACGTTCACCAAAGGCTATGCCTATGTGCGTGGCCTGGCCTCGCCGAAAGTGAGCATGATCAATCTGGCCAGCCTGGGCGAGGCCAGGCCACGCACATAGGCATAGCCTTTGGTGA                                                                    <  2:891818/142‑1                                                                                                                                              AACCGTACCAGCTGACGTTCACCAAAGGCTATGCCTATGTGCGTGGCCTGGCCGCGCCGAAAGTGAGCAGGATCAATCTGGCCAGCATGGGCGAGGCCAGGCC                                                                                              <  2:158224/103‑1                                                                                                                                                                                CTATGTGCGTGGCCTGGCCTCGCCGAAAGTGAGCATGATCAATCTGGCCAGCCTGGGCGAGGCCAGGCCACGCACATAGGCATAGCCTTTGGTGAACGTCAGCTGGTACGGTTCGGCCGCC                                          <  2:900280/121‑1                                                                                                                                                                                                GCCTCGCCGAAAGTGAGCATGATCAATCTGGCCAGCCTGGGCGAGGCCAGGCCACGCACATAGGCATAGCCTTTGGTGAACGTCAGCTGGTACGGTTCGGCCGCCACCGGGATGTGGTGGATCAGCTTGTCGGTACTGGCGTCGA  <  1:409168/145‑1                                                                                                                                                                                                                                                                                                                                                  TCCAGGCTGAACAGGGCCTGGGGCCGATGCGCTTCAGCCATGATGGCCGCTACGGCATCGTGCTCAATACCCTTGAGAACCAGGCACTGGTCATCGACGCCAGTACCGACAAGCTGATCCACCACATCCCGGTGGCGGCCGAACCGTACCAGCTGACGTTCACCAAAGGCTATGCCTATGTGCGTGGCCTGGCCTCGCCGAAAGTGAGCATGATCAATCTGGCCAGCCTGGGCGAGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/3953513‑3953277 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgaggCCAGGCCACGCACATAGGCATAGCCTTTGGTGAACGTCAGCTGGTACGGTTCGGCCGCCACCGGGATGTGGTGGATCAGCTTGTCGGTACTGGCGTCGA  >  AE015451/3953322‑3953420

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