New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 460980943 (1.260)3 (0.090) 3/158 2.7 6.9% coding (390/1224 nt) deoB phosphopentomutase
?NZ_CP009273 4609882 = 41 (1.290)coding (463/1224 nt) deoB phosphopentomutase
Rejected: Frequency below/above cutoff threshold.

GGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGC  >  NZ_CP009273/4609720‑4609897
                                                                                         |                                                                                        
ggTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGAT                                                                                          <  1:919487‑M1/90‑1 (MQ=255)
  tCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATaa                                                                                        >  1:942030/1‑90 (MQ=255)
  tCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATaa                                                                                        >  1:771808/1‑90 (MQ=255)
     cTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACt                                                                                     >  2:662129/1‑90 (MQ=255)
     cTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACt                                                                                     >  1:864687/1‑90 (MQ=255)
          aaaTTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCg                                                                                >  1:257480/1‑90 (MQ=255)
            aTTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGaa                                                                              <  2:231135/90‑1 (MQ=255)
             ttGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAAc                                                                             <  2:743367/90‑1 (MQ=255)
               gCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAAcgc                                                                           <  1:260837/90‑1 (MQ=255)
               gCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAAcgc                                                                           >  1:262109/1‑90 (MQ=255)
                       cccGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCt                                                                   >  1:212820/1‑90 (MQ=255)
                         cGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGc                                                                 <  1:339919/90‑1 (MQ=255)
                         cGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGc                                                                 <  2:178507/90‑1 (MQ=255)
                              cTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGTTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGt                                                            <  2:127564/90‑1 (MQ=255)
                              cTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGt                                                            <  1:746492/90‑1 (MQ=255)
                               tGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggca                                                           >  1:916405‑M1/1‑59 (MQ=255)
                               tGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAAcgc                                                                           <  2:562453/74‑1 (MQ=255)
                               tGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAAcgc                                                                           >  1:562453/1‑74 (MQ=255)
                                        ggggATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTa                                                  <  2:927565/90‑1 (MQ=255)
                                          ggATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGTGTTACCTCGGTAAc                                                >  2:626738/1‑90 (MQ=255)
                                          ggATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAc                                                >  1:670203/1‑90 (MQ=255)
                                           gATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACt                                               <  2:951772/90‑1 (MQ=255)
                                              atTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGcc                                            >  1:318599/1‑90 (MQ=255)
                                              atTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGcc                                            <  2:86485/90‑1 (MQ=255)
                                              atTTCTCCAATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGcc                                            >  2:136709/1‑90 (MQ=255)
                                               tttCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa                                           >  1:467531/1‑90 (MQ=255)
                                                ttCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCAc                                          <  2:937267/90‑1 (MQ=255)
                                                ttCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCAc                                          <  2:262109/90‑1 (MQ=255)
                                                  ctcCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCActc                                        >  1:715803/1‑90 (MQ=255)
                                                    ccGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCtt                                      <  2:917353/90‑1 (MQ=255)
                                                     cGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTc                                     <  2:624358/90‑1 (MQ=255)
                                                      gATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTcc                                    >  2:471254/1‑90 (MQ=255)
                                                         cACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGt                                 <  2:676734/90‑1 (MQ=255)
                                                                  aGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTCCCTCGGTAACTGCCACTCTTCCGGTACGGTCAtt                        <  1:958592/90‑1 (MQ=255)
                                                                  aGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCAtt                        <  2:151438/90‑1 (MQ=255)
                                                                  aGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCAtt                        <  1:865224/90‑1 (MQ=255)
                                                                   gCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGt                                                            <  2:885344/53‑1 (MQ=255)
                                                                   gCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGt                                                            >  1:885344/1‑53 (MQ=255)
                                                                     ttCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgat                                                    >  2:328536‑M1/1‑21 (MQ=255)
                                                                     ttCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgat                                                    <  1:328536‑M1/59‑39 (MQ=255)
                                                                     ttCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGtt                                                           >  2:527480/1‑52 (MQ=255)
                                                                     ttCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGtt                                                           <  1:527480/52‑1 (MQ=255)
                                                                       cccGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGa                   <  2:90121/90‑1 (MQ=255)
                                                                           cAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCaa               >  2:790085/1‑90 (MQ=255)
                                                                                      ggATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctgactccgtgttccagattgcctgccatgaaga    >  2:745961‑M1/1‑4 (MQ=255)
                                                                                        aTAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGc  <  2:393605/90‑1 (MQ=255)
                                                                                        aTAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGc  <  1:328932/90‑1 (MQ=255)
                                                                                        aTAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGc  <  2:159488/90‑1 (MQ=255)
                                                                                         |                                                                                        
GGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGC  >  NZ_CP009273/4609720‑4609897

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