New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AM260480	= 2824906	9 (0.330)	4 (0.160)	4/198	2.0	24.8%	coding (1105/1584 nt)	h16_B2500	sulfite reductase [NADPH] flavoprotein alpha‑component
	?	AM260480	= 2824939	16 (0.630)					coding (1072/1584 nt)	h16_B2500	sulfite reductase [NADPH] flavoprotein alpha‑component

CGCAGCCCGGCCAGCCCGGTGCCGTTGCCGACCAGGATCAGGGGCCGGTCATCGGCCGGCGGGTGGAAGCTGCGGTTGGTGCGGATGCGCAGCGCGATGCGGGTGCCCACGGCCGCGTGTTCGGTCAGCCAGCCGGAGGCCAGGCCCAGGCGGCCGTCGTCGTAGCGGGTCTGCCGCACCAGCAGTTCGAGTTGCCCGTCCCGCGGCAGCGAGGC  >  AM260480/2824791‑2825005                                                                                                                    |                                                                                                    cgcAGCCCGGCCAGCCCGGTGCCGTTGCCGACCAGGATCAGGGGCCGGTCATCGGCCGGCGGGTGGAAGCTGCGGTTGGTGCGGATGCGCAGCGCGATGCGGGTGCCCACGGCCGCGTGTTCGGTCagccagcc                                                                                   >  2:603873/1‑134 (MQ=255)                           gCCGACCAGGATCAGGGGCCGGTCATCGGCCGGCGGGTGGAAGCTGCGGTTGGTGCGGATGCGCAGCGCGATGCGGGTGCCCACGGCCGCGTGTTCGGTCAGCCAGCCGGAGGCCAGGCCCAGGCGGCcgtcg                                                          <  1:603873/133‑1 (MQ=255)                                              cGGTCATCGGCCGGCGGGTGGAAGCTGCGGTTGGTGCGGATGCGCAGCGCGATGCGGGTGCCCACGGCCGCGTGTTCGGTCAGCCAGCCGGAGGCCAGGCCCAGGCGGCCGTCGTCGTAGCGGGTCTGCCGCAc                                      >  1:452610/1‑134 (MQ=255)                                                   aTCGGCCGGCGGGTGGAAGCTGCGGTTGGTGCGGATGCGCAGCGCGATGCGGGTGCCCACGGCCGCctgggcctggcctccggctggctgaccgaacacgcgg                                                                <  2:543981‑M1/103‑38 (MQ=255)                                                   aTCGGCCGGCGGGTGGAAGCTGCGGTTGGTGCGGATGCGCAGCGCGATGCGGGTGCCCACGGCCGCctgggcctggcctccggctggctgaccgaacacgcgg                                                                >  1:543981‑M1/1‑66 (MQ=255)                                                   aTCGGCCGGCGGGTGGAAGCTGCGGTTGGTGCGGATGCGCAGCGCGATGCGGGTGCCCACGGCCGCGTGTTCGGTCAGCCAGCCGGAGGCCAGGCCCAGGCGGCCGTCGTCGTAGCGGGTCt                                             <  1:199730/122‑1 (MQ=255)                                                   aTCGGCCGGCGGGTGGAAGCTGCGGTTGGTGCGGATGCGCAGCGCGATGCGGGTGCCCACGGCCGCGTGTTCGGTCAGCCAGCCGGAGGCCAGGCCCAGGCGGCCGTCGTCGTAGCGGGTCt                                             >  2:199730/1‑122 (MQ=255)                                                           gcggGTGGAAGCTGCGGTTGGTGCGGATGCGCAGCGCGATGCGGGTGCCCACGGCCGCGTGTTCGGTCAGCCAGCCGGAGGCCAGGCCCAGGCGGCCGTCGTCGTAGCGGGTCTGCCGCACCAGCAGTTCGAGt                         <  2:452610/134‑1 (MQ=255)                                                                ttgAAGCTGCGGTTGGTGCGGATGCGCAGCGCGATGCGGGTGCCCACGGCCGCGTGTTCGGTCAGCCAGCCGGAGGCCAGGCCCAGGCGGCCGTCGTCGTAGCGGGTCTGCCGCACCAGCAGTTCGAGTTGccc                    <  2:956529/132‑1 (MQ=255)                                                                             tGGTGCGGATGCGCAGCGCGATGCGGGTGCCCACGGCCGCGTGTTCGGTCAGCCAGCCGGAGGCCAGGCCCAGgcggccgtgg                                                          <  2:453719‑M2/83‑11 (MQ=255)                                                                             tGGTGCGGATGCGCAGCGCGATGCGGGTGCCCACGGCCGCGTGTTCGGTCAGCCAGCCGGAGGCCAGGCCCAGgcggccgtgg                                                          >  1:453719‑M2/1‑73 (MQ=255)                                                                                                           ccACGGCCGCctgggcctggcctccggctggctgaccgaacacgcggccgtgggcacccgcatcgcgctgcgcatccgcacca                            <  2:453719‑M1/83‑74 (MQ=255)                                                                                                           ccACGGCCGCctgggcctggcctccggctggctgaccgaacacgcggccgtgggcacccgcatcgcgctgcgcatccgcacca                            >  1:453719‑M1/1‑10 (MQ=255)                                                                                                           ccACGGCCGCGTGTTCGGTCAGCCAGCCGGAGGCCAGGCCCAGGCGGCCGTCGTCGTAGCGGGTCTGCCGCACCAGAAGTTCGAGTTGCCCGt                  <  2:333634/93‑1 (MQ=255)                                                                                                           ccACGGCCGCGTGTTCGGTCAGCCAGCCGGAGGCCAGGCCCAGGCGGCCGTCGTCGTAGCGGGTCTGCCGCACCAGAAGTTCGAGTTGCCCGt                  >  1:333634/1‑93 (MQ=255)                                                                                                                 ccGCGTGTTCGGTCAGCCAGCCGGAGGCCAGGCCCAGgcggccgtgggcacccgcatcgcgctgcgcatccgcaccaaccgcagcttccacccgccggccgat  >  1:543981‑M2/1‑37 (MQ=255)                                                                                                                 ccGCGTGTTCGGTCAGCCAGCCGGAGGCCAGGCCCAGgcggccgtgggcacccgcatcgcgctgcgcatccgcaccaaccgcagcttccacccgccggccgat  <  2:543981‑M2/103‑67 (MQ=255)                                                                                                                    |                                                                                                    CGCAGCCCGGCCAGCCCGGTGCCGTTGCCGACCAGGATCAGGGGCCGGTCATCGGCCGGCGGGTGGAAGCTGCGGTTGGTGCGGATGCGCAGCGCGATGCGGGTGCCCACGGCCGCGTGTTCGGTCAGCCAGCCGGAGGCCAGGCCCAGGCGGCCGTCGTCGTAGCGGGTCTGCCGCACCAGCAGTTCGAGTTGCCCGTCCCGCGGCAGCGAGGC  >  AM260480/2824791‑2825005

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