使用命令
blastall -i Homo_sapiens.pep.fa -d query_pep.fasta -p blastp -e 1e-5 -b 1 -v 1 -m 8 -o human.out
提示WARNING: ENST00000501597: Could not calculate ungapped Karlin-Altschul parameters due to an invalid query sequence or its translation. Please verify the query sequence(s) and/or filtering options
已去除序列中的终止子符号。
e.g：
>ENST00000620518 ENSG00000081692 JMJD4 1:227732392-227735273(-)
MDRETRALADSHFRGLGVDVPGVGQAPGRVAFVSEPGAFSYADFVRGFLLPNLPCVFSSA
FTQGWGSRRRWVTPAGRPDFDHLLRTYGDVVVPVANCGVQEYNSNPKEHMTLRDYITYWK
EYIQAGYSSPRGCLYLKDWHLCRDFPVEDVFTLPVYFSSDWLNEFWDALDVDDYRFVYAG
e.g.2
>AT1G16710.1 | Symbols: HAC12 | histone acetyltransferase of the CBP family 12 | chr1:5714692-5721782 FORWARD LENGTH=1706
MNVQAHMSGQRSGQVPNQGTVPQNNGNSQMQNLVGSNGAATAVTGAGAATGSGTGVRPSRNIVGAMDHDIMKLRQYMQTLVFNMLQQRQPSPADAASKAK
YMDVARRLEEGLFKMAVTKEDYMNRSTLESRITSLIKGRQINNYNQRHANSSSVGTMIPTPGLSQTAGNPNLMVTSSVDATIVGNTNITSTALNTGNPLI
AGGMHGGNMSNGYQHSSRNFSLGSGGSMTSMGAQRSTAQMIPTPGFVNSVTNNNSGGFSAEPTIVPQSQQQQQRQHTGGQNSHMLSNHMAAGVRPDMQSK
PSGAANSSVNGDVGANEKIVDSGSSYTNASKKLQQGNFSLLSFCPDDLISGQHIESTFHISGEGYSTTNPDPFDGAITSAGTGTKAHNINTASFQPVSRV
应该如何处理？
