New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 401687	33 (0.910)	4 (0.150) +GTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCAT	3/212	4.1	19.4%	intergenic (+214/+215)	PP_0334/PP_0335	transposase/glutathione S‑transferase family protein
	?	NC_002947	= 3969351	12 (0.330)					intergenic (+300/‑15)	PP_3498/PP_3499	transposase/transposase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCAGAAGCAAAAAGCTTCCTCGAATAGGCCTGATACATAGATGCAACCGGGTTCCTATGGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCAACGAAGCTCTTATCCGAAACCGTTCAAAGCCCAGGT  >  NC_002947/3969207‑3969409                                                                                                                                                 |                                                           ccAGAAGCAAAAATGTTCCTCGAATTGGCCTGATACATAGATGCAACCGGGTTCCTATGGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGAGAGCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTaatg                                                         <  1:172347/148‑4 (MQ=21)          aaaaGCTTCCTCGAATAGGCCTGATACATAGATGCAACCGGGTTCCTATGGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGAttt                                                <  2:34205/148‑1 (MQ=255)          aaaaGCTTCCTCGAATAGGCCTGATACATAGATGCAACCGGGTTCCTATGGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGAttt                                                <  5:60616/148‑1 (MQ=255)                ttCCTCGAATAGGCCTGATACATAGATGCAACCGGGTTCCTATGGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATgc                                          <  2:207840/148‑1 (MQ=255)                   cTCGAATAGGCCTGATACATAGATGCAACCGGGTTCCTATGGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATgcg                                         <  6:186136/146‑1 (MQ=255)                    tCGAATAGGCCTGATACATAGATGCAACCGGGTTCCTATGGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCa                                      >  1:111749/1‑148 (MQ=255)                       aaTAGGCCTGATACATAGATGCAACCGGGTTCCTATGGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCaa                                     <  1:132235/146‑1 (MQ=255)                          aGGCCTGATACATAGATGCAACCGGGTTCCTATGGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCAACGAAg                                <  1:23357/148‑1 (MQ=255)                                             aaCCGGGTTCCTATGGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCAACGAAGCTCTTATCCGAAACCGtt              >  5:141805/1‑147 (MQ=255)                                                cGGGTTCCTATGGTAAGCGTCCGGCGTACTCACCTTCAGAACTCCAGCATTATGGGCGATGGTGTCCGCGTATTTTTAAGCGTACGCGATAACCCTTATCGCCAGGATTTCTATGCGCCAACGAAGCTCTTATCCGAAACCGGTCaaa          <  6:37778/148‑1 (MQ=14)                                                cGGGTTCCTATGGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCAACGAAGCTCTTATCCGAAACCGTTCaaa          <  6:103091/148‑1 (MQ=39)                                                cGGGTTCCTATGGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACACGATCACCCTTATCGCCAGGATTTCTATGCGCCAACGAAGCTCTTATCCGAAACCGTTCaaa          >  2:201946/1‑148 (MQ=255)                                                        tATGGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCAACGAAGCTCTTATCCGAAACCGTTCAAAGCCCAGGt  >  2:157941/1‑148 (MQ=32)                                                                                                                                                 |                                                           CCAGAAGCAAAAAGCTTCCTCGAATAGGCCTGATACATAGATGCAACCGGGTTCCTATGGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATCGCCAGGATTTCTATGCGCCAACGAAGCTCTTATCCGAAACCGTTCAAAGCCCAGGT  >  NC_002947/3969207‑3969409

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