New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 46072643 (0.070)49 (1.300) 29/160 0.1 85.4% coding (125/780 nt) deoC deoxyribose‑phosphate aldolase
?NZ_CP009273 4608011 = 14 (0.370)intergenic (+92/‑35) deoC/deoA deoxyribose‑phosphate aldolase/thymidine phosphorylase

GAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCG  >  NZ_CP009273/4607175‑4607350
                                                                                         |                                                                                      
gAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAA                                                                                        >  1:594078‑M1/1‑90 (MQ=255)
    ttGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagc                                                                                    <  2:309114‑M1/90‑5 (MQ=255)
    ttGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagc                                                                                    <  1:548592‑M1/90‑5 (MQ=255)
     tGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagcc                                                                                   >  2:97874‑M1/1‑85 (MQ=255)
     tGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagcc                                                                                   <  2:520874‑M1/90‑6 (MQ=255)
      gATGGACCTGACCACCCTGAATGACGACGACACCGATGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccg                                                                                  <  1:930005‑M1/90‑7 (MQ=255)
      gATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccg                                                                                  <  1:53130‑M1/90‑7 (MQ=255)
         ggcccTGACCCCCCTGAATGACGACGACCCCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACcgccgc                                                                               <  2:1077287/87‑1 (MQ=255)
            ccTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcat                                                                            <  1:360799‑M1/90‑13 (MQ=255)
              tGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatac                                                                          <  1:1027014‑M1/90‑15 (MQ=255)
               gACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatact                                                                         >  2:351589‑M1/1‑75 (MQ=255)
                    cccTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatactctttt                                                                    >  1:1098825‑M1/1‑70 (MQ=255)
                    cccTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatactctttt                                                                    >  2:1091705‑M1/1‑70 (MQ=255)
                    cccTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatactctttt                                                                    >  1:1059630‑M1/1‑70 (MQ=255)
                    cccTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatactctttt                                                                    >  1:357369‑M1/1‑70 (MQ=255)
                    cccTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatactctttt                                                                    <  1:476498‑M1/90‑21 (MQ=255)
                        gAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctc                                                                >  1:102304‑M1/1‑66 (MQ=255)
                          aTGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgg                                                              <  2:669232‑M1/90‑27 (MQ=255)
                          aTGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgg                                                              <  1:814044‑M1/90‑27 (MQ=255)
                              cgacgacACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggagg                                                          >  2:295508‑M1/1‑60 (MQ=255)
                              cgacgacACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggagg                                                          >  1:528897‑M1/1‑60 (MQ=255)
                                      ccGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGa                                                  >  1:666563/1‑90 (MQ=255)
                                           gagaAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCt                                             <  2:52167/90‑1 (MQ=255)
                                             gaAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcg                                           <  1:301597‑M1/90‑46 (MQ=255)
                                              aaaGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgc                                          <  2:317391‑M1/90‑47 (MQ=255)
                                                   gATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcacaag                                     >  2:748811‑M1/1‑39 (MQ=255)
                                                      cGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttt                                                                     >  1:586095‑M1/1‑36 (MQ=255)
                                                      cGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttt                                                                     <  2:586095‑M1/55‑20 (MQ=255)
                                                           tgtgTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcacaagaaattatt                             >  1:493248‑M1/1‑31 (MQ=255)
                                                           tgtgTCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcacaagaaattatt                             >  2:49058‑M1/1‑31 (MQ=255)
                                                            gtgtCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcacaagaaattattc                            <  2:772030‑M1/90‑61 (MQ=255)
                                                            gtgtCATCAGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcacaagaaattattc                            <  2:357369‑M1/90‑61 (MQ=255)
                                                                 atcaGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcacaagaaattattcgtaaa                       >  1:1015387‑M1/1‑25 (MQ=255)
                                                                 atcaGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgg                                                              >  1:162008‑M1/1‑25 (MQ=255)
                                                                 atcaGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgg                                                              <  2:162008‑M1/51‑27 (MQ=255)
                                                                  tcaGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcaca                                       <  1:734533‑M1/73‑50 (MQ=255)
                                                                  tcaGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcaca                                       >  2:734533‑M1/1‑24 (MQ=255)
                                                                  tcaGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggt                                                         >  1:209085‑M1/1‑24 (MQ=255)
                                                                  tcaGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggt                                                         <  2:209085‑M1/55‑32 (MQ=255)
                                                                   caGGCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcacaagaaattattcgtaaaaa                     <  2:307856‑M1/90‑68 (MQ=255)
                                                                      gCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcacaagaaattattcgtaaaaaacg                  >  2:1061655‑M1/1‑20 (MQ=255)
                                                                      gCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgc                                          >  1:127279‑M1/1‑20 (MQ=255)
                                                                      gCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgc                                          <  2:127279‑M1/66‑47 (MQ=255)
                                                                      gCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgc                                          <  2:1053971‑M1/66‑47 (MQ=255)
                                                                      gCCAAAACTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgc                                          >  1:1053971‑M1/1‑20 (MQ=255)
                                                                      gCCAAAACTCCGGTCGGCAAaagacgcctcatactcttttcctcgggaggttaccttgtttctcgcacaagaaattattcgtaaaaaacg                  >  1:300368‑M1/1‑20 (MQ=255)
                                                                          aaaCTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcacaagaaattattcgtaaaagacgtgat              <  1:20520‑M1/90‑75 (MQ=255)
                                                                          aaaCTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcacaagaaattattcgtaaaaaacgtgat              <  1:1061655‑M1/90‑75 (MQ=255)
                                                                          aaaCTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcacaagaaattattcgtaaaaaacgtgat              <  2:132490‑M1/90‑75 (MQ=255)
                                                                          aaaCTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcacaagaaattattcgtaaaaaacgtgat              >  2:143778‑M1/1‑16 (MQ=255)
                                                                          aaaCTCCGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcacaagaaattattcgtaaaaaacgtgat              <  1:972129‑M1/90‑75 (MQ=255)
                                                                               ccGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcacaagaaattattcgtaaaaaacgtgatggtca         >  1:1053827‑M1/1‑11 (MQ=255)
                                                                               ccGGTCGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcacaagaaattattcgtaaaaaacgtgatggtca         <  2:1097138‑M1/90‑80 (MQ=255)
                                                                                    cGGCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcacaagaaattattcgtaaaaaacgtgatggtcatgcgc    <  2:102304‑M1/90‑85 (MQ=255)
                                                                                      gCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcacaagaaattattcgtaaaaaacgtgatggtcatgcgctg  <  2:670713‑M1/90‑87 (MQ=255)
                                                                                      gCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcacaagaaattattcgtaaaaaacgtgatggtcatgcgctg  <  1:687971‑M1/90‑87 (MQ=255)
                                                                                      gCAAaagccgcctcatactcttttcctcgggaggttaccttgtttctcgcacaagaaattattcgtaaaaaacgtgatggtcatgcgctg  <  1:496690‑M1/90‑87 (MQ=255)
                                                                                         |                                                                                      
GAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCG  >  NZ_CP009273/4607175‑4607350

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