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

CTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGA  >  NC_002947/3968874‑3969153                                                                                                                                         |                                                                                                                                                 gttggtaaacctgatcgctctgtaacagggcccatataatccgcaGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGa                                                                                                                                             <  1:565974‑M1/97‑1 (MQ=255) gttggtaaacctgatcgctctgtaacagggcccatataatccgcaGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGa                                                                                                                                             >  2:565974‑M1/46‑142 (MQ=255) gttggtaaacctgatcgctctgtaacagggcccagttaatccgcaGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGATGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATCTgg                                                                                                                                     >  1:192763‑M1/46‑150 (MQ=255)        tCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGt                                                                                                                              <  1:429250/150‑1 (MQ=32)              gcgcgcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACaga                                                                                                                        <  1:94654/150‑1 (MQ=32)                  gcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAG‑CGa                                                                                                                   <  2:570568/150‑1 (MQ=32)                       ccGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGTCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAgg                                                                                                              <  2:32663/150‑1 (MQ=17)                                ccagataatccgcagatttcgattagccagcctgatcgccgcctccttgcgcccaaggcgggatagccatcgcaagaggcggcgatcatctggctgatcggaatcTGCGGATTAACTGGGCCCTGTTACAGAG‑CGATCAGGTTTACCAAc                                                                                                     >  1:192763‑M2/106‑150 (MQ=255)                                 ccAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAGGTTTACCAACg                                                                                                    >  1:355762/1‑150 (MQ=32)                                      tAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAGGTTTACCAACGCAAAc                                                                                               <  2:579358/150‑1 (MQ=32)                                      tAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTCGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATAAGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAGGTTTACCAAAGCAAAc                                                                                               >  1:450854/1‑150 (MQ=255)                                       aaGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGAGCCCTGTTACAGAG‑CGATCAGGTTTACCAACGAAAAcc                                                                                              >  1:88207/1‑150 (MQ=17)                                        tccgcagatttcgattagccagcctgatcgccgcctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcatctggctgatcggaatcTGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAGGTTTACCAAc                                                                                                     <  1:565974‑M2/45‑1 (MQ=255)                                        tccgcagatttcgattagccagcctgatcgccgcctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcatctggctgatcggaatcTGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAGGTTTACCAAc                                                                                                     >  2:565974‑M2/98‑142 (MQ=255)                                                TTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAGGTTTACCAACGCAAACCGAACGaaaa                                                                                     <  2:156107‑M1/150‑1 (MQ=255)                                                        aGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAGGTTTACCAACGCAACCCGAACGAAAAACCGGAgg                                                                             <  2:337336/150‑1 (MQ=17)                                                                     ccgcctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcaggctggctaatcgaaatcTGCGGa                                                                                                                                             <  1:453681‑M2/6‑1 (MQ=255)                                                                     ccgcctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcaggctggctaatcgaaatcTGCGGa                                                                                                                                             >  2:453681‑M2/69‑74 (MQ=255)                                                                     ccgccTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAGGTTTACCAAc                                                                                                     >  1:416983/1‑113 (MQ=33)                                                                     ccgccTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAGGTTTACCAAc                                                                                                     <  2:416983/113‑1 (MQ=33)                                                                        cctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcaggctggctaatcgaaatcTGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAGGTTTACTAACGCTAACCGAACGAAAAACCGGGAGGGGCTTTGAgcctcgg                                                             >  2:124729‑M2/66‑148 (MQ=255)                                                                          tctTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAG‑CGAt                                                                                                                  <  2:417367/95‑1 (MQ=34)                                                                          tctTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAG‑CGAt                                                                                                                  >  1:417367/1‑95 (MQ=34)                                                                                    gCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCTAAATCTGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGa                                                                               <  1:235190/120‑1 (MQ=17)                                                                                    gCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGa                                                                               >  2:235190/1‑120 (MQ=33)                                                                                      tATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAAt                                               <  1:156107/150‑1 (MQ=32)                                                                                                                  gATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAAccc                   <  2:359625/150‑1 (MQ=32)                                                                                                                  gATCAGGCCGGCTAATCGAAATCTGCGGCATGTATCGGCCCCGGT‑CANCGCCGGTCATGGATACCAACTCAAACTGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAAGCGCTCCAAccc                   >  1:172147/1‑150 (MQ=11)                                                                                                                   aTCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGT‑ATGCCACCGGGGCGCTAGGCGCTCCAACCCCt                 <  2:355762/150‑1 (MQ=17)                                                                                                                      aGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCCATGGCCCCGTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGc               <  1:359056/150‑1 (MQ=11)                                                                                                                        gctggctAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTa             >  2:407888/1‑150 (MQ=32)                                                                                                                                   aaaaTCTGTGGCCTATATGGGCGCTGTTACAGAG‑CGATCAAATATCCCAACGCACACCGAGCGAAAAACCGGAGGCTGCTATGAGCCTCTGATCAACCGGGTAATGCCACCGTGGCGCTAGGCGCTCCAGCCCCTGCTAGTGAACATTGa  <  2:172147/149‑1 (MQ=11)                                                                                                                                         |                                                                                                                                                 CTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAG‑CGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGA  >  NC_002947/3968874‑3969153

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.