New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 3782731 =64 (1.280)3 (0.070) 3/90 NT 4.9% coding (721/1494 nt) yidJ predicted sulfatase/phosphatase
?AP009048 3782753 = 61 (1.380)coding (743/1494 nt) yidJ predicted sulfatase/phosphatase
Rejected: Frequency below/above cutoff threshold.

TCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AP009048/3782777‑3782731
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atcaccATCCGCTCTATTTTGCCTGTAATGACTTTGTTGATGACCAAATCGGACGGGTCATCAATGCCTTAAC  >  AP009048/3782753‑3782819
                                                                                                                  
TCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCAC                                                                      >  1:1199804/1‑46
                       ATGGTGATAAAGCCCGTCATCACCATCCGCTCTATTTTGCCTGTAATGA                                            <  1:2148841/49‑1
                          GTGATAAAGCCCGTCATCACCATCCGCTCTATTTTGCCTGTAATGACT                                          <  1:473062/48‑1
                           TGATAAAGCCCGTCATCACCATCCGCTCTATTTTGCCTGTAATGACTTTGTTGATGACCAAATCGGACGGG                  >  1:2374497/1‑71
                                         ATCACCATCCGCTCTATTTTGCCTGTAATGACTTTGTTGATGACCAAATCGGACGGGTCATCAATGCCTT     <  1:4300729/70‑1
                                         ATCACCATCCGCTCTATTTTGCCTGTAATGACTTTGTTGATGACCAAATCGGACGGGTCATCAATGCCTT     <  1:957596/70‑1
                                         ATCACCATCCGCTCTATTTTGCCTGTAATGACTTTGTTGATGACCAAATCGGACGGGTCATCAATGCCTT     <  1:99399/70‑1
                                          TCACCATCCGCTCTATTTTGCCTGTAATGACTTTGTTGATGACCAAATCGGACGGGTCATCAATGCCTTA    <  1:4273959/70‑1
                                          TCACCATCCGCTCTATTTTGCCTGTAATGACTTTGTTGATGACCAAATCG                        <  1:2196006/50‑1
                                          TCACCATCCGCTCTATTTTGCCTGTAATGACTTTGTTGATGACCA                             >  1:3297687/1‑45
                                           CACCATCCGCTCTATTTTGCCTGTAATGACTTTGTTGATGACCAAATCGGACGGGTCATCAATGCCTTAA   >  1:208159/1‑70
                                           CACCATCCGCTCTATTTTGCCTGTAATGACTTTGTTGAT                                  >  1:2322229/1‑39
                                            ACCATCCGCTCTATTTTGCCTGTAATGACTTTGTTGATGACCAAATCGGACGGGTCATCAATGCCTTAAC  >  1:871492/1‑70
                                            ACCATCCGCTCTATTTTGCCTGTAATGACTTTGTTGATGACCAAATCGGAC                     >  1:3240625/1‑51
                                              CATCCGCTCTATTTTGCCTGTAATGACTTTGTTGATGACCAA                            <  1:2119740/42‑1
                                              CATCCGCTCTATTTTGCCTGTAATGACTTTGTTGATGACCAA                            <  1:2125200/42‑1
                                                                                                                  
TCATTACAGGCAAAATAGAGCGGATGGTGATAAAGCCCGTCATCACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AP009048/3782777‑3782731
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atcaccATCCGCTCTATTTTGCCTGTAATGACTTTGTTGATGACCAAATCGGACGGGTCATCAATGCCTTAAC  >  AP009048/3782753‑3782819

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