New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AM260479	= 982724	12 (0.470)	4 (0.170) +ATGTCT	3/192	2.1	28.8%	coding (327/909 nt)	h16_A0889	transcriptional regulator, LysR‑family
	?	AM260479	982748 =	9 (0.360)					coding (303/909 nt)	h16_A0889	transcriptional regulator, LysR‑family

ATCCGGCAGCACGCCCATGCGCAGCGCCACGTCGATGCTGGCATCGACCAGTCCGATCACCTGGTCGGTCAGGCGCAGGTCAAGCGCGATATCGGGGTAGGCGGCAAGGAATTCGGGCAGCGCCGGCGAGACATGCTGGCGCCCGAAGGACGCCGGCGCGGTCACGCGCAACTGCCCGCGCGGGCGCTCGGCGCTGGTGCCGACCGATGCGCTGGCAAGAGCC  >  AM260479/982598‑982820                                                                                                                               |                                                                                                 aTCCGGCAGCACGCCCATGCGCAGCGCCACGTCGATGCTGGCATCGACCAGTCCGATCACCTGGTCGGTCAGGCGCAGGTCAAGCGCGATATCGGGGTAGGCGGCAAGGAATTCGGGCAGCGCCGGCGAGACa                                                                                            <  1:640903/133‑1 (MQ=255) aTCCGGCAGCACGCCCATGCGCAGCGCCACGTCGATGCTGGCATCGACCAGTCCGATCACCTGGTCGGTCAGGCGCAGGTCAAGCGCGATATCGGGGTAGGCGGCAAGGAATTCGGGCAGCGCCGGCGAGACa                                                                                            >  2:640903/1‑133 (MQ=255)                             aCGTCGATGCTGGCATCGACCAGTCCGATCACCTGGTCGGTCAGGCGCAGGTCAAGCGCGATATCGGGGTAGGCGGCAAGGAATTCGGGCAGCGCCGGCGAGACATGCTGGCGCCCGAAGGACGCCGGCGCGGt                                                               >  2:1008531/1‑134 (MQ=255)                               gTCGATGCTGGCATCGACCAGTCCGATCACCTGGTCGGTCAGGCGCAGGTCAAGCGCGATATCGGGGTAGGCGGCAAGGAATTCGGGCAGCGCCGGCGAGACATGCTGGCGCCCGAAGGACGCCGGCGCGGTCa                                                             >  2:345021/1‑134 (MQ=255)                                   aTGCTGGCATCGACCAGTCCGATCACCTGGTCGGTCAGGCGCAGGTCAAGCGCGATATCGGGGTAGGCGGCAAGGAATTCGGGCAGCGCCGGCGAGACATGCTGGCGCCCGAAGGACGCCGGCGCGGTCAcgcg                                                         >  2:351535/1‑134 (MQ=255)                                            tCGACCAGTCCGATCACCTGGTCGGTCAGGCGCAGGTCAAGCGCGATATCGGGGTAGGCGGCAAGGAATTCGGGCAGCGCCGGCGAGACATGCTGGCGCCCGAAGGACGCCGGCGCGGTCACGCGCAACTGccc                                                <  1:345021/134‑1 (MQ=255)                                                      cGATCACCTGGTCGGTCAGGCGCAGGTCAAGCGCGATATCGGGGTAGGCGGCAAGGAATTCGGGCAGCGCCGGCGAGACATGCTGGCGCCCGAAGGACGCCGGCGCGGTCACGCGCAACTGCCCGCGCGGgcg                                       <  1:351535/133‑1 (MQ=255)                                                             cTGGTCGGTCAGGCGCAGGTCAAGCGCGATATCGGGGTAGGCGGCAAGGAATTCGGGCAGCGCCGGCGAGACATGCTGGCGCCCGAAGGACGCCGGCGCGGTCACGCGCAACTGCCCGCGCGGgc                                        <  1:556656/125‑1 (MQ=255)                                                             cTGGTCGGTCAGGCGCAGGTCAAGCGCGATATCGGGGTAGGCGGCAAGGAATTCGGGCAGCGCCGGCGAGACATGCTGGCGCCCGAAGGACGCCGGCGCGGTCACGCGCAACTGCCCGCGCGGgc                                        >  2:556656/1‑125 (MQ=255)                                                                                       cgATATCGGGGTAGGCGGCAAGGAATTCGGGCAGCGCCGGCGAGACATg                                                                                          >  1:611037/1‑49 (MQ=255)                                                                                       cgATATCGGGGTAGGCGGCAAGGAATTCGGGCAGCGCCGGCGAGACATg                                                                                          <  2:611037/49‑1 (MQ=255)                                                                                          tatCGGGGTAGGCGGCAAGGAATTCGGGCAGCGCCGGCGAGACATGCTGGCGCCCGAAGGACGCCGGCGCGGTCACGCGCAACTGCCCGCGCGGGCGCTCGGCGCTGGTGCCGACCGATGCGCTGGCAAGAGcc  <  2:959313/134‑1 (MQ=255)                                                                                             cGGGGTAGGCGGCAAGGAATTCGGGCAGCGCCAGCatgtctcgccggcgcggtcacgcgcaactgcccgcgcgggcgctcg                                                    <  1:1031239‑M1/81‑47 (MQ=255)                                                                                             cGGGGTAGGCGGCAAGGAATTCGGGCAGCGCCAGCatgtctcgccggcgcggtcacgcgcaactgcccgcgcgggcgctcg                                                    >  2:1031239‑M1/1‑35 (MQ=255)                                                                                                                      gggcgccaGCatgtctcgccggcgcggtcacgcgcaactgcccgcgcgggcgctcggcgctggtgccgaccgatgcgctggcaagagccgcggt              <  1:882954‑M1/92‑85 (MQ=255)                                                                                                                      gggcgccaGCatgtctcgccggagcggtcacgcgcaactgcccgcgcgggcgctcggcgctggtgccgaccgatgcgctggcaagagccgcggt              >  2:882954‑M1/3‑10 (MQ=255)                                                                                                                               |                                                                                                 ATCCGGCAGCACGCCCATGCGCAGCGCCACGTCGATGCTGGCATCGACCAGTCCGATCACCTGGTCGGTCAGGCGCAGGTCAAGCGCGATATCGGGGTAGGCGGCAAGGAATTCGGGCAGCGCCGGCGAGACATGCTGGCGCCCGAAGGACGCCGGCGCGGTCACGCGCAACTGCCCGCGCGGGCGCTCGGCGCTGGTGCCGACCGATGCGCTGGCAAGAGCC  >  AM260479/982598‑982820

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.