New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4513364118 (0.690)5 (0.030) 5/242 NT 4.2% coding (46/999 nt) fecC ferric citrate ABC transporter permease
?NC_000913 = 4513381 114 (0.710)coding (29/999 nt) fecC ferric citrate ABC transporter permease
Rejected: Frequency below/above cutoff threshold.

GCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4513210‑4513364
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtgctgcTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCCTGCTGCCTGGACACACGCCAACGCTACCAGAAGCGCTGGTGCAAAAC  <  NC_000913/4513381‑4513236
                                                                                                                                                                                                                                                                                                             
GCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAG                                                                                                                                                          <  3:569411/149‑1
GCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAG                                                                                                                                                          >  3:586673/1‑149
    CCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCT                                                                                                                                                      >  4:133261/1‑149
      AGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGC                                                                                                                                                    >  3:994999/1‑149
                  AAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGC                                                                                                                                                    >  1:460928/1‑137
                  AAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGC                                                                                                                                                    <  2:460928/137‑1
                              TGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCAC                                                                                                                            <  2:1138506/149‑1
                                CACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTT                                                                                                                          <  1:1471218/149‑1
                                CACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGC                                                                                                                                                       >  2:289982/1‑120
                                CACCAGCGCTTCTGGTAGCGTTGGCGGGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGC                                                                                                                                                       <  1:289982/120‑1
                                    AGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTA                                                                                                                      <  1:874035/149‑1
                                                                                      ATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTAT                                                                                                                        <  3:297738/97‑1
                                                                                      ATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCTGTGGGGGCTTCCCGTTGCAGCACTTAT                                                                                                                        >  4:297738/1‑97
                                                                                                                                                       CTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCCTGCTGCCTGGACACACGCCAACGCTACCAGAAGCGCTGGTGCAAAA   <  3:1200657/149‑1
                                                                                                                                                       CTGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCCTGCTGCCTGGACACACGCCAACGCTACCAGAAGCGCTGGTGCAAAA   <  4:1366529/149‑1
                                                                                                                                                        TGCTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCCTGCTGCCTGGACACACGCCAACGCTACCAGAAGCGCTGGTGCAAAAC  <  1:1098590/149‑1
                                                                                                                                                                                                                                                                                                             
GCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4513210‑4513364
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtgctgcTGTGGGGGCTTCCCGTTGCAGCACTTATTATTATTTTCTGGCTGAGTCTGTTTTGCTACTCGGCCATTCCTGTTTCCGGAGCAGATGCAACCCGCGCCCTGCTGCCTGGACACACGCCAACGCTACCAGAAGCGCTGGTGCAAAAC  <  NC_000913/4513381‑4513236

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