New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 34121073 (0.530)6 (0.050) 6/216 NT 7.3% coding (86/1383 nt) yahJ putative metallo‑dependent hydrolase domain deaminase
?NC_000913 = 341232 83 (0.630)coding (108/1383 nt) yahJ putative metallo‑dependent hydrolase domain deaminase

TATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/341100‑341210
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtgccAGCTACCGCCGCATGGGCGAGCGGCACAGAGG  <  NC_000913/341232‑341201
                                                                                                                                               
TATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGT                                     <  3:1257983/108‑1
TATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGT                                     >  4:1257983/1‑108
     TTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCAGCTACCGCCGCATGGGCGAGCGGCACAGAG   <  1:674341/137‑1
     TTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCAGCTACCGCCGCATGGGCGAGCGGCACAGAG   >  2:674341/1‑137
                                AAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGT                                     <  3:1135373/76‑1
                                AAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGT                                     >  4:1135373/1‑76
                                                    TTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCAGCTACCGCCGCATGGGCGAGCGGCACAGAGG  >  3:1089684/1‑91
                                                    TTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCAGCTACCGCCGCATGGGCGAGCGGCACAGAGG  <  4:1089684/91‑1
                                                            CCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCAGCTACCGCCGCATGGGCGAGCG           <  1:1163634/74‑1
                                                            CCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCAGCTACCGCCGCATGGGCGAGCG           >  2:1163634/1‑74
                                                                                                                                               
TATTATTAAAACGTAAGGACGTAAGATGAAAGAAAGCAATAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/341100‑341210
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtgccAGCTACCGCCGCATGGGCGAGCGGCACAGAGG  <  NC_000913/341232‑341201

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