New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4516187	87 (0.820)	3 (0.030)	3/264	NT	3.4%	coding (491/2325 nt)	fecA	TonB‑dependent outer membrane ferric citrate transporter and signal transducer; ferric citrate extracelluar receptor; FecR‑interacting protein
	?	NC_000913	= 4516268	90 (0.890)					coding (410/2325 nt)	fecA	TonB‑dependent outer membrane ferric citrate transporter and signal transducer; ferric citrate extracelluar receptor; FecR‑interacting protein
Rejected: Frequency below/above cutoff threshold.

GAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4516032‑4516187 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgcgcGTGACGTGATCCGCCGTGAGGATTTCGCCAAAACCGGCGCAACCACCATGCGTGAGGTACTTAACCGCATCCCTGGCGTCAGCGCGCCGGAAAACAACGGCACCGGCAGCCACGACCTGGCGATGAACTTTGGCATCCGGGGCCTG  <  NC_000913/4516268‑4516123                                                                                                                                                                                                                                                                                                                GAGCCCGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGACGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCG                                                                                                                                                          <  1:712331/150‑1  cgccccTGAAAGCTGCGGCTGACCGTAAGGCGCGAAGGGGACGGGGATGCCGTCCAGCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGG                                                                                                                                                         <  2:293688/144‑1   GCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCCGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGC                                                                                                                                                        >  1:1473062/1‑150   GCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGC                                                                                                                                                        >  2:1366004/1‑150   GCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGC                                                                                                                                                        >  2:1780137/1‑150     CAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGC                                                                                                                                                      >  1:47671/1‑150                                                  GCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCGTGACGTGATCCGCCGTGAGGATTTCGCCAAAACCGGCGCAAC                                                                                                         >  1:1445185/1‑150                                                                                       GGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCTGGCGCGCGTGACGTGATCCGCCGTGAGGATTTCGCCAAAACCGGCGCAACCACCATGCGTGAGGTACTTAACCGCATCCCTGGCGTC                                                                    <  1:1198619/150‑1                                                                                                  CCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCGTGACGTGATCCGCCGTGAGGATTTCGCCAAAACCGGCGCAACCACCATGCGTGAGGTACTTAACCGCATCCCTGGCGTCAGCGCGCCGGA                                                         <  2:1445185/150‑1                                                                                                                                                       GCGCGCGGGACGTGATCCGCCGTGAGGATTTCGCCAAAACCGGCGCAACCACCATGCGGGAGGTACTTAACCGCATCCCTGGCGGCAGCGCGCCGGAAAACAACGGCACCGGCAGCCACGACCTGGCGATGAACTTTGGCATCCGGGGCC    <  1:1584609/150‑1                                                                                                                                                        CGCGCGTGACGTGATCCGCCGTGAGGATTTCGCCAAAACCGGCGCAACCACCATGCGTGAGGTACTTAACCGCATCCCTGGCGTCAGCGCGCCGGAAAACAACGGCACCGGCCGCCACGCCCTGGCGATGAACTTTGGCATCCGGGGCCT   >  2:481974/1‑150                                                                                                                                                        CGCGCGTGACGTGATCCGCCGTGAGGATTTCGCCAAAACCGGCGCAACCACCATGCGTGAGGTACTTAACCGCATCCCTGGCGTCAGCGCGCCGGAAAACAACGGCACCGGCAGCCACGACCTGGCGATGAACTTTGGCATCCGGGGCCT   >  1:1723990/1‑150                                                                                                                                                        CGCGCGTGACGTGATCCGCCGTGAGGATTTCGCCAAAACCGGCGCAACCACCATGCGTGAGGTACTTAACCGCATCCCTGGCGTCAGCGCGCCGGAAAACAACGGCACCGGCAGCCACGACCTGGCGAGGAACTTTGGCATCCGGGGCCT   <  2:13953/150‑1                                                                                                                                                        CGCGCGTGACGTGATCCGCCGTGAGGATTTCGCAAAAACCGGCGCAAACACCATGCGTGAGGTACTTAACCGCATCCCTGGCGTCAGCGCGCCGGAAAACAACGGCACCGGCAGCCACGACCTGGCGATGAACTTTGGCAACCGGGGCCT   >  1:1360130/1‑150                                                                                                                                                         GCGCGGGACGTGATCCGCCGTGAGGATTTCGCCAAAACCGGCGCAACCCCCATGCGTGAGGTACTTAACCGCATCCCTGGCGTCAGCGCGCCGGAAAACAACGGCACCGGCAGCCACGACCTGGCGATGAACTTTGGCATCCGGGGCCTG  <  1:1101986/150‑1                                                                                                                                                                                                                                                                                                                GAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4516032‑4516187 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgcgcGTGACGTGATCCGCCGTGAGGATTTCGCCAAAACCGGCGCAACCACCATGCGTGAGGTACTTAACCGCATCCCTGGCGTCAGCGCGCCGGAAAACAACGGCACCGGCAGCCACGACCTGGCGATGAACTTTGGCATCCGGGGCCTG  <  NC_000913/4516268‑4516123

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