New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	5639467 =	10 (0.600)	3 (0.190)	3/274	1.8	21.8%	coding (28/1443 nt)	PP_4949	TldD/PmbA family peptidase
	?	NC_002947	5639485 =	12 (0.750)					coding (46/1443 nt)	PP_4949	TldD/PmbA family peptidase

AGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCTCGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCACCAGGCGCTTTCGGTTCGCCGCAACGTCGCCGAACCACCCTTCTTCAGCCAGGATGAGGGCGCCATGCTCACCGTGCGG  >  NC_002947/5639327‑5639610                                                                                                                                             |                                                                                                                                                tccacgccattgacccgcacggtgagcatggcaccctcatcctggctgaagaagggtggttcggcgacgttgcggcgaaccgaaagcgcctggtgcgactgtttcacatggcgcagggaaaacaactcggccgtactgcgCAGCGCTTtg                                                                                                                                        >  2:197325‑M1/141‑150 (MQ=255)                    gAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCTCGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGc                                                                                                                     >  1:341391/1‑150 (MQ=255)                       cTGGCATTGCTCGAGGCACACGACGATACCTGCCTCGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGa                                                                                                                  <  1:31381/150‑1 (MQ=255)                          gCATTGCTCGAGGCACACGACGATACCTGCCTCGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAgtt                                                                                                               >  2:121290/1‑150 (MQ=255)                              tGCTCGAGGCACACGACGATACCTGCCTCGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGttt                                                                                                           <  2:321/150‑1 (MQ=255)                               gCTCGAGGCACACGACGATACCTGCCTCGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACAGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGtttt                                                                                                          <  1:227687/150‑1 (MQ=255)                                 tCGAGGCACACGACGATACCTGCCTCGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTcc                                                                                                        >  1:288376/1‑150 (MQ=255)                                                     tGCCTCGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAgttgt                                                                                                             <  1:218416/125‑1 (MQ=255)                                                     tGCCTCGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAgttgt                                                                                                             >  2:218416/1‑125 (MQ=255)                                                                      ttgcggcgaaccgaaagcgcctggtgcgactgtttcacatggcgcagggaaaacaactcggccgtactgcgCAGCGCTTTGCTGCGCTGCGCAGTACgg                                                                                                                      <  1:85993‑M1/28‑1 (MQ=255)                                                                      ttgcggcgaaccgaaagcgcctggtgcgactgtttcacatggcgcagggaaaacaactcggccgtactgcgCAGCGCTTTGCTGCGCTGCGCAGTACgg                                                                                                                      >  2:85993‑M1/72‑99 (MQ=255)                                                                                                        aGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCACCAGGCGCTTTCGGTTCGCCGCAACGTCGCCGAACCACCCTTCTTCAg                                 >  2:69729/1‑150 (MQ=255)                                                                                                                                       ctgcGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCACCAGGCGCTTTCGGTTCGCCGCAACGTCGCCGAACCACCCTTCTTCAGCCAGGATGAGGGCGCCATGCTCACCGTGCgg  <  2:361876/150‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                AGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCTCGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCACCAGGCGCTTTCGGTTCGCCGCAACGTCGCCGAACCACCCTTCTTCAGCCAGGATGAGGGCGCCATGCTCACCGTGCGG  >  NC_002947/5639327‑5639610

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