New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 550165116 (0.960)3 (0.030) 3/512 NT 2.2% coding (539/816 nt) ylbF putative anaerobic allantoin catabolic oxamate carbamoyltransferase, DUF2877 family protein
?NC_000913 = 551959 150 (1.260)coding (636/1068 nt) purK N5‑carboxyaminoimidazole ribonucleotide synthase
Rejected: Frequency below/above cutoff threshold.

GGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/549876‑550165
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcGATTATGCAGGAGCTGGGCTATGTGGGCGTGATGGCGATGGAGTGTTTTGTCACCCCGCAAGGTCTGTTGATCAACGAACTGGCACCGCGTGTGCATAACAGCGGTCACTGGACACAAAACGGTGCCAGCATCAGCCAGTTTGAGCTGCATCTGCGGGCGATTACCGATCTGCCGTTACCGCAACCAGTGGTGAATAATCCGTCGGTGATGATCAATCTGATTGGTAGCGATGTGAATTATGACTGGCTGAAACTGCCGCTGGTGCATCTGCACTGGTACGACAAA  <  NC_000913/551959‑551674
                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
GGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGCGTCATGATGACACGCTGAGCGGCATGCTGCTGGCGG                                                                                                                                                                                                                                                                                                 >  2:121462/1‑289
                                                                         TTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCGATTATGCAGGAGCTGGGCTATGTGGGCGTGATGGCGATGGAGTGTTTTGTCACCCCGCAAGGTCTGTTGA                                                                                                                                                                                                                         >  2:203713/1‑288
                                                                                                                             GGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGCGATTATGCAGGAGCTGGGCTATGTGGGCGTGATGGCGATGGAGTGTTTTGTCACCCCGCAAGGTCTGTTGATCAACGAACTGGCACCGCGTGTGCATAACAGCGGTCACTGGACACAAAACGGT                                                                                                                                                                    <  1:203713/289‑1
                                                                                                                                                                                                                                                                                       tgGCTGGCGGCCATTATGCAGGGGCTGGGCTATGTGGTGGGGTTGGCGATGGAATGTTTTGTCACCCCCCAAGGTCTGTTGAGCAACGAACTGGCACCGCGGGTGCATAACAGCGGTCACTGGACACAAAACGGTGCCAGCATCAGCCAGTTTGAGCTGCATCTGCAGGCGATTACCGATCTGCCGTTACCGCAACCAGTGGTGAATAATCCGTCGGTGATTATCAATCTGATTGGTAGCGATGTGAATTATGACTGGCTGAAACTGCCGCTGGTGCATCTGCACTGG           <  1:239319/286‑1
                                                                                                                                                                                                                                                                                       aTGCTGGCGGGGTTTATGAAGGTGCTGGGCTATGTGGGCGTAATGGCGATGTAGTGTTTTGTCACCCCGCAAGGTCTGTTGATCAACGAACTGGCACCGCGTGTGCATAACAGCGGTCACTGTACACAAAACGGTGCCAGCATCAGCCAGTTTGAGCTGCATCTGCGGTCGATTACCGATCTGCCGTTACCGCAACCAGTTGTGAATAATCCGTCGGTGATGGTCAATCTGATTGGTAGCGATGTGAATTATGACTGGCTGAAACTGCCGCTGGTGCATCTGCACTGG           <  2:452708/287‑1
                                                                                                                                                                                                                                                                                               GGCGATTATGCAGGAGCTGGGCTATGTGGGCGTGATGGCGATGGAGTGTTTTGTCACCCCGCAAGGTCTGTTGATCAACGAACTGGCACCGCGTGTGCATAACAGCGGTCACTGGACACAAAACGGTGCCAGCATCAGCCAGTTTGAGCTGCATCTGCGGGCGATTACCGATCTGCCGTTACCGCAACCAGTGGTGAATAATCCGTCGGTGATGATCAATCTGATTGGTAGCGATGTGAATTATGACTGGCTGAAACTGCCGCTGGTGCATCTGCACTGGTACGACAAA  >  2:426670/1‑289
                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                
GGGCGTTTCACGGTTAAACAGCCACAGCGTTATTGTTTGCTGCGTATTACGCCGCCTGCGCATCCTCAACCACTTGCAGCTGCATGGATGCAACGCGCGGAGGAAACCGGGCTTTTCGGGCCACTGGCGTTGGCGGCAAGCGATCCGCTGCCTGCTGAGTTACGCCAGTTTCGTCACTGTTTTCAGGCCGCGCTCAATGGCGTTAAGACCGACTGGCGGCACTGGCTGGGTAAAGGCCCCGGATTAACGCCGAGTCATGATGACACGCTGAGCGGAATGCTGCTGGCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/549876‑550165
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcGATTATGCAGGAGCTGGGCTATGTGGGCGTGATGGCGATGGAGTGTTTTGTCACCCCGCAAGGTCTGTTGATCAACGAACTGGCACCGCGTGTGCATAACAGCGGTCACTGGACACAAAACGGTGCCAGCATCAGCCAGTTTGAGCTGCATCTGCGGGCGATTACCGATCTGCCGTTACCGCAACCAGTGGTGAATAATCCGTCGGTGATGATCAATCTGATTGGTAGCGATGTGAATTATGACTGGCTGAAACTGCCGCTGGTGCATCTGCACTGGTACGACAAA  <  NC_000913/551959‑551674

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