New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4597603 =	69 (0.900)	4 (0.060)	3/162	NT	6.1%	coding (1839/2292 nt)	opgB	OPG periplasmic biosynthetic phosphoglycerol transferases I (membrane‑bound) and II (soluble)
	?	NC_000913	4597636 =	60 (0.860)					coding (1806/2292 nt)	opgB	OPG periplasmic biosynthetic phosphoglycerol transferases I (membrane‑bound) and II (soluble)
Rejected: Frequency below cutoff threshold.

TTGGCGGTCAGCAAATTGTTCAGCATGTTGACAAAACAACATGGCAGGGCAAAACGGCATTTAAAGATACGGTCATCGACATGGCGCGTTACAAAGGCAATGTCGAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4597709‑4597603 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atgtcgatGACCGTATCTTTAAATGCCGTTTTGCCCTGCCATGTTGTTTTGTCAACATGCTGAACAATTTGCTGACCGCCAAGCTGCCCTTGCGAGACACACCAGTC  >  NC_000913/4597636‑4597734                                                                                                                                                                                                                TTGGCGGTCAGCAAATTGTTCAGCATGTTGACAAAACAACATGGCAGGGCAAAACGGCATTTAAAGATACGGTCATCGACATGGCGCGTTACAAAGGCAA                                                                                                            <  2:92012/100‑1     CGGTCAGCAAATTGTTCAGCATGTTGACAAAACAACATGGCAGGGCAAAACGGCATTTAAAGATACGGTCATCGACATGGCGCGTTACAAAGGCAATGTC                                                                                                        >  2:696524/1‑100     CGGTCAGCAAATTGTTCAGCATGTTGACAAAACAACATGGCAGGGCAAAACGGCATTTAAAGATACGGTCATCGACATGGCGCGTTACAAAGGCAATGTC                                                                                                        >  1:895005/1‑100      GGTCAGCAAATTGTTCAGCATGTTGACAAAACAACATGGCAGGGCAAAACGGCATTTAAAGATACGGTCATCGACATGGCGCGTTACAAAGGCAATGTCG                                                                                                       >  2:886426/1‑100      GGTCAGCAAATTGTTCAGCATGTTGACAAAACAACATGGCAGGGCAAAACGGCATTTAAAGATACGGTCATCGACATGGCGCGTTACAAAGGCAATGTCG                                                                                                       >  2:886416/1‑100        TCAGCAAATTGTTCAGCATGTTGACAAAACAACATGGCAGGGCAAAACGGCATTTAAAGATACGGTCATCGACATGGCGCGTTACAAAGGCAATGTCGAT                                                                                                     >  2:1152748/1‑100        TCAGCAAATTGTTCAGCATGTTGACAAAACAACATGGCAGGGCAAAACGGCATTTAAAGATACGGTCATCGACATGGCGCGTTACAAAGGCAATGTCGAT                                                                                                     <  1:1938342/100‑1                                                               AAAGATACGGTCATCGACATGGCGCGTTACAAAGGCAATGTCGATGACCGTATCTTTAAATGCCGTTTTGCCCTGCCATGTTGTTTTGTCAACATGCTGA                                              >  1:1014758/1‑100                                                                      CGGTCATCGACATGGCGCGTTACAAAGGCAATGTCGATGACCGTATCTTTAAATGCCGTTTTGCCCTGCCATGTTGTTTTGTCAACATGCTGAACAATTT                                       >  1:1306559/1‑100                                                                            TCGACATGGCGCGTTACAAAGGCAATGTCGATGACCGTATCTTTAAATGCCGTTTTGCCCTGCCATGTTGTTTTGTCAACATGCTGAACAATTTGCTGAC                                 >  2:1594666/1‑100                                                                            TCGACATGGCGCGTTACAAAGGCAATGTCGATGACCGTATCTTTAAATGCCGTTTTGCCCTGCCATGTTGTTTTGTCAACATGCTGAACAATTTGCTGAC                                 >  2:2076863/1‑100                                                                                                      GTCGATGACCGTATCTTTAAATGCCGTTTTGCCCTGCCATGTTGTTTTGTCAACATGCTGAACAATTTGCTGACCGCCAAGCTGCCCTTGCGAGACACAC       <  2:1351537/100‑1                                                                                                      GTCGATGACCGTATCTTTAAATGCCGTTTTGCCCTGCCATGTTGTTTTGTCAACATGCTGAACAATTTGCTGACCGCCAAGCTGCCCTTGCGAGACACAC       <  1:1576257/100‑1                                                                                                      GTCGATGACCGTATCTTTAAATGCCGTTTTGCCCTGCCATGTTGTTTTGTCAACATGCTGAACAATTTGCTGACCGCCAAGCTGCCCTTGCGAGACACAC       <  1:1104690/100‑1                                                                                                       TCGATGACCGTATCTTTAAATGCCGTTTTGCCCTGCCATGTTGTTTTGTCAACATGCTGAACAATTTGCTGACCGCCAAGCTGCCCTTGCGAGACACACC      <  1:1762187/100‑1                                                                                                       TCGATGACCGTATCTTTAAATGCCGTTTTGCCCTGCCATGTTGTTTTGTCAACATGCTGAACAATTTGCTGACCGCCAAGCTGCCCTTGCGAGACACACC      <  2:736624/100‑1                                                                                                           TGACCGTATCTTTAAATGCCGTTTTGCCCTGCCATGTTGTTTTGTCAACATGCTGAACAATTTGCTGACCGCCAAGCTGCCCTTGCGAGACACACCAGTC  >  1:1950982/1‑100                                                                                                           TGACCGTATCTTTAAATGCCGTTTTGCCCTGCCATGTTGTTTTGTCAACATGCTGAACAATTTGCTGACCGCCAAGCTGCCCTTGCGAGACACACCAGTC  <  1:859490/100‑1                                                                                                                                                                                                                TTGGCGGTCAGCAAATTGTTCAGCATGTTGACAAAACAACATGGCAGGGCAAAACGGCATTTAAAGATACGGTCATCGACATGGCGCGTTACAAAGGCAATGTCGAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4597709‑4597603 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atgtcgatGACCGTATCTTTAAATGCCGTTTTGCCCTGCCATGTTGTTTTGTCAACATGCTGAACAATTTGCTGACCGCCAAGCTGCCCTTGCGAGACACACCAGTC  >  NC_000913/4597636‑4597734

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