-- Leo's gemini proxy
-- Connecting to freeshell.de:1965...
-- Connected
-- Sending request
-- Meta line: 20 text/gemini;lang=en-GB
> A swiss-army knife of tools for genomic-analysis tasks.
> Used to intersect, group, convert and count data in BAM, BED, GFF/GTF, VCF format.
Intersect two files regarding the sequences' strand and save the result to the specified file:
bedtools intersect -a {path/to/file_1} -b {path/to/file_2} -s > {path/to/output_file}
Intersect two files with a left outer join, i.e. report each feature from {file_1} and NULL if no overlap with {file_2}:
bedtools intersect -a {path/to/file_1} -b {path/to/file_2} -lof > {path/to/output_file}
Using more efficient algorithm to intersect two pre-sorted files:
bedtools intersect -a {path/to/file_1} -b {path/to/file_2} -sorted > {path/to/output_file}
Group file {`path/to/file`} based on the first three and the fifth column and summarize the sixth column by summing it up:
bedtools groupby -i {path/to/file} -c 1-3,5 -g 6 -o sum
Convert bam-formatted file to a bed-formatted one:
bedtools bamtobed -i {path/to/file}.bam > {path/to/file}.bed
Find for all features in {file_1}.bed the closest one in {file_2}.bed and write their distance in an extra column (input files must be sorted):
bedtools closest -a {path/to/file_1}.bed -b {path/to/file_2}.bed -d
> Copyright © 2014—present the tldr-pages team and contributors.
> This work is licensed under the Creative Commons Attribution 4.0 International License (CC-BY).
-- Response ended
-- Page fetched on Fri May 17 15:57:37 2024