New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4512524	136 (0.970)	4 (0.030)	4/164	NT	3.0%	coding (886/999 nt)	fecC	ferric citrate ABC transporter permease
	?	NC_000913	= 4512556	137 (1.060)					coding (854/999 nt)	fecC	ferric citrate ABC transporter permease
Rejected: Frequency below cutoff threshold.

GCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4512419‑4512524 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgccagTGAGCATGCTGCTGGGGGCCACGCTGATGCTGCTGGCAGATGTACTCGCACGCGCGCTGGCCTTCCCCGGAGATCTGCCCGCAGGCGCAGTGCTGGCG  <  NC_000913/4512556‑4512459                                                                                                                                                                                                              GCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATC                                                                                                          <  1:805529/100‑1  CCTCCACACAAGCCAGAAAAAACAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCT                                                                                                         <  1:3612885/100‑1   CTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTG                                                                                                        <  1:1357967/100‑1   CTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTG                                                                                                        >  1:1464378/1‑100   CTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTG                                                                                                        >  1:683007/1‑100    TCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGC                                                                                                       >  1:111740/1‑100    TCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCATGCGAGTACATCTGC                                                                                                       <  1:3018940/100‑1      CTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGACGGCCAGCGCGCGTGCGAGTACATCTGCC                                                                                                      >  2:1810638/1‑99        CACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTCCATCTGCAAG                                                                                                    >  1:480358/1‑99                                         CGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGTGAGCATGCTGCTGGGGGCCACGCTGATGCTGCT                                                                  <  1:1320869/100‑1                                                    CGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGTGAGCATGCTGCTGGGGGCCACGCTGATGCTGCTGGCAGATGTAC                                                       <  1:3280969/100‑1                                                     GCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGTGAGCATGCTGCTGGGGGCCACGCTGATGCTGCTGGCAGATGTACT                                                      <  1:2513423/100‑1                                                       CTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGTGAGCATGCTGCTGGGGGCCACGCTGATGCTGCTGGCAGATGTACTCG                                                    <  1:1694968/100‑1                                                                                                    CTGCCAGTGAGCATGCTGCTGGGGGCCACGCTGATGCTGCTGGCAGATGTACTCGCACGCGCGCTGGCCTTCCCCGGAGATCTGCCCGCAGGCGCAGTGC       >  2:2282748/1‑100                                                                                                    CTGCCAGTGAGCATGCTGCTGGGGGCCACGCTGATGCTGCTGGCAGATGTACTCGCACGCGCGCTGGCCTTCCCCGGAGATCTGCCCGCAGGCGCAGTGC       <  2:3034769/100‑1                                                                                                     TGCCAGTGAGCATGCGGCTGGGGGCCACGCTGATGCTGCTGGCAGATGTACTCGCACGCGCGCTGGCCTTCCCCGGAGATCTGCCCGCAGGCGCAGTGCT      <  1:1639667/100‑1                                                                                                      GCCAGTGAGCATGCTGCTGGGGGCCACGCTGATGCTGCTGGCAGATGTACTCGCACGCGCGCTGGCCTTCCCCGGAGATCTGCCCGCAGGCGCAGTGCTG     <  2:3546843/100‑1                                                                                                      GCCAGTGAGCATGCTGCTGGGGGCCACGCTGATGCTGCTGGCAGATGTACTCGCACGCGCGCTGGCCTTCCCCGGAGATCTGCCCGCAGGCGCAGTGCTG     >  2:3627897/1‑100                                                                                                      GCCAGTGAGCATGCTGCTGGGGGCCACGCTGATGCTGCTGGCAGATGTACTCGCACGCGCGCTGGCCTTCCCCGGAGATCTGCCCGCAGGCGCAGTGCTG     >  1:3462613/1‑100                                                                                                       CCAGTGAGCATGCTGCTGGGGGCCACGCTGATGCTGCTGGCAGATGTACTCGCACGCGCGCTGGCCTTCCCCGGAGATCTGCCCGCAGGCGCAGTGCTGG    >  1:355023/1‑100                                                                                                         AGTGAGCATGCTGCTGGGGGCCACGCTGATGCTGCTGGCAGATGTACTCGCACGCGCGCTGGCCTTCCCCGGAGATCTGCCCGCAGGCGCAGTGCTGGCG  <  2:3032531/100‑1                                                                                                                                                                                                              GCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4512419‑4512524 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgccagTGAGCATGCTGCTGGGGGCCACGCTGATGCTGCTGGCAGATGTACTCGCACGCGCGCTGGCCTTCCCCGGAGATCTGCCCGCAGGCGCAGTGCTGGCG  <  NC_000913/4512556‑4512459

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