New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 34133576 (0.810)10 (0.110) 7/190 NT 12.5% coding (211/1383 nt) yahJ putative metallo‑dependent hydrolase domain deaminase
?NC_000913 = 341358 68 (0.770)coding (234/1383 nt) yahJ putative metallo‑dependent hydrolase domain deaminase

GCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/341224‑341335
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggcgCGCCGTGCGGGTCTGCACCGCCACGCCATTTT  <  NC_000913/341358‑341327
                                                                                                                                                
GCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGTGCGGGTC                     >  3:322270/1‑125
GCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGTGCGGGTC                     <  4:322270/125‑1
          GCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGCGCCGTGCGGGTCTGCACCGCCACGCCATTTT  <  3:300931/134‑1
          GCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGCGCCGTGCGGGTCTGCACCGCCACGCCATTTT  >  4:300931/1‑134
          GCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGG                                    <  3:591683/100‑1
          GCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGG                                    >  4:591683/1‑100
                         AAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGCGCCGTGCGGGT                      >  4:952467/1‑99
                         AAGCGAACAACACCATGAAAATCACTGACCCGCACTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGCGCCGTGCGGGT                      <  3:952467/99‑1
                                                       CGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGCGCCGTGCGGGTCTGCACCGCCACGCCATTTT  >  3:456970/1‑89
                                                       CGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGCGCCGTGCGGGTCTGCACCGCCACGCCATTTT  <  4:456970/89‑1
                                                                  CTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGCGCCGTGCGGGTCTGCACCGCCA           <  3:839656/69‑1
                                                                  CTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGCGCCGTGCGGGTCTACACCGCCA           >  4:839656/1‑69
                                                                                                                                                
GCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/341224‑341335
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggcgCGCCGTGCGGGTCTGCACCGCCACGCCATTTT  <  NC_000913/341358‑341327

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