Bedtools Slop for Specific Features
I want to use Bedtools Slop on a GTF file, but only on specific features (for example only for those labeled as gene, or exon, ect.)
I have been trying to pipeline the command with awk, but to no avail:
awk '{ if ($3 == "gene") | bedtools slop -i A.bed -g my.genome -l 2 -r 3
Any help would be great!
1 answer
Please give it like below
awk '$3=="gene"{print}' my.gtf | bedtools slop -i - -g my.genome -l 2 -r 3
• 0 views
•
link
written
by
Jeffin Rockey
Log in to answer this question.
More posts like this
-
How do I consider that my reads might map to the 3' UTR (3'-tag RNA seq)?
written by applepeanut •
-
Bedtools fisher error
written by hemr3 •
-
How to get the best overlap
written by g744695539 •
-
Can bedtools subtract support multiple files?
written by msimmer92
-
Mapping Gene Models Back to a Reference Genome of Closely Related Species
written by gtasource
-
Augustus Gene Frequency and Gene Names (Using BRAKER Pipeline)
written by isu2017 •
-
Bedtools Flank Unexpected Behavior
written by gtasource
-
Bedtools Coverage Read Counts
written by gtasource
-
Explanation of ENSEMBL GTF features
written by Benn
-
tophat GTF 2.2 format question
written by Ann
