New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4647603 =	18 (0.570)	4 (0.180)	4/170	NT	29.3%	intergenic (+10/+121)	fusB/PP_4112	elongation factor G 2/sulfate transporter
	?	NC_002947	4647711 =	7 (0.320)					intergenic (+118/+13)	fusB/PP_4112	elongation factor G 2/sulfate transporter

TGAAATGTTCGGCTACGCCACCTCGATGCGCTCGATGACCCAGGGCCGGGCCAGCTTCTCGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGACCGCAACAGTTT  >  NC_002947/4647458‑4647723                                                                                                                                                  |                                                                                                                         gtcgattattcggccgttgcggcgctgaaaaccctgcgcgagcgctatgccaaggcgggcaagcatttgcgcgtggtgcatttgtcggagcggtgcaagaagctgttgaagcgggcgggcgagcagcattgaaaactgttgcggtCtgtg                                                                                                                      >  2:292530‑M1/146‑150 (MQ=255)  gAAATGTTCGGCTACGCCACCTCGATGCGCTCGATGACCCAGGGCCGGGCCAGCTTCTCGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGc                                                                                                                     >  1:378027/1‑150 (MQ=255)     attattcggccgttgcggcgctgaaaaccctgcgcgagcgctatgccaaggcgggcaagcatttgcgcgtggtgcatttgtcggagcggtgcaagaagctgttgaagcgggcgggcgagcagcattgaaaactgttgcggtCTGTGCCgg                                                                                                                  <  1:489059‑M1/9‑1 (MQ=255)           ggCTACGCCACCTCGATGCGCTCGATGACCCAGGGCCGGGCCAGCTTCTCGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCtt                                                                                                            >  2:415554/1‑150 (MQ=255)           ggCTACGCCACCTCGATGCGCTCGATGACCCAGGGCCGGGCCAGCTTCTCGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCtt                                                                                                            <  1:601888/150‑1 (MQ=255)                gttgcggcgctgaaaaccctgcgcgagcgctatgccaaggcgggcaagcatttgcgcgtggtgcatttgtcggagcggtgcaagaagctgttgaagcgggcgggcgagcagcattgaaaactgttgcggtCTGTGCCGGCCTCTTCGCgg                                                                                                       >  1:824771‑M1/131‑150 (MQ=255)                cGCCACCTCGATGCGCTCGATGACCCAGGGCCGGGCCAGCTTCTCGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCgg                                                                                                       <  2:40123/150‑1 (MQ=255)                               ctctATGACCCAGGGCCGGGCCAGCTTCTCGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGAGGGTAAACGCGCTCCTa                                                                                        >  1:90598/1‑150 (MQ=255)                                  gATGACCCAGGGCCGGGCCAGCTTCTCGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAg                                                                                     <  1:61746/150‑1 (MQ=255)                                  gATGACCCAGGGACGGGCACACTTTTCCGGACACTTCTCCCTCTATGCCTAGTCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCTAccg                                                                                     <  2:520576/150‑3 (MQ=25)                                     gACCCAGGGCCGGGCCAGCTTCTCGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAt                                                                                  <  2:378027/150‑1 (MQ=255)                                      aCCCAGGGCCGGGCCAGCTTCTCGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATg                                                                                 <  1:415554/150‑1 (MQ=255)                                          gcgctatgccaaggcgggcaagcatttgcgcgtggtgcatttgtcggagcggtgcaagaagctgttgaagcgggtgggcgagcagcattgaaaactgttgcggtCTGTGCCGGCCTCTTcgc                                                                                                         <  1:376877‑M1/18‑1 (MQ=255)                                          gcgctatgccaaggcgggcaagcatttgcgcgtggtgcatttgtcggagcggtgcaagaagctgttgaagcgggtgggcgagcagcattgaaaactgttgcggtCTGTGCCGGCCTCTTcgc                                                                                                         >  2:376877‑M1/105‑122 (MQ=255)                                             gccgggccAGCTTCTCGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAg                                                                          >  1:18625/1‑150 (MQ=255)                                                           tcGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATtgtg                                                            <  2:18625/150‑1 (MQ=255)                                                              aTGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTgcgc                                                         <  1:356436/150‑1 (MQ=255)                                                                    ttCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCt                                                   >  1:496267/1‑150 (MQ=255)                                                                                  cGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAg                                                                                     >  1:693571/1‑102 (MQ=255)                                                                                  cGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAg                                                                                     <  2:693571/102‑1 (MQ=255)                                                                                                 cATTGCCGATGGCATCGTCAAGACGAGCCGCCGGGACTTAGCGTTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGcc                      <  2:496267/150‑1 (MQ=39)                                                                                                  aTTGCCGAGGGCATGGTCAAGAAGAGCCGCGGGGAGTAAGCCGTTCTGCTGTCCCGGCCTCTTCGCGGCTAAACCCGCTCCTACCGGATTGCCGCAGGCCCAAACAGTGTGTGGCCCACTGTAAGAGCGGGTTTACTCGCGAAGAGGCCg                     <  2:90598/150‑1 (MQ=255)                                                                                                                     agaagaGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCACGGGTAAACCCCCTCCTACAGGCTGGCCGCACGCCCAAACATTGTGCGCCCCCCTGTAGGAGCGGGATTACCGGCGAAGAGGCCGGCACAGACCGCAACAGtct  >  2:213122/1‑148 (MQ=38)                                                                                                                         gagcagcattgaaaactgttgcggtCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCCAACATTGTGCGCCCCCCt                                                   >  2:652658‑M1/26‑97 (MQ=255)                                                                                                                         gagcagcattgaaaactgttgcggtCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCt                                                   <  1:652658‑M1/72‑1 (MQ=255)                                                                                                                            cagcattgaaaactgttgcggtCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGc                                           >  1:498311‑M1/23‑102 (MQ=255)                                                                                                                            cagcattgaaaactgttgcggtCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGc                                           <  2:498311‑M1/80‑1 (MQ=255)                                                                                                                                             ttCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGCGGGt                                       <  2:305943/89‑1 (MQ=255)                                                                                                                                             ttCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGCGGGt                                       >  1:305943/1‑89 (MQ=255)                                                                                                                                                  CTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGagcgggatt                                     <  1:101088‑M1/86‑4 (MQ=255)                                                                                                                                                  CTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGCGGGttt                                     >  2:101088‑M1/1‑86 (MQ=255)                                                                                                                                                  |                                                                                                                         TGAAATGTTCGGCTACGCCACCTCGATGCGCTCGATGACCCAGGGCCGGGCCAGCTTCTCGATGGAGTTCACCCGCTATGCCGAGGCCCCGGCGAGCATTGCCGATGGCATCGTCAAGAAGAGCCGCGGGGAGTAAGCGGTTCTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGATGGCCGCAGGCCCAAACATTGTGCGCCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGACCGCAACAGTTT  >  NC_002947/4647458‑4647723

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.