New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4516891 =	37 (0.800)	5 (0.110)	5/278	3.0	13.0%	coding (827/954 nt)	fecR	anti‑sigma transmembrane signal transducer for ferric citrate transport; periplasmic FecA‑bound ferric citrate sensor and cytoplasmic FecI ECF sigma factor activator
	?	NC_000913	4516925 =	32 (0.730)					coding (793/954 nt)	fecR	anti‑sigma transmembrane signal transducer for ferric citrate transport; periplasmic FecA‑bound ferric citrate sensor and cytoplasmic FecI ECF sigma factor activator
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below cutoff threshold.

ATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACG  >  NC_000913/4516777‑4517073                                                                                                                                                     |                                                                                                                                                     aTGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGt                                                                                                                                                    >  2:642817/1‑151 (MQ=255)  tGTTTATCCAGTACCGCGTACTAGACTGAATTTTAACGGGACGCGTTTGCGCGATAACGTTCAGCATCGCATCGGTATTTTTCAGCGGGAACGGCCCACTCAGGCACAGCCCGGCAACGGCGGGGACGCAGCGCAGCACGCCGTTGCggca                                                                                                                                                   >  1:706292/1‑149 (MQ=255)  tGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTa                                                                                                                                                   <  1:149567/151‑1 (MQ=255)  tGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTa                                                                                                                                                   >  1:476905/1‑151 (MQ=255)       aTCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACggg                                                                                                                                              >  2:577115/1‑151 (MQ=255)            gTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAg                                                                                                                                          >  2:706813/1‑150 (MQ=255)            gTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGGAACGGGTTAg                                                                                                                                          >  1:403284/1‑150 (MQ=255)                       tAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAtcac                                                                                                                              <  1:263635/151‑1 (MQ=255)                              aaTTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCAccc                                                                                                                       <  2:39286/151‑1 (MQ=255)                                ttttAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAg                                                                                                                     <  2:186999/151‑1 (MQ=255)                                        ggAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTAt                                                                                                             >  1:527098/1‑151 (MQ=255)                                          aaGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCg                                                                                                           <  1:654386/151‑1 (MQ=255)                                                       aTAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGa                                                                                              >  2:79566/1‑151 (MQ=255)                                                            gTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCATTCTCATCCCGGCAACGGCGGGATCGCAGGGCAGCAAGCCTTTTCTTTAACTTTTTATCTTTGCTATCACCTCACCCAGCGGTTTATCTCTTAAGCTCAGGATTTc                                                                                          >  2:338969/1‑150 (MQ=255)                                                            gTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTc                                                                                          >  1:694795/1‑150 (MQ=255)                                                            gTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTc                                                                                          >  2:303007/1‑150 (MQ=255)                                                            gTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAAAGGGGTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTc                                                                                          >  2:418600/1‑150 (MQ=255)                                                                 ggATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGt                                                                                    >  2:321441/1‑151 (MQ=255)                                                                            gTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTa                                                                          >  2:698631/1‑150 (MQ=255)                                                                             tATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTAct                                                                        <  1:498608/151‑1 (MQ=255)                                                                             tATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTAc                                                                         <  2:90919/150‑1 (MQ=255)                                                                                       cGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATc                                                               <  1:303005/150‑1 (MQ=255)                                                                                           aaCGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGc                                                           <  2:134/150‑1 (MQ=255)                                                                                            aCGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCgg                                                         <  1:121576/151‑1 (MQ=255)                                                                                              gTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGtt                                                       >  2:447643/1‑151 (MQ=255)                                                                                                        ggCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGGGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTgcgc                                              >  1:47623/1‑150 (MQ=255)                                                                                                              gCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACt                                        >  2:534743/1‑150 (MQ=255)                                                                                                               cccGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAGGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTc                                       <  2:242601/150‑1 (MQ=255)                                                                                                                       cggcggGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGgcg                                >  2:209503/1‑149 (MQ=255)                                                                                                                                    gcgcagcACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCtt                 >  2:135630/1‑151 (MQ=255)                                                                                                                                     cgcagcACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGc                   <  1:83923/148‑1 (MQ=255)                                                                                                                                       cagcaCGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCa              >  2:712277/1‑151 (MQ=255)                                                                                                                                                    cGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACg  <  1:612920/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     ATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACG  >  NC_000913/4516777‑4517073

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