New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3514409	4 (0.250)	3 (0.190)	3/246	NT	37.9%	coding (3034/4158 nt)	PP_3108	rhs family protein
	?	NC_002947	= 3514436	6 (0.390)					coding (3061/4158 nt)	PP_3108	rhs family protein

TGCTGGTCGAGCACGCCTACAACCCGATCCACCGCCGTTACCAGTACGACCCGGCCGGTGAACTGGTGCGTACCCTCGACAAGCTACGCGGCGAGATCAAGTACGAGTACGAAGCCAACGGCCAGTTGCGCAGTCGCGATACCGGCTCGTTGATCGGCAGTGAGGAGTTCCGCTACGACCCGGCAGCCAACCGCCTGGACTTCAATGCGCGGCAATTCGACAAGGTCAAGGACAACCGGATCAAGCAGTGGCGGGACCAGGAGTACCGCTAC  >  NC_002947/3514290‑3514561                                                                                                                                                   |                                                                                                                              tGCTGGTCGAGCACGCCTACAACCCGATCCACCGCCGTTACCAGTACGACCCGGCCGGTGAACTGGTGCGTACCCTCGACAAGCTACGCGGCGAGATCAAGTACGAGTACGAAGCCAACGGCCAGTTGCGCAGTCGCGATACCGGCTcgt                                                                                                                            >  2:168901‑M2/1‑147 (MQ=255)  gCTGGTCGACCACGCCTACAACCCGATCCACCGCCGTTACCAGTACGACCCGGCCGGTGAACTGGTGCGTACCCTCGACAAGCTACGCGGCGAGATCAAGTACGAGTACGAAGCCAACGGCCAGTTGCGCAGTCGCGATACCGGCTcgtt                                                                                                                           >  2:116550‑M2/1‑146 (MQ=255)         gAGCACGCCTACAACCCGATCCACCGCAGTTACCAGTACGACCCGGCCGGTGAACTGTTTCGTACCCTCAACAAGCTACGCGGCGAGATCAAGTACGAGTACGAAGCCAACGGCCAGTTGCGCAGTCGCGATACCGGCTCGTTGATCGGc                                                                                                                    >  1:277782/1‑150 (MQ=255)                                    cgTTACCAGTACGACCCGGCCGGTGAACTGGTGCGTACCCTCGACAAGCTACGCGGCGAGATCAAGTACGAGTACGAAGCCAACGGCCAGTTGCGCAGTCGCGATACCGGCTCGTTGATCGGCAGTGAGGAGTTCCGCTACGACCCGGCa                                                                                         <  2:350736/150‑1 (MQ=255)                                                     ggccggTGAACTGGTGCGTACCCTCGACAAGCTACGCGGCGAGATCAAGTACGAGTACGAAGCCAACGGCCAGTTGCGCAGTCGCGATACCGGCTCGTTGATCGGCAGTGAGGAGTTCCGCTACGACCCGGCAGCCAACCGCCTGGACtt                                                                        <  2:72248/150‑1 (MQ=255)                                                                                                aTCAAGTACGAGTACGAAGCCAACGGCCAGTTGCGCAGTCGCGATACCGGCTCGTTGATCGGCAGTGAGGAGTTCCGCTACGACCCGGCAGCCAACCGCCTGGACTTCAATGCGCGGCAATTCGACAAGGTCAAGGACAACCGGATCAAg                             <  1:380323/150‑1 (MQ=255)                                                                                                    agtacgagtacgAAGCCAACGGCCAGTTGCGCAGTCGCGATACCGGCTcgttggcttcgtactcgtacttgatctcgccgcgtagcttgtcgagggtacgcaccagttcaccggccgggtcgtactggtaacggcggtggatcgggttgt                         <  1:396661‑M2/150‑103 (MQ=255)                                                                                                               gaAGCCAACGGCCAGTTGCGCAGTCGCGATACCGGCTcgttggcttcgtactcgtacttgatctcgccgcgtagcttgtcgagggtactcaccagttcaccggccgggtcgtactg                                                <  1:379686‑M2/116‑80 (MQ=255)                                                                                                               gaAGCCAACGGCCAGTTGCGCAGTCGCGATACCGGCTcgttggcttcgtactcgtacttgatctcgccgcgtagcttgtcgagggtactcaccagttcaccggccgggtcgtactg                                                >  2:379686‑M2/1‑37 (MQ=255)                                                                                                                     aaCGGCCAGTTGCGCAGTCGCGATACCGGCTCGTTGATCGGCAGTGAGGAGTTCCGCTACGACCCGGCAGCCAACCGCCTGGACTTCAATGCGCGGCAATTCGACAAGGTCAAGGACAACCGGATCAAGCAGTGGCGGGACCAGGAGTAc        >  1:420823/1‑150 (MQ=255)                                                                                                                           cAGTTGCGCAGTCGCGATACCGGCTCGTTGATCGGCAGTGAGGAGTTCCGCTACGACCCGGCAGCCAACCGCCTGGACTTCAATGCGCGGCAATTCGACAAGGTCAAGGACAACCGGATCAAGCAGTGGCGGGACCAGGAGTACCGCTAc  <  2:277782/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                              TGCTGGTCGAGCACGCCTACAACCCGATCCACCGCCGTTACCAGTACGACCCGGCCGGTGAACTGGTGCGTACCCTCGACAAGCTACGCGGCGAGATCAAGTACGAGTACGAAGCCAACGGCCAGTTGCGCAGTCGCGATACCGGCTCGTTGATCGGCAGTGAGGAGTTCCGCTACGACCCGGCAGCCAACCGCCTGGACTTCAATGCGCGGCAATTCGACAAGGTCAAGGACAACCGGATCAAGCAGTGGCGGGACCAGGAGTACCGCTAC  >  NC_002947/3514290‑3514561

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