New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	5926821 =	33 (0.830)	4 (0.110)	4/320	NT	11.2%	coding (1074/1623 nt)	PP_5195	putative Iron ABC transporter, permease protein
	?	AE015451	5926859 =	32 (0.840)					coding (1036/1623 nt)	PP_5195	putative Iron ABC transporter, permease protein

CAAGGCGCTGTTGGCCAGTGGCTGGTGCCTGCTGGTGTTCGCCTGCGCTTTCGTCATCCCGCTGCTGCAGTTGCTGGCGTGGTTCTGGCAACGCGGCCGGCATGACCTGGATGAGCGTTATATCGGCCTGGTGCTGCACACCCTGTACCTGGGCAGCATGGCGGCGCTGGTCACGGTATGCGTGGCGCTGCTGCTGGCCTTTGCCCGGCGCCAGGCACCGACGGGCGGTATCCGCGCCGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/5927061‑5926821 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgccgggCAAAGGCCAGCAGCAGCGCCACGCATACCGTGACCAGCGCCGCCATGCTGCCCAGGTACAGGGTGTGCAGCACCAGGCCGATATAACGCTCATCCAGGTCATGCCGGCCGCGTTGCCAGAACCACGCCAGCAACTGCAGCAGCGGGATGACGAAAGCGCAGGCGAACACCAGCAGGCACCAGCCACTGGCCAACAGCGCCTTGAGGCCGCGCAGGTGATACAGCGCCT  >  AE015451/5926859‑5927086                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       CAAGGCGCTGTTGGCCAGTGGCTGGTGCCTGCTGGTGTTCGCCTGCGCTTTCGTCATCCCGCTGCTGCAGTTGCTGGCGTGGTTCTGGCAACGCGGCCGGCATGACCTGGATGAGCGTTATATCGGCCTGGTGCTGCACACCCTGTACCTGGGCAGCATGGCGGCGCTGGTCACGGTATGCGTGGCGCTGCTGCTGGCCTTTGCCCGGCGCCAGGCACCGACGGGCGGTATCCG                                                                                                                                                                                                                                             <  2:657167/234‑1     GCGCTGTTGGCCAGTGGCTGGTGCCTGCTGGTGTTCGCCTGCGCTTTCGTCATCCCGCTGCTGCAGTTGCTGGCGTGGTTCTGGCAACGCGGCCGGCATGACCTGGATGAGCGTTATATCGGCCTGGTGCTGCACACCCTGTACCTGGGCAGCATGGCGGCGCTGGTCACGGTATGCGTGGCGCTGCTGCTGGCCTTTGCCCGGCGCCAGGCACCGACGGGCGGTATCCGCGCCG                                                                                                                                                                                                                                        <  1:679964/235‑1                                                                                                                           TCGGCCTGGTGCTGCACACCCTGTACCTGGGCAGCATGGCGGCGCTGGTCACGGTATGCGTGGCGCTGCTGCTGGCCTTTGCCCGGCGCCAGGCACCGACGGGCGGTATCCGCGCCGGGCAAAGGCCAGCAGCAGCGCCACGC                                                                                                                                                                                                              <  1:601077/143‑1                                                                                                                                                               TGGCGGCGCTGGTCACGGTATGCGTGGCGCTGCTGCTGGCCTTTGCCCGGCGCCAGGCACCGACGGGCGGTATCCGCGCCGGGCAAAGGCCAGCAGCAGCGCCACGCATACC                                                                                                                                                                                                         >  1:648946/1‑112                                                                                                                                                                GGCGGCGCTGGTCACGGTATGCGTGGCGCTGCTGCTGGCCTTTGCCCGGCGCCAGGCACCGACGGGCGGTATCCGCGCCGGGCAAAGGCCAGCAGCAG                                                                                                                                                                                                                      <  2:648945/98‑1                                                                                                                                                                       CTGGTCACGGTATGCGTGGCGCTGCTGCTGGCCTTTGCCCGGCGCCAGGCACCGACGGGCGGTATCCGCGCCGGGCAAAGGCCAGCAGCAGCGCCACGCA                                                                                                                                                                                                             <  1:409646/100‑1                                                                                                                                                                                   TGCGTGGCGCTGCTGCTGGCCTTTGCCCGGCGCCAGGCACCGACGGGCGGTATCCG                                                                                                                                                                                                                                             <  1:675409/56‑1                                                                                                                                                                                   TGCGTGGCGCTGCTGCTGGCCTTTGCCCGGCGCCAGGCACCGACGGGCGGTATCCG                                                                                                                                                                                                                                             >  2:675408/1‑56                                                                                                                                                                                                                                          GCGCCGGGCAAAGGCCAGCAGCAGCGCCACGCATACCGTGACCAGCGCCGCCATGCTGCCCAGGTACAGGGTGTGCAGCACCAGGCCGATATAACGCTCATCCAGGTCATGCCGGCCGCGTTGCCAGAACCACGCCAGCAACTGCAGCAGCGGGATGACGAAAGCGCAGGCGAACACCAGCAGGCACCAGCCACTGGCCAACAGCGCCTTGAGGCCGCGCAGGTGATACAGCGCCT  <  1:457563/236‑1                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       CAAGGCGCTGTTGGCCAGTGGCTGGTGCCTGCTGGTGTTCGCCTGCGCTTTCGTCATCCCGCTGCTGCAGTTGCTGGCGTGGTTCTGGCAACGCGGCCGGCATGACCTGGATGAGCGTTATATCGGCCTGGTGCTGCACACCCTGTACCTGGGCAGCATGGCGGCGCTGGTCACGGTATGCGTGGCGCTGCTGCTGGCCTTTGCCCGGCGCCAGGCACCGACGGGCGGTATCCGCGCCGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/5927061‑5926821 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgccgggCAAAGGCCAGCAGCAGCGCCACGCATACCGTGACCAGCGCCGCCATGCTGCCCAGGTACAGGGTGTGCAGCACCAGGCCGATATAACGCTCATCCAGGTCATGCCGGCCGCGTTGCCAGAACCACGCCAGCAACTGCAGCAGCGGGATGACGAAAGCGCAGGCGAACACCAGCAGGCACCAGCCACTGGCCAACAGCGCCTTGAGGCCGCGCAGGTGATACAGCGCCT  >  AE015451/5926859‑5927086

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