User: sjneph
sjneph • 610
- Reputation:
- 610
- Status:
- Trusted
- Location:
- Last seen:
- 5 years, 8 months ago
- Joined:
- 7 years ago
- Email:
- s*****@gmail.com
about me
Posts by sjneph
<prev
• 28 results •
page 1 of 3 •
next >
0
votes
3
answers
3.4k
views
3
answers
... updated to put MAPQ in field 5 of your-file.bed since that's what you want to average over.
...
written 5.9 years ago by
sjneph • 610
3
votes
3
answers
3.4k
views
3
answers
...
bam2bed < your-file.bam | cut -f1-4,7 > your-file.bed
make-regions 1000 200 regions.bed | bedmap --echo --mean - your-file.bed > answer.bed
Where make-regions is an awk statement, perhaps buried in a shell script
awk -v winsize=$1 -v stepsize=$2 \
'{ \
i = $2; \
while ( i < $3 ) ...
written 5.9 years ago by
sjneph • 610
4
votes
1
answer
3.0k
views
1
answers
... These are available through ftp at EBI.
As the README explains, all 41 cell type footprints are available in one file, and the FOS is shown in column 5 (lower score is a more significant score).
...
written 5.9 years ago by
sjneph • 610
2
votes
3
answers
3.2k
views
3
answers
... Alex is spot on that you will receive a tremendous number of files that would be difficult to manage. Instead, you could consider demarcating every row with an indicator of what file it would be in with your method.
$ awk \
'{ \
regionChromosome = $1; \
regionStart = $2; \
...
written 5.9 years ago by
sjneph • 610
1
vote
2
answers
4.5k
views
2
answers
... Here is a straightforward way to do it with awk:
A: Does bedops have a command similar to the bedtools makewindows?
...
written 5.9 years ago by
sjneph • 610
1
vote
2
answers
2.2k
views
2
answers
... This is pretty straightforward with BEDOPS:
sort-bed file.bed > file.sort.bed
closest-features --no-overlaps --dist --closest file.sort.bed file.sort.bed | awk -F'|' '$3*$3==1'
You basically ask what elements are closest 'to the left' and 'to the right' of every element in the file using an ed ...
written 6.0 years ago by
sjneph • 610
2
votes
2
answers
1.5k
views
2
answers
... sort-bed file1 > file1.sorted
sort-bed file2 > file2.sorted
bedops -e -1 file1.sorted file2.sorted > answer.bed
bedops -e -1 file2.sorted file1.sorted >> answer.bed
The answer.bed file will not be sorted per sort-bed since you catted on the results. If you want to maintain sort ...
written 6.0 years ago by
sjneph • 610
• updated
6.0 years ago by
Alex Reynolds ♦ 29k
1
vote
3
answers
2.9k
views
3
answers
... use
bedops --merge a.bed b.bed > merged.bed
the requirements are that both bed files are sorted, which is easy and smart to do, as it makes all downstream analyses more efficient. The output is already sorted, so you only need to do this once for initial, unsorted files.
sort-bed A.bed > ...
written 6.5 years ago by
sjneph • 610
2
votes
3
answers
5.8k
views
3
answers
... You might consider converting the bam file to starch format, which is a highly compressed version of BED. The starch file will be smaller than the original BAM file, and there is no loss of information. Once in starch format, you can use other BEDOPS utilities to annotate things the way you want. ...
written 6.7 years ago by
sjneph • 610
3
votes
73
answers
32k
views
73
answers
... Yea, until you make the assumption that -w actually works only on whitespace.
printf "foo-choo" | grep -Fw -e "foo"
returns foo-choo. Hate that.
...
written 6.8 years ago by
sjneph • 610
Latest awards to sjneph
Supporter
5.8 years ago,
voted at least 25 times.
Teacher
5.9 years ago,
created an answer with at least 3 up-votes.
For A: Where To Find Encode Dnase I Footprint Data
Appreciated
7.0 years ago,
created a post with more than 5 votes.
For A: Comparing 2 Large Lists (Millions Of Rows) To Identify Shared And Exclusive Elem
Teacher
7.0 years ago,
created an answer with at least 3 up-votes.
For A: Comparing 2 Large Lists (Millions Of Rows) To Identify Shared And Exclusive Elem
Good Answer
7.0 years ago,
created an answer that was upvoted at least 5 times.
For A: Comparing 2 Large Lists (Millions Of Rows) To Identify Shared And Exclusive Elem
Teacher
7.0 years ago,
created an answer with at least 3 up-votes.
For A: Group Reads In Bed File Based On The Distance
Commentator
7.0 years ago,
created a comment with at least 3 up-votes.
For C: Bedtools: Analyzing Genomic Features
Teacher
7.0 years ago,
created an answer with at least 3 up-votes.
For A: Does Bedops Have A Command Similar To The Bedtools Makewindows?
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 2.3.0
Traffic: 825 users visited in the last hour