New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4558976 =45 (1.330)3 (0.090) 3/480 NT 6.6% coding (552/1173 nt) iadA isoaspartyl dipeptidase
?NC_000913 4559022 = 41 (1.260)coding (506/1173 nt) iadA isoaspartyl dipeptidase

TGGCGCTAAGTCGCCTGACGGAAGCGGGCGTCACGTCAGTGGTTGGTCTGCTGGGCACCGACTCTATCTCTCGCCACCCGGAATCCCTGCTCGCCAAGACCCGTGCGCTCAATGAAGAAGGCATCAGCGCCTGGATGCTGACCGGCGCTTATCATGTCCCTTCCCGCACCATTACGGGTTCCGTGGAAAAAGACGTGGCGATTATCGATCGTGTGATTGGCGTGAAATGCGCCATCTCTGATCACCGTTCTGCCGCACCGGACGTTTATCACCTGGCCAATATGGCGGCAGAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4559268‑4558976
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcggcagaaCGGTGATCAGAGATGGCGCATTTCACGCCAATCACA  >  NC_000913/4559022‑4559057
                                                                                                                                                                                                                                                                                                                                         
TGGCGCAAAGTCGCCTGACGGCAGCGGGCGTCAGGGCAGTGGTTGGTCTGCTGGGCACCGACTCTATCTCTCGCCACCCGGAATCCCTGCTCGCCAAGACCCGTGCGCTCAATGAAGAAGGCAGCAGCGCCTGGATGCTGACCGGCGCTTATCATGTCCCTTCCCGCACCATTACGGGTTCCGTGGAAAAAGACGTGGCGATTATCGATCGTGTGATTGGCGTGAAATGCGCCATCTCTGATCACCGTTCTGCCGCACCGGACGTTTATCACCTGGCCAATATGGCGG                                           <  2:188170/288‑1
  GCGCTAAGTCGCCTGACGGAAGCGGGCGTCACGTCAGTGGTTGGTCTGCTGGGCACCGACTCTATCTCTCGCCACCCGGAATCCCTGCTCGCCAAGACCCGTGCGCTCAATGAAGAAGGCATCAGCGCCTGGATGCTGACCGGCGCTTATCATGTCCCTTCCCGCACCATTACGGGTTCCGTGGAAAAAGACGTGGCGATTATCGATCGTGTGATTGGCGTGAAATGCGCCATCTCTGATCACCGTTCTGCCGCACCGGACGTTTATCACCTGGCCAATATGGCGGCAG                                        >  1:33919/1‑289
  GCGCTAAGTCGCCTGACGGAAGCGGGCGTCACGTCAGTGGTTGGTCTGCTGGGCACCGACTCTATCTCTCGCCACCCGGAATCCCTGCTCGCCAAGACCCGTGCGCTCAATGAAGAAGGCATCAGCGCCTGGATGCTGACCGGCGCTTATCATGTCCCTTCCCGCACCATTACGGGTTCCGTGGAAAAAGACGTGGCGATTATCGATCGTGTGATTGGCGTGAAATGCGCCATCTCTGATCACCGTTCTGCCGCACCGGACGTTTATCACCTGGCCAATATGGCGGCAG                                        >  2:258360/1‑289
              CTGACGGAAGCGGGCGGCACTTCAGTGGGTGGTCGGCTGGGAACCGACTCTATCTCTCGCCACCCGGAATCCCTGCTCGCCAAGACCCGTGCGCTCAATGAAGAAGGCATCAGCGCCTGGATGCTGACCGGCGCTTATCATGTCCCTTCCCGCACCATTACGGGTTCCGTGGAAAAAGACGTGGCGATTATCGATCGTGTGATTGGCGTGAAATGCGCCATCTCTGATCACCGTTCTGCCGCACCGGACGTTTATCACCTGGCCAATATGGCGGCAGAACGGTGATCA                             <  2:314162/288‑1
                                                                                                                                                                                     CGTGGAAAAAGACGTGGCGATTATCGATCGTGTGATTGGCGTGAAATGCGCCATCTCTGATCACCGTTCTGCCGCACCGGACGTTTATCACCTGGCCAATATGGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACA  >  1:32112/1‑148
                                                                                                                                                                                     CGTGGAAAAAGACGTGGCGATTATCGATCGTGTGATTGGCGTGAAATGCGCCATCTCTGATCACCGTTCTGCCGCACCGGACGTTTATCACCTGGCCAATATGGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACA  <  2:32112/148‑1
                                                                                                                                                                                                                                                                                                                                         
TGGCGCTAAGTCGCCTGACGGAAGCGGGCGTCACGTCAGTGGTTGGTCTGCTGGGCACCGACTCTATCTCTCGCCACCCGGAATCCCTGCTCGCCAAGACCCGTGCGCTCAATGAAGAAGGCATCAGCGCCTGGATGCTGACCGGCGCTTATCATGTCCCTTCCCGCACCATTACGGGTTCCGTGGAAAAAGACGTGGCGATTATCGATCGTGTGATTGGCGTGAAATGCGCCATCTCTGATCACCGTTCTGCCGCACCGGACGTTTATCACCTGGCCAATATGGCGGCAGAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4559268‑4558976
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcggcagaaCGGTGATCAGAGATGGCGCATTTCACGCCAATCACA  >  NC_000913/4559022‑4559057

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