New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4512063	111 (1.180)	3 (0.030)	3/234	NT	3.1%	coding (352/957 nt)	fecD	ferric citrate ABC transporter permease
	?	NC_000913	= 4512138	87 (1.000)					coding (277/957 nt)	fecD	ferric citrate ABC transporter permease
Rejected: Frequency below/above cutoff threshold.

ACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4511907‑4512063 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccatgccgccAGCCTGGCCTCTGTGGGGGCTCTACTTCTTATGCCGTCACTGCCCGTGATGGTGCTGCCGCTGCTGGCCTTTGCGGGCGGCATGGCGGGGTTGATATTGCTGAAGATGCTGGCAAAGACCCACCAGCCGATGAAGCTGGCGCTCACC  <  NC_000913/4512138‑4511992                                                                                                                                                                                                                                                                                                                  ACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGC                                                                                                                                                              >  8:454894/1‑148      AGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCC                                                                                                                                                        <  5:270288/149‑1                                          AGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGC                                                                                                                                                      >  8:356505/1‑115                                          AGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGC                                                                                                                                                      <  7:356505/115‑1                                                                  CATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCA                                                                                                                                                            >  2:344321/1‑85                                                                  CATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCA                                                                                                                                                            <  1:344321/85‑1                                                                                       GTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCAGCCTGGCCTCTGTGGGGGCTCTACTTCTTATGCCGTCACTGCCCGTGATGGTGCTGCCGCTGCTGGCCTTTGCGGGC                                                                       >  8:245002/1‑149                                                                                             GCCAGCTTCATCGGCGTGTGGGTCTTTGCGAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCAGCCTGGCCTCTGTGGGGGCTCTACTTCTTATGCCGTCACTGCCCTTGATGGTGCTGCCGCTGCTGGCCTTTGCGGGCGGCATG                                                                 <  5:384240/149‑1                                                                                                 GCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCAGCCTGGCCTCTGTGGGGGCTCTACTTCTTATGCCGTCACTGCCCGTGATGGTGCTGCCGCTGCTGGCCTTTGCGGGCGGCATGGCGG                                                             <  7:245002/149‑1                                                                                                                                                     CATGCCGCCAGCCTGGCCTCTGTGGGGGCTCTACTTCTTATGCCGTCACTGCCCGTGATGGTGCTGCCGCTGCTGGCCTTTGCGGGCGGCATGGCGGGGTTGATATTGCTGAAG                                            <  4:27236/114‑1                                                                                                                                                     CATGCCGCCAGCCTGGCCTCTGTGGGGGCTCTACTTCTTATGCCGTCACTGCCCGTGATGGTGCTGCCGCTGCTGGCCTTTGCGGGCGGCATGGCGGGGGTGATATTGCTGAAG                                            >  3:27236/1‑114                                                                                                                                                      ATGCCGCCAGCCTGGCCTCTGTGGGGGCTCTACTTCTTATGCCGTCACTGCCCGTGATGGTGCTGCCGCTGCTGGCCTTTGCGGGCGGCATGGCGGGGTTGATATTGCTGAAGATGCTGGCAAAGACCCACCAGCCGATGAAGCTGGCG        <  5:251019/149‑1                                                                                                                                                        GCCGCCAGCCTGGCCTCTGTGGGGGCTCTACTTCTTATGCCGTCACTGCCCGTGATGGTGCTGCCGCTGCTGGCCTTTGCGGGCGGCATGGCGGGGTTGATATTGCTGAAGATGCTGGCAAAGACCCACCAGCCGATGAAGCTGGCGCT      >  2:459838/1‑149                                                                                                                                                          CGCCAGCCTGGCCTCTGTGGGGGCTCTACTTCTTATGCCGTCACTGCCCGTGATGGTGCTGCCGCTGCTGGCCTTTGCGGGCGGCATGGCGGGGTTGATATTGCTGAAGATGCTGGCAAAGACCCACCAGCCGATGAAGCTGGCGCTCA    >  5:101122/1‑149                                                                                                                                                          CGCCAGCCTGGCCTCTGTGGGGGCTCTACTTCTTATGCCGTCACTGCCCGTGATGGTGCTGCCGCTGCTGGCCTTTGCGGGCGGCATGGCGGGGTTGATATTGCTGAAGATGCTGGCAAAGACCCACCAGCCGATGAAG              <  5:458657/139‑1                                                                                                                                                          CGCCAGCCTGGCCTCTGTGGGGGCTCTACTTCTTATGCCGTCACTGCCCGTGATGGTGCTGCCGCTGCTGGCCTTTGCGGGCGGCATGGCGGGGTTGATATTGCTGAAGATGCTGGCAAAGACCCACCAGCCGATGAAG              >  6:458657/1‑139                                                                                                                                                            CCAGCCTGGCCTCTGTGGGGGCTCTACTTCTTATGCCGTCACTGCCCGTGATGGTGCTGCCGCTGCTGGCCTTTGCGGGCGGCATGGCGGGGTTGATATTGCTGAAGATGCTGGCAAAGACCCACCAGCCGATGAAGCTGGCGCTCACC  <  3:34524/149‑1                                                                                                                                                            CCAGCCTGGCCTCTGTGGGGGCTCTACTTCTTATGCCGTCACTGCCCGTGATGGTGCTGCCGCTGCTGGCCTTTGCGGGCGGCATGGCGGGGTTGATATTGCTGAAGATGCTGGCAAAGACCCACCAGCCGATGAAGCTGGCGCTCACC  >  6:62822/1‑149                                                                                                                                                                                                                                                                                                                  ACAGCAGGGCGTTGTTCACATCCTGCGGGCGCGAGAGCATCAGATAATCCGTCAGGCTGGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4511907‑4512063 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccatgccgccAGCCTGGCCTCTGTGGGGGCTCTACTTCTTATGCCGTCACTGCCCGTGATGGTGCTGCCGCTGCTGGCCTTTGCGGGCGGCATGGCGGGGTTGATATTGCTGAAGATGCTGGCAAAGACCCACCAGCCGATGAAGCTGGCGCTCACC  <  NC_000913/4512138‑4511992

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