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

TCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACACGCT  >  NC_002947/4625433‑4625726                                                                                                                                                    |                                                                                                                                                   tCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGAcc                                                                                                                                                  <  1:112125‑M2/150‑3 (MQ=255)         acaGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGaa                                                                                                                                          <  2:504345/150‑1 (MQ=255)          caGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAAc                                                                                                                                         <  1:1311706/150‑1 (MQ=255)                acgaGGGCAAGAACCACCAACTCGATTACACCTACAAACCCAACGGGCAACTCCTCAGCGAAACCACCAGCTCGGGCCTGCTGCTGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATc                                                                                                                                   <  1:12497/147‑1 (MQ=255)                 ccAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGACTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCa                                                                                                                                  >  1:317523/1‑150 (MQ=255)                   aGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCAcc                                                                                                                                <  2:1254538/150‑1 (MQ=255)                          gaagCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGtacta                                                                                                                         >  1:807127/1‑150 (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                                                                                       <  1:725900/150‑1 (MQ=255)                                                            gggCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAAc                                                                                       <  2:112125/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                                                                                                                  >  1:1373712/1‑100 (MQ=255)                                                                                   cAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGc                                                                                                                  <  2:1373712/100‑1 (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:990621/1‑150 (MQ=255)                                                                                               gctgcAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGt                                                    >  1:1077937/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                                                   <  1:130992/141‑1 (MQ=255)                                                                                                         cGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTg                                                   >  2:130992/1‑141 (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:48299‑M2/150‑118 (MQ=255)                                                                                                                    gcTGGGCAATTTGCAGACCCTGGTGCTGCCCGAccagcgcgtcgtagcggtactgcagcaggcccgcgctgttggtttcgctgagcagttgcccgttggggtcgtaggtgtaatcgagttgctgcttcttgccctggttgtctgtgaagg                               <  1:594529‑M2/150‑118 (MQ=255)                                                                                                                    gcTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCaa                                                                                              >  1:504729/1‑87 (MQ=255)                                                                                                                    gcTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCaa                                                                                              <  2:504729/87‑1 (MQ=255)                                                                                                                        ggCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTCCTGCGCACCCaa                           >  1:782406/1‑150 (MQ=255)                                                                                                                          cAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGTCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTCCTGCGCACCCAAgg                         <  1:1134757/150‑1 (MQ=255)                                                                                                                              ttGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCATATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCaaa                     <  1:555664/150‑1 (MQ=255)                                                                                                                                 cAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTg                  >  1:462977/1‑150 (MQ=255)                                                                                                                                   gACCCTGGTGCTGCCCGACCAGCGCGAACTAAATCACCTGTACTGCGAAACGGGAGACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGt                <  2:1106771/150‑1 (MQ=17)                                                                                                                                    aCCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTc               <  2:807127/150‑1 (MQ=255)                                                                                                                                      ccTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGt                                                                               >  1:334429/1‑84 (MQ=255)                                                                                                                                      ccTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGt                                                                               <  2:334429/84‑1 (MQ=255)                                                                                                                                      ccTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGt                                                    <  2:1255628/111‑1 (MQ=255)                                                                                                                                      ccTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGt                                                    >  1:1255628/1‑111 (MQ=255)                                                                                                                                      ccTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCAc             >  1:230548/1‑150 (MQ=255)                                                                                                                                                 ccGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACACGCt  >  2:407417/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TCACCTTCACAGACAACCAGGGCAAGAAGCAGCAACTCGATTACACCTACGACCCCAACGGGCAACTGCTCAGCGAAACCAACAGCGCGGGCCTGCTGCAGTACCGCTACGACGCGCTGGGCAATTTGCAGACCCTGGTGCTGCCCGACCAGCGCGAACTCAATCACCTGTACTACGGCAGCGGACACCTGCACCAGATCAACTTGAACGGCCGGGTGATCAGCGATTTCGAGCGCGATGCGGTGCACGACGAGGTCCTGCGCACCCAAGGCAAACTGGTCACCCGCACACGCT  >  NC_002947/4625433‑4625726

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