Computing.Net > Forums > Unix > perl script

Computer Problems? Computing.Net has over 1,000,000 posts about all things technology related! Over 90% answered within 24 hours! Click here to start participating now! Also, be sure to check out the New User Guide.

perl script

Reply to Message Icon

Name: rayken1
Date: June 4, 2009 at 23:42:15 Pacific
OS: Macintosh
Product: R / (darwin - unix)
Subcategory: General
Comment:

hi,
I needed to clarify that I wanted a perl script to read annotations of two .embl files and identify the tags that are the same in the two files hence that copy that line from one file to another,

Thanks.



Sponsored Link
Ads by Google

Response Number 1
Name: David Perry
Date: June 5, 2009 at 06:50:44 Pacific
Reply:

Can you post a sample of your data and what portions you would like to extract?


0

Response Number 2
Name: rayken1
Date: June 7, 2009 at 22:48:46 Pacific
Reply:

m.embl

FT CDS complement(7216..17805)
FT /locus_tag="TP01_0003"
FT /codon_start=1
FT /protein_id="XP_765530.1"
FT /db_xref="GI:71031777"
FT /db_xref="GeneID:3502673"
FT gene complement(<7216..>17805)
FT /locus_tag="TP01_0003"
FT /db_xref="GeneID:3502673"
FT mRNA complement(<7216..>17805)
FT /locus_tag="TP01_0003"
FT /product="hypothetical telomeric SfiI fragment 20 protein
FT 3"
FT /transcript_id="XM_760437.1"
FT /db_xref="GI:71031776"
FT /db_xref="GeneID:3502673"
FT CDS complement(join(18028..18116,19351..20668))
FT /locus_tag="TP01_0004"
FT /note="go_function: nutrient reservoir activity [goid
FT 0045735]"
FT /codon_start=1
FT /protein_id="XP_765531.1"
FT /db_xref="GI:71031779"
FT /db_xref="GeneID:3503550"

n.embl

with
FT T.parva_chr1.gb"
FT misc_feature 1..1941
FT /note="Automatically generated region of difference with
FT T.parva_chr1.gb"
FT misc_feature 1..1941
FT /note="Automatically generated region of difference with
FT T.parva_chr1.gb"
FT misc_feature 3544..16471
FT /note="Automatically generated region of difference with
FT T.parva_chr1.gb"
FT misc_feature 3544..16471
FT /note="Automatically generated region of difference with
FT T.parva_chr1.gb"
FT misc_feature 3544..16471
FT /note="Automatically generated region of difference with
FT T.parva_chr1.gb"
FT CDS complement(join(17937..18025,19259..20576))
FT /note="none"
FT /locus_tag="TP01_0004"
FT CDS complement(20859..21875)
FT /note="none"
FT /locus_tag="TP01_0005"


Hi David, the above are parts of the embl files am working on . If you check on the indented areas, You will see that they are both CDS features with the same locus tag but one has an empty note.

I need to read the same files as arrays and locate the same locus tags then copy the notes from the one with a note to the one that has no notes.

THANKS ALOT.


0

Sponsored Link
Ads by Google
Reply to Message Icon

Related Posts

See More






Use following form to reply to current message:

Login or Register to Reply
LoginRegister


Sponsored links

Ads by Google


Results for: perl script

Perl Scripts..search & display www.computing.net/answers/unix/perl-scriptssearch-amp-display/6939.html

Perl Script www.computing.net/answers/unix/perl-script/6893.html

perl script - record altering www.computing.net/answers/unix/perl-script-record-altering/3546.html