New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3760851 =	202 (3.050)	6 (0.090)	5/262	NT	3.1%	coding (399/1392 nt)	selA	selenocysteine synthase
	?	NC_000913	3760879 =	184 (2.870)					coding (371/1392 nt)	selA	selenocysteine synthase
Rejected: Frequency below/above cutoff threshold.

TGATCAACCTGACGGGAACCGTGCTGCATACCAACCTTGGGCGAGCTTTACAGGCGGAAGCCGCGGTGGAAGCCGTTGCGCAGGCTATGCGTTCGCCAGTGACCCTCGAGTATGATCTGGACGACGCCGGA‑CGCGGACATCGCGATCGGGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3761002‑3760851 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGC  >  NC_000913/3760879‑3761024                                                                                                                                                                                                                                                                                                             TGATCAACCTGACGGGAACCGTGCTGCATACCAACCTTGGGCGAGCTTTACAGGCGGAAGCCGCGGTGGAAGCCGTTGCGCAGGCTATGCGTTCGCCAGTGACCCTCGAGTATGATCTGGACGACGCCGGA‑CGCGGACATCGCGATCGG                                                                                                                                                       <  2:479498/149‑1  GATCAACCTGACGGGAACCGTGCTGCATACCAACCTTGGGCGAGCTTTACAGGCGGAAGCCGCGGTGGAAGCCGTTGCGCAGGCTATGCGTTCGCCAGTGACCCTCGAGTATGATCTGGACGACGCCGGA‑CGCGGACATCGCGATCGGG                                                                                                                                                      >  3:562143/1‑149  GATCAACCTGACGGGAACCGTGCTGCATACCAACCTTGGGCGAGCTTTACAGGCGGAAGCCGCGGTGGAAGCCGTTGCGCAGGCTATGCGTTCGCCAGTGACCCTCGAGTATGATCTGGACGACGCCGGA‑CGCGGACATCGCGATCGGG                                                                                                                                                      >  3:421293/1‑149   ATCAACCTGACGGGAACCGTGCTGCCTACCAACCTTGGGCGAGCTTTACAGGCGGAAGCCGTTGTGGAAGCCGTTGCGCAGGCTATGCGTTCGCCAGTGACCCTCGAGTATGATCTGGACGACGCCGGA‑CGCGGACATCGCGATCGGGC                                                                                                                                                     >  3:284945/1‑149   ATCAACCTGACGGGAACCGTGCTGCATACCAACCTTGGGCGAGCTTTACAGGCGGAAGCCGCGGTGGAAGCCGTTGCGCAGGCTATGCGTTCGCCATTGACCCTCGAGTATGATCTGGACGACGCCGGA‑CGCGGACATCGCGATCGGGC                                                                                                                                                     <  4:85712/149‑1   ATCAACCTGACGGGAACCGTGCTGCATACCAACCTTGGGCGAGCTTTACAGGCGGAAGCCGCGGTGGAAGCCGTTGCGCAGGCTATGCGTTCGCCAGTGACCCTCGAGTATGATCTGGACGACGCCGGA‑CGCGGACATCGCGATCGGGC                                                                                                                                                     <  2:143352/149‑1   ATCAACCTGACGGGAACCGTGCTGCATACCAACCTTGGGCGAGCTTTACAGGCGGAAGCCGCGGTGGAAGCCGTTGCGCAGGCTATGCGTTCGCCAGTGACCCTCGAGTATGATCTGGACGACGCCGGA‑CGCGGACATCGCGATCGGGC                                                                                                                                                     >  3:194837/1‑149   ATCAACCTGACGGGAACCGTGCTGCATACCAACCTTGGGCGAGCTTTACAGGCGGAAGCCGCGGTGGAAGCCGTTGCGCAGGCTATGCGTTCGCCAGTGACCCTCGAGTATGATCTGGACGACGCCGGA‑CGCGGACATCGCGATCGGGC                                                                                                                                                     >  3:501473/1‑149    TCAACCTGACGGGAACCGTGCTGCATACCAACCTTGGGCGAGCTTTACAGGCGGAAGCCGCGGTGGAAGCCGTTGCGCAGGCTATGCGTTCGCCAGTGACCCTCGAGTATGATCTGGACGACGCCGGA‑CGCGGACATCGCGATCGGGCG                                                                                                                                                    <  4:55393/149‑1    TCAACCTGACGGGAACCGTGCTGCATACCAACCTTGGGCGAGCTTTACAGGCGGAAGCCGCGGTGGAAGCCGTTGCGCAGGCTATGCGTTCGCCAGTGACCCTCGAGTATGATCTGGACGACGCCGGA‑CGCGGACATCGCGATCGGGCG                                                                                                                                                    <  2:11412/149‑1                                 AACCTTGGGCGAGCTTTACAGGCGGAAGCCGCGGTGGAAGCCGTTGCGCAGGCTATGCGTTCGCCAGTGACCCTCGAGTGTGATGTGGACGACGCCCGATCGCAATGTCCGCG‑TCCGGCGTCGTCCAGATCATACTCGAGG                                                                                                                               >  3:411042/1‑141                                 AACCTTGGGCGAGCTTTACAGGCGGAAGCCGCGGTGGAAGCCGTTGCGCAGGCTATGCGTTCGCCAGTGACCCTCGAGTATGATCTGGACGACGCCCGATCGCAATGTCCGCG‑TCCGGCGTCGTCCAGATCATACTCGAGG                                                                                                                               <  4:411042/141‑1                                                                                                                           GACGCCGGA‑CGCGGACATCGCGATCGGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGT                                                     <  3:408759/125‑1                                                                                                                           GACGCCGGA‑CGCGGACATCGCGATCGGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGT                                                     >  4:408759/1‑125                                                                                                                                 GGA‑CGCGGACATCGCGATC‑GGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGTTGATCACC                      >  1:348301/1‑149                                                                                                                                 GGA‑CGCGGACATCGCGATC‑GGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGTTGATCACC                      >  1:533517/1‑149                                                                                                                                                      GGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCG   >  1:576730/1‑149                                                                                                                                                       GCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGC  >  4:480065/1‑149                                                                                                                                                                                                                                                                                                             TGATCAACCTGACGGGAACCGTGCTGCATACCAACCTTGGGCGAGCTTTACAGGCGGAAGCCGCGGTGGAAGCCGTTGCGCAGGCTATGCGTTCGCCAGTGACCCTCGAGTATGATCTGGACGACGCCGGA‑CGCGGACATCGCGATCGGGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3761002‑3760851 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGC  >  NC_000913/3760879‑3761024

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