New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 460525479 (1.030)3 (0.040) 3/168 NT 4.0% coding (410/789 nt) fhuF ferric iron reductase involved in ferric hydroximate transport
?NC_000913 = 4605268 71 (0.980)coding (396/789 nt) fhuF ferric iron reductase involved in ferric hydroximate transport
Rejected: Frequency below cutoff threshold.

TGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4605151‑4605254
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gaaacCGGACGCGTCGCCTGTTTCTGGGTCGATGTGTGTGAAGATAAAAACGCAACACCACATTCGCCGCAGCATCGAATGGAAACGTTAATCAGCCAGGC  <  NC_000913/4605268‑4605173
                                                                                                                                                                                                        
TGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAG                                                                                                      <  2:1733454/100‑1
TGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAG                                                                                                      <  2:281172/100‑1
 nCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGA                                                                                                     >  2:1056026/2‑100
 GCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGA                                                                                                     >  1:1343789/1‑100
    TGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAAC                                                                                                  <  2:545880/100‑1
                           GCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACCGGACGCGTCGCCTGTTTCTGGG                                                                           <  1:560482/100‑1
                            CTGATTAGCGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACGTCGACCCAGAAACCGGACGCGTCGCCTGTTTCTGGGT                                                                          <  2:1109695/100‑1
                                          CATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACCGGACGCGTCGCCTGTTTCTGGGTCGATGTGTGTGAAG                                                            <  2:424178/100‑1
                                                                                                   GAAACCGGACGCGTCGCCTGTTTCTGGGTCGATGTGTGTGAAGATAAAAACGCAACACCACATTCGCCGCAGCATCGAATGGAAACGTTAATCAGCCAGG   >  2:449003/1‑100
                                                                                                     AACCGGACGCGTCGCCTGTTTCTGGGTCGATGTGTGTGAAGATAAAAACGCAACACCACATTCGCCGCAGCATCGAATGGAAACGTTAATCAGCCAGGC  >  2:529717/1‑99
                                                                                                                                                                                                        
TGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4605151‑4605254
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gaaacCGGACGCGTCGCCTGTTTCTGGGTCGATGTGTGTGAAGATAAAAACGCAACACCACATTCGCCGCAGCATCGAATGGAAACGTTAATCAGCCAGGC  <  NC_000913/4605268‑4605173

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