New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	850440 =	43 (0.490)	4 (0.050)	5/238	NT	8.6%	intergenic (‑20/‑785)	pcs/hemA	phosphatidylcholine synthase/glutamyl‑tRNA reductase
	?	NC_002947	850461 =	44 (0.520)					intergenic (‑41/‑764)	pcs/hemA	phosphatidylcholine synthase/glutamyl‑tRNA reductase

TTGATGAGGTTTAGCAGCGGTAATGATGAACTGGAATGATGCACACTTACTTGGTGACACACCAACAAGAAGCATTTTTATATTTCAGTTTTGTTAAGCTTATCTTAAAAACGATTTCCAATGCAGGGTATGCTTCGTG‑CATGGCAGGAAGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/850591‑850440 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gaagcATACCCTGCATTGGAAATCGTTT  >  NC_002947/850461‑850483                                                                                                                                                                                  TTGATGAGGTTTAGCAGCGGTAATGATGAACTGGAATGATGCACACTTACTTGGTGACACACCAACAAGAAGCATTTCTCTACCTCCGTACTGTTATTCTTTTCTTAAAAACGATTTCCAATGCAGGGTATGCTTCGTT‑CATGGCAGGAA                           >  1:862081/1‑150   GATGAGGTTTAGCAGCGGTAATGATGAACTGGAATGATGCACACTTACTTGGTGACACACCAACAAGAAGCATTTTCATATTTCAGTTTTGTTAAGCTTATCTTAAAAACGATTTCCAATGCAGGGTATGCTTCGTG‑CATGGCAGGAAGC                         >  1:152096/1‑150       AGGTTTAGCAGCGGTAATGATGAACTGGAATGATGCACACTTACTTGGTGACACACCAACAAGAAGCATTTTTATATTTCAGTTTTGTTAAGCTTATCTTAAAAACGATTTCCAATGCAGGGTATGCTTCCTGCCAT‑GCACGAAGCATAC                     <  2:663758/150‑1              GCAGCGGTAATGATGAACTGGAATGATGCACACTTACTTGGTGACACACCAACAAGAAGCATTTTTATATTTCAGTTTTGTTAAGCTTATCTTAAAAACGATTTCCAATGCAGGGTATGCTTCCTGCCAT‑GCACGAAGCATACCCTGCAT              >  2:891846/1‑150                          ATGAACTGGAATGATGCACACTTACTTGGTGACACACCAACAAGAAGCATTTTTATATTTCAGTTTTGTTAAGCTTATCTTAAAAACGATTTCCAATGCAGGGTATGCTTCCTGCCAT‑GCACGAAGCATACCCTGCATTGGAAATCGTTT  <  1:891846/150‑1                                  GAATGATGCACACTTACTTGGTGACACACCAACAAGAAGCATTTTTATATTTCAGTTTTGTTACGCTTATCTTAAAAACGATTTCCAATGCAGGGTATGCTTCCTGCCAT‑GCACGAAGCATACCCTGCATTGG           <  1:754693/133‑1                                  GAATGATGCACACTTACTTGGTGACACACCAACAAGAAGCATTTTTATATTTCAGTTTTGTTAAGCTTATCTTAAAAACGATTTCCAATGCAGGGTATGCTTCCTGCCAT‑GCACGAAGCATACCCTGCATTGG           >  2:754693/1‑133                                                                                                                                                                                  TTGATGAGGTTTAGCAGCGGTAATGATGAACTGGAATGATGCACACTTACTTGGTGACACACCAACAAGAAGCATTTTTATATTTCAGTTTTGTTAAGCTTATCTTAAAAACGATTTCCAATGCAGGGTATGCTTCGTG‑CATGGCAGGAAGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/850591‑850440 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gaagcATACCCTGCATTGGAAATCGTTT  >  NC_002947/850461‑850483

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