New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	687851 =	NA (NA)	30 (0.630)	23/282	1.5	96.8%	intergenic (‑104/+146)	gltI/insH1	glutamate/aspartate periplasmic binding protein/IS5 transposase and trans‑activator
	?	NC_000913_3_bme_pgi	869874 =	1 (0.020)					coding (464/1539 nt)	gsiB	glutathione ABC transporter periplasmic binding protein

AGCGCCAGGGCGAGGATTGCTGTGGCAGGTTTACGTAATTGCATATCCAACTCCTTTATCTTCTGCGTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAG  >  NC_000913_3_bme_pgi/687704‑687993                                                                                                                                                    |                                                                                                                                               aGCGCCAGGGCGAGGATTGCTGTGGCAGGTTTACGTAATTGCATATCCAACTCCTTTATCTTCTGCGTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAggg                                                                                                                                               >  1:511925/1‑149 (MQ=255)  gcgcCAGGGCGAGGATGGCTGTGGCAGGTTTACGTAATTGCATATCCAACTCCTTTATCTTCTGCGTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGa                                                                                                                                              <  1:673140/149‑1 (MQ=255)  gccccAGGGCGAGGATGGCTGTGGCAGGTTTATGTAATTGCATATCCAACTCCTTTATCTTCTGCGTTAAGACCGCATCGACACTCGTTGTGATTGTTTGTTTTATTATCGTCTGCAACTTTATTGTGCATTGTTGTGCCTGTTAGGGa                                                                                                                                              <  2:349813/146‑1 (MQ=255)           cGAGGATTGCTGTGGCAGGTTTACGTAATTGCATATCCAATTCCTTTTTCGTCTGCGTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGaa                                                                                                                                     <  1:113132/149‑1 (MQ=255)                  tGCTGTGGCAGGTTTACGTAATTGCATATCCAACTCCTTTATCTTCTGCGTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCgg                                                                                                                              <  1:110663/149‑1 (MQ=255)                   gCTGTGGCAGGTTTACGTAATTGCATATCCAACTCCTTTATCTTCTGCGTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCggg                                                                                                                             >  2:670215/1‑149 (MQ=255)                    cTGTGGCAGGTTTACGTAATTGCATATCCAACTCCTTTATCTTCTGCGTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCgggg                                                                                                                            >  2:440588/1‑149 (MQ=255)                     tgtgGCAGGTTTACGTAATTGCATATCCAACTCCTTTATCTTCTGCGTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGa                                                                                                                           >  1:171634/1‑149 (MQ=255)                         gCAGGTTTACGTAATTGCATATCCAACTCCTTTATCTTCTGCGTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAAtt                                                                                                                       >  2:641772/1‑149 (MQ=255)                           aGGTTTACGTAATTGCATATCCAACTCCTTTATCTTCTGCGTTAAGAACGCATGGATACCCGTTGAGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAAttct                                                                                                                     >  1:699364/1‑149 (MQ=255)                              tttACGTAATTGCATATCCAACTCCTTTATCTTCTGCGTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATtcttct                                                                                                                  >  1:603583/1‑149 (MQ=255)                                           agatCCAGCTGCTTTATCTTCTGCGTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGt                                                                                                     <  1:812287/147‑1 (MQ=255)                                            tatCCAACTCCTTTATCTTCTGCGTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtc                                                                                                    >  2:246161/1‑149 (MQ=255)                                            tatCCAACTCCTTTATCTTCTGCGTTAAGAACGCATGGATACCCGTTGCGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtc                                                                                                    <  2:734644/149‑1 (MQ=255)                                              tCCAACGTCTTGATTTTGGGCGTTAAGAACGTATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTAGTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcat                                                                                                  <  2:695879/149‑1 (MQ=37)                                                cAACTCCTTCATCGTATGCGTTAAGAACGCACGGATACCCGTTGTAAGTGTTTGTGTTGTTATCGTCTGCAACTTTACTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcattt                                                                                                <  2:667587/149‑1 (MQ=37)                                                        ttATCTTCTGCGTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTAGGCTGACTCAGTCATTTCATttctt                                                                                        >  2:554059/1‑149 (MQ=255)                                                          aTCTTCTGCGTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAAtt                                                                                                                       <  1:338629/116‑1 (MQ=255)                                                          aTCTTCTGCGTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAAtt                                                                                                                       >  2:338629/1‑116 (MQ=255)                                                              tctGCGTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTGTCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGtt                                                                                  <  1:368428/149‑1 (MQ=255)                                                               ctGCGTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGttt                                                                                 <  1:474217/149‑1 (MQ=255)                                                                   gTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTTTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGc                                                                             <  1:366331/149‑1 (MQ=255)                                                                   gTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGc                                                                             >  2:385950/1‑149 (MQ=255)                                                                    ttAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGcc                                                                            >  2:641170/1‑149 (MQ=255)                                                                            gCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAtttttt                                                                    >  2:819198/1‑149 (MQ=255)                                                                                ggATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTcc                                                                >  1:217308/1‑149 (MQ=255)                                                                                   tACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAGTTTTTCTACCGt                                                             >  1:61015/1‑149 (MQ=255)                                                                                     cccGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTaa                                                           >  1:249658/1‑149 (MQ=255)                                                                                      ccGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGt                                                                                                                                          <  2:829468/69‑1 (MQ=255)                                                                                      ccGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGt                                                                                                                                          >  1:829468/1‑69 (MQ=255)                                                                                      ccGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTaaa                                                          >  2:312685/1‑149 (MQ=255)                                                                                      ccGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTaaa                                                          >  2:645163/1‑149 (MQ=255)                                                                                       cGTTGTGAGTGTTTGTGTTGTTATCGTCGGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGCAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAAt                                                         <  2:743078/149‑1 (MQ=255)                                                                                        gTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCggg                                                                                                                             >  1:511369/1‑80 (MQ=255)                                                                                        gTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCggg                                                                                                                             <  2:511369/80‑1 (MQ=255)                                                                                           gtgGGTGTTTGGGTTGGTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCt                                                     <  1:210918/149‑1 (MQ=255)                                                                                           gtgAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCt                                                     <  1:312685/149‑1 (MQ=255)                                                                                            tgAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCtt                                                    <  2:786483/149‑1 (MQ=255)                                                                                             gAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTg                                                   >  1:437049/1‑149 (MQ=255)                                                                                               gtgtTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCa                                                                                      <  1:454320/112‑1 (MQ=255)                                                                                               gtgtTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCa                                                                                      >  2:454320/1‑112 (MQ=255)                                                                                                  tttGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAGTAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCa                                              >  1:387567/1‑149 (MQ=255)                                                                                                       gttgttATCGTCTGCAACTTTATTGTGCAGTGTTGTGCGTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCCGCCTa                                         >  1:790077/1‑149 (MQ=255)                                                                                                       gttgttATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGGCTTGAATCAGCCTa                                         >  1:824102/1‑149 (MQ=255)                                                                                                       gttgttATCGTCTGCAACTTTATTGTGCAGTGTTATGCCTGTTAGGGAAGGTGCGAATAAGCGGGAAAATTCTTCTCGGCTGACTCAGTCATTTCATTCCTTCACGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTa                                         >  2:663439/1‑149 (MQ=255)                                                                                                                gTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTTAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCg                                >  1:193052/1‑149 (MQ=255)                                                                                                                      aaCTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttctt                          >  1:537486/1‑149 (MQ=255)                                                                                                                       aCTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttcttc                         <  1:246161/149‑1 (MQ=255)                                                                                                                       aCTTTATTGTGCAGTGGTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGTAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttcttc                         <  2:47435/149‑1 (MQ=255)                                                                                                                         gttATTGGGCAGTGTTGGGTCTCTTCGGGAAGGTGCGAATAAGCGGGGAACTTCTTCTCGGCCAACTCAGTCCTTTCATTTCTTCATGTTTGAACCGATTTTTTCTCCCCTAAATGCCTGGCATCAGCCTATTGAGACCGTTTCTTCGc                       <  2:819848/148‑1 (MQ=11)                                                                                                                                   aGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGc             >  1:192297/1‑149 (MQ=39)                                                                                                                                       ttgtGACTGGTAGGGAAGGTGCGAATAAGCGGGTAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTa         <  2:249658/149‑1 (MQ=12)                                                                                                                                         gggCCTGTTAGGGAAGGTGCGAGTAAGCGGGGAAATTCTTCTCGGCTGACTCAGGCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTACATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATc       <  1:397179/147‑1 (MQ=255)                                                                                                                                           gCCTGTTAGGGAAGGTGCGAATAGGCGGGGAAATTCTTCTCGGCTGCCTCAGTCATTTCATTTTTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATccc     <  2:335332/149‑1 (MQ=12)                                                                                                                                            ccTGTTAGGGAAGGTGCGACTAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCAGATTTTTTCTCCCGTCAGTGCCTTGAGTCAGCCCATTTAGACCCATTATGCGCCGTTCAAAGCGTTATcccc    >  2:706014/1‑149 (MQ=11)                                                                                                                                            ccTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATATTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATATAGACCGTTTCTTCGCCATTTAAGGCGTTATcccc    >  1:121608/1‑149 (MQ=17)                                                                                                                                              tGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAg  >  1:695816/1‑149 (MQ=11)                                                                                                                                                    |                                                                                                                                               AGCGCCAGGGCGAGGATTGCTGTGGCAGGTTTACGTAATTGCATATCCAACTCCTTTATCTTCTGCGTTAAGAACGCATGGATACCCGTTGTGAGTGTTTGTGTTGTTATCGTCTGCAACTTTATTGTGCAGTGTTGTGCCTGTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAG  >  NC_000913_3_bme_pgi/687704‑687993

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