New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3160836 =34 (1.090)3 (0.100) 3/208 3.1 7.7% coding (311/2220 nt) ygiQ Radical SAM superfamily protein
?NC_000913 3160864 = 39 (1.300)coding (283/2220 nt) ygiQ Radical SAM superfamily protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ATCACCCAAGCTTCGGGATGGCGATTTGCGGTCGTATGCTGGAAGCGCAGGGCTTTCGCGTCGGGATCATCGCCCAGCCGGACTGGAGCAGCAAAGACGACTTTATGCGTCTGGGTAAACCGAATCTGTTTTTCGGTGTTAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3160977‑3160836
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttacCCAGACGCATAAAGTCGTCTTTGCTGCTCCAGTCCGGCTGGGCGATGATCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCGCAAATCGCCATCCCGAAGCTTGGGTGATCGACATA  >  NC_000913/3160864‑3160984
                                                                                                                                                                                                                                                                       
ATCACCCAAGCTTCGGGATGGCGATTTGCGGTCGTATGCTGGAAGCGCAGGGCTTTCGCGTCGGGATCATCGCCCAGCCGGACTGGAGCAGCAAAGACGACTTTATGCGTCTGGGTAAACCGAATCTGTTTTTCGGTGT                                                                                                                              >  2:538795/1‑139
       AAGCTTCGGGATGGCGATTTGCGGTCGTATGCTGGAAGCGCAGGGCTTTCGCGTCGGGATCATCGCCCAGCCGGACTGGAGCAGCAAAGACGACTTTATGCGTCTGGGTAAACCGAATCTGTTTTTCGGTGTTA                                                                                                                            >  1:321244/1‑134
       AAGCTTCGGGATGGCGATTTGCGGTCGTATGCTGGAAGCGCAGGGCTTTCGCGTCGGGATCATCGCCCAGCCGGACTGGAGCAGCAAAGACGACTTTATGCGTCTGGGTAAACCGAATCTGTTTTTCGGTGTTA                                                                                                                            <  2:321244/134‑1
                                                                                                                      ACCGAATCTGTTTTTCGGTGTTACCCAGACGCATAAAGTCGTCTTTGCTGCTC                                                                                              <  1:410472/53‑1
                                                                                                                      ACCGAATCTGTTTTTCGGTGTTACCCAGACGCATAAAGTCGTCTTTGCTGCTC                                                                                              >  2:410472/1‑53
                                                                                                                           ATCTGTTTTTCGGTTTACCCAGACGCATAAAGTCGTCTTTGCTGCTCCAGTCCGGCTGGGCGATGATCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCGCAAATCGCCATCCCGAAGCTTGGGTGATCGACATA  >  2:683773/1‑139
                                                                                                                                                                                                                                                                       
ATCACCCAAGCTTCGGGATGGCGATTTGCGGTCGTATGCTGGAAGCGCAGGGCTTTCGCGTCGGGATCATCGCCCAGCCGGACTGGAGCAGCAAAGACGACTTTATGCGTCTGGGTAAACCGAATCTGTTTTTCGGTGTTAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3160977‑3160836
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttacCCAGACGCATAAAGTCGTCTTTGCTGCTCCAGTCCGGCTGGGCGATGATCCCGACGCGAAAGCCCTGCGCTTCCAGCATACGACCGCAAATCGCCATCCCGAAGCTTGGGTGATCGACATA  >  NC_000913/3160864‑3160984

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