This is a test version of Biostars. For the public version, visit https://www.biostars.org.
add strand info to a bed file

Is there a script/tool that can add strand (+/-) information to the fourth column of a bed file? I have a bed file (first few lines below) of 45,000 targets each chromosome, start, and stop. Thank you :).

Target Bed

chr1    955542    955763
chr1    957570    957852
chr1    976034    976270
picard
• 4,049 views • link
updated by Ram Ram • written by bioguy24 bioguy24

1 answer 

sed 's/$/\t+/' < input.bed > output.bed
• 0 views • link
updated by Ram Ram • written by Pierre Lindenbaum Pierre Lindenbaum

Works great... thank you :)

As I am learning these techniques do you mind briefly explaining how this works? Thank you very much.

• 0 views • link
updated by Ram Ram • written by bioguy24 bioguy24

Log in to answer this question.

More posts like this
  • Get gVCF file from genome coordinates
    written by Paul Paul
    Dear all, I have bed file like: Chr \t start \t stop And I woud like to add to fourth column information about reference allele …
  • list variants that map to a genomic location
    written by bioguy24 bioguy24
    The result of my awk script combines and formats all regions of a bed file containing 11 million bases, into targets less than a number …
  • How to sort chromosome name in bed file
    written by raya.girish • raya.girish
    Hello Friends! I have one stupid question. I have a bed file chr1 955541 955763 AGRN.chr1.955552.955753 chr1 957569 957852 AGRN.chr1.957580.95784 in my bed file my …
  • add gene name from bed file to vcf in vcftools
    written by bioguy24 bioguy24
    I am trying to add gene information to a vcf file using vcftools: vcftools annotate \ -a genes.bed.gz \ -c CHROM,FROM,TO,GENE \ -h <(echo '##INFO=<ID=GENE,Number=1,Type=String,Description="Gene …
  • query using mysql on ucsc table browser
    written by bioguy24 bioguy24
    I have a folder of 57 text files each with 1000 locations in them. I was going to use the mySQL browser at mysql --user=genome …
  • picard tools error in strand field
    written by bioguy24 bioguy24
    I have a bed file formatted for picard-tools: <p>BI=</p> <pre> @HD VN:1.4 GO:none SO:coordinate @SQ SN:chr1 LN:249250621 @SQ SN:chr2 LN:243199373 @SQ SN:chr3 LN:198022430 @SQ SN:chr4 …
  • coverage bed error but produces output file
    written by bioguy24 bioguy24
    The below command produces a file (674 MB). ``` coverageBed \ -d \ -sorted \ -a /home/cmccabe/Desktop/NGS/bed/sorted_unix_5column_xgen_targets.bed \ -b /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_newheader.bam \ > /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_output.txt ``` Output …
  • bedtools 2.25 vs bedtools2.20 coverageBed output not consistent
    written by bioguy24 bioguy24
    Hopefully, this edit is more clear in explaining the issue and will help led to a solution. As I am lost as to why there …
  • bedmap output on one line
    written by bioguy24 bioguy24
    I am using: bedmap --echo --skip-unmapped --delim '\t' --echo-map-id-uniq sort_3column_xgen_targets.bed sort_gene.bed > answer.bed `answer.bed` ``` chr1 957570 957852 NOC2L chr1 976034 976270 PERM1 chr1 976542 …
  • interval file for picard
    written by bioguy24 bioguy24
    Appears not to like the interval files, but not sure why and no matter what I try I seem to get another error? Thanks. Code: …