How do I run this perl script proprerly to get GC percent from multiple fasta files?
1
0
Entering edit mode
3.3 years ago
natt1621 • 0

Hi. I just found this script on Researchgate (https://www.researchgate.net/post/How-can-I-estimate-the-GC-content-of-a-genome) to get the GC content from fasta files (all the rights of this script go to Jennifer Meneghin, July 23, 2009), but I don't know clearly how to run it. I use Windows 10 x64 and I downloaded the Pearl interpreter (strawberry-perl-5.32.1.1-64bit). I run it through cmd on windows, executing the .pl file but I don't know how to input the secuences. I'm kind of new in this bioinformatics and scripts world, so any help would be appreciated. This the header of the script.

http://ibb.co/kcffTBx

genome sequencing alignment perl • 976 views
ADD COMMENT
0
Entering edit mode

If your using windows, it's better to use the windows subsystem for linux (WSL) in general so that you can run linux (almost natively) on windows.

I would also recommend seqtk comp or seqkit fx2tab as alternatives.

ADD REPLY
0
Entering edit mode

Thank you so much, I will take a look to it.

ADD REPLY
0
Entering edit mode
3.3 years ago
Joe 21k

The script tells you its usage:

perl scriptname.pl <fasta file>

You simply need to write the filepath/name of the fasta file in afterwards. You need a space after the scriptname, but do not include the < or > as these are just used to denote 'placeholder' values.

ADD COMMENT
0
Entering edit mode

Thank you so much. I can write the name of the fasta file or its path so the script can recognize the file, right?

ADD REPLY
0
Entering edit mode

Yes, but if your working directory isn't the directory where the file is, you will need to provide the absolute or relative filepath, rather than just the filename.

ADD REPLY

Login before adding your answer.

Traffic: 1635 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