New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AM260480	2564349 =	11 (0.400)	3 (0.130)	3/208	1.4	25.3%	coding (1361/3738 nt)	narG2	respiratory nitrate reductase alpha chain (NarG)
	?	AY305378	288727 =	7 (0.340)					coding (1384/3762 nt)	narG	respiratory nitrate reductase catalytic subunit

TCATCGTCGTACGAGGCGGCGACGTTCGGCCCGCCGAGTCCCTGATCCACGCCGTAGTTGGCCATCTGCAGGTCATACACGGTCGCCACCAGCGCAGTGCTGCCATCGGCCAGCGTCAGCCGGCGTGCCGGCACGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AM260480/2564484‑2564349 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgcCGGCTGACGCTGGCCGACGGCAGCACCGCGCTGGTGGCCACCGTGTACGACCTGCA  >  AY305378/288727‑288782                                                                                                                                                                                                  TCATCGTCGTACGAGGCGGCGACGTTCGGCCCGCCGAGTCCCTGATCCACGCCGTAGTTGGCCATCTGCAGGTCATACACGGTCGCCACCAGCGCAGTGCTGCCATCGGCCAGCGTCAGCCGGCGTGCCGGCAC                                                            >  2:257341/1‑134              AGGCGGCGACGTTCGGCCCGCCGAGTCCCTGATCCACGCCGTAGTTGGCCATCTGCAGGTCATACACGGTCGCCACCAGCGCAGTGCTGCCATCGGCCAGCGTCAGCCGGCGTGCCGGCACGCCGGCTGACGCT                                               >  2:239491/1‑134                         TTCGGCCCGCCGAGTCCCTGATCCACGCCGTAGTTGGCCATCTGCAGGTCATACACGGTCGCCACCAGCGCAGTGCTGCCATCGGCCAGCGTCAGCCGGCGTGCCGGCACGCCGGCTGACGCTGGCCGATGGCA                                    <  2:151487/134‑1                                                            GGCCATCTGCAGGTCATACACGGTCGCCACCAGCGCAGTGCTGCCATCGGCCAGCGTCAGCCGGCGTGCCGGCACGCCGGCTGACGCTGGCCGATGGCAGCACTGCGCTGGTGGCGACCGTGTATGACCTGCA  <  1:239491/133‑1                                                                                                                                                                                                  TCATCGTCGTACGAGGCGGCGACGTTCGGCCCGCCGAGTCCCTGATCCACGCCGTAGTTGGCCATCTGCAGGTCATACACGGTCGCCACCAGCGCAGTGCTGCCATCGGCCAGCGTCAGCCGGCGTGCCGGCACGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AM260480/2564484‑2564349 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgcCGGCTGACGCTGGCCGACGGCAGCACCGCGCTGGTGGCCACCGTGTACGACCTGCA  >  AY305378/288727‑288782

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