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.

CGCCTCTGAGTTTGGCGCAGTGAAACCGCTGGATGACGAGAGTACAAGCTGGACGAAGGACATCCTGAGCTTCAGCGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCGTTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4517043‑4516891
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccgttgcGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACG  >  NC_000913/4516925‑4517073
                                                                                                                                                                                                                                                                                                              
CGCCTCTGAGTTTGGCGCAGTGAAACCGCTGGATGACGAGAGTACAAGCTGGACGAAGGACATCCTGAGCTTCAGCGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCG                                                                                                                                                           <  2:209503/149‑1
                            CTGGATGACGAGAGTACAAGCTGGACGAAGGACATCCTGAGCTTCAGCGATAACCCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACC                                                                                                                             <  1:336493/151‑1
                            CTGGATGACGAGAGTACAAGCTGGACGAAGGACATCCTGAGCTTCAGCGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCGTTGCGGTAACGGGTTAGCGTGGCTATCAC                                                                                                                              <  2:720828/150‑1
                                                                           CGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCT                                                                              <  2:252099/151‑1
                                                                                acCCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTAC                                                                         <  1:728115/149‑1
                                                                                                                CGTTACCGCAACGGCGTGCTGCGCGGCAATCCCGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAAC                                         <  1:282575/151‑1
                                                                                                                                                        CGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACG  <  1:612920/150‑1
                                                                                                                                                                                                                                                                                                              
CGCCTCTGAGTTTGGCGCAGTGAAACCGCTGGATGACGAGAGTACAAGCTGGACGAAGGACATCCTGAGCTTCAGCGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCGTTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4517043‑4516891
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccgttgcGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACG  >  NC_000913/4516925‑4517073

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