New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4467477 =	137 (1.070)	8 (0.070)	7/162	NT	6.1%	coding (47/54 nt)	mgtL	regulatory leader peptide for mgtA
	?	NC_000913	4467511 =	122 (1.050)					intergenic (+27/‑114)	mgtL/mgtA	regulatory leader peptide for mgtA/magnesium transporter
Rejected: Frequency below cutoff threshold.

TCAGGCGTGGCGAAGCCTTACCTGAAGGTAATAAAGCTGAGCAGGGACGTTAAAACGTCACGGCCTTACGCATCCGGTAACGCCGTGAAAAGACGGGCTTACCGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4467582‑4467477 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttaccggaTGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGT  >  NC_000913/4467511‑4467609                                                                                                                                                                                                               TCAGGCGTGGCGAAGCCTTACCTGAAGGTAATAAAGCTGAGCAGGGACGTTAAAACGTCACGGCCTTACGCATCCGGTAACGCCGTGAAAAGACGGGCTT                                                                                                           >  2:2215491/1‑100   AGGCGTGGCGAAGCCTTACCTGAAGGTAATAAAGCTGAGCAGGGACGTTAAAACGTCACGGCCTTACGCATCCGGTAACGCCGTGAAAAGACGGGCTTAC                                                                                                         <  1:1616319/100‑1   AGGCGTGGCGAAGCCTTACCTGAAGGTAATAAAGCTGAGCAGGGACGTTAAAACGTCACGGCCTTACGCATCCGGTAACGCCGTGAAAAGACGGGCTTAC                                                                                                         <  1:3399799/100‑1    GGCGTGGCGAAGCCTTACCTGAAGGTAATAAAGCTGAGCAGGGACGTTAAAACGTCACGGCCTTACGCATCCGGTAACGCCGTGAAAAGACGGGCTTACC                                                                                                        >  2:962039/1‑100    GGCGTGGCGAAGCCTTACCTGAAGGTAATAAAGCTGAGCAGGGACGTTAAAACGTCACGGCCTTACGCATCCGGTAACGCCGTGAAAAGACGGGCTTACC                                                                                                        >  2:1970405/1‑100      CGTGGCGAAGCCTTACCTGAAGGTAATAAAGCTGAGCAGGGACGTTAAAACGTCACGGCCTTACGCATCCGGTAACGCCGTGAAAAGACGGGCTTACCGG                                                                                                      <  1:1444948/100‑1                                                            ACGGCCTTACGCATCCGGTAACGCCGTGAAAAGACGGGCTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGT                                                >  1:742114/1‑100                                                                 CTTACGCATCCGGTAACGCCGTGAAAAGACGGGCTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGC                                           >  1:2405288/1‑100                                                                   TACGCATCCGGTAACGCCGTGAAAAGACGGGCTTACCGGATGCGTAAGGCCGTGACGTCTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTT                                         >  1:353544/1‑100                                                                       CATCCGGTAACGCCGTGAAAAGACGGGCTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCC                                     >  2:2810046/1‑100                                                                         TCCGGTAACGCCGTGAAAAGACGGGCTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCAC                                   >  2:2855933/1‑100                                                                           CGGTAACGCCGTGAAAAGACGGGCTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGC                                 >  2:877377/1‑100                                                                           CGGTAACGCCGTGAAAAGACGGGCTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGC                                 >  1:1772476/1‑100                                                                                 CGCCGTGAAAAGACGGGCTTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCGACGCCTGAA                            >  1:364203/1‑99                                                                                                   TTACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTC         <  2:2907771/100‑1                                                                                                    TACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCA        <  2:3435841/100‑1                                                                                                    TACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCA        <  2:2032896/100‑1                                                                                                     ACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAG       >  2:458121/1‑100                                                                                                     ACCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAG       <  2:1391515/100‑1                                                                                                      CCGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCATGAAGACATTTCTGTACTGTTTCAG       >  1:3138917/1‑99                                                                                                       CGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGAC     <  1:2269719/100‑1                                                                                                       CGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGAC     >  2:1146712/1‑100                                                                                                       CGGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGAC     <  1:3416130/100‑1                                                                                                        GGATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACA    >  1:140606/1‑100                                                                                                          ATGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGT  <  2:1758073/100‑1                                                                                                                                                                                                               TCAGGCGTGGCGAAGCCTTACCTGAAGGTAATAAAGCTGAGCAGGGACGTTAAAACGTCACGGCCTTACGCATCCGGTAACGCCGTGAAAAGACGGGCTTACCGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4467582‑4467477 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttaccggaTGCGTAAGGCCGTGACGTTTTAACGTCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGT  >  NC_000913/4467511‑4467609

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