New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 1536866	40 (1.340)	4 (0.140)	4/266	3.6	9.9%	coding (444/696 nt)	narW	nitrate reductase 2 (NRZ), delta subunit (assembly subunit)
	?	NC_000913_3_bme_pgi	= 1536879	36 (1.290)					coding (431/696 nt)	narW	nitrate reductase 2 (NRZ), delta subunit (assembly subunit)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_bme_pgi/1536709‑1536866 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaacaacgcgAGGCACCGTGGTACGCGTTGTTTGATGCTCTGCTGCAACTGGCGGGAAGCTCCCTTTCAAGTGACAGTGTCACAAAACAAGTAAACAGCGAAGAGCGCGATGACACCCGCCAGGCGCTTGATGCAGTGTGGGAAGAGGAACAGGTGAA  <  NC_000913_3_bme_pgi/1536879‑1536732                                                                                                                                                                                                                                                                                                                    CCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCA                                                                                                                                                               >  2:411990/1‑149              TTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCCCGCTCTTCGCTGTTTCCTTGTTTTGTGACACTGTCCCTTGCAAGGGCGCTTCCCGCCAGTTGCAGCATAGCATCAAACCACCCGtccc                                                                                                                                                  >  2:130539/1‑145                                         CCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGAGGCACCGTGGTACGCGTTGT                                                                                                                                 >  1:227198/1‑139                                         CCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGAGGCACCGTGGTACGCGTTGT                                                                                                                                 <  2:227198/139‑1                                                                       ATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGAGGCACCGTGGTACGCGTTGTT                                                                                                                                >  1:265310/1‑110                                                                       ATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGAGGCACCGTGGTACGCGTTGTT                                                                                                                                <  2:265310/110‑1                                                                                                                                                          ACGCGAGGCACCGTGGTACGCGTTGTTTGATGCTCTGCTGCAACTGGCGGGAAGCTCCCTTTCAAGTGACAGTGTCACGAAACAAGTAAACAGCGAAGAGCGCGATGACACTCGCCAGGCGCTTGATGCAGTGTGGGAAGAGGAACAGG      >  1:487770/1‑149                                                                                                                                                            GCGAGGCACCGTGGTACGCGTTGTTTGATGCTCTGCTGCAACTGGCGGGAAGCTCCCTTTCAAGTGACAGTGTCACAAAACAAGTAAACAGCGAAGAGCGCGATGACACCCGCCAGGCGCTTGATGCAGTGTGGGAAGAGGAACAGGTG    >  1:194431/1‑149                                                                                                                                                              GAGGCACCGTGGTACGCGTTGTTTGATGCTCTGCTGCAACTGGCGGGAAGCTCCCTTTCAAGTGACAGTGTCACAAAACAAGTAAACAGCGAAGAGCGCGATGACACCCGCCAGGCGATTGATGCAGTGTGGGAAGAGGAACAGGGGAA  >  2:1880/1‑149                                                                                                                                                                                                                                                                                                                    CCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_bme_pgi/1536709‑1536866 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaacaacgcgAGGCACCGTGGTACGCGTTGTTTGATGCTCTGCTGCAACTGGCGGGAAGCTCCCTTTCAAGTGACAGTGTCACAAAACAAGTAAACAGCGAAGAGCGCGATGACACCCGCCAGGCGCTTGATGCAGTGTGGGAAGAGGAACAGGTGAA  <  NC_000913_3_bme_pgi/1536879‑1536732

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