New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4515778135 (0.980)9 (0.070) 6/166 NT 6.4% coding (900/2325 nt) fecA TonB‑dependent outer membrane ferric citrate transporter and signal transducer; ferric citrate extracelluar receptor; FecR‑interacting protein
?NC_000913 = 4515794 139 (1.080)coding (884/2325 nt) fecA TonB‑dependent outer membrane ferric citrate transporter and signal transducer; ferric citrate extracelluar receptor; FecR‑interacting protein
Rejected: Frequency below cutoff threshold.

CAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4515676‑4515778
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccgcgcTGCTCTACTCCGGCACGCGCGGCAGTGACTG‑‑GCGCGAGCACAGCGCCACCCGCATCGACGACCTGATGCTGAAAAGCAAATATGCGCCGGATGAGGTGC  <  NC_000913/4515794‑4515696
                                                                                                                                                                                                            
CAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCCGTCACTGCCG                                                                                                          >  2:1308614/1‑100
CAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCG                                                                                                          >  1:1949294/1‑100
CAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCG                                                                                                          >  2:2454684/1‑100
  GCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCG                                                                                                        >  1:2580602/1‑100
  GCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCG                                                                                                        <  2:594095/100‑1
   CAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGC                                                                                                            <  1:149721/95‑1
      GCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCTGC                                                                                                    >  2:3769105/1‑100
                GTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCTGCTCTACTCCGG                                                                                          <  1:3769079/100‑1
                        CTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCTGCTCTACTCCGGCACGCGCG                                                                                  <  1:108096/100‑1
                          gatctGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCTGCTCTACTCCGGCACGCGCGGC                                                                                <  2:396779/95‑1
                          CATCCGGCGCATATTTGTTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCTGCTCTACTCCGGCACGCGCGGC                                                                                <  2:290225/100‑1
                          CATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCTGCTCTACTCCGGCACGCGCGGC                                                                                <  2:1383569/100‑1
                            TCCTGCGCATATTTGCTTTTCAGCAGCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCTGCTCTACTCCGGCACGCGCGGCAG                                                                              <  1:3248955/100‑1
                            TCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCTGCTCTACTCCGGCACGCGCGGCAG                                                                              <  2:1674455/100‑1
                             CCGGCGCATATTTGCTTTTCAGCATCAGGCCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCTGCTCTACCCCGGCACGCGCGGCAGT                                                                             <  1:951277/100‑1
                               GGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCTGCTCTACTCCGGCACGCGCGGCAGatc                                                                           >  1:396773/1‑97
                                                                                                  CGCGCTGCTCTACTCCGGCACGCGCGGCAGTGACTG‑‑GCGCGAGCACAGCGCCACCCGCATCGACGACCTGATGCTGAAAAGCAAATATGCGCCGGATGAG      >  1:132357/1‑100
                                                                                                  CGCGCTGCTCTACTCCGGCACGCGCGGCAGTGACTG‑‑GCGCGAGCACAGCGCCACCCGCATCGACGACCTGATGCTGAAAAGCAAATATGCGCCGGATGAG      >  2:215100/1‑100
                                                                                                  CGCGCTGCTCTACTCCGGCACGCGCGGCAGTGACTG‑‑GCGCGAGCACAGCGCCACCCGCATCGACGACCTGATGCTGAAAAGCAAATATGCGCCGGATGAG      >  1:3774032/1‑100
                                                                                                    CGCTGCTCTACCCCGGCACGCGCGGCAGTGACTG‑‑GCGCGAGCACAGCGCCACCCGCATCGACGACCTGATGCTGAAAAGCAAATATGCGCCGGATGAGGT    <  1:3382501/100‑1
                                                                                                     GCTGCTCTACTCCGGCACGCGCGGCAGTGACTG‑‑GCGCGAGCACAGCGCCACCCGCATCGACGACCTGATGCTGAAAAGCAAATATGCGCCGGATGAGGTG   <  1:2745171/100‑1
                                                                                                      CTGCTCTACTCCGGCACGCGCGGCAGTGACTG‑‑GCGCGAGCACAGCGCTACCCGCATCGACGACCTGATGCTGAAAAGCAAATATGCGCCGGATGAGGTGC  <  1:596426/100‑1
                                                                                                      CTGCTCTACTCCGGCACGCGCGGCAGTGACTG‑‑GCGCGAGCACAGCGCCACCCGCATCGACGACCTGATGCTGAAAAGCAAATATGCGCCGGATGAGGTGC  <  1:596869/100‑1
                                                                                                      CTGCTCTACTCCGGCACGCGCGCCAGTCACTGCCGCGCGTGCACAGCGCCACCCGCATCGACGACCTGATGCTGAAAAGCAAATATGCGCCGGATGAGGT    >  2:2064080/1‑100
                                                                                                                                                                                                            
CAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4515676‑4515778
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccgcgcTGCTCTACTCCGGCACGCGCGGCAGTGACTG‑‑GCGCGAGCACAGCGCCACCCGCATCGACGACCTGATGCTGAAAAGCAAATATGCGCCGGATGAGGTGC  <  NC_000913/4515794‑4515696

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