New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 = 203339844 (0.750)3 (0.060) 3/92 NT 7.2% coding (1057/1419 nt) dcm DNA cytosine methylase
?AP009048 = 2033414 38 (0.720)coding (1041/1419 nt) dcm DNA cytosine methylase

CATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGT  >  AP009048/2033344‑2033482
                                                                      |                                                                    
cATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATaa                                                                      >  1:1785953/1‑71 (MQ=255)
 aTCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAAc                                                                     <  1:4875977‑M2/71‑2 (MQ=255)
 aTCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAAc                                                                     <  1:2277160‑M2/71‑2 (MQ=255)
 aTCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAAc                                                                     <  1:2283295‑M2/71‑2 (MQ=255)
   cTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAAcc                                                                    <  1:2645015‑M2/70‑3 (MQ=255)
      tGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAAccg                                                                   >  1:263484‑M2/1‑65 (MQ=255)
       gTAATAACGCGCAGACAGCTTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAAccg                                                                   >  1:4090979‑M2/1‑64 (MQ=255)
       gTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGc                                                               >  1:1212567/1‑71 (MQ=255)
           taaCGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGt                                                            >  1:3695268/1‑70 (MQ=255)
             aCGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTa                                                          <  1:4168105/70‑1 (MQ=255)
              cgcgCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTAcc                                                        <  1:4226135/71‑1 (MQ=255)
               gcgcAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTa                                                          >  1:3485494/1‑68 (MQ=255)
                 gcAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACcgcg                                                     >  1:667192/1‑71 (MQ=255)
                  cAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACcgcgc                                                    <  1:4047234/71‑1 (MQ=255)
                    gACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACcgcgcg                                                   >  1:2690089/1‑70 (MQ=255)
                      cAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCt                                                <  1:4994625/71‑1 (MQ=255)
                         cGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGAt                                             >  1:3194368/1‑71 (MQ=255)
                           tACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTa                                                          <  1:439208/56‑1 (MQ=255)
                            aCGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAAccgaa                                                                 >  1:2134548‑M2/1‑43 (MQ=255)
                             cgcgTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTAcc                                                        >  1:2882733/1‑56 (MQ=255)
                              gcgTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGcc                                                              <  1:1804185/49‑1 (MQ=255)
                              gcgTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCt                                                <  1:3609582/63‑1 (MQ=255)
                              gcgTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGTTGtttt                                        <  1:3945543/71‑1 (MQ=255)
                              gcgTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGtttt                                        <  1:2153250/71‑1 (MQ=255)
                              gcgTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGtttt                                        <  1:1463279/71‑1 (MQ=255)
                               cgTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGttttt                                       <  1:3595254/71‑1 (MQ=255)
                               cgTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGttttt                                       <  1:1053153/71‑1 (MQ=255)
                                 tGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCg                                                             <  1:4146290/47‑1 (MQ=255)
                                  gACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTg                                               >  1:3101301/1‑60 (MQ=255)
                                  gACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGtttttt                                      >  1:2747786/1‑69 (MQ=255)
                                   aCGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGcc                                                              >  1:2359413/1‑44 (MQ=255)
                                   aCGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACcgc                                                      >  1:3420456/1‑52 (MQ=255)
                                     gCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGttt                                         >  1:4767319/1‑63 (MQ=255)
                                       tttGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACcg                                                       >  1:1018180/1‑47 (MQ=255)
                                       tttGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCatat                               >  1:934599/1‑71 (MQ=255)
                                       tttGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCata                                <  1:4700116/70‑1 (MQ=255)
                                         tGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATc                              >  1:4533892/1‑70 (MQ=255)
                                         tGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATc                              <  1:4275740/70‑1 (MQ=255)
                                          gCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCatat                               >  1:173197/1‑68 (MQ=255)
                                          gCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGa                            >  1:3863248/1‑71 (MQ=255)
                                             gATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATag                         <  1:1086540/71‑1 (MQ=255)
                                                  ttCGGATAAACCATTCCATAAccgaacaatccgcaaagcgtcacgcgtacgctgtctgcgcgttat                         >  1:2218021‑M2/1‑21 (MQ=255)
                                                  ttCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGt                     <  1:2433370/70‑1 (MQ=255)
                                                    cGGATAAACCATTCCATAAccgaacaatccgcaaagcgtcacgcgt                                           <  1:4388668‑M2/46‑28 (MQ=255)
                                                    cGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGttt                                         >  1:4485172/1‑48 (MQ=255)
                                                     ggATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACt                  <  1:5068581/70‑1 (MQ=255)
                                                      gATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATag                         >  1:2101360/1‑62 (MQ=255)
                                                      gATAAAACATTCCATAAccgaacaatccgcaaagcgtcacgcgtacgctgtctgcgcgttattacaaagat                <  1:685635‑M2/71‑55 (MQ=255)
                                                          aaCCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACTTCcaca            <  1:4751859/71‑1 (MQ=255)
                                                               ttCCATAACCGAAGCCGTTGCCGCGCGCCTGATGttttt                                       >  1:3491484/1‑39 (MQ=255)
                                                                 ccATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGa                            >  1:349369/1‑48 (MQ=255)
                                                                 ccATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGa                            >  1:4770501/1‑48 (MQ=255)
                                                                 ccATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCgg     >  1:3454860/1‑71 (MQ=255)
                                                                 ccAAAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCg      >  1:3726615/1‑70 (MQ=255)
                                                                  cATAACCGAAGCCGTTACCGCGCGCCTGATGtttttt                                      >  1:4376828/1‑37 (MQ=255)
                                                                   aTAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGt                     <  1:311690/53‑1 (MQ=255)
                                                                    tAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGt  >  1:3550157/1‑71 (MQ=255)
                                                                      |                                                                    
CATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGT  >  AP009048/2033344‑2033482

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.