New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	3968919 =	8 (0.320)	6 (0.250)	5/256	2.4	28.7%	coding (1003/1179 nt)	PP_3498	transposase
	?	NC_002947	3969010 =	22 (0.910)					coding (1094/1179 nt)	PP_3498	transposase

GGGCAAGATGAGCAAGCGCGGTGACGGCTACATACGCAGCCTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACC  >  NC_002947/3968834‑3969061                                                                                      |                                                                                                                                                ccgaggctcaaagcccctcccggtttttcgttcggttagcgttagtaaacctgatcgctctgtaacagggcccatataatccgcaGATTTCGATTAGCCAGCCTGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAaggagg                                                                                 >  2:124729‑M1/86‑150 (MQ=255)                                         gttggtaaacctgatcgctctgtaacagggcccatataatccgcaGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGa                                                <  1:565974‑M1/97‑1 (MQ=255)                                         gttggtaaacctgatcgctctgtaacagggcccatataatccgcaGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGa                                                >  2:565974‑M1/46‑142 (MQ=255)                                         gttggtaaacctgatcgctctgtaacagggcccagttaatccgcaGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGATGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTgg                                        >  1:192763‑M1/46‑150 (MQ=255)                                                tCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGt                                 <  1:429250/150‑1 (MQ=32)                                                      gcgcgcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACaga                           <  1:94654/150‑1 (MQ=32)                                                          gcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGa                       <  2:570568/150‑1 (MQ=32)                                                               ccGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGTCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAgg                  <  2:32663/150‑1 (MQ=17)                                                                         ccAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACg        >  1:355762/1‑150 (MQ=32)                                                                              tAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAAc   <  2:579358/150‑1 (MQ=32)                                                                              tAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTCGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATAAGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAAAGCAAAc   >  1:450854/1‑150 (MQ=255)                                                                               aaGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGAGCCCTGTTACAGAGCGATCAGGTTTACCAACGAAAAcc  >  1:88207/1‑150 (MQ=17)                                                                                tccgcaGATTTCGATTAGCCAGCCTGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAggcgg                                                                              <  1:453681‑M1/68‑1 (MQ=255)                                                                                tccgcaGATTTCGATTAGCCAGCCTGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAggcgg                                                                              >  2:453681‑M1/7‑74 (MQ=255)                                                                                      |                                                                                                                                                GGGCAAGATGAGCAAGCGCGGTGACGGCTACATACGCAGCCTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGA‑TGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACC  >  NC_002947/3968834‑3969061

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