New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 3918584 =89 (1.520)4 (0.080) 4/92 NT 5.0% coding (634/852 nt) insK IS150 conserved protein InsB
?AP009048 3918607 = 72 (1.360)coding (611/852 nt) insK IS150 conserved protein InsB
Rejected: Frequency below/above cutoff threshold.

TGACCAGGGATGGCAGTATCGTATGAGAAGATATCAAAATATCCTTAAAGAACATGGTATTAAACAAAGCATGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AP009048/3918657‑3918584
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑catgtTCTTTAAGGATATTTTGATATCTTCTCATACGATACTGCCATCCCTGGTCAGAGTGCAGAACAGGA  >  AP009048/3918607‑3918672
                                                                                                                                            
TGACCAGGGATGGCAGTATCGTATGAGAAGATATCAAAATATCCTTAAAGAACATGGTATTAAACAAAGCA                                                                       >  1:1005287/1‑71
TGACCAGGGATGGCAGTATCGTATGAGAAGATATCAAAATATCCTTAAAGAACATGGTATTAAACAAAGCA                                                                       >  1:2988412/1‑71
  ACCAGGGATGGCAGTATCGTATGAGAAGATATCAAAATATCCTTAAAGAACATGGTATTAAACAAAGCATG                                                                     <  1:4517169/71‑1
  ACCAGGGATGGCAGTATCGTATGAGAAGATATCAAAATATCCTTAAAGAACATGGTATTAAACAAAGCATG                                                                     <  1:5047326/71‑1
  ACCAGGGATGGCAGTATCGTATGAGAAGATATCAAAATATCCTTAAAGAACATGGTATTAAACAAAGCATG                                                                     >  1:5130358/1‑71
  ACCAGGGATGGCAGTATCGTATGAGAAGATATCAAAATATCCTTAAAGAACATGGTATTAAACAAAGC                                                                        <  1:4021258/68‑1
   CCAGGGATGGCAGTATCGTATGAGAAGATATCAAAATATCCTTAAAGAACATGGTATTAAACAAAGCATGT                                                                    >  1:2815370/1‑71
   CCAGGGATGGCAGTATCGTATGAGAAGATATCAAAATATCCTTAAAGAACATGGTATTAAACAAAGCATGT                                                                    <  1:4199776/71‑1
                               TATCAAAATATCCTTAAAGAACATGGTATTAAACAAAGCA                                                                       >  1:2251971/1‑40
                                 TCAAAATATCCTTAAAGAACATGGTATTAAACAAAGCATGT                                                                    <  1:5296088/41‑1
                                                 GAACATGGTATTAAACAAAGCATGTTCTTTAAGGATATTTTGATATCTTCTCATACGATACTGCCATCC                        <  1:4024362/69‑1
                                                   ACATGGTATTAAACAAAGCATGTTCTTTAAGGATATTTTGATATCTTCTCA                                        >  1:5337222/1‑51
                                                    CATGGTATTAAACAAAGCATGTTCTTTAAGGATATTTTGATATCTTCTCA                                        >  1:2688685/1‑50
                                                            TAAACAAAGCATGTTCTTTAAGGATATTTTGATATCTTCTCATACGATACTGCCATCCCTGGTCAGAGTG            >  1:987698/1‑70
                                                                     CATGTTCTTTAAGGATATTTTGATATCTTCTCATACGATACTGCCATCCCTGGTCAGAGTGCAGAACAGGA  <  1:2992752/71‑1
                                                                       TGTTCTTTAAGGATATTTTGATATCTTCTCATACGATACTGCCATCCCTGGTCAG                <  1:1806638/55‑1
                                                                        GTTCTTTAAGGATATTTTGATATCTTCTCATACGATACTGCCATCCCTGGTCAGAGTGCAGAA       >  1:3981142/1‑63
                                                                         TTCTTTAAGGATATTTTGATATCTTCTCATACGATACTGCCATCCCTGGTCAGAG              >  1:691022/1‑55
                                                                                                                                            
TGACCAGGGATGGCAGTATCGTATGAGAAGATATCAAAATATCCTTAAAGAACATGGTATTAAACAAAGCATGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AP009048/3918657‑3918584
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑catgtTCTTTAAGGATATTTTGATATCTTCTCATACGATACTGCCATCCCTGGTCAGAGTGCAGAACAGGA  >  AP009048/3918607‑3918672

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