New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	3919152 =	18 (0.960)	3 (0.170)	3/238	NT	16.9%	coding (2890/3129 nt)	mexB	multidrug resistance protein MexB
	?	NC_002947	3919185 =	12 (0.660)					coding (2923/3129 nt)	mexB	multidrug resistance protein MexB

TGATCGGCCTGTCGGCGAAGAACGCCATCCTGATCATCGAGTTCGCCAAGCATCTGGTTGACCAGGGCGTGGACGCCGCCGACGCCGCCGTGCAAGCGGCCCGCCTGCGCCTGCGGCCGATCGTGATGACGTCGCTGGCGTTCATCCTCGGCGTGGTGCCCCTGGCCATCGCCAGCGGCGCAAGCTCTGCCAGCCAGCAGGCGATTGGCACCGGGGTGATCGGCGGCATGCTCAGCGCCACCCTGGCGGT  >  NC_002947/3919045‑3919294                                                                                                                                             |                                                                                                              tGATCGGCCTGTCGGCGAAGAACGCCATCCTGATCATCGAGTTCGCCAAGCATCTGGTTGACCAGGGCGTGGACGCCGCCGACGCCGCCGTGCAAGCGGCCCGCCTGCGCCTGCGGCCGATCGTGATGACGTCGCTGGCGTTCATc                                                                                                          <  1:132540/146‑1 (MQ=255) tGATCGGCCTGTCGGCGAAGAACGCCATCCTGATCATCGAGTTCGCCAAGCATCTGGTTGACCAGGGCGTGGACGCCGCCGACGCCGCCGTGCAAGCGGCCCGCCTGCGCCTGCGGCCGATCGTGATGACGTCGCTGGCGTTCATc                                                                                                          >  2:132540/1‑146 (MQ=255) gccagggtggcgctgagcatgccgccgatcaccccggtgccaatcgcctgctggctggcagagcttgcgccgctggcgatggccaggggcaccacgccgaggatgaaCGCCTGCGGCCGATCGTGATGACGTCGCTGGCGTTCATCCTCg                                                                                                      >  1:98982‑M1/108‑150 (MQ=255)               gCGAAGAACGCCATCCTGATCATCGAGTTCGCCAAGCATCTGGTTGACCAGGGCGTGGACGCCGCCGACGCCGCCGTGCAAGCGGCCCGCCTGCGCCTGCGGCCGATCGTGATGACGTCGCTGGCGTTCATCCTCGGCGTGGTGCCCCTg                                                                                        >  2:194566/1‑150 (MQ=255)                                          ttCGCCAAGCATCTGGTTGACCAGGGCGTGGACGCCGCCGACGCCGCCGTGCAAGCGGCCCGCCTGCGCCTGCGGCCGATCGTGATGACGTCGCTGGCGTTCATCCTCGGCGTGGTGCCCCTg                                                                                        <  2:269328/123‑1 (MQ=255)                                          ttCGCCAAGCATCTGGTTGACCAGGGCGTGGACGCCGCCGACGCCGCCGTGCAAGCGGCCCGCCTGCGCCTGCGGCCGATCGTGATGACGTCGCTGGCGTTCATCCTCGGCGTGGTGCCCCTg                                                                                        >  1:269328/1‑123 (MQ=255)                                                aaGCATCTGGTTGACCAGGGCGTGGACGCCGCCGACGCCGCCGTGCAAGCGGCCCGCCTGCGCCTGCGGCCGATCGTGATGACGTCGCTGGCGTTCATCCTCGGCGTGGTGCCCCTGGCCATCGCCAGCGGCGCAAGCTCTGccagccag                                                       >  2:89873/1‑150 (MQ=255)                                                      gctggcagagcttgcgccgctggcgatggccaggggcaccacgccgaggatgaaCGCCTGCGGCCGATCGTGATGACGTCGCTGGCGTTCATCCTCGGCGTGGTGCCCCTg                                                                                        >  2:26417‑M1/55‑111 (MQ=255)                                                      gctggcagagcttgcgccgctggcgatggccaggggcaccacgccgaggatgaaCGCCTGCGGCCGATCGTGATGACGTCGCTGGCGTTCATCCTCGGCGTGGTGCCCCTg                                                                                        <  1:26417‑M1/57‑1 (MQ=255)                                                         gTTGACCAGGGCGTGGACGCCGCCGACGCCGCCGTGCAAGCGGCCCGCCTGCGCCTGCGGCCGATCGTGATGACGTCGCTGGCGTTCATCCTCGGCGTGGTGCCCCTGGCCATCCCCAGCGGCGCAAGCTCTGCCAGCCAGCAGGCGAtt                                              <  1:194566/150‑1 (MQ=255)                                                                              gccgacgccgccgTGCAAGCGGCCCGCCTGCGCCTGCGGCCGATCGTGATGACGTCGCTGGCGTTCATCCTCGGCGTGGTGCCCCTGGCCATCGCCAgcggc                                                                         <  1:68257/102‑1 (MQ=255)                                                                              gccgacgccgccgTGCAAGCGGCCCGCCTGCGCCTGCGGCCGATCGTGATGACGTCGCTGGCGTTCATCCTCGGCGTGGTGCCCCTGGCCATCGCCAgcggc                                                                         >  2:68257/1‑102 (MQ=255)                                                                                    caggggcaccacgccgaggatgaacgccagcgacgtcatcacgatcggccgcaggcgTTCATCCTCGGCGTGGTGCCCCTGGCCATCGCCAGCGGCGCAAGCTCTgccagc                                                          <  1:26417‑M2/54‑1 (MQ=255)                                                                                    caggggcaccacgccgaggatgaacgccagcgacgtcatcacgatcggccgcaggcgTTCATCCTCGGCGTGGTGCCCCTGGCCATCGCCAGCGGCGCAAGCTCTgccagc                                                          >  2:26417‑M2/58‑111 (MQ=255)                                                                                               aGCGGCCCGCCTGCGCCTGCGGCCGATCGTGATGACGTCGCTGGCGTTCATCCTCGGCGTGGTGCCCCTGGCCATCGCCAGCGGCGCAAGCTCTGCCAGCCAGCAGGCGATTGGCACCGGGGTGATCGGCGGCATGCTCAGCGCCACCCt        >  2:434221/1‑150 (MQ=255)                                                                                                 cGGCCCGCCTGCGCCTGCGGCCGATCGTGATGACGTCGCTGGCGTTCATCCTCGGCGTGGTGCCCCTGGCCATCGCCAGCGGCGCAAGCTCTGCCAGCCAGCAGGCGATTGGCACCGGGGTGATCGGCGGCATGCTCAGCGCCACCCTgg      >  1:82392/1‑150 (MQ=255)                                                                                                  cgaggatgaacgccagcgacgtcatcacgatcggccgcaggcgTTCATCCTCGGCGTGGTGCCCCTGGCCATCGCCAGCGGCGCAAGCTCTGCCAGCCAGCAGGCGATTGGCACCGGGGTGATCGGCGGCATGCTCAGCGCCACCCTGGc     >  1:98982‑M2/44‑150 (MQ=255)                                                                                                     ccGCCTGCGCCTGCGGCCGATCGTGATGACGTCGCTGGCGTTCATCCTCGGCGTGGTGCCCCTGGCCATCGCCAGCGGCGCAAGCTCTGCCAGCCAGCAGGCGATTGGCACCGGGGTGATCGGCGGCATGCTCAGCGCCACCCTGGCggt  >  2:357782/1‑150 (MQ=255)                                                                                                                                             |                                                                                                              TGATCGGCCTGTCGGCGAAGAACGCCATCCTGATCATCGAGTTCGCCAAGCATCTGGTTGACCAGGGCGTGGACGCCGCCGACGCCGCCGTGCAAGCGGCCCGCCTGCGCCTGCGGCCGATCGTGATGACGTCGCTGGCGTTCATCCTCGGCGTGGTGCCCCTGGCCATCGCCAGCGGCGCAAGCTCTGCCAGCCAGCAGGCGATTGGCACCGGGGTGATCGGCGGCATGCTCAGCGCCACCCTGGCGGT  >  NC_002947/3919045‑3919294

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