New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 4781905	54 (0.940)	5 (0.090)	5/324	NT	9.3%	coding (5409/7887 nt)	pvdJ	non‑ribosomal peptide synthetase (subunit of ferribactin synthase)
	?	AE015451	= 4781937	46 (0.840)					coding (5377/7887 nt)	pvdJ	non‑ribosomal peptide synthetase (subunit of ferribactin synthase)

GCTCGTGCTCGCGCAGGTTGAGGTTGTGGTCAAGCAGTGACTTGAGCCACGCCTTGGCCGTGACGCCGGCTTCGGCGCCCGGCAACCGGACCCGCAGCGGGATGCTGTTGATGAACAGCCCGACGGTAGCCTGCATCTCGGGGAATTCCACCGGCCGACCGGCCACGGTCACCCCGAACACCACGTCACGGTCGCCGCTGTAGCGTTGCAGTAACAGTGCCCAGGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/4781679‑4781905 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcccaggcgGCACAGTTGACGGTCAACACCCTGGTGCAAGGCGCCTGGGCACTGTTACTGCAACGCTACAGCGGCGACCGTGACGTGGTGTTCGGGGTGACCGTGGCCGGTCGGCCGGTGGAATTCCCCGAGATGCAGGCTACCGTCGGGCTGTTCATCAACAGCATCCCGCTGCGGGTCCGGTTGCCGGGCGCCGAAGCCGGCGTCACGGCCAAGGCGTGGCTCAAGTCACTGCTTGA  <  AE015451/4781937‑4781708                                                                                                                                                                                                                                                                                                                                                                                                                                                                           GCTCGGGCTCGCGCGGGTTGAGGTGGTGGTCAAGCAGTGACTTGAGCCACGCCTTGGCCGTGACGCCGGCTTCGGCGCCCGGCAACCGGACCCGCAGCGGGATGCTGTTGATGAACAGCCCGACGGTAGCCTGCATCTCGGGGAATTCCACCGGCCGACCGGCCACGGTCACCCCGAACACCACGTCACGGTCGCCGCTGTAGCGTTGCAGTAACAGTGCCCAGGC                                                                                                                                                                                                                                         <  2:59225/226‑1                                                                                                                                                                           ACCCCGAACACCACGTCACGGTCGCCGCTGTGGCGATGCCGTAGCAGTGCCCAGGCGGCACAGCTGACGGTCGACACCCc                                                                                                                                                                                                                 >  1:662873/1‑79                                                                                                                                                                                      CACGTCACGGTCGCCGCTGTAGCGTTGCAGTAACAGTGCCCAGGCGGCACAGTTGACGGTCAACACCCTGGTGCAAGGCGCCTGGGCACTGTTACTGCAACGCTACAGCGG                                                                                                                                                                       >  2:754194/1‑111                                                                                                                                                                                              GGTCGCCGCTGTAGCGTTGCAGTAACAGGGCCCAGGCGGCACAGTTGACGGTCAACACCCTGGTGCAAGGCGCCTG                                                                                                                                                                                                  >  2:218422/1‑76                                                                                                                                                                                                 CGCCGCTGTAGCGTTGCAGTAACAGTGCCCAGGCGGCACAGTTGACGGTCAACACCCTGGTGCAAGGCGCCTGGGCACTGTTACTGCAACGCTACAGCGGCGACCGTGACGTGGTGTTCGGGGTGACCGTGGCCGGTCGGCCGGTGGAATTCCCCGAGATGCCGGCTA                                                                                                   >  1:515716/1‑168                                                                                                                                                                                                                   GTAACAGTGCCCAGGCGGCACAGTTGACGGTCAACACCCTGGTGCAAGGCGCCTGGGCACTGTTACTGCAACGCTACAGCGGCGACCGTGACGTGGTGTTCGGGGTGACCGTGGCCGGTCGGCCGGTGGAATTCCCCGAGATGCAGGCTACCGTCGGGCTGTTCATCAACAGCATC                                                                         >  1:335067/1‑176                                                                                                                                                                                                                            CCCAGGCGGCACAGTTGACGGTCAACACCCTGGTGCAAGGCGCCTGGGCACTGTTACTGCAACGCTACAGCGGCGACCGTGACGTGGTGTTCGGGGTGACCGTGGCCGGTCGGCCGGTGGAATTCCCCGAGATGCAGGCTACCGTCGGGCTGTTCATCAACAGCATCCCGCTGCGGGTCCGGTTGCCGGGCGCCGAAGCCGGCGTCACGGCCAAGGCGTGGCTCAAGTCACTGCT     >  2:646188/1‑235                                                                                                                                                                                                                             CCAGGCGGCACAGTTGACGGTCAACACCCTGGTGCAAGGCGCCTGGGCACTGTTACTGCAACGCTACAGCGGCGACCGTGACGTGGTGTTCGGGGTGACCGTGGCCGGTCGGCCGGTGGAATTCCCCGAGATGCAGGCTA                                                                                                   <  2:863100/140‑1                                                                                                                                                                                                                              CAGGCGGCACAGTTGACGGTCAACACCCTGGTGCAAGGCGCCTGGGCACTGTTACTGCAACGCTACAGCGGCGACCGTGACGTGGTGTTCGGGGTGACCGTGGCCGGTCGGCCGGTGGAATTCCCCGAGATGCAGGCTACCGTCGGGCTGTTCATCAACAGCATCCCGCTGCGGGTCCGGTTGCCGGGCGCCGAAGCCGGCGTCACGGCCAAGGCGTGGCTCAAGTCACTGCTTGA  <  1:331933/236‑1                                                                                                                                                                                                                                                                                                                                                                                                                                                                           GCTCGTGCTCGCGCAGGTTGAGGTTGTGGTCAAGCAGTGACTTGAGCCACGCCTTGGCCGTGACGCCGGCTTCGGCGCCCGGCAACCGGACCCGCAGCGGGATGCTGTTGATGAACAGCCCGACGGTAGCCTGCATCTCGGGGAATTCCACCGGCCGACCGGCCACGGTCACCCCGAACACCACGTCACGGTCGCCGCTGTAGCGTTGCAGTAACAGTGCCCAGGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/4781679‑4781905 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcccaggcgGCACAGTTGACGGTCAACACCCTGGTGCAAGGCGCCTGGGCACTGTTACTGCAACGCTACAGCGGCGACCGTGACGTGGTGTTCGGGGTGACCGTGGCCGGTCGGCCGGTGGAATTCCCCGAGATGCAGGCTACCGTCGGGCTGTTCATCAACAGCATCCCGCTGCGGGTCCGGTTGCCGGGCGCCGAAGCCGGCGTCACGGCCAAGGCGTGGCTCAAGTCACTGCTTGA  <  AE015451/4781937‑4781708

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