New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 255386 =34 (1.140)7 (0.250) 5/280 2.0 16.5% coding (331/1458 nt) pepD aminoacyl‑histidine dipeptidase (peptidase D)
?NC_000913 255482 = 39 (1.390)coding (235/1458 nt) pepD aminoacyl‑histidine dipeptidase (peptidase D)
Rejected: Frequency below cutoff threshold.

GTATGGAAAATCGTAAACCGGTCGTCTTACAGGCCCACCTCGATATGGTGCCGCAGAAAAATAACGACACCGTGCATGACTTCACGAAAGATCCTATCCAGCCTTATATTGATGGCGAATGGGTTAAAGCGCGCGGCACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/255526‑255386
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcggcaccaTATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTT  >  NC_000913/255482‑255631
                                                                                                                                                                                                                                                                                                   
GTATGGAAAATCGTAAACCGGTCGTCTTACAGGCCCACCTCGATATGGTGCCGCAGAAAAATAACGACACCGTGCATGACTTCACGAAAGATCCTATCCAGCCTTATATTGATGGCGAATGGGTTAAAGCGCGCGGCACCATATCGAGG                                                                                                                                                <  1:162279/149‑1
GTATGGAAAATCGTAAACCGGTCGTCTTACAGGCCCACCTCGATATGGTGCCGCAGAAAAATAACGACACCGTGCATGACTTCACGAAAGATCCTATCCAGCCTTATATTGATGGCGAATGGGTTAAAGCGCGCGGCACCATATCGAGG                                                                                                                                                <  2:426881/149‑1
GTATGGAAAATCGTAAACCGGTCGTCTTACAGGCCCACCTCGATATGGTGCCGCAGAAAAATAACGACACCGTGCATGACTTCACAAAAGATCCTATCCAGCCTTATATTGATGGCGAATGGGTTAAAGCGCGCGGCACCATATCGAGG                                                                                                                                                <  2:168696/149‑1
                                                        AAAAATAACGACACCGTGCATGACTTCACGAAAGATCCTATCCAGCCTTATATTGATGGCGAATGGGTTAAAGCGCGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGA                                                                                       <  2:414316/150‑1
                                                                    ACCGTGCATGACTTCACGAAAGATCCTATCCAGCCTTATATTGATGGCGAATGGGTTAAAGCGCGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTAC                                                                          <  2:368253/151‑1
                                                                     CCGTGCATGACTTCACGAAAGATCCTATCCAGCCTTATATTGATGGCGAATGGGTTAAAGCGCGCGGCACCATATCGAGGTGGGCCTATAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACC                                                                         <  2:153150/151‑1
                                                                                                             TGATGGCGAATGGGTTAAAGCGCGCGGCGCCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCC                                 <  1:403125/151‑1
                                                                                                                                       GCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTG       >  1:267485/1‑151
                                                                                                                                         cCCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTCTACGAATCATGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTG       >  2:157931/2‑149
                                                                                                                                            ATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTT  >  2:12144/1‑151
                                                                                                                                            ATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTT  >  2:56469/1‑151
                                                                                                                                                                                                                                                                                                   
GTATGGAAAATCGTAAACCGGTCGTCTTACAGGCCCACCTCGATATGGTGCCGCAGAAAAATAACGACACCGTGCATGACTTCACGAAAGATCCTATCCAGCCTTATATTGATGGCGAATGGGTTAAAGCGCGCGGCACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/255526‑255386
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcggcaccaTATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTT  >  NC_000913/255482‑255631

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