New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 462550643 (0.850)3 (0.060) 3/242 NT 6.6% intergenic (+4099/+465) PP_5654/PP_4091 phosphotriesterase/transposase
?NC_002947 = 4625534 45 (0.950)intergenic (+4127/+437) PP_5654/PP_4091 phosphotriesterase/transposase

GACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTCCTGCGCACCCA  >  NC_002947/4625362‑4625654
                                                                                                                                                 |                                                                                                                                                    
gacgCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCACCAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCaa                                                                                                                                                  <  7:139722/149‑1 (MQ=255)
   gcgGCCGACAATCTGCTGTCCATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGAACCGCTACGACGCGCTGGGCCAttt                                                                                                                                               >  8:182811/1‑149 (MQ=255)
    cgGCCGACAATCTGCTGTCCATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGTAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTg                                                                                                                                              <  1:258651/149‑1 (MQ=255)
      gCCGACAATCTGCTGTCCATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCa                                                                                                                                            <  4:110778/149‑1 (MQ=255)
              tCTGCTGTCCATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTgg                                                                                                                                    >  3:34902/1‑149 (MQ=255)
                  ctgTCCATCGCCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACACCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCGGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGtgct                                                                                                                                >  8:126213/1‑149 (MQ=255)
                       cATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGAccc                                                                                                                                       >  7:230455/1‑137 (MQ=255)
                       cATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGAccc                                                                                                                                       <  8:230455/137‑1 (MQ=255)
                        aTCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGCCCCTGGTGCTGCCCgc                                                                                                                          >  8:150871/1‑148 (MQ=255)
                         tCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGAc                                                                                                                         >  1:18784‑M2/1‑148 (MQ=255)
                             cTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGc                                                                                                                                                    >  8:65236/1‑118 (MQ=255)
                             cTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGc                                                                                                                                                    <  7:65236/118‑1 (MQ=255)
                                cacaGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCGGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGAccagcgc                                                                                                                   >  7:126213‑M2/1‑141 (MQ=255)
                                cacaGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGAccagcgc                                                                                                                   >  5:123746‑M2/1‑141 (MQ=255)
                                  caGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGCAACCAACAGCGCGGGCCTGCTGCAGTACAGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGaa                                                                                                                >  2:284383/1‑149 (MQ=255)
                                       aaCCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCAGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTc                                                                                                             >  1:254024/1‑148 (MQ=255)
                                       aaCCAGGGCAAGAAGCAGCAACTCGATTACACATACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCaa                                                                                                           >  4:197938/1‑149 (MQ=255)
                                          cAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGtcgggcagcaccagggtctgcaaattgcccagcgc                                                                                                                   <  8:112312‑M1/138‑36 (MQ=255)
                                          cAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGtcgggcagcaccagggtctgcaaattgcccagcgc                                                                                                                   >  7:112312‑M1/1‑103 (MQ=255)
                                           aGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCAc                                                                                                       <  8:124238/149‑1 (MQ=255)
                                            gggCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCAcc                                                                                                      >  4:25206/1‑149 (MQ=255)
                                             ggCAAGATCCATCAAATCGATTACTCATACGACCCCAACGGGCAACTGCTCAGCGAAACCTACACTGCGGGCCTGCTGCTGTACCGCTACGACGCGCTGGtcgggcagcaccagggtctgcaaattgcccagcgcgtcgtagcggtact                                                                                                     <  8:115462‑M1/149‑50 (MQ=255)
                                               cAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTg                                                                                                    >  3:97901/1‑148 (MQ=255)
                                                    agcagcAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTGCCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGtactac                                                                                              <  7:186395/149‑1 (MQ=255)
                                                      cagcaACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACgg                                                                                            <  3:97711/149‑1 (MQ=255)
                                                               gATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGaca                                                                                   <  3:25206/149‑1 (MQ=255)
                                                                   acacCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTg                                                                               >  1:28577/1‑149 (MQ=255)
                                                                     acCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCa                                                                             >  4:135424/1‑149 (MQ=255)
                                                                         aCGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCAc                                                                            >  4:295689/1‑146 (MQ=255)
                                                                          cGACCCCACCGGCCACCTGCTCCGCGACACCCACAGCGCGGGCCTGCTGCACTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGa                                                                        <  8:240387/149‑1 (MQ=17)
                                                                          cGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGa                                                                        <  4:97901/149‑1 (MQ=255)
                                                                                          cTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCg                                                        >  2:16958/1‑149 (MQ=255)
                                                                                             cTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCggg                                                      <  2:254024/148‑1 (MQ=255)
                                                                                              tCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGt                                                     >  2:156894/1‑148 (MQ=255)
                                                                                                          aaCAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGc                                                          <  6:331844/132‑1 (MQ=255)
                                                                                                          aaCAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGc                                                          >  5:331844/1‑132 (MQ=255)
                                                                                                           aCAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTc                                        <  1:156894/149‑1 (MQ=255)
                                                                                                           aCAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTc                                        <  1:16958/149‑1 (MQ=255)
                                                                                                             aGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGa                                      <  1:266280/149‑1 (MQ=255)
                                                                                                               cgcgGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAgc                                    >  6:272246/1‑149 (MQ=255)
                                                                                                                         ctgcAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGt                          >  6:282125/1‑149 (MQ=255)
                                                                                                                           gcAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCAc                                                                            >  6:188784/1‑97 (MQ=255)
                                                                                                                           gcAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCAc                                                                            <  5:188784/97‑1 (MQ=255)
                                                                                                                           gcAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTg                         >  1:177600/1‑148 (MQ=255)
                                                                                                                            cAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTg                         <  3:295689/147‑1 (MQ=255)
                                                                                                                             aGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGAccagcgcgtcgtagcggtacagcagcaggcccgcagtgtaggtttcgctgagcagttgcccgttggggtcgtatgagtaatcgatttgatggatcttgcc                      <  8:115462‑M2/149‑101 (MQ=255)
                                                                                                                                  cGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCAcgacga                 <  3:135424/149‑1 (MQ=255)
                                                                                                                                    cTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAgg               <  7:189810/149‑1 (MQ=255)
                                                                                                                                      acgacgCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGAt                                           <  8:83941/119‑1 (MQ=255)
                                                                                                                                      acgacgCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGAt                                           >  7:83941/1‑119 (MQ=255)
                                                                                                                                       cgacgCGCGGGGCAATTTGCAGACCCTGCGGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTc             <  1:284383/148‑1 (MQ=255)
                                                                                                                                          cgcgCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTCCt           <  2:91552/147‑1 (MQ=255)
                                                                                                                                          cgcgCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTCCt           >  1:91552/1‑147 (MQ=255)
                                                                                                                                           gcgcTGGtcgggcagcaccagggtctgcaaattgcccagcgcgtcgtagcggtactgcagcaggcccgcgctgttggtttcgctgagcagttgcccgttggggtcgtaggtgtaatcgagttgctgcttcttgccctggttgtctgtg         >  5:123746‑M1/1‑7 (MQ=255)
                                                                                                                                           gcgcTGGtcgggcagcaccagggtctgcaaattgcccagcgcgtcgtagcggtaccgcagcaggcccgcgctgttggtttcgctgagcagttgcccgttggggtcgtaggtgtaatcgagttgctgcttcttgccctggttgtctgtg         >  7:126213‑M1/1‑7 (MQ=255)
                                                                                                                                           gcgcTGGGCAATTTGCAGACCCTGGTGCTGCCCGAccagcgcgtcgtagcggtactgcagcaggcccgcgctgttggtttcgctgagcagttgcccgttggggtcgtaggtgtaatcgagttgctgcttcttgccctg                   <  8:112312‑M2/138‑104 (MQ=255)
                                                                                                                                           gcgcTGGGCAATTTGCAGACCCTGGTGCTGCCCGAccagcgcgtcgtagcggtactgcagcaggcccgcgctgttggtttcgctgagcagttgcccgttggggtcgtaggtgtaatcgagttgctgcttcttgccctg                   >  7:112312‑M2/1‑35 (MQ=255)
                                                                                                                                               tGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAg                                                                         <  4:170035/80‑1 (MQ=255)
                                                                                                                                               tGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAg                                                                         >  3:170035/1‑80 (MQ=255)
                                                                                                                                                 gtcgggcagcaccagggtctgcaaattgcccagcgcgtcgtagcggtactgcagcaggcccgcgctgttggtttcgctgagcagttgcccgttggggtcgtaggtgtaatcgagttgctgcttcttgccctggttgtctgtgaaggtga  >  1:18784‑M1/1‑1 (MQ=255)
                                                                                                                                                 ggCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTCCTGCGCAccc   >  5:282462/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTCCTGCGCACCCA  >  NC_002947/4625362‑4625654

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