Account Options

  1. Sign in
The old Google Groups will be going away soon, but your browser is incompatible with the new version.
Google Groups Home
« Groups Home
feature request checkbed
There are currently too many topics in this group that display first. To make this topic appear first, remove this option from another topic.
There was an error processing your request. Please try again.
flag
  2 messages - Collapse all  -  Translate all to Translated (View all originals)
The group you are posting to is a Usenet group. Messages posted to this group will make your email address visible to anyone on the Internet.
Your reply message has not been sent.
Your post was successful
 
From:
To:
Cc:
Followup To:
Add Cc | Add Followup-to | Edit Subject
Subject:
Validation:
For verification purposes please type the characters you see in the picture below or the numbers you hear by clicking the accessibility icon. Listen and type the numbers you hear
 
Stéphane Plaisance  
View profile  
 More options Nov 16 2012, 4:51 am
From: Stéphane Plaisance <Stephane.Plaisa...@vib.be>
Date: Fri, 16 Nov 2012 10:51:43 +0100
Local: Fri, Nov 16 2012 4:51 am
Subject: feature request checkbed
Hi Aaron,

Would it be possible to add such a command to bedtools to diagnose all possible erroneous records and put them in a file or on screen.
A curative mode would even be nicer, cleaning the file at the same time as storing the error triggering lines.
I from time to time get errors while running bedtools commands that relate to wrong delimiters (not anymore ;-) ) or 'end less than start'.
I made a primitive and quite inefficient bash script to isolate the later after generating start=0 lines by liftover.

Thanks
Stephane

###
checkbed ()
{
# does not check for other IFS than 'tab'
if [[ $1 == *.bed.gz ]]; then
   name=$(basename $1 ".bed.gz");
   zcat $1 | gawk –v name=${name} 'BEGIN{FS="\t"; OFS="\t"}
        {if ($2==0 || $3<$2) print $0>name"-errors.bed"; else print $0; }'
else
    name=$(basename name ".bed");
    cat $1 | gawk –v name=${name} 'BEGIN{FS="\t"; OFS="\t"}
        {if ($2==0 || $3<$2) print $0>name"-errors.bed"; else print $0; }'
fi


 
You must Sign in before you can post messages.
To post a message you must first join this group.
Please update your nickname on the subscription settings page before posting.
You do not have the permission required to post.
Aaron Quinlan  
View profile  
 More options Nov 19 2012, 8:27 am
From: Aaron Quinlan <aaronquin...@gmail.com>
Date: Mon, 19 Nov 2012 08:27:06 -0500
Local: Mon, Nov 19 2012 8:27 am
Subject: Re: [bedtools-discuss] feature request checkbed

Hi Stephane,

I recognize the general utility of having validation tool.  However, I feel this is somewhat low priority, as the tools themselves will identify the exact cause of many common problems.  As such, we may consider something like this for a future release, but I would have to consider it somewhat low priority given our current development goals.

Thanks for the suggestion,
- Aaron
quinlanlab.org

On Nov 16, 2012, at 4:51 AM, Stéphane Plaisance <stephane.plaisa...@vib.be> wrote:


 
You must Sign in before you can post messages.
To post a message you must first join this group.
Please update your nickname on the subscription settings page before posting.
You do not have the permission required to post.
End of messages
« Back to Discussions « Newer topic     Older topic »