problem with bedtools getfasta
0
0
Entering edit mode
12 months ago
anasjamshed ▴ 120

I am trying this code from bedtools getfasta manual:

 $ cat test.fa
    >chr1
    AAAAAAAACCCCCCCCCCCCCGCTACTGGGGGGGGGGGGGGGGGG

    $ cat test.bed
    chr1 5 10

    $ bedtools getfasta -fi test.fa -bed test.bed
    >chr1:5-10
    AAACC

    # optionally write to an output file
    $ $ cat test.fa
>chr1
AAAAAAAACCCCCCCCCCCCCGCTACTGGGGGGGGGGGGGGGGGG

$ cat test.bed
chr1 5 10

$ bedtools getfasta -fi test.fa -bed test.bed
>chr1:5-10
AAACC

# optionally write to an output file
$ bedtools getfasta -fi test.fa -bed test.bed -fo test.fa.out

But I got this error:

It looks as though you have less than 3 columns at line 1 in file test.bed.  Are you sure your files are tab-delimited?

How can I solve this?

getfasta bedtools • 853 views
ADD COMMENT
1
Entering edit mode

What does cat -vet test.bed show?

ADD REPLY
0
Entering edit mode

Its showing like this: pybed

ADD REPLY
1
Entering edit mode

Looks like you have an extra tab at the end of the line. Perhaps that is why you get that error.

ADD REPLY
0
Entering edit mode

how can I Correct it?

ADD REPLY
0
Entering edit mode

Well, is the bed file tab delimited?

ADD REPLY
0
Entering edit mode

I dont know what is tab delimited? How can I make that file?

ADD REPLY

Login before adding your answer.

Traffic: 1957 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6