New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 2535433 =46 (0.820)3 (0.060) 3/88 NT 7.6% coding (190/2016 nt) ligA DNA ligase, NAD(+)‑dependent
?AP009048 2535454 = 33 (0.680)coding (169/2016 nt) ligA DNA ligase, NAD(+)‑dependent

TCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTGGATGTTTGGTTTCCAGCTCGCGCAGTTCG  >  AP009048/2535363‑2535503
                                                                      |                                                                      
tcaAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAgc                                                                        <  1:4072582/71‑1 (MQ=255)
tcaAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAgc                                                                        <  1:1570358/71‑1 (MQ=255)
   tacgacaggctgatgcgcgaactgcgcgagctggaaaccaaacatccagaactgattacgcctgattCGcc                                                                     <  1:3282746‑M1/4‑1 (MQ=255)
    aaaCGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAgcgc                                                                      >  1:266697/1‑69 (MQ=255)
     cgacaggctgatgcgcgaactgcgcgagctggaaaccaaacatccagaactgattacgcctgattCGCCTa                                                                   <  1:2396879‑M1/6‑1 (MQ=255)
     aaCGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTa                                                                   >  1:1737877/1‑71 (MQ=255)
        caggctgatgcgcgaactgcgcgagctggaaaccaaacatccagaactgattacgcctgattCGCCTacac                                                                <  1:315777‑M1/9‑1 (MQ=255)
         ttATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGcc                                                                     <  1:2104041/65‑1 (MQ=255)
         ttATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACg                                                               <  1:3371729/71‑1 (MQ=255)
          tATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGt                                                              >  1:807842/1‑71 (MQ=255)
           aTCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGt                                                              <  1:4593892/70‑1 (MQ=255)
            ctgatgcgcgaactgcgcgagctggaaaccaaacatccagaactgattacgcctgattCGCCTACACGTTg                                                            >  1:2820352‑M1/59‑71 (MQ=255)
              cAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGtt                                                             <  1:805488/68‑1 (MQ=255)
              cAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAg                                                          >  1:1423525/1‑71 (MQ=255)
              cAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAg                                                          <  1:114661/71‑1 (MQ=255)
              cAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAg                                                          >  1:1111034/1‑71 (MQ=255)
                tgcgcgaactgcgcgagctggaaaccaaacatccagaactgattacgcctgattCGc                                                                      <  1:1955934‑M1/3‑1 (MQ=255)
                gTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAg                                                          >  1:4296671/1‑69 (MQ=255)
                gTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTa                                                        <  1:4171560/71‑1 (MQ=255)
                  gACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAgg                                                      >  1:1061155/1‑71 (MQ=255)
                   aCAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTa                                                        <  1:170702/68‑1 (MQ=255)
                   aCAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGc                                                     >  1:1837660/1‑71 (MQ=255)
                       actgcgcgagctggaaaccaaacatccagaactgattacgcctgattc                                                                        <  1:3088026‑M1/1‑1 (MQ=255)
                        aTTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAAt                                                 >  1:3989839/1‑70 (MQ=255)
                          tGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAg                                              >  1:582127/1‑71 (MQ=255)
                           ggTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAg                                              <  1:2513707/70‑1 (MQ=255)
                           ggTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAg                                              <  1:884780/70‑1 (MQ=255)
                              aCTTCATGGCGTATCTGGCTGAAAGCCGCAAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCg                                           >  1:2809812/1‑70 (MQ=255)
                               cTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATc                                                <  1:870278/64‑1 (MQ=255)
                               cTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGt                                          <  1:4809260/70‑1 (MQ=255)
                                ttCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAAt                                                 >  1:1163797/1‑62 (MQ=255)
                                 tCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAAt                                       >  1:2601293/1‑71 (MQ=255)
                                 tCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAAt                                       >  1:2622122/1‑71 (MQ=255)
                                  cATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAAt                                       >  1:43401/1‑70 (MQ=255)
                                    tGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCa                                     >  1:377246/1‑70 (MQ=255)
                                     ggCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAgg                                                      >  1:2722447/1‑52 (MQ=255)
                                     aaaccaaacatacagaactgattacgcctgattCGCCTACACGTTGAGTa                                                        >  1:48902‑M1/34‑50 (MQ=255)
                                       cGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTc                                 <  1:2101886/71‑1 (MQ=255)
                                        ccaaacatccagaactgattacgcctgattCGCCTACACGTTGAGTAGGc                                                     >  1:1152184‑M1/31‑50 (MQ=255)
                                         tATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGtt                                  <  1:805283/68‑1 (MQ=255)
                                         tATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTg                               <  1:3932971/71‑1 (MQ=255)
                                           tCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTg                                                            <  1:4827360/40‑1 (MQ=255)
                                           tCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAg                                              >  1:4574911/1‑54 (MQ=255)
                                              ggCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTGGATGt                          >  1:2052367/1‑71 (MQ=255)
                                               gCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTGGATGt                          <  1:858315/70‑1 (MQ=255)
                                                   aaaGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCa                                               <  1:690592/45‑1 (MQ=255)
                                                   aaaGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCa                                               <  1:1075487/45‑1 (MQ=255)
                                                   aaaGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTGGATGTTTGGt                     <  1:3724536/71‑1 (MQ=255)
                                                   aaaGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTGGATGTTTGGt                     >  1:5151108/1‑71 (MQ=255)
                                                   aaaGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTGGATGTTTGGt                     >  1:3100892/1‑71 (MQ=255)
                                                        cgccAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTGGATGTTTGGTTTc                  >  1:3633868/1‑69 (MQ=255)
                                                        cgccAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTGGATGTTTGGTTTc                  >  1:2502742/1‑69 (MQ=255)
                                                          ccAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTGGATGTTTGGTTTCCAGc              >  1:3158120/1‑71 (MQ=255)
                                                           cAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTGGATGTTTGGTTTCCAGCt             >  1:4733550/1‑71 (MQ=255)
                                                            aGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTGGATGTTTGGTTTCCAGCTc            >  1:189555/1‑71 (MQ=255)
                                                              cggcgCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTGGATGTTTGGTTTCCAGCTcgc          <  1:2353147/71‑1 (MQ=255)
                                                              cggcgCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTGGATGTTTGGTTTCCAGCTcgc          >  1:670179/1‑71 (MQ=255)
                                                              cggcgCAGCGCCTACACGTTGAGGAGGCGAATCAGGCGTAATCAGTTCTgg                              <  1:366316/51‑1 (MQ=255)
                                                                 cgcAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTGGATGTTTGGTTTCCAGCTCGCGCa       <  1:4849693/71‑1 (MQ=255)
                                                                  gcAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTGGATGTTTGGTTTCCAGCTCGCGCAg      >  1:334232/1‑71 (MQ=255)
                                                                   cAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCt                                >  1:1710214/1‑44 (MQ=255)
                                                                   cAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTGGATGt                          <  1:428348/50‑1 (MQ=255)
                                                                   cAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTGGATGTTTGGTTTCCAGCTCGCGCAGt     <  1:3369794/71‑1 (MQ=255)
                                                                      CGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTGGATGTTTGGTTTCCAGCTCGCGCAGTTcg  <  1:3507927‑M1/71‑1 (MQ=255)
                                                                      |                                                                      
TCAAAAACGTTATCCAGTGACAGCATTGGTACTTCATGGCGTATCTGGCTGAAAGCCGCCAGCGGCGCAGCGCCTACACGTTGAGTAGGCGAATCAGGCGTAATCAGTTCTGGATGTTTGGTTTCCAGCTCGCGCAGTTCG  >  AP009048/2535363‑2535503

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.