New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AM260480	2577612 =	53 (1.800)	6 (0.230)	5/192	1.8	13.4%	coding (827/1665 nt)	nirS	cytochrome cd1 nitrite reductase (NirS)
	?	AM260480	2577663 =	30 (1.130)					coding (878/1665 nt)	nirS	cytochrome cd1 nitrite reductase (NirS)
Rejected: Frequency below/above cutoff threshold.

GACAAGTCCGAGTAGTTGACCATCAGGATCTTGCCGGTCTCCTTGGCGTTGATCACGAACTCCGGATGGAAATGGCTGGCGACGATGGCCGCCACGCGCGGCTCGGGGTGGTATTCGTTGTCCACGGTCATGCCGCGCGTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AM260480/2577753‑2577612 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccgcgcgtggCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCAACCTGAAGACCACCACCATCGATTCGGCCAAGTTCCTG  >  AM260480/2577663‑2577793                                                                                                                                                                                                                                                                                   GACAAGTCCGAGTAGTTGACCATCAGGATCTTGCCGGTCTCCTTGGCGTTGATCACGAACTCCGGATGGAAATGGCTGGCGACGATGGCCGCCACGCGCGGCTCGGGGTGGTATTCGTTGTCCACGGTCATGCC                                                                                                                                             <  1:287594/134‑1   CAAGTCCGAGTAGTTGACCATCAGGATCTTGCCGGTCTCCTTGGCGTTGATCACGAACTCCGGATGGAAATGGCTGGCGACGATGGCCGCCACGCGCGGCTCGGGGTGGTATTCGTTGTCCACGGTCATGCCGC                                                                                                                                           <  2:282037/134‑1      GTCCGAGTAGTTGACCATCAGGATCTTGCCGGTCTCCTTGGCGTTGATCACGAACTCCGGATGGAAATGGCTGGCGACGATGGCCGCCACGCGCGGCTCGGGGTGGTATTCGTTGTCCACGGTCATGCCGCGCG                                                                                                                                        >  1:948884/1‑134         CGAGTAGTTGACCATCAGGATCTTGCCGGTCTCCTTGGCGTTGATCACGAACTCCGGATGGAAATGGCTGGCGACGATGGCCGCCACGCGCGGCTCGGGGTGGTATTCGTTGTCCACGGTCATGCCGCGCGTG                                                                                                                                      <  2:919800/133‑1                                               CGTTGATCACGAACTCCGGATGGAAATGGCTGGCGACGATGGCCGCCACGCGCGGCTCGGGGTGGTATTCGTTGTCCACGGTCATGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTG                                                                                               <  1:551996/134‑1                                                        CGAACTCCGGATGGAAATGGCTGGCGACGATGGCCGCCACGCGCGGCTCGGGGTGGTATTCGTTGTCCACGGTCATGCCGCGCGTG                                                                                                                                      >  1:493141/1‑86                                                        CGAACTCCGGATGGAAATGGCTGGCGACGATGGCCGCCACGCGCGGCTCGGGGTGGTATTCGTTGTCCACGGTCATGCCGCGCGTG                                                                                                                                      <  2:493141/86‑1                                                           atTCCGGATGGAAATGGCTGGTGACGATGGGCGACACGCGCGGCTCGGGGTGGTATTCGTTGTCCACGGTCATGCCGCGCGTGGCGGCTATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAG                                                                                   <  1:784315/132‑1                                                                                                       TCGGGGTGGTATTCGTTGTCCACGGTCATGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCAA                                       >  2:699486/1‑134                                                                                                       TCGGGGTGGTATTCGTTGTCCACGGTCATGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGG                                                             <  1:909964/112‑1                                                                                                       TCGGGGTGGTATTCGTTGTCCACGGTCATGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGG                                                             >  2:909964/1‑112                                                                                                                     GTTGTCCACGGTCATGCCGCGCGTGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCAACCTGAAGACCACCA                         <  1:699486/134‑1                                                                                                                                            TGGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCAACCTGAAGACCACCACCATCGATTCGGCCAAGTTCCTG  >  1:714410/1‑134                                                                                                                                             GGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCAAC                                      <  1:588545/97‑1                                                                                                                                             GGCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCAAC                                      >  2:588545/1‑97                                                                                                                                                                                                                                                                                   GACAAGTCCGAGTAGTTGACCATCAGGATCTTGCCGGTCTCCTTGGCGTTGATCACGAACTCCGGATGGAAATGGCTGGCGACGATGGCCGCCACGCGCGGCTCGGGGTGGTATTCGTTGTCCACGGTCATGCCGCGCGTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AM260480/2577753‑2577612 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccgcgcgtggCGGCCATCGTCGCCAGCCATTTCCATCCGGAGTTCGTGATCAACGCCAAGGAGACCGGCAAGATCCTGATGGTCAACTACTCGGACTTGTCCAACCTGAAGACCACCACCATCGATTCGGCCAAGTTCCTG  >  AM260480/2577663‑2577793

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