New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4289024	125 (1.220)	3 (0.030)	3/288	12.3	2.3%	coding (1261/1437 nt)	nrfA	nitrite reductase, formate‑dependent, cytochrome
	?	NC_000913	= 4289082	139 (1.440)					coding (1319/1437 nt)	nrfA	nitrite reductase, formate‑dependent, cytochrome
Rejected: Position hash score below cutoff.

CGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4288824‑4289024 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcctgctggGCTTTCTCTTTGGTTGAGATATCCGGGATCTGGATTTCATGGGTGATGCCTTTGGTCGCCAGCAGGCGCGCCAGTTTGGTGCGTGCATCCGCCGCTTTATCCATCGCCGTACCGAGCATCCGTAAACCTTCTTCCGGTGCGTGCATATGAATGCCGTGGGAAGCGATCGCCAGATCC  <  NC_000913/4289082‑4288906                                                                                                                                                                                                                                                                                                                                                                                            CGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTG                                                                                                                                                                                    >  4:324000/1‑200                                                                     GCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAACCTGGCGCGCCTGCTG                                                                                                                                                                                    >  2:548112/1‑132                                                                          GTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGG                                                                                                                                                                                   <  3:671261/128‑1                                                                          GTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGG                                                                                                                                                                                   >  4:671305/1‑128                                                                                                                                                  CGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGGCTTTCTCTTTGGTTGAGATATCCGGGATCTGGATTTCATGGGTGATGCCTTTGGTCGCCAGCAGGCGCGCCAGTTTGGTGCGTGCATCCGCCG                                                                                     >  3:849522/1‑150                                                                                                                                                              GGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGGCTTTCTCTTTGGTTGAGATGTCCGGGATCTGGATTTCATGGGTGATGCCTTTGGTCGCCAGCAGGCGCGCCAGTTTGGTGCGTGCATCCGCCGCTTTATCCATCG                                                                         >  3:133327/1‑150                                                                                                                                                                                   CACCAAACTGGCGCGCCTGCTGGGCTTTCTCTTTGGTTGAGATATCCGGGATCTGGATTTCATGGGTGATGCCTTTGGTCGCCAGCAGGCGCGCCAGTTTGGTGCGTGCATCCGCCGCTTTATCCATCGCCGTACCGAGCATCCGTAAACCTTCTTCCGGTGCGTGCATATGAATGCCGTGGGAAGCGATCGCCAGATCC  >  4:809556/1‑200                                                                                                                                                                                                   CTGCTGGGCTTTCTCTTTGGTTGAGATATCCGGGATCTGGATTTCATGGGTGATGCCTTTGGTCGCCAGCAGGCGCGCCAGTTTGGTGCGTGCATCCGCCGCTTTATCCATCGTCGTACCGAGCATCCGTAA                                                      >  2:166512/1‑132                                                                                                                                                                                                   CTGCTGGGCTTTCTCTTTGGTTGAGATATCCGGGATCTGGATTTCATGGGTGATGCCTTTGGTCGCCAGCAGGCGCGCCAGTTTGGTGCGTGCATCCGCCGCTTTATCCATCGCCGTACCGAGCATCCGTAAACCTTCTTCCGGTGCGTG                                    >  3:498751/1‑150                                                                                                                                                                                                   CTGCTGGGCTTTCTCTTTGGTTGAGATATCCGGGATCTGGATTTCATGGGTGATGCCTTTGGTCGCCAGCAGGCGCGCCAGTTTGGTGCGTGCATCCGCCGCTTTATCCATCGCCGTACCGAGCATCCGTAA                                                      >  2:126757/1‑132                                                                                                                                                                                                      cgGGGCTTTCTCTTTGGTTGAGATATCCGGGATCTGGATTTCATGGGTGATGCCTTTGGTCGCCAGCAGGCGCGCCAGTTTGGTGCGTGCATCCGCCGCTTTATCCATCGCCGTACCGAGCATCCGTAAACCTTCTTCCGGTGCGTGCAT                                 <  3:117347/148‑1                                                                                                                                                                                                         GGCTTTCTCTTTGGTTGAGATATCCGGGATCTGGATTTCATGGGTGATGCCTTTGGTCGCCAGCAGGCGCGCCAGTTTGGTGCGTGCATCCGCCGCTTTATCCATCGCCGTACCGAGCATCCGTAAACCTTCTTCCGGTGCGTGCATATG                              >  3:134030/1‑150                                                                                                                                                                                                         GGCTTTCTCTTTGGTTGAGATATCCGGGATCTGGATTTCATGGGTGATGCCTTTGGTCGCCAGCAGGCGCGCCAGTTTGGTGCGTGCATCCGCCGCTTTATCCATCGCCGTACCGAGCATCCGTAAACCTTCTTCCGGTGCGTGCATATG                              >  3:235537/1‑150                                                                                                                                                                                                         GGCTTTCTCTTTGGTTGAGATATCCGGGATCTGGATTTCATGGGTGATGCCTTTGGTCGCCAGCAGGCGCGCCAGTTTGGTGCGTGCATCCGCCGCTTTATCCATCGCCGTACCGAGCATCCGTAAACCTTCTTCCGGTGCGTGCATATG                              >  3:798259/1‑150                                                                                                                                                                                                         GGCTTTCTCTTTGGTTGAGATATCCGGGATCTGGATTTCATGGGTGATGCCTGTGGTCGCCAGCAGGCGCGCCAGTTTGGTGCGTGCATCCGCCGCTTTATCCATCGCCGTACCGAGCATCCGTAAACCTTCTTCCGGTGCGTGCATATG                              <  3:565744/150‑1                                                                                                                                                                                                                                                                                                                                                                                            CGAAAGCAGCGCTGGATGCAGGCGCGACGGAAGCTGAAATGAAGCCAATTCAGGACGATATCCGTCATGCCCAGTGGCGCTGGGATCTGGCGATCGCTTCCCACGGCATTCATATGCACGCACCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4288824‑4289024 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcctgctggGCTTTCTCTTTGGTTGAGATATCCGGGATCTGGATTTCATGGGTGATGCCTTTGGTCGCCAGCAGGCGCGCCAGTTTGGTGCGTGCATCCGCCGCTTTATCCATCGCCGTACCGAGCATCCGTAAACCTTCTTCCGGTGCGTGCATATGAATGCCGTGGGAAGCGATCGCCAGATCC  <  NC_000913/4289082‑4288906

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