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

GTGCACAACACCGGCATCAATACCTCGTTGCTGGTCGAGCACGCCTACAACCCGATCCACCGCCGTTACCAGTACGACCCGGCCGGTGAACTGGTGCGTACCCTCGACAAGCTACGCGGCGAGATCAAGTACGAGTACGAAGCCAACGGCCAGTTGCGCAGTCGCGATACCGGCTCGTTGATCGGCAGTGAGGAGTTCCGCTACGACCCGGCAGCCAACCGCCTGGACTTCAATGCGCGGCAATTCGACAAGGTCAAGGACAACCGGATCAAGCAGTGGCGGGACCAGGAGTACC  >  NC_002947/3514262‑3514556                                                                                                                                                    |                                                                                                                                                    gTGCACAACACCGGCATCAATACCTCGTTGCTGGTCGAGCACGCCTACAACCCGATCCACCGCCGTTACCAGTACGACCCGGCCGGTGAACTGGTGCGTACCCTCGACAAGCTACGCGGCGAGATCAAGTACGAGTACGAAGCCAACGGc                                                                                                                                                   <  1:186456/150‑1 (MQ=255)                             tGCTGGTCGAGCACGCCTACAACCCGATCCACCGCCGTTACCAGTACGACCCGGCCGGTGAACTGGTGCGTACCCTCGACAAGCTACGCGGCGAGATCAAGTACGAGTACGAAGCCAACGGCCAGTTGCGCAGTCGCGATACCGGCTcgt                                                                                                                       >  2:168901‑M2/1‑147 (MQ=255)                              gCTGGTCGACCACGCCTACAACCCGATCCACCGCCGTTACCAGTACGACCCGGCCGGTGAACTGGTGCGTACCCTCGACAAGCTACGCGGCGAGATCAAGTACGAGTACGAAGCCAACGGCCAGTTGCGCAGTCGCGATACCGGCTcgtt                                                                                                                      >  2:116550‑M2/1‑146 (MQ=255)                                     gAGCACGCCTACAACCCGATCCACCGCAGTTACCAGTACGACCCGGCCGGTGAACTGTTTCGTACCCTCAACAAGCTACGCGGCGAGATCAAGTACGAGTACGAAGCCAACGGCCAGTTGCGCAGTCGCGATACCGGCTCGTTGATCGGc                                                                                                               >  1:277782/1‑150 (MQ=255)                                               aCAACCCGATCCACCGCCGTTACCAGTACGACCCGGCCGGTGAACTGGTGCGTACCCTCGACAAGCTACGCGGCGAGATCAAGTACGAGTACGAAGCCAACGagccggtatcgcgactgcgcaactggccgttggcttcgtactcgtact                                                                                                     <  1:396661‑M1/150‑49 (MQ=255)                                                                cgTTACCAGTACGACCCGGCCGGTGAACTGGTGCGTACCCTCGACAAGCTACGCGGCGAGATCAAGTACGAGTACGAAGCCAACGGCCAGTTGCGCAGTCGCGATACCGGCTCGTTGATCGGCAGTGAGGAGTTCCGCTACGACCCGGCa                                                                                    <  2:350736/150‑1 (MQ=255)                                                                      cAGTACGACCCGGCCGGTGAACTGGTGAGTACCCTCGACAAGCTACGCGGCGAGATCAAGTACGAGTACGAAGCCAACGagccggtatcgcgactgcgcaactggccgttggcttc                                                                                                                <  1:379686‑M1/116‑38 (MQ=255)                                                                      cAGTACGACCCGGCCGGTGAACTGGTGAGTACCCTCGACAAGCTACGCGGCGAGATCAAGTACGAGTACGAAGCCAACGagccggtatcgcgactgcgcaactggccgttggcttc                                                                                                                >  2:379686‑M1/1‑79 (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)                                                                                                                                                 aaCGagccggtatcgcgactgcgcaactggccgttggcttcgtactcgtacttgatctcgccgcgtagcttgtcgagggtacgcaccagttcaccggccgggtcgtactggtaacggcggtggatcgggttgtaggcgtggtcgaccagc   >  2:116550‑M1/1‑4 (MQ=255)                                                                                                                                                 aaCGGCCAGTTGCGCAGTCGCGATACCGGCTCGTTGATCGGCAGTGAGGAGTTCCGCTACGACCCGGCAGCCAACCGCCTGGACTTCAATGCGCGGCAATTCGACAAGGTCAAGGACAACCGGATCAAGCAGTGGCGGGACCAGGAGTAc   >  1:420823/1‑150 (MQ=255)                                                                                                                                                  aCGagccggtatcgcgactgcgcaactggccgttggcttcgtactcgtacttgatctcgccgcgtagcttgtcgagggtacgcaccagttcaccggccgggtcgtactggtaacggcggtggatcgggttgtaggcgtgctcgaccagca  >  2:168901‑M1/1‑3 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GTGCACAACACCGGCATCAATACCTCGTTGCTGGTCGAGCACGCCTACAACCCGATCCACCGCCGTTACCAGTACGACCCGGCCGGTGAACTGGTGCGTACCCTCGACAAGCTACGCGGCGAGATCAAGTACGAGTACGAAGCCAACGGCCAGTTGCGCAGTCGCGATACCGGCTCGTTGATCGGCAGTGAGGAGTTCCGCTACGACCCGGCAGCCAACCGCCTGGACTTCAATGCGCGGCAATTCGACAAGGTCAAGGACAACCGGATCAAGCAGTGGCGGGACCAGGAGTACC  >  NC_002947/3514262‑3514556

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.