New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AM260480	= 2485169	7 (0.250)	3 (0.120)	3/196	2.3	25.6%	coding (1897/2424 nt)	h16_B2195	signal transduction histidine kinase containing CHASE2 sensor domain
	?	AM260480	= 2485189	11 (0.440)					coding (1877/2424 nt)	h16_B2195	signal transduction histidine kinase containing CHASE2 sensor domain

GTCGCATCCAGCACCAGGCCACCCAATTCCACTTCCATGAACGCCAGGGTCTTCGACTCCGCGCGCGCCACGCTGATGAAGTCGTCGGCCAGCGTCATCGTGCGGCGTGCATGGGCGGCGATGCGTGCCAGCACGCCGCTGTCCGTGGCGACGTCGCCGGCGGCTGCCCCATCGCGCCGCGCGCGGCGCTCATGCAGCTCGATCAAGGCCAGGATCGATGCCTGCGGCGAGCGCATGTCG  >  AM260480/2485060‑2485299                                                                                                              |                                                                                                                                   gTCGCATCCAGCACCAGGCCACCCAATTCCACTTCCATGAACGCCAGGGTCTTCGACTCCGCGCGCGCCACGCTGATGAAGTCGTCGGCCAGCGTCATCGTGCGGCGTGCATGGGCGGCGATGCGTGCCAgcac                                                                                                            >  2:918046‑M2/1‑130 (MQ=255)                  gCCACCCAATTCCACTTCCATGAACGCCAGGGTCTTCGACTCCGCGCGCGCCACGCTGATGAAGTCGTCGGCCAGCGTCATCGTGCGGCGTGCATGGGCGGCGATGCGTGCCAgcacgccgcacgatgacgctg                                                                                           >  2:18955‑M2/1‑113 (MQ=255)                                          cGCCAGGGTCTTCGACTCCGCGCGCGCCACGCTGATGAAGTCGTCGGCCAGCGTCATCGTGCGGCGTGCATGGGCGGCGATGCGTGCCAGCACGCCGCTGTCCGTGGCGACGTCGCCGGCGGCTGCCCCATcgc                                                                   >  1:247438/1‑134 (MQ=255)                                                 gTCTTCGACTCCGCGCGCGCCACGCTGATGAAGTCGTCGGCCAGCGTCATCGTGCGGCGTGCtggcacgcatcgccgcccatgcacgccg                                                                                                        >  1:48232‑M1/1‑62 (MQ=255)                                                 gTCTTCGACTCCGCGCGCGCCACGCTGATGAAGTCGTCGGCCAGCGTCATCGTGCGGCGTGCtggcacgcatcgccgcccatgcacgccg                                                                                                        <  2:48232‑M1/90‑29 (MQ=255)                                                 gTCTTCGACTCCGCGCGCGCCACGCTGATGAAGTCGTCGGCCAGCGTCATCGTGCGGCGTGCATGGGCGGCGATGCGTGCCAGCACGCCGCTGTCCGTGGCGACGTCGCCGGCGGCTGCCCCATcgcgccgcgc                                                            >  1:482982/1‑134 (MQ=255)                                                                            atgaAGTCGTCGGCCAGCGTCATCGTGCGGCGTGCATGGGCGGCGATGCGTGCCAGCACGCCGCTGTCCGTGGCGACGTCGCCGGCGGCTGCCCCATCGCGCCGCGCGCGGCGCTCATGCAGCTCGATCAAg                                   <  2:823804/132‑1 (MQ=255)                                                                               aaGTCGTCGGCCAGCGTCATCGTGCGGCGTGCATGGGCGGCGATGCGTGCCAGCACGCCGCTGTCCGTGGCGACGTCGCCGGCGGCTGCCCCATCGCGCCGCGCGCGGCGCTCATGCAGCTCGATCAAGGCCAg                              <  2:247438/134‑1 (MQ=255)                                                                                        gCCAGCGTCATCGTGCGGCGTGCATGGGCGGCGATGCGTGCCAGCACGCCGCTGTCCGTGGCGACGTCGCCGGCGGCTGCCCCATCGCGCCGCGCGCGGCGCTCATGCAGCTCGATCAAGGCCAGGATCGATGc                     <  2:230740/134‑1 (MQ=255)                                                                                          cAGCGTCATCGTGCGGCGTGCtggcacgcatcgccgcccatgcacgccgcacgatgacgctggccgacgacttcatcagcgtggcgcgcgcggagtcgaagaccctggcgttcatggaagtggaattgggtggc                   >  2:18955‑M1/1‑21 (MQ=255)                                                                                                       cggcgTGCATGGGCGGCGATGCGTGCCAgcacgccgcacgatgacgctggccgacgacttcatcagcgtggcgcgcgcggagtcgaagac                                                  >  1:48232‑M2/1‑28 (MQ=255)                                                                                                       cggcgTGCATGGGCGGCGATGCGTGCCAgcacgccgcacgatgacgctggccgacgacttcatcagcgtggcgcgcgcggagtcgaagac                                                  <  2:48232‑M2/90‑63 (MQ=255)                                                                                                           gTGCtggcacgcatcgccgcccatgcacgccgcacgatgacgctggccgacgacttcatcagcgtggcgcgcgcggagtcgaagaccctggcgttcatggaagtggaattgggtggcctggtgctggatgcgac  >  2:918046‑M1/1‑4 (MQ=255)                                                                                                           gTGCATGGGCGGCGATGCGTGCCAGCACGCCGCTGTCCGTGGCGACGTCGCCGGCGGCTGCCCCATCGCGCCGCGCGCGGCGCTCATGCAGCTc                                          >  1:321808/1‑94 (MQ=255)                                                                                                           gTGCATGGGCGGCGATGCGTGCCAGCACGCCGCTGTCCGTGGCGACGTCGCCGGCGGCTGCCCCATCGCGCCGCGCGCGGCGCTCATGCAGCTc                                          <  2:321808/94‑1 (MQ=255)                                                                                                              |                                                                                                                                   GTCGCATCCAGCACCAGGCCACCCAATTCCACTTCCATGAACGCCAGGGTCTTCGACTCCGCGCGCGCCACGCTGATGAAGTCGTCGGCCAGCGTCATCGTGCGGCGTGCATGGGCGGCGATGCGTGCCAGCACGCCGCTGTCCGTGGCGACGTCGCCGGCGGCTGCCCCATCGCGCCGCGCGCGGCGCTCATGCAGCTCGATCAAGGCCAGGATCGATGCCTGCGGCGAGCGCATGTCG  >  AM260480/2485060‑2485299

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