New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 4625552	43 (0.820)	4 (0.080)	3/230	NT	8.2%	intergenic (+4099/+465)	PP_5654/PP_4091	phosphotriesterase/transposase
	?	NC_002947	= 4625580	49 (1.000)					intergenic (+4127/+437)	PP_5654/PP_4091	phosphotriesterase/transposase

GCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTCCTGCGCACCCAA  >  NC_002947/4625403‑4625701                                                                                                                                                      |                                                                                                                                                      gCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGG                                                                                                                                                       <  2:78023‑M1/150‑1 (MQ=255) gCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGG                                                                                                                                                       <  1:1285081‑M1/150‑1 (MQ=255)     ccacgCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCaa                                                                                                                                                   <  1:143679/148‑1 (MQ=255)          cgGCCGACAATCTGCTGTCCATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGc                                                                                                                                              <  1:740477/150‑1 (MQ=255)                    tCTGCTGTCCATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTTGGCAATTTGCAGACCCTGGt                                                                                                                                    <  2:164443/150‑1 (MQ=255)                            ccATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCg                                                                                                                            >  2:1037859/1‑150 (MQ=255)                               tCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGAcc                                                                                                                         <  1:112125‑M2/150‑3 (MQ=255)                                  ccTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGtcgggcagcaccagggtctgcaaattgcccagc                                                                                                                      <  1:594529‑M1/150‑34 (MQ=255)                                  ccTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGtcgggcagcaccagggtctgcaaattgcccagc                                                                                                                      <  1:48299‑M1/150‑34 (MQ=255)                                       acaGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGaa                                                                                                                 <  2:504345/150‑1 (MQ=255)                                        caGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAAc                                                                                                                <  1:1311706/150‑1 (MQ=255)                                          gACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGtcgggcagcaccagggtctgcaaattgcccagcgcgtcgta                                                                                                              >  1:756571‑M1/1‑109 (MQ=255)                                              acgaGGGCAAGAACCACCAACTCGATTACACCTACAAACCCAACGGGCAACTCCTCAGCGAAACCACCAGCTCGGGCCTGCTGCTGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATc                                                                                                          <  1:12497/147‑1 (MQ=255)                                               ccAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGtcgggcagcaccagggtctgcaaattgcccagcgcgtcgtagcggt                                                                                                         <  2:756571‑M1/150‑47 (MQ=255)                                               ccAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGACTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCa                                                                                                         >  1:317523/1‑150 (MQ=255)                                                 aGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCAcc                                                                                                       <  2:1254538/150‑1 (MQ=255)                                                        gaagCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGtacta                                                                                                >  1:807127/1‑150 (MQ=255)                                                        gaagCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGG                                                                                                                                                       <  1:543680‑M1/95‑1 (MQ=255)                                                        gaagCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGG                                                                                                                                                       >  2:543680‑M1/1‑95 (MQ=255)                                                           gcagcaACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCa                                                                                                                                             <  1:224723/102‑1 (MQ=255)                                                           gcagcaACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCACGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCa                                                                                                                                             >  2:224723/1‑102 (MQ=255)                                                               caACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGc                                                                                         <  1:729837/150‑1 (MQ=255)                                                                               ccgcccccAACGGGCAACTGCTCAGCGAAATCAACAGCGCGGGCCCGCTGCAGTACCGCTACGACGCGCTGTGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGa                                                                         <  2:1009250/146‑1 (MQ=255)                                                                                 gACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCCGACACCTGCACCAGATc                                                                       <  2:729353/150‑1 (MQ=255)                                                                                    cccAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCATATCAAc                                                                    >  2:555664/1‑150 (MQ=255)                                                                                      cAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCTGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCg                                                                                        <  1:359290/128‑1 (MQ=255)                                                                                      cAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCg                                                                                        >  2:359290/1‑128 (MQ=255)                                                                                          gggCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAAc                                                              <  2:112125/150‑1 (MQ=255)                                                                                          gggCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAAc                                                              <  1:725900/150‑1 (MQ=255)                                                                                               aCTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGacac                                                                                   <  1:681420/124‑1 (MQ=255)                                                                                               aCTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGacac                                                                                   >  2:681420/1‑124 (MQ=255)                                                                                                   cTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTg                                                     <  1:1221066/150‑1 (MQ=255)                                                                                                         gAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGc                                               <  2:1157488/150‑1 (MQ=255)                                                                                                           aaCCAACAGCGCGGGCCTGCTGCAGAACCGTAACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGa                                             >  1:1106771/1‑150 (MQ=255)                                                                                                            aCCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGAt                                            <  2:1037597/150‑1 (MQ=255)                                                                                                              cAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTTCACCAGATCAACTTGAACGGCCGGGTGATCAGCGAttt                                          >  2:1302257/1‑150 (MQ=255)                                                                                                               aaCAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGTCCGGGTGATCAGCGATTTc                                         >  2:1134757/1‑150 (MQ=255)                                                                                                                 cAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGc                                                                                         <  2:1373712/100‑1 (MQ=255)                                                                                                                 cAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGc                                                                                         >  1:1373712/1‑100 (MQ=255)                                                                                                                       gggCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAg                                                                          >  2:689400/1‑109 (MQ=255)                                                                                                                       gggCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAg                                                                          <  1:689400/109‑1 (MQ=255)                                                                                                                           ctgctgCAGTACCGCTACGACGCGCTGGGCAATTTGCAGATCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACGTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCg                             >  1:1357328/1‑150 (MQ=255)                                                                                                                            tgctgcAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCgg                            <  2:317523/150‑1 (MQ=255)                                                                                                                             gctgcAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGt                           >  1:1077937/1‑150 (MQ=255)                                                                                                                             gctgcAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGt                           >  1:990621/1‑150 (MQ=255)                                                                                                                                gcAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCa                        <  1:1037859/150‑1 (MQ=255)                                                                                                                                 cAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCac                       >  1:517477/1‑150 (MQ=255)                                                                                                                                 cAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCac                       <  1:235704/150‑1 (MQ=255)                                                                                                                                     aCCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGAccagcgcgtcgtagcggtactgcagcaggcccgcgctgttggtttcgctgagcagttgcccgttggggtcgtaggtgtaatcgagttgctgcttcttgccctgg                   <  2:756571‑M2/150‑105 (MQ=255)                                                                                                                                       cGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTg                          >  2:130992/1‑141 (MQ=255)                                                                                                                                       cGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTg                          <  1:130992/141‑1 (MQ=255)                                                                                                                                       cGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAg                 <  2:1077937/150‑1 (MQ=255)                                                                                                                                        gCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAgg                >  1:190038/1‑150 (MQ=255)                                                                                                                                          tACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGAccagcgcgtcgtagcggtactgcagcaggcccgcgctgttggtttcgctgagcagttgcccgttggggtcgtaggtgtaatcgagttgctgcttcttgccctggttgtc              >  1:756571‑M2/1‑41 (MQ=255)                                                                                                                                             gacgCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTCCTg           >  1:162253/1‑150 (MQ=255)                                                                                                                                                  gcTGGGCAATTTGCAGACCCTGGTGCTGCCCGAccagcgcgtcgtagcggtactgcagcaggcccgcgctgttggtttcgctgagcagttgcccgttggggtcgtaggtgtaatcgagttgctgcttcttgccctggttgtctgtgaagg      <  1:594529‑M2/150‑118 (MQ=255)                                                                                                                                                  gcTGGGCAATTTGCAGACCCTGGTGCTGCCCGAccagcgcgtcgtagcggtactgcagcaggcccgcgctgttggtttcgctgagcagttgcccgttggggtcgtaggtgtaatcgagttgctgcttcttgccctggttgtctgtgaagg      <  1:48299‑M2/150‑118 (MQ=255)                                                                                                                                                  gcTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCaa                                                                     <  2:504729/87‑1 (MQ=255)                                                                                                                                                  gcTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCaa                                                                     >  1:504729/1‑87 (MQ=255)                                                                                                                                                     ggtcgggcagcaccagggtctgcaaattgcccagcgcgtcgtagcggtactgcagcaggcccgcgctgttggtttcgctgagcagttgcccgttggggtcgtaggtgtaatcgagttgctgcttcttgccctggttgtctgtgaaggtga   <  1:112125‑M1/150‑149 (MQ=255)                                                                                                                                                      ggCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTCCTGCGCACCCaa  >  1:782406/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GCTACGACGCGGCCGACAATCTGCTGTCCATCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTCCTGCGCACCCAA  >  NC_002947/4625403‑4625701

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