New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4511780	185 (0.990)	4 (0.020)	4/230	NT	2.1%	coding (635/957 nt)	fecD	ferric citrate ABC transporter permease
	?	NC_000913	= 4511807	189 (1.050)					coding (608/957 nt)	fecD	ferric citrate ABC transporter permease
Rejected: Frequency below/above cutoff threshold.

ACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4511629‑4511780 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgcgATCTCGACCTCCTTGCACTCGGCGATGCGCGCGCCACCACGCTCGGTGTGTCGGTGCCCCATACCCGATTCTGGGCTTTGTTACTAGCTGTCGCCATGACATCTACCGGCGTGGCCGCCTGCGGCCCGATTAGCTTTATTGGTCTCGT  <  NC_000913/4511807‑4511660                                                                                                                                                                                                                                                                                                              ACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGC                                                                                                                                                         >  3:1629211/1‑149                             ACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGC                                                                                                                                                         <  4:231568/121‑1                             ACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGC                                                                                                                                                         >  3:231568/1‑121                                                         CGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGATCTCGACCTCCTTGCACTCGGCGATGCGCG                                                                                                                       >  4:1186014/1‑127                                                         CGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGATCTCGACCTCCTTGCACTCGGCGATGCGCG                                                                                                                       <  3:1186014/127‑1                                                                    ACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGATCTCGACCTCCTTGCACTCGGCGATGCGCGC                                                                                                                      >  1:314272/1‑117                                                                    ACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCGATCTCGACCTCCTTGCACTCGGCGATGCGCGC                                                                                                                      <  2:314272/117‑1                                                                                                                                                     CGCGATCTCGACCTCCTTGCACTCGGCGATGCGCGCGCCACCACGCTCGGTGTGTCGGTGCCCCATACCCGATTCTGGGCTTTGTTACTAGCTGTCGCCATGACATCTACCGGCGTGGCCGCCTG                             >  1:951347/1‑125                                                                                                                                                     CGCGATCTCGACCTCCTTGCACTCGGCGATGCGCGCGCCACCACGCTCGGTGTGTCGGTGCCCCATACCCGATTCTGGGCTTTGTTACTAGCTGTCGCCATGACATCTACCGGCGTGGCCGCCTG                             <  2:951347/125‑1                                                                                                                                                     CGCGATCTCGACCTCCTTGCACTCGGCGATGCGCGCGCCACCACGCTCGGTGTGTCGGTGCCCCATACCCGATTCTGGGCTTTGTTACTAGCTGT                                                           <  4:577345/95‑1                                                                                                                                                     CGCGATCTCGACCTCCTTGCACTCGGCGATGCGCGCGCCACCACGCTCGGTGTGTCGGTGCCCCATACCCGATTCTGGGCTTTGTTACTAGCTGT                                                           >  3:577345/1‑95                                                                                                                                                      GCGATCTCGACCTCCTTGCACTCGGCGATGCGCGCGCCACCACGCTCGGTGTGTCGGTGCCCCATACCCGATTCTGGGCTTTGTTACTAGCTGTCGCCATGACATCTACCGGCGTGGCCGCAGGCGGCCac                      <  2:293268/131‑3                                                                                                                                                      GCGATCTCGACCTCCTTGCACTCGGCGATGCGCGCGCCACCACGCTCGGTGTGTCGGTGCCCCATACCCGATTCTGGGCTTTGTTACTAGCTGTCGCCATGACATCTACCGGCGTGGCCGCAGGCGGCCac                      >  1:293268/1‑129                                                                                                                                                       CGATCTCGACCTCCTTGCACTCGGCGATGCGCGCGCCACCACGCTCGGTGTGTCGGTGCCCCATACCCGATTCTGGGCTTTGTTACTAGCTGTCGCCATGACATCTACCGGCGTGGCCGCCTGCGG                          >  2:24185/1‑126                                                                                                                                                       CGATCTCGACCTCCTTGCACTCGGCGATGCGCGCGCCACCACGCTCGGTGTGTCGGTGCCCCATACCCGATTCTGGGCTTTGTTACTAGCTGTCGCCATGACATCTACCGGCGTGGCCGCCTGCGG                          <  1:24185/126‑1                                                                                                                                                        GATCTCGACCTCCTTGCACTCGGCGATGCGCGCGCCACCACGCTCGGTGTGTCGGTGCCCCATACCCGATTCTGGGCTTTGTTACTAGCTGTCGCCATGACATCTACCGGCGTGGCCGCCTGCGGCCCGATTAGCTTTATTGGTCTCGT  <  2:1686897/149‑1                                                                                                                                                        GATCTCGACCTCCTTGCACTCGGCGATGCGCGCGCCACCACGCTCGGTGTGTCGGTGCCCCATACCCGATTCTGGGCTTTGTTACTAGCTGTCGCCATGACATCTACCGGCGTGGCCGCCTGCGGCCCGATTAGCTTTATTGGTCTCGT  <  4:1537702/149‑1                                                                                                                                                                                                                                                                                                              ACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACACCGAGCGTGGTGGCGCGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4511629‑4511780 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgcgATCTCGACCTCCTTGCACTCGGCGATGCGCGCGCCACCACGCTCGGTGTGTCGGTGCCCCATACCCGATTCTGGGCTTTGTTACTAGCTGTCGCCATGACATCTACCGGCGTGGCCGCCTGCGGCCCGATTAGCTTTATTGGTCTCGT  <  NC_000913/4511807‑4511660

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