New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3948771 =88 (1.070)3 (0.040) 3/250 NT 3.7% coding (1199/1521 nt) yifB magnesium chelatase family protein and putative transcriptional regulator
?NC_000913 3948803 = 73 (0.940)coding (1167/1521 nt) yifB magnesium chelatase family protein and putative transcriptional regulator
Rejected: Frequency below/above cutoff threshold.

GACATTATCAGGGAAACCATAACCGCTGCACGCCAGAACAGACATTACGTTATCTCAACCGGCTCTCGGGGCCCTTTCTCGACCGCTTCGATCTCTCACTGGAGATCCCATTACCACCCCCCGGCATTTTGAGTAAAACGGTAGTGCCGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3948921‑3948771
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgccgggGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCG  >  NC_000913/3948803‑3948945
                                                                                                                                                                                                                                                                                                      
GACTTTTTCAGGGAACCGATAACCGCTGCACGCCAGAACAGACATTACGTTATCTCAACCGGCTCTCGGGGCCCTTTCTCGACCGCTTCGATCTCTCACTGGAGATCCCATTACCACCCCCCGGCATTTTGAGTAAAACGGTAGTGCCG                                                                                                                                                   <  6:213182/149‑1
 ACATTATCAGGGAAACCATAACCGCTGCACGCCAGAACAGACATTACGTTATCTCAACCGGCTCTCGGGGCCCTTTCTCGACCGCTTCGATCTCTCACTGGAGATCCCATTACCACCCCCCGGCATTTTGAGTAAAACGGTAGTGCCGG                                                                                                                                                  <  6:132292/149‑1
                                            TTACGTTATCTCAACCGGCTCTCGGGGCCCTTTCTCGACCGCTTCGATCTCTCACTGGAGATCCCATTACCACCCCCCGGCATTTTGAGTAAAACGGTAGTGCCGG                                                                                                                                                  >  5:55570/1‑106
                                            TTACGTTATCTCAACCGGCTCTCGGGGCCCTTTCTCGACCGCTTCGATCTCTCACTGGAGATCCCATTACCACCCCCCGGCATTTTGAGTAAAACGGTAGTGCCGG                                                                                                                                                  <  6:55570/106‑1
                                                                                                                    CCCCCCGGCATTTTGAGTAAAACGGTAGTGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAA                                                                           >  3:245950/1‑105
                                                                                                                    CCCCCCGGCATTTTGAGTAAAACGGTAGTGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAA                                                                           <  4:245950/105‑1
                                                                                                                     CCCCCGGCATTTTGAGTAAAACGGTAGTGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAA                              >  1:222765/1‑149
                                                                                                                                                TGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCC   <  8:17376/149‑1
                                                                                                                                                 GCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTTAGATAACGTAATGTCTGTTCTGGCGTGCAGAGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCG  >  6:218375/1‑149
                                                                                                                                                                                                                                                                                                      
GACATTATCAGGGAAACCATAACCGCTGCACGCCAGAACAGACATTACGTTATCTCAACCGGCTCTCGGGGCCCTTTCTCGACCGCTTCGATCTCTCACTGGAGATCCCATTACCACCCCCCGGCATTTTGAGTAAAACGGTAGTGCCGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3948921‑3948771
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgccgggGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGCTGGGATTCATCGCCG  >  NC_000913/3948803‑3948945

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