How to print the header of bed file and the lines that belong to chromosome 2
I have a bed file consists of many lines and I want to extract the header and the lines that belong to chromosome 2 only in another file
• 2,410 views
•
link
updated
by
Ram
•
written
by
FadyNabil •
•
written
by
FadyNabil •
2 answers
awk 'NR == 1 || ($1=="chr2")' lamina.bed
• 0 views
•
link
updated
by
Ram
•
written
by
FadyNabil •
•
written
by
FadyNabil •
If the BED file is sorted, here's a faster way:
$ cat <(head -1 regions.bed) <(bedextract chr2 regions.bed) > answer.bed
Generically, using awk to solve this will require reading through the entire file. This can be expensive for whole-genome scale files.
• 0 views
•
link
written
by
Alex Reynolds
Log in to answer this question.
More posts like this
-
Calculating polygenic score using plink
written by m.shoaib •
-
samtools view extract from .bed list to .bam file
written by matthews.alix •
-
Parsing header of FASTA File
written by Shelle •
-
bedtools getfasta chromosome not found
written by from the mountains
-
identify the type of repeat genomic region of a BED file
written by jing.mengrabbit •
-
How to extract specific lines from a MCL dump file
written by nancydong20
-
Get cds using bedtools
written by User000
-
Changing the bam header
written by hpapoli
-
Aligning fasta sequences against a reference genome with scaffolds
written by thjnant
-
Intersectbed Overlap
written by Nicolas Rosewick
