New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 219634131 (0.900)5 (0.140) 5/266 3.6 13.5% intergenic (+8/‑133) metG/yehH methionyl‑tRNA synthetase/DUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function
?NC_000913 = 2196367 33 (0.960)intergenic (+34/‑107) metG/yehH methionyl‑tRNA synthetase/DUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGGTGGCTAACCTGGCACCACGTAAAATGCGCTTCGGTATCTCTGAAGGCATGGTGATGGCTGCCGGTCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2196197‑2196341
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AAAGGCGCTGTCGATGCAGCGCCTTGAAGGGGGATTATTTCACCTGATGACCCGGTTTAGCACCGGC  <  NC_000913/2196367‑2196301
                                                                                                                                                                                                                    
TGGTGGCTAACCTGGCACCACGTAAAATGCGCTTCGGTATCTCTGAAGGCATGGTGATGGCTGCCGGTCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTAAAG                                                                 <  3:135146/149‑1
                                                              GCCGGTCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTAAAGGCGCTGTCGATGCAGCGCCTTGAAGGGGGATTATTTCACCTGATGACCCGGTTTAGCACCGG   >  1:27836/1‑149
                                                               CCGGTCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGTCGATGCAGCGCCTTGAAGGGGGATTATTTCACCTGATGACCCGGTTTAGCACCGGC  <  2:27836/149‑1
                                                                                                            GTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTAAAGGCGCTGTCGATGCAGCGCCTTGAAGGGGGA                                   <  1:77253/71‑1
                                                                                                            GTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTAAAGGCGCTGTCGATGCAGCGCCTTGAAGGGGGA                                   >  2:77253/1‑71
                                                                                                                                                                                                                    
TGGTGGCTAACCTGGCACCACGTAAAATGCGCTTCGGTATCTCTGAAGGCATGGTGATGGCTGCCGGTCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2196197‑2196341
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AAAGGCGCTGTCGATGCAGCGCCTTGAAGGGGGATTATTTCACCTGATGACCCGGTTTAGCACCGGC  <  NC_000913/2196367‑2196301

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.