New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4096861 =35 (0.940)8 (0.230) 5/248 4.0 18.4% coding (588/1470 nt) rhaB rhamnulokinase
?NC_000913 4096886 = 38 (1.080)coding (563/1470 nt) rhaB rhamnulokinase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CAACAACCTGAACTTATTCCACACATTGCTCACGCTCTGCTGATGCCGGATTACTTCAGTTATCGCCTGACCGGCAAGATGAACTGGGAATATACCAACGCCACGACCACGCAACTGGTCAATATCAATAGCGACGACTGGGACGAGTCGCTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4097013‑4096861
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtcgctaTTGATATTGACCAGTTGCGTGGTCGTGGCGTTGGTATATTCCCAGTTCATCTTGCCGGTCAGGCGATAACTGAAGTAATCCGGCATCAGCAGAGCGTGAGCAATGTG  >  NC_000913/4096886‑4096992
                                                                                                                                                                                                                                                                    
CAACAACCTGAACTTATTCCACACATTGCTCACGCTCTGCTGATGCCGGATTACTTCAGTTATCGCCTGACCGGCAAGATGAACTGGGAATATACCAACGCCACGACCACGCAACTGGTCAATATCAATAGCGACGACTGGGACGAG                                                                                                                   >  5:20308/1‑147
CAACAACCTGAACTTATTCCACACATTGCTCACGCTCTGCTGATGCCGGATTACTTCAGTTATCGCCTGACCGGCAAGATGAACTGGGAATATACCAACGCCACGACCACGCAACTGGTCAATATCAATAGCGACGACTGGGACGAG                                                                                                                   <  6:20308/147‑1
                                                                                                                          TATCAATAGCGACGACTGGGACGAGTCGCTATTGATATTGACCAGTTGCGTGGTCGTGGCGTTGGTATATTCCCAGTTCATCTTGCCGGTCAGGCGATAACTGAAGTAATCCGGCATCAGCAGAGCGTGAGCAATGTG  <  5:155850/138‑1
                                                                                                                          TATCAATAGCGACGACTGGGACGAGTCGCTATTGATATTGACCAGTTGCGTGGTCGTGGCGTTGGTATATTCCCAGTTCATCTTGCCGGTCAGGCGATAACTGAAGTAATCCGGCATCAGCAGAGCGTGAGCAATGTG  >  6:155850/1‑138
                                                                                                                          TATCAATAGCGACGACTGGGACGAGTCGCTATTGATATTGACCAGTTGCGTGGTCGTGGCGTTGGTATATTCCCAGTTCATCTTGCCGGTCAGGCGATAAC                                       >  1:167747/1‑101
                                                                                                                          TATCAATAGCGACGACTGGGACGAGTCGCTATTGATATTGACCAGTTGCGTGGTCGTGGCGTTGGTATATTCCCAGTTCATCTTGCCGGTCAGGCGATAAC                                       <  2:167747/101‑1
                                                                                                                               ATAGCGACGACTGGGACGAGTCGCTATTGATATTGACCAGTTGCGTGGTCGTGGCGTTGGTATATTCCCAGTTCATCTTGCC                                                     <  1:177465/82‑1
                                                                                                                               ATAGCGACGACTGGGACGAGTCGCTATTGATATTGACCAGTTGCGTGGTCGTGGCGTTGGTATATTCCCAGTTCATCTTGCC                                                     >  2:177465/1‑82
                                                                                                                               ATAGCGACGACTGGGACGAGTCGCTATTGATATTGACCAGTTGCGTGGTCGTGGCGTTGGTATATTCCCAGTTCATCTTGCC                                                     <  7:64660/82‑1
                                                                                                                               ATAGCGACGACTGGGACGAGTCGCTATTGATATTGACCAGTTGCGTGGTCGTGGCGTTGGTATATTCCCAGTTCATCTTGCC                                                     >  8:64660/1‑82
                                                                                                                                                     GCTATTGATATTGACCAGTTGCGTGGTCGTGGCGTTGGTATATTCCCAGTTCATCTTGCCGGTCAGGCGATAACTG                                     <  2:78847/76‑1
                                                                                                                                                                                                                                                                    
CAACAACCTGAACTTATTCCACACATTGCTCACGCTCTGCTGATGCCGGATTACTTCAGTTATCGCCTGACCGGCAAGATGAACTGGGAATATACCAACGCCACGACCACGCAACTGGTCAATATCAATAGCGACGACTGGGACGAGTCGCTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4097013‑4096861
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtcgctaTTGATATTGACCAGTTGCGTGGTCGTGGCGTTGGTATATTCCCAGTTCATCTTGCCGGTCAGGCGATAACTGAAGTAATCCGGCATCAGCAGAGCGTGAGCAATGTG  >  NC_000913/4096886‑4096992

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