New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AM260479	4024124 =	18 (0.660)	4 (0.160)	4/194	2.0	19.4%	coding (1065/1455 nt)	h16_A3725	signal transduction histidine kinase
	?	AM260479	4024274 =	17 (0.690)					coding (915/1455 nt)	h16_A3725	signal transduction histidine kinase
Rejected: Frequency below/above cutoff threshold.

GGATCGCGTTGTGGACCAGGTTGGCCAGCGCCTCGCGCAGCAGCGCGGGGTCGGCGCGCACCGGCAGCTTCAGCCCCGGCACCGGCTCCCAGCCGAAGTCCTGCTGCTTGCCGCGCGCCAGCGGCAGGTAGTCCAGCGCCACCTGCTCGGCCACGGTGACGGCGTCGATGATGTCGATGCCATCGGCTTCCTGGGCGGCTTCGCCGTGGTGCTGGCGTACCCG  >  AM260479/4024015‑4024237                                                                                                              |                                                                                                                  ggATCGCGTTGTGGACCAGGTTGGCCAGCGCCTCGCGCAGCAGCGCGGGGTCGGCGCGCACCGGCAGCTTCAGCCCCGGCACCGGCTCCCAGCCGAAGTCCTGCTGCTTGCCGCGCGCCAGCGGCAGGTAGTc                                                                                            >  2:57457/1‑133 (MQ=255)      gcgTTGTGGACCAGGTTGGCCAGCGCCTCGCGCAGCAGCGCGGGGTCGGCGCGCACCGGCAGCTTCAGCCCCGGCACCGGCTCCCAGCCGAAGTCCTGCTGCTTGCCGCGCGCCAGCGGCAGGTAGTCCAgcgc                                                                                      >  2:298268/1‑134 (MQ=255)           gtgGACCAGGTTGGCCAGCGCCTCGCGCAGCAGCGCGGGGTCGGCGCGCACCGGCAGCTTCAGCCCCGGCACCGGCTCCCAGCCGAAGTCCTGCTGCTTGCCGCGCGCCAGCGGCAGGTAGTCCAGCGCCAcc                                                                                  <  1:57457/133‑1 (MQ=255)           gtgGACCAGGTTGGCCAGCGCCTCGCGCAGCAGCGCGGGGTCGGCGCGCACCGGCAGCTTCAGCCCCGGCACCGGCTCCCAGCCGAAGTCCTGCTGCTTGCCGCGCGCCAGCGGCAGGTAGTCCAGCGCCAcc                                                                                  <  1:811607/133‑1 (MQ=255)           gtgGACCAGGTTGGCCAGCGCCTCGCGCAGCAGCGCGGGGTCGGCGCGCACCGGCAGCTTCAGCCCCGGCACCGGCTCCCAGCCGAAGTCCTGCTGCTTGCCGCGCGCCAGCGGCAGGTAGTCCAGCGCCACCt                                                                                 >  2:480186/1‑134 (MQ=255)                     cccagcgccagctgcgcgcccctgaccccgagaccgcgcgccaggccgccggcgaagccgtgggcggcatcgtcacccagctggaacagGCCGCGCGCCAGCGGCAGGTAGTCCAGCGCCACCTGCTCGGCCAc                                                                       <  1:17052‑M1/45‑1 (MQ=255)                      ccagcgccagctgcgcgcccctgaccccgagaccgcgcgccaggccgccggcgaagccgtgggcggcatcgtcacccagctggaacagGCCGCGCGCCAGCGGCAGGTAGTCCAGCGCCACCTGCTCGGCCACg                                                                      >  2:17052‑M1/89‑134 (MQ=255)                        gCCAGCGCCTCGCGCAGCAGCGCGGGGTCGGCGCGCACCGGCAGCTTCAGCCCCGGCACCGGCTCCCAGCCGAAGTCCTGCTGCTTGCcgc                                                                                                               >  1:1145924/1‑91 (MQ=255)                        gCCAGCGCCTCGCGCAGCAGCGCGGGGTCGGCGCGCACCGGCAGCTTCAGCCCCGGCACCGGCTCCCAGCCGAAGTCCTGCTGCTTGCcgc                                                                                                               <  2:1145924/91‑1 (MQ=255)                            gcgcCTCGCGCAGCAGCGCGGGGTCGGCGCGCACCGGCAGCTTCAGCCCCGGCACCGGCTCCCAGCCGAAGTCCTGCTGCTTGCCGCGCGCCAGCGGCAGGTAGTCCAGCGCCACCTGCTCGGCCACGGTGACg                                                                >  2:1044362/1‑134 (MQ=255)                                cTCGCGCAGCAGCGCGGGGTCGGCGCGCACCGGCAGCTTCAGCCCCGGCACCGGCTCCCAGCCGAAGTCCTGCTGCTTGCCGCGCGCCAGCGGCAGGTAGTCCAGCGCCACCTGCTCGGCCACGGTGACGGCGt                                                            <  1:1044362/134‑1 (MQ=255)                                                gggTCGGCGCGCACCGGCAGCTTCAGCCCCGGCACCGGCTCCCAGCCGAAGTCCTGCTGCTTGCCGCGCGcc                                                                                                          <  2:175245/72‑1 (MQ=255)                                                gggTCGGCGCGCACCGGCAGCTTCAGCCCCGGCACCGGCTCCCAGCCGAAGTCCTGCTGCTTGCCGCGCGcc                                                                                                          >  1:175245/1‑72 (MQ=255)                                                  gTCGGCGCGCACCGGCAGCTTCAGCCCCGGCACCGGCTCCCAGCCGAAGTCCTGCTGCTTGCCGCGCGCCAGCGGCAGGTAGTCCAGCGCCACCTGCTCGGCCACGGTGACGGCGTCGATGATGTCGATGCCAt                                          >  1:693590/1‑134 (MQ=255)                                                                          ccccGGCACCGGCTCCCAGCCGAAGTCCTGCTGCTTGCCGCGCGCCAGCgg                                                                                                     <  1:568807/51‑1 (MQ=255)                                                                          ccccGGCACCGGCTCCCAGCCGAAGTCCTGCTGCTTGCCGCGCGCCAGCgg                                                                                                     >  2:568807/1‑51 (MQ=255)                                                                          ccccGGCACCGGCTCCCAGCCGAAGTCCTGCTGCTTGCCGCGCGCCAGCGGCAGGTAGTCCAGCGCCACCTGCTCGGCCACGGTGACGGCGTCGATGATGTCGATGCCATCGGCTTCCTGGGCGGCTTCGCCgt                  <  2:693590/134‑1 (MQ=255)                                                                                   cGGCTCCCAGCCGAAGTCCTGCTGCTTGCCGCGCGCCAGCGGCAGGTAGTCCAGCGCCACCTGCTCGGCCACGGTGACGGCGTCGATGATGTCGATGCCATCGGCTTCCTGGGCGGCTTCGCCGTGGTGCTGGc         <  1:298268/134‑1 (MQ=255)                                                                                     gCTCCCAGCCGAAGTCCTGCTGCTTGCCGCGCGCCAGCGGCAGGTAGTCCAGCGCCACCTGCTCGGCCACGGTGACGGCGTCGATGATGTCGATGCCATCGGCTTCCTGGGCGGCTTCGCCGTGGTGCTGGCGt       <  1:474411/134‑1 (MQ=255)                                                                                          cAGCCGAAGTCCTGCTGCTTGCCGCGCGCCAGCGGCAGGTAGTCCAGCGCCACCTGCTCGGCCACGGTGACGGCGTCGATGATGTCGATGCCATCGGCTTCCTGGGCGGCTTCGCCGTGGTGCTGGCGTACCcg  >  2:1065343/1‑134 (MQ=255)                                                                                            gtcacccagctggaacagGCCGCGCGCCAGCGGCAGGTAGTCCAGCGCCACCTGCTCGGCCACGGTGACGGCGt                                                            <  1:781554‑M1/56‑1 (MQ=255)                                                                                            gtcacccagctggaacagGCCGCGCGCCAGCGGCAGGTAGTCCAGCGCCACCTGCTCGGCCACGGTGACGGCGt                                                            >  2:781554‑M1/19‑74 (MQ=255)                                                                                                              |                                                                                                                  GGATCGCGTTGTGGACCAGGTTGGCCAGCGCCTCGCGCAGCAGCGCGGGGTCGGCGCGCACCGGCAGCTTCAGCCCCGGCACCGGCTCCCAGCCGAAGTCCTGCTGCTTGCCGCGCGCCAGCGGCAGGTAGTCCAGCGCCACCTGCTCGGCCACGGTGACGGCGTCGATGATGTCGATGCCATCGGCTTCCTGGGCGGCTTCGCCGTGGTGCTGGCGTACCCG  >  AM260479/4024015‑4024237

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