New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 3140629	42 (0.910)	4 (0.090)	4/268	6.2	8.6%	coding (171/489 nt)	hybE	hydrogenase 2‑specific chaperone
	?	NC_000913_3_bme_pgi	= 3140635	45 (1.040)					coding (165/489 nt)	hybE	hydrogenase 2‑specific chaperone
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_bme_pgi/3140480‑3140629 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑acggggtgTGTGATCACCCCGTGGATGCTGAGTGCAGTTATCTTCCCCGGCCCGGATCAACTCTGGCCGCTGCGCAAAGTGAGTGAAAAAATTGGTCTGCAACTGCCGTATGGCACTATGACCTTTACCGTTGGCGAACTGGACGGTGTTTCGCA  <  NC_000913_3_bme_pgi/3140635‑3140489                                                                                                                                                                                                                                                                                                           GAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCACACCAATTTTTTCACTCACTTTGCGCAGCGGCCCGAGTTGATCCGGGCCGGGGAAGAACACTGCACTCAGCATCCACGGGGT                                                                                                                                                      >  1:196774/1‑149              AAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGTGTGATCACCCC                                                                                                                                         >  1:316466/1‑149                ACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTGTGTGATCACCCCGT                                                                                                                                       <  2:316466/149‑1                                          GGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGG                                                                                                                                                        >  1:38206/1‑106                                          GGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGG                                                                                                                                                        <  2:38206/106‑1                                                                                                                           GATAACTGCACTCAGCATCCACGGGGTGTGTGATCACCCCGTGGATGCT                                                                                                                                <  1:780435/49‑1                                                                                                                           GATAACTGCACTCAGCATCCACGGGGTGTGTGATCACCCCGTGGATGCT                                                                                                                                >  2:780435/1‑49                                                                                                                                               ACGGGGTGTGTGATCACCCCGTGGATGCTGAGTGCAGTTATCTTCCCCGGCCCGGATCAACTCTGGCCGCTGCGCAAAGTGAGTGAAAAAATTGGTCTGCAACTGCCGTATGGCACTATGACCTTTACCGTTGGCGAACTGGACGGTGT        <  2:536723/149‑1                                                                                                                                                     TGTGTGATCACCCCGTGGATGCTGAGTGCAGTTATCTTCGCCGGCCCGGATGAACTCTGGCCGCTGCGCAAAGTGAGTGACAACATTGGTCTACACCTGCCGTATGGCACTATGACCTTTACCGTTACCGTCCTAGACGGTATTTCGCA  >  2:533007/1‑149                                                                                                                                                                                                                                                                                                           GAGAGATATTGCGAAACACCGTCCAGTTCGCCAACGGTAAAGGTCATAGTGCCATACGGCAGTTGCAGACCAATTTTTTCACTCACTTTGCGCAGCGGCCAGAGTTGATCCGGGCCGGGGAAGATAACTGCACTCAGCATCCACGGGGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_bme_pgi/3140480‑3140629 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑acggggtgTGTGATCACCCCGTGGATGCTGAGTGCAGTTATCTTCCCCGGCCCGGATCAACTCTGGCCGCTGCGCAAAGTGAGTGAAAAAATTGGTCTGCAACTGCCGTATGGCACTATGACCTTTACCGTTGGCGAACTGGACGGTGTTTCGCA  <  NC_000913_3_bme_pgi/3140635‑3140489

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