New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 = 405918671 (1.260)3 (0.060) 3/94 NT 4.2% coding (92/489 nt) yhhY predicted acetyltransferase
?AP009048 = 4059197 71 (1.370)coding (81/489 nt) yhhY predicted acetyltransferase
Rejected: Frequency below/above cutoff threshold.

GACGAGTTGCTTGATGCCGGGACGATCGGCGAGTCGCTCCTGCCACATATGATCGGAAGGATGAGGCACCTGTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AP009048/4059113‑4059186
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgtaTTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGG  <  AP009048/4059197‑4059129
                                                                                                                                               
GACGAGTTGCTTGATGCCGGGACGATCGGCGAGTCGCTCCTGCCACATATGATCGGAAGGATGAGGCACCT                                                                          >  1:4456334/1‑71
   GAGTTGCTTGATGCCGGGACGATCGGCGAGTCGCTCCTGCCACATATGATCGGAAGGATGAGGCACCTGTA                                                                       <  1:2563420/71‑1
   GAGTTGCTTGATGCCGGGACGATCGGCGAGTCGCTCCTGCCACATATGATCGGAAGGATGAGGCACCTGTA                                                                       <  1:3370167/71‑1
   GAGTTGCTTGATGCCGGGACGATCGGCGAGTCGCTCCTGCCACATATGATCGGAAGGATGAGGCACCTG                                                                         >  1:862946/1‑69
                          CGGCGAGTCGCTCCTGCCACATATGATCGGAAGGATGAGGCACCTGTATTGCAACACACTA                                                          <  1:2218589/61‑1
                                    CTCCTGCCACATATGATCGGAAGGATGAGGCACCTG                                                                         <  1:4191785/36‑1
                                      CCTGCCACATATGATCGGAAGGATGAGGCACCTGTATTGCAACACACTACA                                                        >  1:4086122/1‑51
                                              ATATGATCGGAAGGATGAGGCACCTGTATTGCAACACACTACAG                                                       <  1:3195630/44‑1
                                                                      TGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCC    >  1:3119415/1‑71
                                                                        TATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGG  <  1:4397352/71‑1
                                                                         ATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGG  >  1:1997744/1‑70
                                                                                                                                               
GACGAGTTGCTTGATGCCGGGACGATCGGCGAGTCGCTCCTGCCACATATGATCGGAAGGATGAGGCACCTGTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AP009048/4059113‑4059186
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgtaTTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGG  <  AP009048/4059197‑4059129

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.