New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4513364	91 (1.000)	3 (0.030)	3/262	NT	3.2%	coding (46/999 nt)	fecC	ferric citrate ABC transporter permease
	?	NC_000913	= 4513381	97 (1.130)					coding (29/999 nt)	fecC	ferric citrate ABC transporter permease
Rejected: Frequency below/above cutoff threshold.

CAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4513215‑4513364 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtgctgcTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCCTGCTGCCTGGACACACGCCAACGCTACCAGAAGCGCTGGTGCAAAACC  <  NC_000913/4513381‑4513235                                                                                                                                                                                                                                                                                                           CAGGCTGCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTG                                                                                                                                                      <  1:1489905/149‑1  AGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGC                                                                                                                                                     >  2:147387/1‑149                                     TTCTGGTAGCGTTGGAGTGTGTCCAGTAACCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATT                                                                                                                  <  2:1146183/149‑1                                        TGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAAAAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTC                                                                                                               <  2:88449/149‑1                                           TAGCGTTGGCGTGTGTCACGGCAGCAGGGCGCGGGTGGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTGAGCCAGAAAATAATAAAATGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGG                                                                                                            <  2:1248937/149‑1                                                                                                                                                GTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCCTGCTGCCTGGACACACGCCAACGCTACCAGAAGCGCTGGTGCA       >  2:1281015/1‑149                                                                                                                                                    TGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCCTGCTGCCTGGACACACGCCAACGCTACCAGAAGCGCTGGTGCAAAAC   >  2:119311/1‑149                                                                                                                                                     GCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCCTGCTGCCTGGACACACGCCAACGCTACCAGAAGCGCTGGTGCAAAACC  >  1:644119/1‑149                                                                                                                                                     GCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCCTGCTGCCTGGACACACGCCAACGCTACCAGAAGCGCTGGTGCAAAACC  <  2:471766/149‑1                                                                                                                                                                                                                                                                                                           CAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4513215‑4513364 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtgctgcTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCCTGCTGCCTGGACACACGCCAACGCTACCAGAAGCGCTGGTGCAAAACC  <  NC_000913/4513381‑4513235

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