New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	9898 =	29 (0.700)	7 (0.180)	7/252	3.9	20.6%	intergenic (+5/+30)	mog/satP	molybdochelatase incorporating molybdenum into molybdopterin/succinate‑acetate transporter
	?	NC_000913	9918 =	27 (0.700)					intergenic (+25/+10)	mog/satP	molybdochelatase incorporating molybdenum into molybdopterin/succinate‑acetate transporter
Rejected: Coverage evenness skew score above cutoff.

GCAATCATCCACTTTGCCGGCTGGATTGGGCTGATCTGCGGTGCCAGCGCAATCTATCTGGCGATGGGTGAAGTACTGAACGAGCAGTTTGGTCGCACCGTTCTGCCGATTGGTGAATCCCACTAATTGTTTTGAGATCCCCCCGCTCGGGGGGAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/10053‑9898 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggggggatCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGC  >  NC_000913/9918‑10056                                                                                                                                                                                                                                                                                                         GCAATCATCCACTTTGCCGGCTGGATTGGGCTGATCTGCGGTGCCAGCGCAATCTATCTGGCGATGGGTGAAGTACTGAACGAGCAGTTTGGTCGCACCGTTCTGCCGATTGGTGAATCCCACTAATTGTTTTGAGATCCCCCCGCTCG                                                                                                                                                    >  2:31998/1‑149  CAATCATCCACTTTGCCGGCTGGATTGGGCTGATCTGCGGTGCCAGCGCAATCTATCTGGCGATGGGTGAAGTACTGAACGAGCAGTTTGGTCGCACCGTTCTGCCGATTGGTGAATCCCACTAATTGTTTTGAGATCCCCCCGCTCGG                                                                                                                                                   <  3:24653/149‑1                                                       TATCTGGCGATGGGTGAAGTACTGAACGAGCAGTTTGGTCGCACCGTTCTGCCGATTGGTGAATCCCACTAATTGTTTTGAGATCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGG                                                                                                              >  5:123197/1‑133                                                       TATCTGGCGATGGGTGAAGTACTGAACGAGCAGTTTGGTCGCACCGTTCTGCCGATTGGTGAATCCCACTAATTGTTTTGAGATCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGG                                                                                                              <  6:123197/133‑1                                                                                                         GCCGATTGGTGAATCCCACTAATTGTTTTGAGATCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACC                                            >  6:185131/1‑149                                                                                                                                   TTTGAGATCCCCCCGCTCGGGGGGATCTCACAACAATTAGTGGGATTCACCAATCGGCAGCCCGGTGCGACCAAACGGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCCGCCCACTCCAGCCGGC                    >  3:154229/1‑147                                                                                                                                        GATCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGG             >  3:168743/1‑149                                                                                                                                         ATCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGC                                  >  1:172650/1‑127                                                                                                                                         ATCCCCCCGCTCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGC                                  <  2:172650/127‑1                                                                                                                                                     GGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGC  <  5:185131/147‑1                                                                                                                                                                                                                                                                                                         GCAATCATCCACTTTGCCGGCTGGATTGGGCTGATCTGCGGTGCCAGCGCAATCTATCTGGCGATGGGTGAAGTACTGAACGAGCAGTTTGGTCGCACCGTTCTGCCGATTGGTGAATCCCACTAATTGTTTTGAGATCCCCCCGCTCGGGGGGAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/10053‑9898 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggggggatCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGC  >  NC_000913/9918‑10056

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