New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	16219 =	10 (0.440)	3 (0.140)	3/258	1.3	18.3%	intergenic (+469/‑370)	gyrB/PP_5421	DNA gyrase subunit B/hypothetical protein
	?	NC_002947	= 16362	17 (0.770)					intergenic (+612/‑227)	gyrB/PP_5421	DNA gyrase subunit B/hypothetical protein
Rejected: Frequency below/above cutoff threshold.

TCCAAATTGCTCCAGCGATCACGAAGTGGCCGTTTAGAAATAAAGAGAGCAAGCCGATTGCTTTAGATTCTGATCAGTCATGAGAGTCATGACCAATGAGAGAAAAGGAGACCTGAAAATGTAGAAAAACGCGATCAATCCCGAGAGGACTGACCGCGCTTAGCGAAGCAGACTTGCGTCCTCAGAACCCGCAATGTCTGCTTTTTCGAAGCGACTGTCAATCCTCTCCCATTTTTTCGACCTTAGCGAGGTTACTCGACG  >  NC_002947/16229‑16489                                                                                                                                      |                                                                                                                                accaAATTGCTCCAGCGATCACGAAGTGGCCGTTTAGAAATAAAGAGAGCAAGCCGATTGCTTTAGATTCTGATCAGTCATGAGAGTCATGACCAATGAGAGAAAAGGAGACCTGAAAATGTAGAAAAACGCGAtc                                                                                                                               >  2:565136‑M2/2‑134 (MQ=255)   cAAATTGCTCCAGCGATCACGAAGTGGCCGTTTAGAAATAAAGAGAGGAAGCCGATTGCTTTAGATTCTGATCAGTCATGAGAGTCATGACCAATGAGAGAAAAGGAGACCTGAAAATGTAGAAAAACGCGATCaa                                                                                                                             >  2:501574/1‑136 (MQ=255)   cAAATTGCTCCAGCGATCACGAAGTGGCCGTTTAGAAATAAAGAGAGCAAGCCGATTGCTTTAGATTCTGATCAGTCATGAGAGTCATGACCAATGAGAGAAAAGGAGACCTGAAAATGTAGAAAAACGCGATCaa                                                                                                                             >  1:326572/1‑136 (MQ=255)   cAAATTGCTCCAGCGATCACGAAGTGGCCGTTTAGAAATAAAGAGAGCAAGCCGATTGCTTTAGATTCTGATCAGTCATGAGAGTCATGACCAATGAGAGAAAAGGAGACCTGAAAATGTAGAAAAACGCGATCaa                                                                                                                             >  2:368317/1‑136 (MQ=255)        tGCTCCAGCGATCACGAAGTGGCCGTTTAGAAATAAAGAGAGCAAGCCGATTGCTTTAGATTCTGATCAGTCATGAGAGTCATGACCAATGAGAGAAAAGGAGACCTGAAAATGTAGAAAAACGCGAtcagcgaaa                                                                                                                        <  1:338653‑M2/136‑10 (MQ=255)             cAGCGATCACGAAGTGGCCGTTTAGAAATAAAGAGAGCAAGCCGATTGCTTTAGATTCTGATCAGTCATGAGAGTCATGACCAATGAGAGAAAAGGAGACCTGAAAATGTAGAAAAACGCGATCAATCCCGAGAgg                                                                                                                   >  2:572599/1‑136 (MQ=255)                    cACGAAGTGGCCGTTTAGAAATAAAGAGAGCAAGCCGATTGCTTTAGATTCTGATCAGTCATGAGAGTCATGACCAATGAGAGAAAAGGAGACCTGAAAATGTAGAAAAACGCGATCAATCCCGAGAGGACTGAcc                                                                                                            <  2:101938/136‑1 (MQ=255)                    cACGAAGTGGCCGTTTAGAAATAAAGAGAGCAAGCCGATTGCTTTAGATTCTGATCAGTCATGAGAGTCATGACCAATGAGAGAAAAGGAGACCTGAAAATGTAGAAAAACGCGATCAATCCCGAGAGGACTGAcc                                                                                                            <  2:220971/136‑1 (MQ=255)                    cACGAAGTGGCCGTTTAGAAATAAAGAGAGCAAGCCGATTGCTTTAGATTCTGATCAGTCATGAGAGTCATGACCAATGAGAGAAAAGGAGACCTGAAAATGTAGAAAAACGCGATCAATCCCGAGAGGACTGAcc                                                                                                            <  2:2534/136‑1 (MQ=255)                    cACGAAGTGGCCGTTTAGAAATAAAGAGAGCAAGCCGATTGCTTTAGATTCTGATCAGTCATGAGAGTCATGACCAATGAGAGAAAAGGAGACCTGAAAATGTAGAAAAACGCGATCAATCCCGAGAGGACTGAcc                                                                                                            <  2:481709/136‑1 (MQ=255)                                                                          tCAGTCATGAGAGTCATGACCAATGAGAGAAAAGGAGACCTGAAAATGTAGAAAAACGCGATCAATCCCGAGAGGACTGAcc                                                                                                            <  2:417536/82‑1 (MQ=255)                                                                          tCAGTCATGAGAGTCATGACCAATGAGAGAAAAGGAGACCTGAAAATGTAGAAAAACGCGATCAATCCCGAGAGGACTGAcc                                                                                                            >  1:417536/1‑82 (MQ=255)                                                                                aTGAGAGTCATGACCAATGAGAGAAAAGGAGACCTGAAAATGTAGAAAAACGCGATCAATCCCGAGAGGACTGACCGCGCTTAGCGAAGCAGACTTGCGTCCTCAGAACCCGCAATGTCTGCTTTTTCGAAGCGAc                                                <  1:572599/136‑1 (MQ=255)                                                                                           gACCAATGAGAGAAAAGGAGACCTGAAAATGTAGAAAAACGCGATCAATCCCGAGAGGACTGACCGCGCTTAGCGAAGCAGACTTGCGTCCTCAGAACCCGCAATGTCTGCTTTTTCGAAGCGACTGTCAATCctc                                     >  2:447222/1‑136 (MQ=255)                                                                                                          aGGAGACCTGAAAATGTAGAAAAACGCGAtcagacctgatccaaattgctccagcgatcacgaagtggccgtttagaaataaagagagcaagccgattgctttagattctgatc                                            <  1:50084‑M2/114‑86 (MQ=255)                                                                                                          aGGAGACCTGAAAATGTAGAAAAACGCGAtcagacctgatccaaattgctccagcgatcacgaagtggccgtttagaaataaagagagcaagccgattgctttagattctgatc                                            >  2:50084‑M2/1‑29 (MQ=255)                                                                                                          aGGAGACCTGAAAATGTAGAAAAACGCGATCAATCCCGAGAGGACTGACCGCGCTTAGCGAAGCAGACTTGCGTCCTCAGAACCCGCAATGTCTGCTTTTTCGAAGCGACTGTCAATCCTCTCCCATTTTTTCGAc                      >  1:119202/1‑136 (MQ=255)                                                                                                           ggAGACCTGAATATGTAGAAAAACGCGTTCAATCCCGAGAGGACTGACCGCGCTTAGCGAAGCAGACTTGCGTCCTCAGAACCCGCAATGTCTGCTTTTTCGAAGCGACTGTCAATCCTCTCCCATTTTTTCGAcc                     <  2:477663/136‑1 (MQ=255)                                                                                                                      aaTGTAGAAAAACGCGATCAATCCCGAGAGGACTGACCGCGCTTAGCGAAGCAGACTTGCGTCCTCAGAACCCGCAATGTCTGCTTTTTCGAAGCGACTGTCAATCCTCTCCCATTTTTTCGACCTTAGCGAGGtt          >  2:263646/1‑136 (MQ=255)                                                                                                                         gTAGAAAAACGCGATCAATCCCGAGAGGACTGACCGCGCTTAGCGAAGCAGACTTGCGTCCTCAGAACCCGCAATGTCTGCTTTTTCGAAGCGACTGTCAATCCTCTCCCATTTTTTCGACCTTAGCGAGGTTACt       <  1:447222/136‑1 (MQ=255)                                                                                                                              aaaaCGCGATCAATCCCGAGAGGACTGACCGCGCTTAGCGAAGCAGACTTGCGTCCTCAGAACCCGCAATGTCTGCTTTTTCGAAGCGACTGTCAATCCTCTCCCATTTTTTCGACCTTAGCGAGGTTACTCGACg  <  2:498998/136‑1 (MQ=255)                                                                                                                                      |                                                                                                                                TCCAAATTGCTCCAGCGATCACGAAGTGGCCGTTTAGAAATAAAGAGAGCAAGCCGATTGCTTTAGATTCTGATCAGTCATGAGAGTCATGACCAATGAGAGAAAAGGAGACCTGAAAATGTAGAAAAACGCGATCAATCCCGAGAGGACTGACCGCGCTTAGCGAAGCAGACTTGCGTCCTCAGAACCCGCAATGTCTGCTTTTTCGAAGCGACTGTCAATCCTCTCCCATTTTTTCGACCTTAGCGAGGTTACTCGACG  >  NC_002947/16229‑16489

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