New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 5396823 (0.620)7 (0.190) 5/276 4.0 23.2% coding (735/1287 nt) surA peptidyl‑prolyl cis‑trans isomerase (PPIase)
?NC_000913 = 53983 24 (0.660)coding (720/1287 nt) surA peptidyl‑prolyl cis‑trans isomerase (PPIase)
Rejected: Coverage evenness skew score above cutoff.

GGTTTCAGCAGAATATGGCGAGCATGAACTTCGGTCACCGAGATATTTTTGCTTTCGCCGCGCAGGTCGTTAACTTTCAGAATATGGAAGCCAACGCCGGAACGAATCGGGCCAACAATGTCGCCTTTCTTCGCGGTGCTTAATGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/53822‑53968
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgccCGGGATCTTCGCCCAGGCATTAAGCACC  <  NC_000913/53983‑53956
                                                                                                                                                                               
GGTTTCAGCAGAATATGGCGAGCATGAACTTCGGTCACCGAGATATTTTTGCTTTCGCCGCGCAGGTCGTTAACTTTCAGAATATGGAAGCCAACGCCGGAACGAATCGGGCCAACAATGTCGCCTTTCTTCGCGGTGCTTAATGCCCG                            <  5:56038/149‑1
        CAGAATATGGCGAGCATGAACTTCGGTCACCGAGATATTTTTGCTTTCGCCGCGCAGGTCGTTAACTTTCAGAATATGGAAGCCAACGCCGGAACGAATCGGGCCAACAATGTCGCCTTTCTTCGCGGTGCTTAATGCCCGGGATCTTC                    >  4:65928/1‑149
                        TGAACTTCGGTCACCGAGATATTTTTGCTTTCGCCGCGCAGGTCGTTAACTTTCAGAATATGGAAGCCAACGCCGGAACGAATCGGGCCAACAATGTCGCCTTTCTTCGCGGTGCTTAATGCCCGGGATCTTCGCCCAGGCATTAAGCA    <  3:65928/149‑1
                                                               AGGTCGTTAACTTTCAGAATATGGAAGCCACCGCCGGAACGAAGCGGGCCAACAATGTCGCCGTTCTTCGCGGTGCTTAATGCCCGGGATCTTCGCCCAGGCATTAAGCACC  <  8:10451/112‑1
                                                               AGGTCGTTAACTTTCAGAATATGGAAGCCAACGCCGGAACGAATCGGGCCAACAATGTCGCCTTTCTTCGCGGTGCTTAATGCCCGGGATCTTCGCCCAGGCATTAAGCACC  >  7:8366/1‑112
                                                               AGGTCGTTAACTTTCAGAATATGGAAGCCAACGCCGGAACGAATCGGGCCAACAATGTCGCCTTTCTTCGCGGTGCTTAATGCCCGGGATCTTCGCCCAGGCATTAAGCACC  <  8:8366/112‑1
                                                               AGGTCGTTAACGTTCAGAATATGGAAGCCAACGCCGGAACGAATCGGGCCAACAATGTCGCCTTTCTTCGCGGTGCTTAATGCCCGGGATCTTCGCCCAGGCATTAAGCACC  >  7:10451/1‑112
                                                                                                                                                                               
GGTTTCAGCAGAATATGGCGAGCATGAACTTCGGTCACCGAGATATTTTTGCTTTCGCCGCGCAGGTCGTTAACTTTCAGAATATGGAAGCCAACGCCGGAACGAATCGGGCCAACAATGTCGCCTTTCTTCGCGGTGCTTAATGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/53822‑53968
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgccCGGGATCTTCGCCCAGGCATTAAGCACC  <  NC_000913/53983‑53956

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.