New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? minE = 195307098 (1.420)4 (0.070) 3/92 NT 4.4% coding (1039/1326 nt) pepP proline aminopeptidase P II
?minE = 1953107 88 (1.450)coding (1002/1326 nt) pepP proline aminopeptidase P II
Rejected: Frequency below/above cutoff threshold.

GCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  minE/1952996‑1953070
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gaaaggTGATGTTGATGAACTGATCGCTCAGAACGCCCATCGTCCTTTCTTTATGCATGGCCTTAGCCACTGGTTA  <  minE/1953107‑1953038
                                                                                                                                                 
GCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAG                                                                             >  1:2118252/1‑70
   ATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGA                                                                            <  1:78812/68‑1
   ATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGA                                                                            <  1:1932306/68‑1
    TCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAG                                                                         >  1:2502695/1‑70
            ATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAG                                                                             <  1:1038290/58‑1
               AACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGTGATGTTGAT                                                              >  1:2679646/1‑70
               AACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGTGATGTTGAT                                                              >  1:2869768/1‑70
                      ACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGG                                                                        >  1:1765289/1‑53
                                   CCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGA                                                                            <  1:1876074/36‑1
                                              CAGTGGCTAAGGCCATGCATAAAGAAAGGTGATGTTGATGAACTGATCGCTCAGAACGCCCATCGTCCTTT                              <  1:1741563/71‑1
                                                   GCTAAGGCCATGCATAAAGAAAGGTGATGTTGATGAACTGATCGCTCAGAACGCCCATCGTCCTTTCTTT                          <  1:3777776/70‑1
                                                                     GAAAGGTGATGTTGATGAACTGATCGCTCAGAACGCCCATCGTCCTTTCTTTATGCATGGCCTTAGCCACT       >  1:2196360/1‑71
                                                                      AAAGGTGATGTTGATGAACTGATCGCTCAGAACGCCCATCGTCCTTTCTTTATGCATGGCCTTAGCCAC        <  1:1546644/69‑1
                                                                      AAAGGTGATGTTGATGAACTGATCGCTCAGAACGCCCATCGTCCTTTCTTTATGCATGGCCTTAGCCAC        <  1:625825/69‑1
                                                                        AGGTGATGTTGATGAACTGATCGCTCAGAACGCCCATCGTCCTTTCTTTATGCATGGCCTTAGCCACTGGT    <  1:1770534/71‑1
                                                                         GGTGATGTTGATGAACTGATCGCTCAGAACGCCCATCGTCCTTTCTTTATGCATGGCCTTAGCCACTGGTT   <  1:3078805/71‑1
                                                                         GGTGATGTTGATGAACTGATCGCTCAGAACGCCCATCGTCCTTTCTTTATGCATGGCCTTAGCCACTGG     <  1:248282/69‑1
                                                                         GGTGATGTTGATGAACTGATCGCTCAGAACGCCCATCGTCCTTTCTTTATGCATGGCCTTAGCCACTGG     <  1:3824778/69‑1
                                                                          GTGATGTTGATGAACTGATCGCTCAGAACGCCCATCGTCCTTTCTTTATGCATGGCCTTAGCCACTGGTTA  <  1:599296/71‑1
                                                                          GTGATGTTGATGAACTGATCGCTCAGAACGCCCATCGTCCTTTCTTTATGCATGGCCTTAGCCACTGGTTA  <  1:1969854/71‑1
                                                                                                                                                 
GCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  minE/1952996‑1953070
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gaaaggTGATGTTGATGAACTGATCGCTCAGAACGCCCATCGTCCTTTCTTTATGCATGGCCTTAGCCACTGGTTA  <  minE/1953107‑1953038

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