New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 = 42784156 (1.000)3 (0.060) 3/90 NT 5.5% coding (971/1848 nt) secD SecYEG protein translocase auxillary subunit
?AP009048 = 427846 54 (1.090)coding (976/1848 nt) secD SecYEG protein translocase auxillary subunit

CGGGCGACTCTGAAGTGAAACAGACCCGCGAAGGTCAGCCAGTTGTGCTGTACAAACGCGTAATTCTGACCGGTGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AP009048/427766‑427841
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggtgaTATGGTCACCGGTCAGAATTACGCGTTTGTACAGCACAACTGGCTGACCTTCGCGGGTCTGTTTCACTT  <  AP009048/427846‑427778
                                                                                                                                                 
CGGGCGACTCTGAAGTGAAACAGACCCGCGAAGGTCAGCCAGTTGTGCTGTACAAACGCGTAATTCTGACC                                                                            >  1:279720/1‑71
 GGGCGACTCTGAAGTGAAACAGACCCGCGAAGGTCAGCCAGTTGTGCTGTACAAACGCGTAATTCTGACCG                                                                           <  1:4350094/71‑1
   GCGACTCTGAAGTGAAACAGACCCGCGAAGGTCAGCCAGTTGTGCTGTACAAACGCGTAATTCTGACCGGT                                                                         <  1:2777087/71‑1
    CGACTCTGAAGTGAAACAGACCCGCGAAGGTCAGCCAGTTGTGCTGTACAAACGCGTAATTCTGACCGGTG                                                                        <  1:1013538/71‑1
    CGACTCTGAAGTGAAACAGACCCGCGAAGGTCAGCCAGTTGTGCTGTACAAACGCGTAATTCTGACCGGT                                                                         >  1:1579272/1‑70
                 AAACAGACCCGCGAAGGTCAGCCAGTTGTGCTGTACAAACGCGTAATTCTGACCGGTGATATGGTCACCGG                                                           <  1:4079482/71‑1
                   ACAGACCCGCGAAGGTCAGCCAGTTGTGCTGTACAAACGCGTAATTCTGACCG                                                                           <  1:3191282/53‑1
                                            GTGCTGTACAAACGCGTAATTCTGACCGGTGATATGGTCACCG                                                            >  1:65399/1‑43
                                             TGCTGTACAAACGCGTAATTCTGACCGGTGATATGGTCACCG                                                            <  1:4886482/42‑1
                                                                      CGGTGATATGGTCACCGGTCAGAATTACGCGTTTGTACAGCACAACTGGCTGACCTTCGCGGGTCTGTTT       >  1:5125110/1‑70
                                                                       GGTGATATGGTCACCGGTCAGAATTACGCGTTTGTACAGCACAACTGGCTGACCTTCGCGGGTCTGTTTCA     >  1:682859/1‑71
                                                                           ATATGGTCACCGGTCAGAATTACGCGTTTGTACAGCACAACTGGCTGACCTTCGCGGGTCTGTTTCACTT  >  1:1109373/1‑70
                                                                                                                                                 
CGGGCGACTCTGAAGTGAAACAGACCCGCGAAGGTCAGCCAGTTGTGCTGTACAAACGCGTAATTCTGACCGGTGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AP009048/427766‑427841
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggtgaTATGGTCACCGGTCAGAATTACGCGTTTGTACAGCACAACTGGCTGACCTTCGCGGGTCTGTTTCACTT  <  AP009048/427846‑427778

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.