New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 16647 =5 (0.150)3 (0.120)
+GATGGTTCAACCCTTAAGTTAGCGCTTAT
3/186 3.5 53.0% noncoding (1261/1345 nt) IS186 repeat region
?NC_000913 = 4543616 2 (0.060)coding (502/549 nt) fimA major type 1 subunit fimbrin (pilin)
Rejected: Coverage evenness skew score above cutoff.

ATCTGAAGAATGCGACGCCTCGGTGCCTCGTTAAGACGATGCCTCGCGTTCTTCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/16701‑16647
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTCTGAAC  <  NC_000913/4543616‑4543539
                                                       |||||||||||||||||||||||||||||                                                                              
ATCTGAAGAATGCGACGCCTCGGTGCCTCGTTAAGACGATGCCTCGCGTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGG               >  2:144939/1‑149
 TCTGAAGAATGCGACGCCTCGGTGCCTCGTTAAGACGATGCCTCGCGTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGG              >  5:59010/1‑149
             GACGCCTCGGTGCCTCGTTAAGACGATGCCTCGCGTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTCTGAAC  <  1:144939/149‑1
                                                       |||||||||||||||||||||||||||||                                                                              
ATCTGAAGAATGCGACGCCTCGGTGCCTCGTTAAGACGATGCCTCGCGTTCTTCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/16701‑16647
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTCTGAAC  <  NC_000913/4543616‑4543539

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