bedSeq === USAGE ==== Usage: ./bedSeq <seqDir> <bedfilename> <bedfiletype:bed|ebed> [options] > <ofilename> Description:	append sequence to bed files (bed format for the feature coordinate, ebed format for feature block coordinates)	seqDir contains per chromosome <chr>.seq which is one line raw sequence without header. convert fa to seq file using convertFaToPureSeq.py Options: --output-fasta	output fasta file instead of appending seq to bed file --use-coord-as-name	use coordinate as name of sequence instead of the name field of bed file. Only valid when --output-fasta option is specified --use-block-coord-as-name	use block coordinate as name. Only when --output-fasta is specified and format is ebed === INSTALLATION === compile using bash make.sh install into path by changing, e.g., ~/.bashrc export PATH=${PATH}:/path/to/bedSeq/folder === TIPS for making seq files from fasta files (one fasta per chromosome) === Say you have <chr>.fa for each chromosome, e.g., chr1.fa You can use the convertFaToPureSeq.py python script to do conversion. If it is in your path, then for i in *.fa; do convertFaToPureSeq.py $i ${i/.fa/}.seq; done this will make <chr>.seq for each chromosome. e.g., chr1.fa => chr1.seq === EXAMPLES === example file is top 10 lines from mm9 acembly (http://www.ncbi.nlm.nih.gov/IEB/Research/Acembly/) bed file downlaoded from UCSC genome browser table -> mouse genome mm9 -> gene prediction track -> AceView -> acembly -> format bed and renamed as acembly.ebed, then head acembly.ebed > acembly.top10.ebed --- to append sequences to the end of the bed file (assuming that you installed the seq files in the seq/ subfolder: bedSeq seq/ acembly.top10.ebed ebed > acembly.top10.ebedseq to get fasta files of the sequences: bedSeq seq/ acembly.top10.ebed ebed --output-fasta > acembly.top10.fasta