New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4096011 =166 (0.970)3 (0.020) 3/266 NT 2.0% coding (1438/1470 nt) rhaB rhamnulokinase
?NC_000913 4096045 = 140 (0.850)coding (1404/1470 nt) rhaB rhamnulokinase
Rejected: Frequency below/above cutoff threshold.

CTCGACGCTCGGCAATATCGGCATCCAGTTAATGACGCTGGATGAACTCAACAATGTGGATGATTTCCGTCAGGTCGTCAGCACCACCGCGAATCTGACCACCTTTACCCCTAATCCTGACAGTGAAATTGCCCACTATGTGGCGCAGATTCACT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4096165‑4096011
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttcactGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAA  >  NC_000913/4096045‑4096193
                                                                                                                                                                                                                                                                                                                
CTCGACGCTCGGCAATATCGGCATCCAGTTAATGACGCTGGATGAACTCAACAATGTGGATGATTTCCGTCAGGTCGTCAGCACCACCGCGAATCTGACCACCTTTACCCCTAATCCTGACAGTGAAATTGCCCACTATGTGGCGCAGAT                                                                                                                                                            <  1:2501890/150‑1
 TCGACGCTCGGCAATATCGGCATCCAGTTAATGACGCTGGATGAACTCAACAATGTGGATGATTTCCGTCAGGTCGTCAGCACCACCGCGAATCTGACCACCTTTACCCCTAATCCTGACAGTGAAATTGCCCACTATGTGGCGCAGATT                                                                                                                                                           >  1:2143992/1‑150
 TCGACGCTCGGCAATATCGGCATCCAGTTAATGACGCTGGATGAACTCAACAATGTGGATGATTTCCGTCAGGTCGTCAGCACCACCGCGAATCTGACCACCTTTACCCCTAATCCTGACAGTGAAATTGCCCACTATGTGGCGCAGATT                                                                                                                                                           >  1:921155/1‑150
  CGACGCTCGGCAATATCGGCATCCAGTTAATGACGCTGGATGAACTCAACAATGTGGATGATTTCCGTCAGGTCGTCAGCACCACCGCGAATCTGACCACCTTTACCCCTAATCCTGACAGTGAAATTGCCCACTATGTGGCGCAGATTC                                                                                                                                                          <  2:2424710/150‑1
    ACGCTCGGCAATATCGGCATCCAGTTAATGACGCTGGATGAACTCAACAATGTGGATGATTTCCGTCAGGTCGTCAGCACCACCGCGAATCTGACCACCTTTACCCCTAATCCTGACAGTGAAATTGCCCACTATGTGGCGCAGATTCAC                                                                                                                                                        <  2:62541/150‑1
                                                                                                                     TGACAGTGAAATTGCCCACTATGTGGCGCAGATTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATT                                           >  1:1315160/1‑146
                                                                                                                     TGACAGTGAAATTGCCCACTATGTGGCGCAGATTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATT                                           <  2:1315160/146‑1
                                                                                                                             AAATTGCCCACTATGTGGCGCAATTTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGA                               >  2:383991/1‑150
                                                                                                                                                     TTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCGACAGGCCCGGCGATCAC       <  2:1593352/150‑1
                                                                                                                                                     TTCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCAC       <  2:2694781/150‑1
                                                                                                                                                      TCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACC      <  2:2185939/150‑1
                                                                                                                                                      TCACTGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACC      <  2:2763671/150‑1
                                                                                                                                                          TGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAA  <  2:1136262/150‑1
                                                                                                                                                          TGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAA  <  2:372876/150‑1
                                                                                                                                                                                                                                                                                                                
CTCGACGCTCGGCAATATCGGCATCCAGTTAATGACGCTGGATGAACTCAACAATGTGGATGATTTCCGTCAGGTCGTCAGCACCACCGCGAATCTGACCACCTTTACCCCTAATCCTGACAGTGAAATTGCCCACTATGTGGCGCAGATTCACT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4096165‑4096011
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttcactGTCAGGATTAGGGGTAAAGGTGGTCAGATTCGCGGTGGTGCTGACGACCTGACGGAAATCATCCACATTGTTGAGTTCATCCAGCGTCATTAACTGGATGCCGATATTGCCGAGCGTCGAGGCTTCAACAGGCCCGGCGATCACCCGAA  >  NC_000913/4096045‑4096193

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