New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3968879	NA (NA)	8 (0.140)	6/248	NT	14.7%	coding (1007/1179 nt)	PP_3498	transposase
	?	NC_002947	3968942 =	49 (0.830)					coding (1070/1179 nt)	PP_3498	transposase

TGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCT‑‑GGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAAC  >  NC_002947/3968831‑3969021                                                 |                                                                                                                                                 tgatgaTCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCT‑‑GGCTAATCGAAATCTGCGGATTATATgg                                            >  6:181363/1‑149 (MQ=32) tgatgaTCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCT‑‑GGCTAATCGAAATCTGCGGATTATATgg                                            >  3:331550/1‑149 (MQ=32)   atgatCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCATGATTCCGATCAGCCAGATGAACGCCGCCTATTGCGATGGCTATCCCGCGTTGGGCGCAAGGAGGCGGCGAGCAGGAT‑‑GGGTAATCGAAATATGCGGATTATATGGGc                                          >  3:92869/1‑149 (MQ=11)       tCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCT‑‑GGCTAATCGAAATCTGCGGATTATATGGGCCCTg                                      <  2:13097/149‑1 (MQ=32)       tCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCT‑‑GGCTAATCGAAATCTGCGGATTATATGGGCCCTg                                      >  7:42536/1‑149 (MQ=32)          aGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCT‑‑GGCTAATCGAAATCTGCGGATTATATGGGCCATGtt                                    >  5:51998/1‑148 (MQ=17)            ggCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCT‑‑GGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTa                                   >  4:172484/1‑147 (MQ=32)            ggCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCT‑‑GGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACa                                 <  8:341820/149‑1 (MQ=32)             gcgcgcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCT‑‑GGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTa                                   >  5:142263/1‑146 (MQ=32)               gcgcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCT‑‑GGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACagag                              <  3:144125/149‑1 (MQ=32)                 gcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGAGGATAATAGAAATGAGGCGATTAAAAGGGGCCAGTTACagag                              >  1:123707/1‑149 (MQ=11)                       cttCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATTATCGCCGCCTCTTGCGCTGGCTATCCGGCCTTGGGCGCAAGGAGGCGGCGATCAGGCT‑‑GGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAg                       <  6:331344/146‑1 (MQ=12)                              aGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCT‑‑GGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTAcc                >  3:335408/1‑148 (MQ=32)                                ccAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCT‑‑GGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAAc             <  6:51998/149‑1 (MQ=32)                                     tAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCT‑‑GGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCaaa        >  7:106836/1‑149 (MQ=32)                                      aaGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCT‑‑GGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAAc       >  4:228065/1‑149 (MQ=32)                                         gCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCT‑‑GGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcg     >  5:345182/1‑148 (MQ=32)                                           cTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCT‑‑GGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaac  <  8:199276/149‑1 (MQ=32)                                            tGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCT‑‑GGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaac  >  3:119591/1‑148 (MQ=32)                                                 |                                                                                                                                                 TGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCT‑‑GGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAAC  >  NC_002947/3968831‑3969021

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