Hello, I am trying to convert the FASTA sequence in to two-bit format by faToTwoBit program downloaded from UCSC genome browser. I have converted it to executable bit but am unable to run the program on my FASTA sequence. What is the command line to run it?Please help.
If you run the program on its own without any options it will show you how to use it
$ ./faToTwoBit
faToTwoBit - Convert DNA from fasta to 2bit format
usage:
faToTwoBit in.fa [in2.fa in3.fa ...] out.2bit
So use one or more input fasta files and provide an output name out.2bit.
I have found this usage but did not understand it, please clarify.
I have put this command faToTwoBit Chr1.fa out.2bit but the terminal shows faToTwoBit: command not found....why is this happening?
I did as you mentioned. I put the command ./faToTwoBit Chr1.fa.2bit This appeared again:
faToTwoBit - Convert DNA from fasta to 2bit format usage: faToTwoBit in.fa [in2.fa in3.fa ...] out.2bit options: -long use 64-bit offsets for index. Allow for twoBit to contain more than 4Gb of sequence. NOT COMPATIBLE WITH OLDER CODE. -noMask Ignore lower-case masking in fa file. -stripVersion Strip off version number after '.' for GenBank accessions. -ignoreDups Convert first sequence only if there are duplicate sequence names. Use 'twoBitDup' to find duplicate sequences.
What is wrong in the command?
I have used the command you mentioned, actually nothing happened. It returned to its original $ position again
That means the command worked and made the Chr1.2bit file. You can check that by listing the files.
Sounds like you need to get familiar with basic command line skills. Please spend some time on that first. I recommend: https://angus.readthedocs.io/en/2019/shell_intro/index.html
Yes, I am completely new to Linux and still learning. Thank You! It worked. Thanks for the link too.
Log in to answer this question.
