New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 421610052 (0.790)3 (0.050) 3/464 5.4 4.5% coding (623/1602 nt) aceB malate synthase A
?NC_000913 4216861 = 77 (1.220)coding (1384/1602 nt) aceB malate synthase A
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCACCATCAGTTACACCAATGAAGCAGGCAAAATTTACCAGCTCAAGCCCAATCCAGCGGTTTTGATTTGTCGGGTACGCGGTCTGCACTTGCCGGAAAAACATGTCACCTGGCGTGGTGAGGCAATCCCCGGCAGCCTGTTTGATTTTGCGCTCTATTTCTTCCACAACTATCAGGCACTGTTGGCAAAGGGCAGTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4215858‑4216100
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcagtggATCCATCATCAAAAAACGTTGAGCAATGGCAAACCGGTGA  >  NC_000913/4216861‑4216900
                                                                                                                                                                                                                                                                                           
CAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCACCATCAGTTACACCAATGAAGCAGGCAAAATTTACCAGCTCAAGCCCACTCCAGCGGTTTTGATTTGTCGGGTACGCGGTCTGCACTTGCCGGAAAAACATGTCCCCTGGCGTGGTGAGGCAATCCCCGGCAGCCTGTTTGATTTTGCGCTCTATTTCTTCCACAACTATCAGGCCCTGTTGGCAACGGGCAGTAGATCggca                                   >  2:165790/1‑246
                                                                                                                     GGTACGCGGTCTGCACTTGCCGGAAAAACATGTCACCTGGCGTGGTGAGGCAATCCCCGGCAGCCTGTTTGATTTTGCGCTCTATTTCTTCCACAACTATCAGGTACTGTTGGCAAAGGGCAGTGGATCCATCATCAAAAAACGTTGAGCAATGGCAAACCGGTGA  <  1:164153/166‑1
                                                                                                                     GGTACGCGGTCTGCACTTGCCGGAAAAACATGTCACCTGGCGTGGTGAGGCAATCCCCGGCAGCCTGTTTGATTTTGCGCTCTATTTCTTCCACAACTATCAGGTACTGTTGGCAAAGGGCAGTGGATCCATCATCAAAAAACGTTGAGCAATGGCAAACCGGTGA  >  2:164153/1‑166
                                                                                                                                                                                                                                                                                           
CAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCACCATCAGTTACACCAATGAAGCAGGCAAAATTTACCAGCTCAAGCCCAATCCAGCGGTTTTGATTTGTCGGGTACGCGGTCTGCACTTGCCGGAAAAACATGTCACCTGGCGTGGTGAGGCAATCCCCGGCAGCCTGTTTGATTTTGCGCTCTATTTCTTCCACAACTATCAGGCACTGTTGGCAAAGGGCAGTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4215858‑4216100
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcagtggATCCATCATCAAAAAACGTTGAGCAATGGCAAACCGGTGA  >  NC_000913/4216861‑4216900

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