You are right. I use fasta files that | downloaded from "https://www.uniprot.org/uniprot/" and use the address of downloaded files in NcbiblastpCommandline() directly.
For example A0A0E1R3H3.fasta file (downloaded form https://www.uniprot.org/uniprot/A0A0E1R3H3.fasta):
>tr|A0A0E1R3H3|A0A0E1R3H3_LISMN Penicillin-binding protein 3
OS=Listeria monocytogenes serotype 4b str. LL195 OX=1230340 GN=pbpC PE=3 SV=1
MKGCIMASYGGKKRNNKKAIIIGSIAAVAVLAGIAIYFFIQNQHKDEKNALAAAETFTSN
IAKEKYDKLGSSVSSESLKKVEFTKKEMEDKYQAVYDGIGAKDIKVKNLKSVYDDKENKF
NLTYELEMRTSLGKLATQKYKTTISKQDDDWKIDWKPALIFPGMVKTDKVRITEDEAERG
QIVDRNGNPLATTGQFAEAGVVPSKLGEGDEKVKNIADISKKLEVSTDYINKQLDQKWVQ
ADSFVPLVTLNKDNLPEATGLTYAQKEMRTYPLNEATSHLIGYVGEVSAEDIEKNPKLSV
GDVIGKSGLERYYDKQLRGKNGGAIKIINDQTKQEDTLQKIDKKDGEEIKLTIDAAVQKK
AFDSLGSETGAVTMINPTNGELLALVSTPSYDANQMVLGITSEDYAKYNDDKRLPFLARY
ANRYAPGSTFKTITATIGLDTGVTKPDKVREISGLQWQKDASWGKYFVTRVHDVPKVNMT
DALVHSDNIYFAQEGLEIGKDKLTAGLKKFDFDKEYNLPFTMKPAQISNDGLNSEILLAD
TSYGQGELLMSPIQQAIAYSAIASDGKMPYPKLDTKEKAGKETQATEAASANQVKAALVK
TVSDPAGTAHALQIQGHSIAAKTGTAELKEKQGEDGLENGFVYAFDADNPNYLMVGMIEN
VKGRGGSGLVIDKLKPVIESMYK
A0A0H2UQB5.fasta file (downloaded from https://www.uniprot.org/uniprot/A0A0H2UQB5.fasta):
>tr|A0A0H2UQB5|A0A0H2UQB5_STRPN LicC protein
OS=Streptococcus pneumoniae serotype 4 (strain ATCC BAA-334 / TIGR4)
OX=170187 GN=licC PE=1 SV=1
MKAIILAAGLGTRLRPMTENTPKALVQVNQKPLIEYQIEFLKEKGINDIIIIVGYLKEQF
DYLKEKYGVRLVFNDKYADYNNFYSLYLVKEELANSYVIDADNYLFKNMFRNDLTRSTYF
SVYREDCTNEWFLVYGDDYKVQDIIVDSKAGRILSGVSFWDAPTAEKIVSFIDKAYVSGE
FVDLYWDNMVKDNIKELDVYVEELEGNSIYEIDSVQDYRKLEEILKNEN
And code is:
output = NcbiblastpCommandline(query=A0A0E1R3H3_FileAdress, subject=A0A0H2UQB5_fileAdress, outfmt=0)()[0]
and I tried this :
output =NcbiblastpCommandline(query=A0A0E1R3H3_FileAdress, subject=A0A0H2UQB5_fileAdress, ungapped=False)()
However the results with code is:
Query= tr|A0A023W3H0|A0A023W3H0_KLEPN Class D OXA-48 carbapenemase (Fragment)
OS=Klebsiella pneumoniae OX=573 GN=OXA-48 PE=4 SV=1
Length=145
Score (Bits) E Value
Sequences producing significant alignments: 13.5 6.0
tr|O08355|O08355_PSEFL Mannitol dehydrogenase OS=Pseudomonas fluo...
> tr|O08355|O08355_PSEFL Mannitol dehydrogenase OS=Pseudomonas fluorescens
OX=294 GN=mtlD PE=1 SV=1 Length=493
Score = 13.5 bits (23), Expect = 6.0, Method: Composition-based stats.
Identities = 11/31 (35%), Positives = 12/31 (39%), Gaps = 3/31 (10%)
Query 105 RIVKQAMLTEANGDYIIRAKTGYSTRIEPKI 135
RIV LT G Y I G P+I
Sbjct 125 RIVS---LTITEGGYCIDDSNGEFMAHLPQI 152
Score = 13.1 bits (22), Expect = 6.6, Method: Composition-based stats.
Identities = 3/6 (50%), Positives = 4/6 (67%), Gaps = 0/6 (0%)
Query 132 EPKIGW 137
EP + W Sbjct 259 EPFVQW 264
Lambda K H a alpha
0.320 0.135 0.421 0.792 4.96
Gapped Lambda K H a alpha sigma
0.267 0.0410 0.140 1.90 42.6 43.6
Effective search space used: 56160
Any another suggestion? Anybody can help me?
So this issue is not fixed? Why did you accept the answer from Mensur? Accepting an answer indicates that your problem has been fixed.
You will need to provide the source files (fasta files you are using) for people to try and reproduce the problem. Use pastebin.com to post sequences.
Unfortunately the problem still exists. Excuse me. I though accepting answer means I read and tried it! However, I posted the link of two of fasta file to easily download them from uniprot site. The links are:
https://www.uniprot.org/uniprot/A0A0E1R3H3.fasta and https://www.uniprot.org/uniprot/A0A0H2UQB5.fasta
I tried all codes with my sequences and no changes happened in the results. Even I tried to change format of sequences in order to fix the problem but I didn't succeed. What should I do? Thanks
No. Accepting an answer means it actually solved the problem described in original post.