Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

ERROR: Cannot find seed sequence in genome of reference species for Dracaena_cambodiana_DN10000_c0_g1_i1.p1! #31

Closed
majssssa opened this issue Mar 28, 2023 · 4 comments

Comments

@majssssa
Copy link

I encountered the following error while running fDOG
ERROR: Cannot find seed sequence in genome of reference species for Dracaena_cambodiana_DN10000_c0_g1_i1.p1! This is my code fdog.run --seqFile Dracaena_cambodiana.pep --jobName test --refspec ARATH@3702@3`.Also, nothing appears in the window when I run this command.Finally, the above error occurs.
This is my seqFile

Dracaena_cambodiana_DN10000_c0_g1_i1.p1
MVAVFNKELLSCYLITLKLKQTVEAGLAKSQPNSPTPKPPLRLTQGQPPQPPVPRADSEWVISIREKLDEARQEQAACPWAKLSIYRVPKSLREGDSKAYVPQVVSLGPYHRGKHHLRGMDHHKWRALHKVLRRTGHDVELYLDSIKMLEERARACYEGSLSLTSNEFVETMVLDGIFILELFRGAAGEGFKQLGYSHNDPVFAMRATMHSIQRDMIMLENQIPLFVLDRLLALQICKPEQSGLVATLAVQFFDPLIPTDEPLRKIDRNKLQSTSPSLKATAAVFDP
Dracaena_cambodiana_DN10004_c0_g1_i1.p1
VKEEKLSKDSNSKNGAKSGGELKRASKRPDLKAVSKMIENHQVLYPEKMIGHLPGIDVGDQFLSRAEMCVLGLHNHWLNGIDYMGQSFAKLERYKNYTFPLAICIVLSGMYEDDLDNSEDIVYTGQGGHDLLGSKHQVSDQKMERGNLALKNNCDLGVPVRVVRGHELKSSYSGKVYTYDGLYRVIKTWAEKGVSGFTVYKYKLRRVEGQPVLRTNQVQFTRAT
Dracaena_cambodiana_DN10014_c0_g1_i1.p1
VRLLLLVGPAQVQIQQLLLEKLPEHFDAVSPSRNLNDDVARLIVNQFRWLDFLVDSRGFAEKLMEVLSIAPLLLKKEIIGSLPEIIGDQNSSMVIASLERLLQEDSEVIVPVLDSFSNLNLDDQLQEQAVTIALSCIRTVDAEHMPHLLRFLLLSATSVNVGRIISQIREQLKFVGVSDPCVAKQKMLKGKYLADSTEASILDALRSSLRFKNTTCEAIFKELKSLDHPRDHKVIDVWFLMLIYTNGGSLQKNAEKIL
Dracaena_cambodiana_DN10015_c0_g1_i1.p1
MTPNANTLPPPQTPTPNSLIMAGKSPKPSILTDGLLFFGGAAVALLLFLTLYSFLSPLAPSSPSLPFRHSLPLLANPHSPRHHHNSCDSNLHIDPSSPNLYDDPSVSYTLSYSGSLITDWDRKRTNWLHSHPTFAAKRNRILMVSGSQPGPCRNPIGDHLLLRFYKNKADYCRLHDIDLFYNTALLQPDMFSFWAKIPVVRAAMVAHPEAEWIWWVDSDAAITDMEFELPLDRYEGYNLVLHGWRDLVYKAKSWTSLNAGVFLMRNCQWSLDFIDTWSRMGPQSPEYERWGRIQRAIFKDKLYAASDDQTGIAYLLLEEKRKWADKIYLESAFYFEGYWVEIVNKLERIEAEYLEMDRRVKGLRRRRAEKVAVAYGRMREELLEKQGVRRGMEGWRRPFITHFTGCQPCSGDHNKEYSGENCFGGMQRALNFADNQVLRSYGFRHERLVNSTGVRPLPFDFPAAG*
Dracaena_cambodiana_DN10025_c0_g1_i1.p1
MSPCGGVGHEPAVLQLQKWGYLQFQLEPSEFRLASISPTRDLLLLLSYQCEALLLPLLLGKFQSENICEPNSSEQVITCIPDDVDSAQCSKNDEESVKGAALPILESSSASKSYPVISGVKSLAWGHS
Dracaena_cambodiana_DN10026_c0_g1_i1.p1
VICLGIFFRSWPSPRPSFIRDISQQNSSESYGLVHVCALGGEFLCRLAVFAILERLGSQILYSMFASFCLMAAIFLRRNIVETKGKTLQEIEISLRQPE*
Dracaena_cambodiana_DN10027_c0_g1_i1.p4
MSPTDRGGPTSSAVGSLHLLPSELLHSILLRLPLPDLLRLRSLSRLLLSSISAPEFRRTYQSTSPWLFLFQKRPPRSSLLRGFDPRSARWFTFPSLSALIAAPALPPGDDLYLLA
Dracaena_cambodiana_DN10028_c0_g1_i1.p1
MGFTAAELERLISSNARILIVAPAIPRLEFWQAFIGGDNRKDLVSVLTRNRGLITHDIVSGIAPKMLLLKEHGLSQRDIVGLVKRGHGFITRSSKTIEAVLNSAKELGLDTKSPMFGHTLSSLVSFSSDSFKAKMEIFRGFGWSEEELLAAFKKAPSFLHLSEENIREKMEFLVGRAGCKQSYIALNPLLLTFSLEKRLRPRHYVMEVLKSKGIMGRARFSKIMCLTEKKFVESVILPCKEQVPNIHELYIA
Dracaena_cambodiana_DN10030_c0_g1_i1.p1
LFRNPLSLSLLRKFFPKQISNLMGYLSCRAESSISTCRSISAISPSQSPTSTTNTKKKHLKPYKSLQNQENDDDDDDKSPCRRIEHFTYAELESATNNFSNSSLLGRGS

@trvinh
Copy link
Member

trvinh commented Mar 28, 2023

Please prefer to this wiki to make sure you have the correct input data for fDOG https://github.com/BIONF/fDOG/wiki/Input-and-Output-Files (and https://github.com/BIONF/fDOG/wiki/Use-standalone-version)

  • For fdog.run, it must be a SINGLE sequence in FASTA format
  • The reference species (specified with --refspec) must be present in the blast database directory

Also check this page https://github.com/BIONF/fDOG/wiki/How-fDOG-works to know the differences between core taxa and search taxa. I think you are mistaken between these, since you said that your reference species is Dracaena_cambodiana_DN10000_c0_g1_i1.p1 but you have set --refspec ARATH@3702@3.

@trvinh
Copy link
Member

trvinh commented Mar 30, 2023

@majssssa This issue is solved, right?

@majssssa
Copy link
Author

I'm sorry to see your reply just now. This problem has been successfully solved

@trvinh
Copy link
Member

trvinh commented Mar 30, 2023

Nice! Then I will close it :)

@trvinh trvinh closed this as completed Mar 30, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants