Reverse complement of a bacterial gene

[seema12345]

Hello everyone,

I am trying to create reverse complement of a bacteria gene in minus strand. Though I have read a lot on its mechanism and even got a perl script done that generates reverse complement of a sequence in Fasta format.
But, I need to generate it for bases in rows, with separate columns of genome position, reference, mutation, reference supporting reads and variant supporting read. Will appreciate suggestion regarding suitable language for it.

Regards,
Seema

[TiborNagy]

I think you would like to create a VCF-like file (from a VCF file?). If this is the case Awk/Perl/Python is the suitable language.

[bi_maniac]

If you post an example, I could help you. I am new in Bioinformatics but seasoned in computers

[seema12345]

If you post an example, I could help you. I am new in Bioinformatics but seasoned in computers

Hi,

Attaching a sample data...column 1, 2, 3, 4, 5 are genome position of bases, reference, mutation, reference supporting read and variant supporting read, respectively.
I need to get reverse complement of both reference and mutation column.
The output will be exactly like this format...just the bases in reference and mutation column will be reverse complement ..

data:
2153889 T 12 0
2153890 C +G 11 1
2153891 A -G 11 1
2153892 G +T 10 1
2153893 C +A 11 1

I am trying to write the code in perl-shell script..

[bi_maniac]

I see here two issues: computer programming and bioinformatics ...

... I mean: I have no problems with perl script.
I do not know the meaning of those "+" and "-" signs.
Please correct me if I am wrong.