changeset 10:e73da4ffaa6a draft

Uploaded v0.0.20, preview 1. Adds unit tests for BLASTN and TBLASTX.
author peterjc
date Mon, 29 Apr 2013 09:19:50 -0400
parents b2eb4f48122f
children 579acb6fb8a8
files test-data/blastn_rhodopsin_vs_three_human.tabular test-data/tblastx_rhodopsin_vs_three_human_ext.tabular test-data/three_human_mRNA.fasta tools/ncbi_blast_plus/ncbi_blast_plus.txt tools/ncbi_blast_plus/ncbi_blastn_wrapper.xml tools/ncbi_blast_plus/ncbi_tblastx_wrapper.xml
diffstat 6 files changed, 334 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/blastn_rhodopsin_vs_three_human.tabular	Mon Apr 29 09:19:50 2013 -0400
@@ -0,0 +1,7 @@
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	92.10	1050	77	6	1	1047	88	1134	0.0	1474
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	91.59	333	28	0	1	333	118	450	9e-133	 460
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	91.36	243	19	2	3127	3368	782	1023	7e-94	 331
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	94.22	173	10	0	1410	1582	448	620	8e-74	 265
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	92.98	171	10	2	2854	3023	615	784	8e-69	 248
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	91.58	962	75	6	1	959	118	1076	0.0	1323
+gi|18148870|dbj|AB062417.1|	ENA|BC112106|BC112106.1	87.55	1052	121	10	1	1047	88	1134	0.0	1208
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/tblastx_rhodopsin_vs_three_human_ext.tabular	Mon Apr 29 09:19:50 2013 -0400
@@ -0,0 +1,117 @@
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	97.39	230	6	0	1	690	88	777	0.0	 559	ENA|BC112106|BC112106.1	1215	224	227	0	98.70	1	1	MNGTEGPNFYVPFSNKTGVVRSPFEYPQYYLAEPWQFSMLAAYMFLLIVLGFPINFLTLYVTVQHKKLRTPLNYILLNLAVADLFMVFGGFTTTLYTSLHGYFVFGPTGCNLEGFFATLGGEIALWSLVVLAIERYVVVCKPMSNFRFGENHAIMGVAFTWVMALACAAPPLVGWSRYIPEGMQCSCGIDYYTLKPEVNNESFVIYMFVVHFTIPMIVIFFCYGQLVFTV	MNGTEGPNFYVPFSNATGVVRSPFEYPQYYLAEPWQFSMLAAYMFLLIVLGFPINFLTLYVTVQHKKLRTPLNYILLNLAVADLFMVLGGFTSTLYTSLHGYFVFGPTGCNLEGFFATLGGEIALWSLVVLAIERYVVVCKPMSNFRFGENHAIMGVAFTWVMALACAAPPLAGWSRYIPEGLQCSCGIDYYTLKPEVNNESFVIYMFVVHFTIPMIIIFFCYGQLVFTV	1047	1213
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	94.12	102	6	0	742	1047	829	1134	0.0	 236	ENA|BC112106|BC112106.1	510	96	100	0	98.04	1	1	KEVTRMVIIMVIAFLICWVPYASVAFYIFTHQGSNFGPIFMTLPAFFAKSSSIYNPVIYIMMNKQFRNCMLTTLCCGKNPLGDDEASTTGSKTETSQVAPA*	KEVTRMVIIMVIAFLICWVPYASVAFYIFTHQGSNFGPIFMTIPAFFAKSAAIYNPVIYIMMNKQFRNCMLTTICCGKNPLGDDEASATVSKTETSQVAPA*	1047	1213
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	91.22	148	13	0	1046	603	1133	690	0.0	 308	ENA|BC112106|BC112106.1	667	135	138	0	93.24	-2	-3	*AGATWLVSVLEPVVEASSSPSGFLPQQRVVSMQFRNCLFIMM*MTGL*MEEDFAKNAGSVMKMGPKLDPWWVKM*NATLA*GTQQIRKAMTMMMTMRVTSFSAFWVVADSCCWAAASLTVKTSCP*QKKMTIMGMVKWTTNM*MTKD	*AGATWLVSVLDTVAEASSSPSGFLPQQMVVSMQFRNCLFIMI*MTGL*MAALLAKNAGMVMKMGPKLEPWWVKM*NATLA*GTQQIRKAMTMMMTMRVTSFSAFCVVADSCCWAAASLTVKTSCP*QKKMIIMGMVKWTTNM*MTKD	1047	1213
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	94.32	88	5	0	566	303	653	390	0.0	 207	ENA|BC112106|BC112106.1	447	83	84	0	95.45	-2	-3	IPHEHCMPSGMYLDQPTRGGAAQASAMTQVKATPIMAWFSPKRKLLMGLHTTTYRSMARTTKDHRAISPPSVAKKPSKLHPVGPKTKY	IPHEHCRPSGMYLDQPASGGAAQASAMTQVKATPMMAWFSPKRKLLMGLHTTTYRSMARTTKDHRAISPPRVAKKPSKLHPVGPKTKY	1047	1213
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	96.34	82	3	0	248	3	335	90	0.0	 182	ENA|BC112106|BC112106.1	393	79	80	0	97.56	-2	-3	SATARLSRM*LRGVRSFLCWTVT*SVRKLMGKPSTMSRNM*AASMENCHGSAR*YCGYSKGLRTTPVLLEKGT*KFGPSVPF	SATARLSRM*LRGVRSFLCWTVT*SVRKLMGKPSTISRNM*AASMENCHGSAR*YCGYSKGLRTTPVALEKGT*KLGPSVPF	1047	1213
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	83.33	204	34	0	18	629	105	716	4e-158	 404	ENA|BC112106|BC112106.1	877	170	174	0	85.29	3	3	PELLRALLQQNGCGTQPLRVPTVLPG*AMAVLHAGRLHVPAHRAWLPHQLPHALRHGPAQEAAHASQLHPAQPGRG*PLHGLRWLHHHPLHLSAWILCLWAHRMQFGGLLCHTGR*NCPVVFGGPGH*AVRGGV*AHEQLPLWGEPCHNGRRFHLGHGTGLRCTPPRWLVQVHP*RHAVFMRDRLLHTQARSQQRVLCHLHVRG	P*LLRALLQCDGCGTQPLRVPTVLPG*AMAVLHAGRLHVSADRAGLPHQLPHALRHRPAQEAAHASQLHPAQPSRG*PLHGPRWLHQHPLHLSAWILRLRAHRMQFGGLLCHPGR*NCPVVLGGPGHRAVRGGV*AHEQLPLRGEPCHHGRCLHLGHGAGLRRTPTRRLVQVHPRGPAVLVWNRLLHAQAGGQQRVFCHLHVRG	1047	1213
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	75.28	89	22	0	780	1046	867	1133	4e-158	 161	ENA|BC112106|BC112106.1	346	67	72	0	80.90	3	3	FPDLLGALRQRGILHLHPPGVQLWPHLHDTPGVLRKVLLHLQPCHLHHDEQAVPELHAHYPLLWQEPTG**RGFHNRFQDGDQPGGTGL	FPDLLGALRQRGILHLHPPGLQLRSHLHDHPSVLCQERRHLQPCHLYHDEQAVPELHAHHHLLRQEPTG*R*GLCYRVQDGDEPGGPGL	1047	1213
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	81.28	203	38	0	609	1	696	88	5e-153	 360	ENA|BC112106|BC112106.1	782	165	175	0	86.21	-1	-2	KGLVVDFWLECVVVDPA*TLHAFRDVPGPTNEGGCSAGQCHDPGESDAHYGMVLPKAEVAHGLTHHHVPLNGQDHQRPQGNFTAQCGKEALQIASCGPKDKVSMQRGVEGGGEATEDHEEVSHGQVEQDVVERRAQLLVLDRDVEREEVDGEAKHDEQEHVGGQHGELPWLSQVVLWVLEGAAYHTRFVGEGHVEVRALRPVH	KRLVVDLRLERVVVDSTRALQALGDVPGPAGEWGCGAGQRHDPGEGNAHDGMVLPEAEVAHGLTHHHVPLDGQDHQGPQGNFTAQGGKEALQIASCGPEDEVSMQRGVEGAGEAT*DHEEVSHG*VEQDVVERRAQLLVLDGDVEREEVDGEAQHDQQKHVGGQHGELPWLSQVVLWVLEGAAYHTRRIGEGHVEVRAFCAIH	1047	1213
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	80.60	67	13	0	916	716	1003	803	5e-153	 135	ENA|BC112106|BC112106.1	290	54	60	0	89.55	-3	-1	R*QGCRWRRTLRRTPGVS*RWGQSWTPGG*RCRMPRWRRAPNRSGKQ*P***PCE*PPSQPSGWWLT	R*QGCRWRRSWQRTLGWS*RWDRSWSPGG*RCRMPRWRRAPSRSGKR*P***PCG*PPSLPSVWWLT	1047	1213
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	70.27	37	11	0	1047	937	1134	1024	5e-153	64.2	ENA|BC112106|BC112106.1	134	26	29	0	78.38	-1	-2	LGRCHLAGLRLGTGCGSLVITQWVLATAEGSEHAVPE	LGRGHLARLRLGHGSRGLIVTQWVLAAADGGEHAVPE	1047	1213
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	100.00	7	0	0	646	626	733	713	5e-153	24.0	ENA|BC112106|BC112106.1	46	7	7	0	100.00	-3	-1	WGW*SGP	WGW*SGP	1047	1213
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	89.23	65	7	0	460	266	547	353	4e-105	 167	ENA|BC112106|BC112106.1	360	58	63	0	96.92	-3	-1	WHGSPQSGSCSWAYTPPRTAQWPGPPKTTGQFHRPVWQRSPPNCILWAQRQSIHAERCRGWW*SH	WHGSPRSGSCSWAYTPPRTARWPGPPRTTGQFHRPGWQRSPPNCILWARRRSIHAERCRGCW*SH	1047	1213
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	89.58	48	5	0	184	41	271	128	4e-105	 104	ENA|BC112106|BC112106.1	222	43	43	0	89.58	-3	-1	*RRA*GS*WGSQAR*AGTCRRPAWRTAMAQPGSTVGTRRGCVPHPFCW	*RRA*GS*WGSPARSAETCRRPAWRTAMAQPGSTVGTRRGCVPHPSHW	1047	1213
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	77.78	45	10	0	882	748	969	835	4e-105	93.9	ENA|BC112106|BC112106.1	199	35	39	0	86.67	-1	-2	EERRECHEDGAKVGPLVGEDVECHAGVGHPTDQESNDHDDDHASD	KERWDGHEDGTEVGALVGEDVECHAGVGHPADQESDDHDDDHAGD	1047	1213
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	67.86	28	9	0	1045	962	1132	1049	4e-105	51.9	ENA|BC112106|BC112106.1	107	19	20	0	71.43	-3	-1	RPVPPGWSPSWNRLWKPRHHPVGSCHSR	RPGPPGSSPSWTR*QRPHRHPVGSCRSR	1047	1213
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	59.09	22	9	0	586	521	673	608	4e-105	33.1	ENA|BC112106|BC112106.1	66	13	16	0	72.73	-3	-1	A*VCSSRSRMNTACLQGCTWTN	A*ACSSRFHTSTAGPRGCTWTS	1047	1213
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	81.40	86	16	0	296	553	383	640	2e-87	 185	ENA|BC112106|BC112106.1	399	70	71	0	82.56	2	2	CMDTLSLGPQDAIWRASLPHWAVKLPCGLWWSWPLSGTWWCVSP*ATSALGRTMP*WASLSPGSWHWPALHPPSLVGPGTSLKACS	CMDTSSSGPQDAIWRASLPPWAVKLPCGPWWSWPSSGTWWCVSP*ATSASGRTMPSWALPSPGSWRWPAPHPHSPAGPGTSPRACS	1047	1213
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	84.38	32	5	0	11	106	98	193	2e-87	74.8	ENA|BC112106|BC112106.1	157	27	30	0	93.75	2	2	RRARTSTCPSPTKRVWYAAPSSTHSTTWLSHG	QKALTSTCPSPMRRVWYAAPSSTHSTTWLSHG	1047	1213
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	71.43	35	10	0	941	1045	1028	1132	2e-87	61.6	ENA|BC112106|BC112106.1	128	25	25	0	71.43	2	2	GTACSLPSAVARTHWVMTRLPQPVPRRRPARWHRP	GTACSPPSAAARTHWVTMRPLLPCPRRRRARWPRP	1047	1213
+gi|57163782|ref|NM_001009242.1|	ENA|BC112106|BC112106.1	94.44	18	1	0	794	847	881	934	2e-87	50.1	ENA|BC112106|BC112106.1	103	17	17	0	94.44	2	2	GCPTPAWHSTSSPTRGPT	GCPTPAWHSTSSPTRAPT	1047	1213
+gi|2734705|gb|U59921.1|BBU59921	ENA|BC112106|BC112106.1	83.61	238	39	0	18	731	64	777	0.0	 507	ENA|BC112106|BC112106.1	1102	199	219	0	92.02	3	1	GGL*RAANMNGTEGPNFYIPMSNKTGVVRSPFEYPQYYLAEPWQYSILCAYMFLLILLGFPINFMTLYVTIQHKKLRTPLNYILLNLAFANHFMVLCGFTVTMYSSMNGYFILGATGCYVEGFFATLGGEIALWSLVVLAIERYVVVCKPMSNFRFSENHAVMGVAFTWIMALSCAVPPLLGWSRYIPEGMQCSCGVDYYTLKPEVNNESFVIYMFVVHFTIPLIIIFFCYGRLVCTV	GSATRATAMNGTEGPNFYVPFSNATGVVRSPFEYPQYYLAEPWQFSMLAAYMFLLIVLGFPINFLTLYVTVQHKKLRTPLNYILLNLAVADLFMVLGGFTSTLYTSLHGYFVFGPTGCNLEGFFATLGGEIALWSLVVLAIERYVVVCKPMSNFRFGENHAIMGVAFTWVMALACAAPPLAGWSRYIPEGLQCSCGIDYYTLKPEVNNESFVIYMFVVHFTIPMIIIFFCYGQLVFTV	1574	1213
+gi|2734705|gb|U59921.1|BBU59921	ENA|BC112106|BC112106.1	82.35	85	15	0	783	1037	829	1083	0.0	 188	ENA|BC112106|BC112106.1	404	70	82	0	96.47	3	1	KEVTRMVIIMVVFFLICWVPYASVAFFIFSNQGSEFGPIFMTVPAFFAKSSSIYNPVIYIMLNKQFRNCMITTLCCGKNPFGEDD	KEVTRMVIIMVIAFLICWVPYASVAFYIFTHQGSNFGPIFMTIPAFFAKSAAIYNPVIYIMMNKQFRNCMLTTICCGKNPLGDDE	1574	1213
+gi|2734705|gb|U59921.1|BBU59921	ENA|BC112106|BC112106.1	70.96	303	88	0	925	17	971	63	2e-130	 435	ENA|BC112106|BC112106.1	945	215	230	0	75.91	-2	-3	AKKAGTVMKMGPNSEP*LLKMKKATEA*GTQQIRKKTTMMITILVTSFSAFWVVADSCC*AAASFTVHTRRP*QKKMIIRGMVKWTTNM*MTKDSLLTSGFRV**STPQEHCMPSGMYLDHPRSGGTAQDKAIIQVNATPMTAWFSLNRKLLMGLQTTTYRSMAKTTRDQRAISPPRVAKKPST*QPVAPRMKYPFIEEYIVTVNPHRTMKWLAKAKFSRI*FKGVRSFLCWMVTYKVMKLIGNPSRMSRNMYAHRMEYCHGSAR*Y*GYSKGLRTTPVLLDMGM*KFGPSVPFILAAL*SPP	AKNAGMVMKMGPKLEPWWVKM*NATLA*GTQQIRKAMTMMMTMRVTSFSAFCVVADSCCWAAASLTVKTSCP*QKKMIIMGMVKWTTNM*MTKDSLLTSGLSV**SIPHEHCRPSGMYLDQPASGGAAQASAMTQVKATPMMAWFSPKRKLLMGLHTTTYRSMARTTKDHRAISPPRVAKKPSKLHPVGPKTKYPCREV*RVLVKPPRTMKRSATARLSRM*LRGVRSFLCWTVT*SVRKLMGKPSTISRNM*AASMENCHGSAR*YCGYSKGLRTTPVALEKGT*KLGPSVPFMAVALVADP	1574	1213
+gi|2734705|gb|U59921.1|BBU59921	ENA|BC112106|BC112106.1	72.22	18	5	0	1027	974	1073	1020	2e-130	35.0	ENA|BC112106|BC112106.1	70	13	13	0	72.22	-2	-3	PKGFLPQHRVVIMQLRNC	PSGFLPQQMVVSMQFRNC	1574	1213
+gi|2734705|gb|U59921.1|BBU59921	ENA|BC112106|BC112106.1	55.32	188	84	0	605	42	651	88	7e-89	 245	ENA|BC112106|BC112106.1	529	104	134	0	71.28	-1	-2	SAGALHALGDVPGPSKEWRNSTGQSHYPGERYAHDGMVLTKSEVAHGLTDDHVSFNGQDHQGPKGDFTTKGSEEAFNIATGGSEDEVSVH*GVHCHRESTQDHEVVGKGQIQQDIV*GCPELLVLDGDVQGHEVDWEPK*NEQEHVRAQNGILPWLCQVILRVFKGASHHPSLVGHGYVKVWAFCSVH	STRALQALGDVPGPAGEWGCGAGQRHDPGEGNAHDGMVLPEAEVAHGLTHHHVPLDGQDHQGPQGNFTAQGGKEALQIASCGPEDEVSMQRGVEGAGEAT*DHEEVSHG*VEQDVVERRAQLLVLDGDVEREEVDGEAQHDQQKHVGGQHGELPWLSQVVLWVLEGAAYHTRRIGEGHVEVRAFCAIH	1574	1213
+gi|2734705|gb|U59921.1|BBU59921	ENA|BC112106|BC112106.1	61.11	72	28	0	1037	822	1083	868	7e-89	91.3	ENA|BC112106|BC112106.1	193	44	54	0	75.00	-1	-2	IVFSKGILAAAQGGDHAVTELLVEHDVDDGVVDGRTLGKESWDGHEDGAELRALIAEDEESDRGVGDPTDKE	LIVTQWVLAAADGGEHAVPELLVHHDIDDRVVDGGALGKERWDGHEDGTEVGALVGEDVECHAGVGHPADQE	1574	1213
+gi|2734705|gb|U59921.1|BBU59921	ENA|BC112106|BC112106.1	49.02	204	104	0	29	640	75	686	4e-78	 197	ENA|BC112106|BC112106.1	424	100	132	0	64.71	2	3	KSRQYERNRRPKLLHTHVQQDWGGAKPL*IPSVLPGRAMAIFHSVRVHVPAHSTWVPNQLHDLVRHHPAQEAPDTLKLYPAEFGLCQPLHGPVWIHGDNVLLNERILHPRSHRLLC*RLLRYPWW*NRPLVPGGLGH*TIRGRL*AHEQLPI**EPCRHGRSVHLDNGFVLCCSSTPWMVQVHPRGHAVLLRSRLLHPEARGQQ	KGHSHEWHRRP*LLRALLQCDGCGTQPLRVPTVLPG*AMAVLHAGRLHVSADRAGLPHQLPHALRHRPAQEAAHASQLHPAQPSRG*PLHGPRWLHQHPLHLSAWILRLRAHRMQFGGLLCHPGR*NCPVVLGGPGHRAVRGGV*AHEQLPLRGEPCHHGRCLHLGHGAGLRRTPTRRLVQVHPRGPAVLVWNRLLHAQAGGQQ	1574	1213
+gi|2734705|gb|U59921.1|BBU59921	ENA|BC112106|BC112106.1	66.04	53	18	0	860	1018	906	1064	4e-78	85.8	ENA|BC112106|BC112106.1	181	35	39	0	73.58	2	3	LHLQQSGL*VRPHLHDRPSFLCQEFFHLQPRHLHHAQQAVP*LHDHHPVLRQE	LHLHPPGLQLRSHLHDHPSVLCQERRHLQPCHLYHDEQAVPELHAHHHLLRQE	1574	1213
+gi|2734705|gb|U59921.1|BBU59921	ENA|BC112106|BC112106.1	44.44	27	15	0	689	769	735	815	4e-78	32.2	ENA|BC112106|BC112106.1	64	12	14	0	51.85	2	3	DYHFLLLWPPGVHCERGCSSTARVRHH	DYHLFLLWAARLHRQGGRCPAAGVSHH	1574	1213
+gi|2734705|gb|U59921.1|BBU59921	ENA|BC112106|BC112106.1	47.47	198	104	0	633	40	679	86	4e-65	 177	ENA|BC112106|BC112106.1	380	94	106	0	53.54	-3	-1	PRASGCSSRLRRSTACPRGCTWTIQGVEEQHRTKPLSR*TLRP*RHGSH*IGSCSWAYRRPRIVQWPRPPGTKGRFHHQG*RRSLQHSNRWLRG*SIRSLRSTLSP*IHTGP*SGWQRPNSAGYSLRVSGASCAGW*RTRS*S*LGTQVE*AGTCTRTEWNIAMALPGNTEGIQRGFAPPQSCWTWVCKSLGLLFRSY	PPA*ACSSRFHTSTAGPRGCTWTSRRVGVRRRPAP*PR*RQRP*WHGSPRSGSCSWAYTPPRTARWPGPPRTTGQFHRPGWQRSPPNCILWARRRSIHAERCRGCW*SHLGP*RGQPRLG*AGCS*EACAASCAGR*RRA*GS*WGSPARSAETCRRPAWRTAMAQPGSTVGTRRGCVPHPSHWRRARRS*GLLCHSW	1574	1213
+gi|2734705|gb|U59921.1|BBU59921	ENA|BC112106|BC112106.1	68.09	47	15	0	1017	877	1063	923	4e-65	80.3	ENA|BC112106|BC112106.1	169	32	37	0	78.72	-3	-1	SCRSTGW*SCSYGTAC*A*CR*RGCRWKNSWQRKLGRS*RWGRTQSP	SCRSRWW*ACSSGTACSS*YR*QGCRWRRSWQRTLGWS*RWDRSWSP	1574	1213
+gi|2734705|gb|U59921.1|BBU59921	ENA|BC112106|BC112106.1	57.89	114	48	0	265	606	311	652	3e-46	 137	ENA|BC112106|BC112106.1	294	66	70	0	61.40	1	2	SC*IWPLPTTSWSCVDSR*QCTPQ*TDTSSSEPPVAMLKASSLPLVVKSPFGPWWSWPLNDTWSSVSP*ATSDLVRTMPSWA*RSPG*WLCPVLFLHSLDGPGTSPRACSAPAE	SCST*PWLTSSWS*VASPAPSTPLCMDTSSSGPQDAIWRASLPPWAVKLPCGPWWSWPSSGTWWCVSP*ATSASGRTMPSWALPSPGSWRWPAPHPHSPAGPGTSPRACSARVE	1574	1213
+gi|2734705|gb|U59921.1|BBU59921	ENA|BC112106|BC112106.1	46.30	54	29	0	19	180	65	226	3e-46	52.4	ENA|BC112106|BC112106.1	108	25	28	0	51.85	1	2	GDFKEPPI*TEQKAQTFTYPCPTRLGWCEAPLNTLSITWQSHGNIPFCARTCSC	GQPQGPQP*MAQKALTSTCPSPMRRVWYAAPSSTHSTTWLSHGSSPCWPPTCFC	1574	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	96.40	111	4	0	1	333	118	450	0.0	 264	ENA|BC112106|BC112106.1	570	107	109	0	98.20	1	1	VPFSNKTGVVRSPFEHPQYYLAEPWQFSMLAAYMFLLIVLGFPINFLTLYVTVQHKKLRTPLNYILLNLAVADLFMVFGGFTTTLYTSLHGYFVFGPTGCNLEGFFATLGG	VPFSNATGVVRSPFEYPQYYLAEPWQFSMLAAYMFLLIVLGFPINFLTLYVTVQHKKLRTPLNYILLNLAVADLFMVLGGFTSTLYTSLHGYFVFGPTGCNLEGFFATLGG	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	92.31	65	5	0	3174	3368	829	1023	0.0	 151	ENA|BC112106|BC112106.1	324	60	64	0	98.46	3	1	KEVTRMVIIMVIAFLICWLPYAGVAFYIFTHQGSNFGPIFMTLPAFFAKSSSIYNPVIYIMMNKQ	KEVTRMVIIMVIAFLICWVPYASVAFYIFTHQGSNFGPIFMTIPAFFAKSAAIYNPVIYIMMNKQ	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	96.43	56	2	0	2855	3022	616	783	0.0	 141	ENA|BC112106|BC112106.1	302	54	56	0	100.00	2	1	RYIPEGMQCSCGIDYYTLKPEVNNESFVIYMFVVHFTIPMIVIFFCYGQLVFTVKE	RYIPEGLQCSCGIDYYTLKPEVNNESFVIYMFVVHFTIPMIIIFFCYGQLVFTVKE	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	93.22	59	4	0	1404	1580	442	618	0.0	 138	ENA|BC112106|BC112106.1	295	55	56	0	94.92	3	1	LAGEIALWSLVVLAIERYVVVCKPMSNFRFGENHAIMGLALTWVMALACAAPPLVGWSR	LGGEIALWSLVVLAIERYVVVCKPMSNFRFGENHAIMGVAFTWVMALACAAPPLAGWSR	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	92.00	25	2	0	4222	4296	1021	1095	0.0	64.3	ENA|BC112106|BC112106.1	134	23	24	0	96.00	1	1	QFRNCMLTTLCCGKNPLGDDEASTT	QFRNCMLTTICCGKNPLGDDEASAT	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	88.89	9	1	0	3128	3154	783	809	0.0	22.6	ENA|BC112106|BC112106.1	43	8	8	0	88.89	2	3	GSCPAAGVS	GRCPAAGVS	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	78.38	111	24	0	333	1	450	118	7e-171	 212	ENA|BC112106|BC112106.1	458	87	97	0	87.39	-3	-2	TSQGGKKALQIASRRPEDKVSMQGGVEGGGEASEDHEEVGHGQVEQDVVERRTQLLVLNRDIEREEVDGEAEHDQQKHVGGQHGELPWLGQVVLWMLEGTAHHACLVGEGH	TAQGGKEALQIASCGPEDEVSMQRGVEGAGEAT*DHEEVSHG*VEQDVVERRAQLLVLDGDVEREEVDGEAQHDQQKHVGGQHGELPWLSQVVLWVLEGAAYHTRRIGEGH	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	88.75	80	9	0	3367	3128	1022	783	7e-171	 161	ENA|BC112106|BC112106.1	347	71	73	0	91.25	-2	-3	CLFIMI*MTGL*MEDDLAKNAGRVMKMGPKLEPWWVKM*NATPAYGSQQIRKAMTMMMTIRVTSFSAFWVVADSCCWAAA	CLFIMI*MTGL*MAALLAKNAGMVMKMGPKLEPWWVKM*NATLA*GTQQIRKAMTMMMTMRVTSFSAFCVVADSCCWAAA	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	93.33	60	4	0	1582	1403	620	441	7e-171	 136	ENA|BC112106|BC112106.1	291	56	57	0	95.00	-2	-3	YLDQPTSGGAAQASAMTQVRARPMMAWFSPKRKLLMGLHTTTYRSMARTTKDHRAISPAR	YLDQPASGGAAQASAMTQVKATPMMAWFSPKRKLLMGLHTTTYRSMARTTKDHRAISPPR	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	91.07	56	5	0	3021	2854	782	615	7e-171	 119	ENA|BC112106|BC112106.1	254	51	52	0	92.86	-3	-3	SLTVKTSCP*QKNMTIIGMVKWTTNM*MTKDSLLTSGLRV**SIPHEHCMPSGMYL	SLTVKTSCP*QKKMIIMGMVKWTTNM*MTKDSLLTSGLSV**SIPHEHCRPSGMYL	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	84.62	26	4	0	4301	4224	1100	1023	7e-171	52.8	ENA|BC112106|BC112106.1	109	22	23	0	88.46	-1	-3	EAVVEASSSPSGFLPQQRVVSMQFRN	DTVAEASSSPSGFLPQQMVVSMQFRN	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	95.83	72	3	0	218	3	335	120	8e-142	 152	ENA|BC112106|BC112106.1	327	69	69	0	95.83	-1	-3	SATARLSRM*LRGVRSFLC*TVT*SVRKLMGKPSTISRNM*AASMENCHGSAR*YCGCSKGLRTTPVLLEKG	SATARLSRM*LRGVRSFLCWTVT*SVRKLMGKPSTISRNM*AASMENCHGSAR*YCGYSKGLRTTPVALEKG	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	77.78	63	14	0	3368	3180	1023	835	8e-142	 125	ENA|BC112106|BC112106.1	267	49	53	0	84.13	-1	-2	LLVHHDIDDRVVDGGRLGKECREGHEDGAKVRALVGEDVECHAGIRQPTD*ESNDHDDDHTGD	LLVHHDIDDRVVDGGALGKERWDGHEDGTEVGALVGEDVECHAGVGHPADQESDDHDDDHAGD	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	85.11	47	7	0	1544	1404	582	442	8e-142	 108	ENA|BC112106|BC112106.1	231	40	43	0	91.49	-1	-2	QCHDPGEGKAHDGVVLPEAEVAHGLAYHHVPLDGQDHQGPQGNFTCK	QRHDPGEGNAHDGMVLPEAEVAHGLTHHHVPLDGQDHQGPQGNFTAQ	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	80.36	56	11	0	3022	2855	783	616	8e-142	 101	ENA|BC112106|BC112106.1	216	45	48	0	85.71	-2	-2	LLDGEDQLSIAKEYDNHRDGEVDHEHVDDKGLVVDLRLEGVVVDST*ALHALWDVP	LLDGEDELPIAEKDDNHGDGEVDHEHVDDKRLVVDLRLERVVVDSTRALQALGDVP	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	75.86	29	7	0	325	239	442	356	8e-142	58.3	ENA|BC112106|BC112106.1	121	22	25	0	86.21	-2	-1	GWQKSPPDCIP*ARRQSIHAGRCRGWW*S	GWQRSPPNCILWARRRSIHAERCRGCW*S	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	86.36	22	3	0	4287	4222	1086	1021	8e-142	48.7	ENA|BC112106|BC112106.1	100	19	20	0	90.91	-3	-2	GLIVTQWVLATAEGSEHAVPEL	GLIVTQWVLAAADGGEHAVPEL	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	90.91	11	1	0	3159	3127	814	782	8e-142	31.3	ENA|BC112106|BC112106.1	62	10	10	0	90.91	-3	-1	WWLTPAAGQLP	WWLTPAAGQRP	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	81.03	58	11	0	2854	3027	615	788	2e-122	 128	ENA|BC112106|BC112106.1	274	47	51	0	87.93	1	3	QVHPRGHAVLMWNRLLHPQAGGQQRVLCHLHVRGPLHHPYDCHILLLWTAGLHRQGGK	QVHPRGPAVLVWNRLLHAQAGGQQRVFCHLHVRGPLHHPHDYHLFLLWAARLHRQGGR	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	83.33	60	10	0	1403	1582	441	620	2e-122	 125	ENA|BC112106|BC112106.1	267	50	51	0	85.00	2	3	PCR*NCPVVLGGPGHRAVRGGMQAHEQLPLRGEPRHHGPCPHLGHGTGLRRAPASRLVQV	PGR*NCPVVLGGPGHRAVRGGV*AHEQLPLRGEPCHHGRCLHLGHGAGLRRTPTRRLVQV	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	80.60	67	13	0	3	203	120	320	2e-122	 119	ENA|BC112106|BC112106.1	254	54	55	0	82.09	3	3	ALLQQDRRGAQSLRASTVLPGRAMAVLHAGRLHVSADRARLPHQLPHALCHGSAQEAAYASQLHPAQ	ALLQCDGCGTQPLRVPTVLPG*AMAVLHAGRLHVSADRAGLPHQLPHALRHRPAQEAAHASQLHPAQ	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	91.30	23	2	0	4220	4288	1019	1087	2e-122	53.8	ENA|BC112106|BC112106.1	111	21	21	0	91.30	2	2	SSSGTACSLPSAVARTHWVTMRP	SSSGTACSPPSAAARTHWVTMRP	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	78.26	23	5	0	266	334	383	451	2e-122	48.3	ENA|BC112106|BC112106.1	99	18	19	0	82.61	2	2	CMDTLSSGLRDAIWRAFLPPWEV	CMDTSSSGPQDAIWRASLPPWAV	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	80.00	20	4	0	3308	3367	963	1022	2e-122	46.0	ENA|BC112106|BC112106.1	94	16	17	0	85.00	2	3	ILCQVVLHLQPCHLYHDEQA	VLCQERRHLQPCHLYHDEQA	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	83.33	18	3	0	3226	3279	881	934	2e-122	40.5	ENA|BC112106|BC112106.1	82	15	15	0	83.33	1	2	GCRMPAWHSTSSPTRALT	GCPTPAWHSTSSPTRAPT	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	80.39	51	10	0	154	2	271	119	1e-92	97.3	ENA|BC112106|BC112106.1	206	41	42	0	82.35	-2	-1	*HRA*GS*WGSRARSAETCRRPAWRTAMARPGSTVDARRDCAPRLSCWRRA	*RRA*GS*WGSPARSAETCRRPAWRTAMAQPGSTVGTRRGCVPHPSHWRRA	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	73.68	57	15	0	3366	3196	1021	851	1e-92	95.9	ENA|BC112106|BC112106.1	203	42	46	0	80.70	-3	-1	ACSS*YR*QGCRWRTTWQRMPGGS*RWGQS*SPGG*RCRMPRRHTAANRLGKQ*P**	ACSS*YR*QGCRWRRSWQRTLGWS*RWDRSWSPGG*RCRMPRWRRAPSRSGKR*P**	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	81.82	44	8	0	1521	1390	559	428	1e-92	95.0	ENA|BC112106|BC112106.1	201	36	36	0	81.82	-3	-1	QGP*WRGSPRSGSCSWACIPPRTARWPGPPRTTGQFHLQGSQLS	QRP*WHGSPRSGSCSWAYTPPRTARWPGPPRTTGQFHRPGWQRS	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	95.00	20	1	0	332	273	449	390	1e-92	52.4	ENA|BC112106|BC112106.1	108	19	20	0	100.00	-1	-3	PPRVAKKPSRLHPVGPKTKY	PPRVAKKPSKLHPVGPKTKY	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	86.36	22	3	0	4288	4223	1087	1022	1e-92	48.3	ENA|BC112106|BC112106.1	99	19	19	0	86.36	-2	-1	RPHRHPVGSCHSRG**ACSSGT	RPHRHPVGSCRSRWW*ACSSGT	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	61.54	26	10	0	2912	2835	673	596	1e-92	37.7	ENA|BC112106|BC112106.1	76	16	17	0	65.38	-1	-1	A*GCSSRFHMSTACPLGCT*GPAEIG	A*ACSSRFHTSTAGPRGCTWTSRRVG	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	69.57	23	7	0	3023	2955	784	716	1e-92	33.1	ENA|BC112106|BC112106.1	66	16	18	0	78.26	-1	-1	PP*R*RPAVHSKRI*QS*GW*SG	PP*R*RRAAHSRKR**SWGW*SG	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	76.67	30	7	0	2	91	119	208	6e-43	68.0	ENA|BC112106|BC112106.1	142	23	24	0	80.00	2	2	CPSPTRQAWCAVPSSIHSTTWPSHGSSPCW	CPSPMRRVWYAAPSSTHSTTWLSHGSSPCW	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	91.11	45	4	0	1411	1545	449	583	6e-43	68.0	ENA|BC112106|BC112106.1	142	41	41	0	91.11	1	2	VKLPCGPWWSWPSSGTWWYASP*ATSASGRTTPSWALPSPGSWHW	VKLPCGPWWSWPSSGTWWCVSP*ATSASGRTMPSWALPSPGSWRW	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	77.27	22	5	0	267	332	384	449	6e-43	45.1	ENA|BC112106|BC112106.1	92	17	18	0	81.82	3	3	AWILCLRAYGMQSGGLFCHPGR	AWILRLRAHRMQFGGLLCHPGR	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	86.36	22	3	0	4224	4289	1023	1088	6e-43	44.1	ENA|BC112106|BC112106.1	90	19	20	0	90.91	3	3	VPELHAHYPLLWQEPTG*R*GL	VPELHAHHHLLRQEPTG*R*GL	4301	1213
+gi|283855845|gb|GQ290303.1|	ENA|BC112106|BC112106.1	83.33	12	2	0	2856	2891	617	652	6e-43	25.4	ENA|BC112106|BC112106.1	49	10	10	0	83.33	3	2	GTSQRACSAHVE	GTSPRACSARVE	4301	1213
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	95.91	220	9	0	1	660	118	777	0.0	 526	ENA|BC112106|BC112106.1	1142	211	217	0	98.64	1	1	VPFSNKTGVVRSPFEYPQYYLAEPWQFSMLAAYMFLLIVLGFPINFLTLYVTVQHKKLRTPLNYILLNLAVANLFMVFGGFTTTLYTSMHGYFVFGATGCNLEGFFATLGGEIALWSLVVLAIERYVVVCKPMSNFRFGENHAIMGLAFTWVMALACAAPPLAGWSRYIPEGMQCSCGIDYYTLKPEVNNESFVIYMFVVHFTIPMIVIFFCYGQLVFTV	VPFSNATGVVRSPFEYPQYYLAEPWQFSMLAAYMFLLIVLGFPINFLTLYVTVQHKKLRTPLNYILLNLAVADLFMVLGGFTSTLYTSLHGYFVFGPTGCNLEGFFATLGGEIALWSLVVLAIERYVVVCKPMSNFRFGENHAIMGVAFTWVMALACAAPPLAGWSRYIPEGLQCSCGIDYYTLKPEVNNESFVIYMFVVHFTIPMIIIFFCYGQLVFTV	983	1213
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	92.13	89	7	0	712	978	829	1095	0.0	 212	ENA|BC112106|BC112106.1	457	82	88	0	98.88	1	1	KEVTRMVIIMVVAFLICWLPYASVAFYIFTHQGSNFGPVFMTIPAFFAKSSSIYNPVIYIMMNKQFRNCMLTTLCCGKNPLGDDEASTT	KEVTRMVIIMVIAFLICWVPYASVAFYIFTHQGSNFGPIFMTIPAFFAKSAAIYNPVIYIMMNKQFRNCMLTTICCGKNPLGDDEASAT	983	1213
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	91.01	178	16	0	536	3	653	120	1e-178	 353	ENA|BC112106|BC112106.1	765	162	166	0	93.26	-1	-3	IPHEHCMPSGMYLDQPASGGAAQASAMTHVKARPMMAWFSPKRKLLMGLQTTTYRSMARTTRDHRAISPPSVAKKPSRLHPVAPKTKYPCIEVYRVVVKPPKTMKRLATARLSRM*LRGVRSFLCWTVT*SVRKLMGNPSTISRNM*AASMENCQGSAR*YCGYSKGLRTTPVLLEKG	IPHEHCRPSGMYLDQPASGGAAQASAMTQVKATPMMAWFSPKRKLLMGLHTTTYRSMARTTKDHRAISPPRVAKKPSKLHPVGPKTKYPCREV*RVLVKPPRTMKRSATARLSRM*LRGVRSFLCWTVT*SVRKLMGKPSTISRNM*AASMENCHGSAR*YCGYSKGLRTTPVALEKG	983	1213
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	88.32	137	16	0	983	573	1100	690	1e-178	 277	ENA|BC112106|BC112106.1	599	121	124	0	90.51	-1	-3	EAVVDASSSPSGFLPQQRVVSMQFRNCLFIMI*MTGL*MEDDLAKNAGMVMKTGPKLEPWWVKM*NATLA*GSQQIRNATTMMMTMRVTSFSAFWVVADSCCWAAASFTVNTSCP*QKKMTIMGMVKWTTNM*MTKD	DTVAEASSSPSGFLPQQMVVSMQFRNCLFIMI*MTGL*MAALLAKNAGMVMKMGPKLEPWWVKM*NATLA*GTQQIRKAMTMMMTMRVTSFSAFCVVADSCCWAAASLTVKTSCP*QKKMIIMGMVKWTTNM*MTKD	983	1213
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	77.43	319	72	0	3	959	120	1076	4e-174	 593	ENA|BC112106|BC112106.1	1290	247	259	0	81.19	3	3	ALLQQDGCGAQPLRVPAVLPG*ALAVLHAGCLHVSADRARIPHQLPHALRHRPAQEAAHASQLHPAQPGCGQPLHGLWRLHHHPVYLYAWILRLRGHGMQSGGLLCHAGR*NRPVVPGGPGHRAVCGGLQAHEQLPLWGEPRHHGPRLHVGHGTGLRCTPTSRLVQVHPRGHAVLVWD*LLHAQTGGQQRVLRHLHVRGPLHHPHDCHFLLLRTAGVHSEGGGCPAAGVSHHPEGREGSHAHGHHHGRCVPNLLAALRQRGILHLYPPGL*LWPCLHDHPGILRQVILHLQPGHLYHDEQAVPELHAHHPLLWQEPTG*	ALLQCDGCGTQPLRVPTVLPG*AMAVLHAGRLHVSADRAGLPHQLPHALRHRPAQEAAHASQLHPAQPSRG*PLHGPRWLHQHPLHLSAWILRLRAHRMQFGGLLCHPGR*NCPVVLGGPGHRAVRGGV*AHEQLPLRGEPCHHGRCLHLGHGAGLRRTPTRRLVQVHPRGPAVLVWNRLLHAQAGGQQRVFCHLHVRGPLHHPHDYHLFLLWAARLHRQGGRCPAAGVSHHTEGREGGHPHGHHHGHRFPDLLGALRQRGILHLHPPGLQLRSHLHDHPSVLCQERRHLQPCHLYHDEQAVPELHAHHHLLRQEPTG*	983	1213
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	79.07	129	27	0	558	172	675	289	2e-133	 248	ENA|BC112106|BC112106.1	537	102	110	0	85.27	-3	-2	RFERVVVNPTRALHALWDVPGPAG*WGCSAGQCHDPREGEAHDGVVLPKAEVAHGLADHHIPLDGQDHQGPQGDFTAQRGKEALQIASRGPEDEVSMHRGIQGGGEASKDHEEVGHSQVEQDVVERRAQ	RLERVVVDSTRALQALGDVPGPAGEWGCGAGQRHDPGEGNAHDGMVLPEAEVAHGLTHHHVPLDGQDHQGPQGNFTAQGGKEALQIASCGPEDEVSMQRGVEGAGEAT*DHEEVSHG*VEQDVVERRAQ	983	1213
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	76.83	82	19	0	963	718	1080	835	2e-133	 159	ENA|BC112106|BC112106.1	342	63	70	0	85.37	-3	-2	VITQWVLATAEGGEHAVPELLVHHDIDDRVVDGG*LGEECRDGHEDRAKVRALVGKDVECHAGVGQPTD*ERNDHDDDHARD	IVTQWVLAAADGGEHAVPELLVHHDIDDRVVDGGALGKERWDGHEDGTEVGALVGEDVECHAGVGHPADQESDDHDDDHAGD	983	1213
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	84.09	44	7	0	133	2	250	119	2e-133	97.3	ENA|BC112106|BC112106.1	206	37	38	0	86.36	-2	-1	*WGIRARSAETCRQPAWRTARAQPGSTAGTRRGCAPHPSCWRRA	*WGSPARSAETCRRPAWRTAMAQPGSTVGTRRGCVPHPSHWRRA	983	1213
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	73.08	78	21	0	433	200	550	317	6e-102	 145	ENA|BC112106|BC112106.1	311	57	61	0	78.21	-2	-1	*WRGSPQSGSCSWACRPPHTARWPGPPGTTGRFHRPAWQRSPPDCIPWPRRRSIHA*RYTGWW*SLQRP*RGWPQPG*	*WHGSPRSGSCSWAYTPPRTARWPGPPRTTGQFHRPGWQRSPPNCILWARRRSIHAERCRGCW*SHLGP*RGQPRLG*	983	1213
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	70.15	67	20	0	799	599	916	716	6e-102	 106	ENA|BC112106|BC112106.1	226	47	51	0	76.12	-2	-1	*RCRMPRWRRAANRLGTQRP***PCA*LPSRPSGWWLTPAAGQPPPSL*TPAVRSRRK*QSWGW*SG	*RCRMPRWRRAPSRSGKR*P***PCG*PPSLPSVWWLTPAAGQRPP*R*RRAAHSRKR**SWGW*SG	983	1213
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	80.49	41	8	0	123	1	240	118	6e-102	84.5	ENA|BC112106|BC112106.1	178	33	37	0	90.24	-3	-2	SEHDQQKHVGSQHGELPGLSQVVLRVLEGAAHHTRLVGEGH	AQHDQQKHVGGQHGELPWLSQVVLWVLEGAAYHTRRIGEGH	983	1213
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	77.78	27	6	0	553	473	670	590	6e-102	51.9	ENA|BC112106|BC112106.1	107	21	23	0	85.19	-2	-1	*ACSSQSHTSTACPLGCTWTSRLVGVQ	*ACSSRFHTSTAGPRGCTWTSRRVGVR	983	1213
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	71.43	14	4	0	889	848	1006	965	6e-102	32.7	ENA|BC112106|BC112106.1	65	10	12	0	85.71	-2	-1	YR*PGCRWRMTWRR	YR*QGCRWRRSWQR	983	1213
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	81.82	11	2	0	958	926	1075	1043	6e-102	28.6	ENA|BC112106|BC112106.1	56	9	9	0	81.82	-2	-1	HPVGSCHSRGW	HPVGSCRSRWW	983	1213
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	84.85	33	5	0	239	337	356	454	4e-48	72.5	ENA|BC112106|BC112106.1	152	28	29	0	87.88	2	2	ASPPPCIPLCMDTSSSGPRDAIWRASLPRWAVK	ASPAPSTPLCMDTSSSGPQDAIWRASLPPWAVK	983	1213
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	86.67	30	4	0	2	91	119	208	4e-48	71.2	ENA|BC112106|BC112106.1	149	26	26	0	86.67	2	2	CPSPTRRVWCAAPSSTRSTTWLSPGSSPCW	CPSPMRRVWYAAPSSTHSTTWLSHGSSPCW	983	1213
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	81.40	43	8	0	404	532	521	649	4e-48	47.3	ENA|BC112106|BC112106.1	97	35	35	0	81.40	2	2	TSALGRTTPSWASPSRGSWHWPALHPH*PAGPGTSQRACSARV	TSASGRTMPSWALPSPGSWRWPAPHPHSPAGPGTSPRACSARV	983	1213
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	88.89	18	2	0	764	817	881	934	4e-48	44.6	ENA|BC112106|BC112106.1	91	16	16	0	88.89	2	2	GCPTPAWHSTSLPTRALT	GCPTPAWHSTSSPTRAPT	983	1213
+gi|283855822|gb|GQ290312.1|	ENA|BC112106|BC112106.1	87.50	8	1	0	935	958	1052	1075	4e-48	21.7	ENA|BC112106|BC112106.1	41	7	7	0	87.50	2	2	AVARTHWV	AAARTHWV	983	1213
+gi|18148870|dbj|AB062417.1|	ENA|BC112106|BC112106.1	93.91	230	14	0	1	690	88	777	0.0	 538	ENA|BC112106|BC112106.1	1169	216	221	0	96.09	1	1	MNGTEGPNFYVPFSNKTGVVRSPFEAPQYYLAEPWQFSMLAAYMFLLIMLGFPINFLTLYVTVQHKKLRTPLNYILLNLAVADLFMVFGGFTTTLYTSLHGYFVFGPTGCNLEGFFATLGGEIALWSLVVLAIERYVVVCKPMSNFRFGENHAIMGVAFTWVMALACAAPPLVGWSRYIPEGMQCSCGIDYYTPHEETNNESFVIYMFVVHFIIPLIVIFFCYGQLVFTV	MNGTEGPNFYVPFSNATGVVRSPFEYPQYYLAEPWQFSMLAAYMFLLIVLGFPINFLTLYVTVQHKKLRTPLNYILLNLAVADLFMVLGGFTSTLYTSLHGYFVFGPTGCNLEGFFATLGGEIALWSLVVLAIERYVVVCKPMSNFRFGENHAIMGVAFTWVMALACAAPPLAGWSRYIPEGLQCSCGIDYYTLKPEVNNESFVIYMFVVHFTIPMIIIFFCYGQLVFTV	1047	1213
+gi|18148870|dbj|AB062417.1|	ENA|BC112106|BC112106.1	91.18	102	9	0	742	1047	829	1134	0.0	 233	ENA|BC112106|BC112106.1	503	93	100	0	98.04	1	1	KEVTRMVIIMVIAFLICWLPYAGVAFYIFTHQGSDFGPIFMTIPAFFAKTSAVYNPVIYIMMNKQFRNCMVTTLCCGKNPLGDDEASTTVSKTETSQVAPA*	KEVTRMVIIMVIAFLICWVPYASVAFYIFTHQGSNFGPIFMTIPAFFAKSAAIYNPVIYIMMNKQFRNCMLTTICCGKNPLGDDEASATVSKTETSQVAPA*	1047	1213
+gi|18148870|dbj|AB062417.1|	ENA|BC112106|BC112106.1	88.83	188	21	0	566	3	653	90	0.0	 394	ENA|BC112106|BC112106.1	854	167	170	0	90.43	-2	-3	IPHEHCIPSGMYLDQPTRGGAAQARAMTQVNATPMMAWFSPKRKLLMGLHTTTYRSMASTTKDHSAISPPKVAKKPSRLQPVGPNTKYPCREVYRVVVKPPNTMKRSATARLSRM*LRGVLSFLCWTVTYSVRKLMGKPSMISRNM*AASMLNCHGSAR*YCGASKGLRTTPVLLEKGT*KFGPSVPF	IPHEHCRPSGMYLDQPASGGAAQASAMTQVKATPMMAWFSPKRKLLMGLHTTTYRSMARTTKDHRAISPPRVAKKPSKLHPVGPKTKYPCREV*RVLVKPPRTMKRSATARLSRM*LRGVRSFLCWTVT*SVRKLMGKPSTISRNM*AASMENCHGSAR*YCGYSKGLRTTPVALEKGT*KLGPSVPF	1047	1213
+gi|18148870|dbj|AB062417.1|	ENA|BC112106|BC112106.1	84.06	138	22	0	1046	633	1133	720	0.0	 260	ENA|BC112106|BC112106.1	562	116	124	0	89.86	-2	-3	*AGATWLVSVLDTVVEASSSPRGFLPQQRVVTMQFRNCLFIMM*MTGL*TADVLAKKAGMVMKMGPKSDP*WVKM*NATPA*GSQHMRKAMTMMMTIRVTSFSAF*VVADSCCWAAASLTVKTSCP*QKNMTIRGMMK	*AGATWLVSVLDTVAEASSSPSGFLPQQMVVSMQFRNCLFIMI*MTGL*MAALLAKNAGMVMKMGPKLEPWWVKM*NATLA*GTQQIRKAMTMMMTMRVTSFSAFCVVADSCCWAAASLTVKTSCP*QKKMIIMGMVK	1047	1213
+gi|18148870|dbj|AB062417.1|	ENA|BC112106|BC112106.1	67.11	228	75	0	684	1	771	88	7e-132	 333	ENA|BC112106|BC112106.1	725	153	176	0	77.19	-1	-2	EDQLPVAEEYDNQGDDEVYNEHVDDERLIVGLLVGRVVIDPARALHPLRDVSRPADEGGGRTGQSHDPGECDAHDGVVLPEAEVAHGLAHHYVPLDG*HHQGPQCNFTA*GGKEALEVAARRTEHEVPVQRGIQGGGEAPEHHEEICHGQVEQDVVEGGPQLLVLDCDVQREEVDGEAKHDQQEHVGGQHAELPWLSQVVLRRLEGAAYDARLVGERHVEVWALGPVH	EDELPIAEKDDNHGDGEVDHEHVDDKRLVVDLRLERVVVDSTRALQALGDVPGPAGEWGCGAGQRHDPGEGNAHDGMVLPEAEVAHGLTHHHVPLDGQDHQGPQGNFTAQGGKEALQIASCGPEDEVSMQRGVEGAGEAT*DHEEVSHG*VEQDVVERRAQLLVLDGDVEREEVDGEAQHDQQKHVGGQHGELPWLSQVVLWVLEGAAYHTRRIGEGHVEVRAFCAIH	1047	1213
+gi|18148870|dbj|AB062417.1|	ENA|BC112106|BC112106.1	67.27	110	36	0	1045	716	1132	803	7e-132	 141	ENA|BC112106|BC112106.1	302	74	77	0	70.00	-3	-1	RQGPLGWSLSWTRSWRPRRHLGGSCHSREW*PCSSGTACSS*CR*RGYRRQTSWQRKPGWS*RWGQSRIPDG*RCRTPHLRRAASI*GKR*P***PYA*PPSRPSEWWPT	RPGPPGSSPSWTR*QRPHRHPVGSCRSRWW*ACSSGTACSS*YR*QGCRWRRSWQRTLGWS*RWDRSWSPGG*RCRMPRWRRAPSRSGKR*P***PCG*PPSLPSVWWLT	1047	1213
+gi|18148870|dbj|AB062417.1|	ENA|BC112106|BC112106.1	70.20	151	45	0	3	455	90	542	1e-128	 236	ENA|BC112106|BC112106.1	509	106	113	0	74.83	3	3	ERDRGPKLLRAFLQQDGRRTQPLRGAAVLPG*AMAVQHAGRLHVPADHAWLPHQLPHAVRHSPAQEAEDPPQLHPAQPGRGRSLHGVRGLHHHPVYLSARVLRVRSDGLQPRGLLCHLRR*NCTVVLGGASHRAVRSGVQAHEQLPLRGEP	EWHRRP*LLRALLQCDGCGTQPLRVPTVLPG*AMAVLHAGRLHVSADRAGLPHQLPHALRHRPAQEAAHASQLHPAQPSRG*PLHGPRWLHQHPLHLSAWILRLRAHRMQFGGLLCHPGR*NCPVVLGGPGHRAVRGGV*AHEQLPLRGEP	1047	1213
+gi|18148870|dbj|AB062417.1|	ENA|BC112106|BC112106.1	64.04	89	32	0	780	1046	867	1133	1e-128	 136	ENA|BC112106|BC112106.1	292	57	65	0	73.03	3	3	FPHMLAALRRCGVLHLHPSGIRLWPHLHDHPGFLCQDVCRL*PRHLHHDEQAVPELHGHHSLLWQEPPR*RRGLHDRVQDRDQPSGPCL	FPDLLGALRQRGILHLHPPGLQLRSHLHDHPSVLCQERRHLQPCHLYHDEQAVPELHAHHHLLRQEPTG*R*GLCYRVQDGDEPGGPGL	1047	1213
+gi|18148870|dbj|AB062417.1|	ENA|BC112106|BC112106.1	66.22	74	25	0	510	731	597	818	1e-128	 111	ENA|BC112106|BC112106.1	236	49	53	0	71.62	3	3	PPRRLV*IHPGGDAVLVRDRLLHAPRGDQQ*VVRHLHVRCTLHHPPDCHILLLRAAGLHRQGGCSPAAGVGHHS	PTRRLVQVHPRGPAVLVWNRLLHAQAGGQQRVFCHLHVRGPLHHPHDYHLFLLWAARLHRQGGRCPAAGVSHHT	1047	1213
+gi|18148870|dbj|AB062417.1|	ENA|BC112106|BC112106.1	66.04	106	36	0	242	559	329	646	2e-58	 161	ENA|BC112106|BC112106.1	347	70	71	0	66.98	2	2	WQISSWCSGASPPPCIPLCTGTSCSVRRAATSRASLPP*AVKLHCGPWWC*PSSGT*WCASP*ATSASGRTTPSWASHSPGSWLWPVRPPPSSAGLDTSRRGCSAR	WLTSSWS*VASPAPSTPLCMDTSSSGPQDAIWRASLPPWAVKLPCGPWWSWPSSGTWWCVSP*ATSASGRTMPSWALPSPGSWRWPAPHPHSPAGPGTSPRACSAR	1047	1213
+gi|18148870|dbj|AB062417.1|	ENA|BC112106|BC112106.1	85.71	21	3	0	92	154	179	241	2e-58	53.8	ENA|BC112106|BC112106.1	111	18	18	0	85.71	2	2	WLSHGSSACWPPTCSC*SCLA	WLSHGSSPCWPPTCFC*SCWA	1047	1213
+gi|18148870|dbj|AB062417.1|	ENA|BC112106|BC112106.1	73.68	19	5	0	791	847	878	934	2e-58	39.1	ENA|BC112106|BC112106.1	79	14	14	0	73.68	2	2	AGCPTQVWRSTSSPIRDPT	AGCPTPAWHSTSSPTRAPT	1047	1213
+gi|18148870|dbj|AB062417.1|	ENA|BC112106|BC112106.1	61.29	62	24	0	424	239	511	326	4e-55	81.3	ENA|BC112106|BC112106.1	171	38	38	0	61.29	-3	-1	ACTPLRTARWLAPPRTTVQFHRLRWQRSPRGCSPSDRTRSTRAERYTGWW*SPRTP*RDLPR	AYTPPRTARWPGPPRTTGQFHRPGWQRSPPNCILWARRRSIHAERCRGCW*SHLGP*RGQPR	1047	1213
+gi|18148870|dbj|AB062417.1|	ENA|BC112106|BC112106.1	74.36	39	10	0	133	17	220	104	4e-55	69.8	ENA|BC112106|BC112106.1	146	29	30	0	76.92	-3	-1	TCRRPAC*TAMAQPGSTAAPRRGCVRRPSCWRKARRSLG	TCRRPAWRTAMAQPGSTVGTRRGCVPHPSHWRRARRS*G	1047	1213
+gi|18148870|dbj|AB062417.1|	ENA|BC112106|BC112106.1	65.71	35	12	0	882	778	969	865	4e-55	56.3	ENA|BC112106|BC112106.1	117	23	27	0	77.14	-1	-2	KESRDGHEDGAKVGSLMGEDVERHTCVGQPAYEES	KERWDGHEDGTEVGALVGEDVECHAGVGHPADQES	1047	1213
+gi|18148870|dbj|AB062417.1|	ENA|BC112106|BC112106.1	58.14	43	18	0	649	521	736	608	4e-55	50.6	ENA|BC112106|BC112106.1	104	25	27	0	62.79	-3	-1	SGG**SVQRTCR*RTTHCWSPRGACSNRSRTSTASPPGCI*TS	SWGW*SGPRTCR*QKTRC*PPA*ACSSRFHTSTAGPRGCTWTS	1047	1213
+gi|18148870|dbj|AB062417.1|	ENA|BC112106|BC112106.1	66.67	12	4	0	972	937	1059	1024	4e-55	23.9	ENA|BC112106|BC112106.1	46	8	10	0	83.33	-1	-2	ATAESGDHAVPE	AAADGGEHAVPE	1047	1213
+gi|12583664|dbj|AB043817.1|	ENA|BC112106|BC112106.1	82.13	235	42	0	11	715	76	780	0.0	 498	ENA|BC112106|BC112106.1	1081	193	212	0	90.21	2	1	RTATMNGTEGPNFYIPMSNATGVVRSPFEYPQYYLAEPWAFSALSAYMFFLIIAGFPINFLTLYVTIEHKKLRTPLNYILLNLAVADLFMVFGGFTTTMYTSMHGYFVFGPTGCNIEGFFATLGGEIALWCLVVLAIERWMVVCKPVTNFRFGESHAIMGVMVTWTMALACALPPLFGWSRYIPEGLQCSCGIDYYTRAPGINNESFVIYMFTCHFSIPLAVISFCYGRLVCTVK	RATAMNGTEGPNFYVPFSNATGVVRSPFEYPQYYLAEPWQFSMLAAYMFLLIVLGFPINFLTLYVTVQHKKLRTPLNYILLNLAVADLFMVLGGFTSTLYTSLHGYFVFGPTGCNLEGFFATLGGEIALWSLVVLAIERYVVVCKPMSNFRFGENHAIMGVAFTWVMALACAAPPLAGWSRYIPEGLQCSCGIDYYTLKPEVNNESFVIYMFVVHFTIPMIIIFFCYGQLVFTVK	1344	1213
+gi|12583664|dbj|AB043817.1|	ENA|BC112106|BC112106.1	78.31	83	18	0	770	1018	835	1083	0.0	 177	ENA|BC112106|BC112106.1	382	65	79	0	95.18	2	1	VTRMVVIMVISFLVCWVPYASVAWYIFTHQGSTFGPIFMTIPSFFAKSSALYNPMIYICMNKQFRHCMITTLCCGKNPFEEED	VTRMVIIMVIAFLICWVPYASVAFYIFTHQGSNFGPIFMTIPAFFAKSAAIYNPVIYIMMNKQFRNCMLTTICCGKNPLGDDE	1344	1213
+gi|12583664|dbj|AB043817.1|	ENA|BC112106|BC112106.1	72.29	332	92	0	1017	22	1082	87	1e-150	 516	ENA|BC112106|BC112106.1	1120	240	256	0	77.11	-1	-3	SSSSKGFFPQQRVVIMQWRNCLFMQM*IMGL*RAELLAKKDGMVMKMGPKVLPWWVKMYQATLA*GTQQTRKEMTMMTTMRVTSRSALWVVSDSCCWAAASLTVHTSRP*QKEMTASGMEKWQVNM*ITKDSLLIPGARV**SIPHEHCRPSGMYRDQPKRGGRAQANAMVQVTMTPMMAWLSPKRKLVTGLQTTIHLSMARTTRHQRAISPPRVAKNPSMLQPVGPKTK*PCMEVYIVVVNPPNTMKRSATARFSRM*FSGVLSFLCSMVTYRVRKLMGNPAIIRKNM*ADRAEKAHGSAR*YCGYSNGLLTTPVAFDMGM*KLGPSVPFM	SSSPSGFLPQQMVVSMQFRNCLFIMI*MTGL*MAALLAKNAGMVMKMGPKLEPWWVKM*NATLA*GTQQIRKAMTMMMTMRVTSFSAFCVVADSCCWAAASLTVKTSCP*QKKMIIMGMVKWTTNM*MTKDSLLTSGLSV**SIPHEHCRPSGMYLDQPASGGAAQASAMTQVKATPMMAWFSPKRKLLMGLHTTTYRSMARTTKDHRAISPPRVAKKPSKLHPVGPKTKYPCREV*RVLVKPPRTMKRSATARLSRM*LRGVRSFLCWTVT*SVRKLMGKPSTISRNM*AASMENCHGSAR*YCGYSKGLRTTPVALEKGT*KLGPSVPFM	1344	1213
+gi|12583664|dbj|AB043817.1|	ENA|BC112106|BC112106.1	48.30	147	76	0	712	272	777	337	2e-98	 169	ENA|BC112106|BC112106.1	363	71	92	0	62.59	-3	-2	DGAHQSAVAERDDGEWDGEVAGKHVDHKGLIVDPRRPGIVVDPARALQTFRNVPRPAEERGEGTGQCHGPGHHDTHDGMALTEAEIGHWLADDHPPFNGQDNEAPEGNLAAEGGEEPFDVAAGGAEDKVAVHGGVHRGGESAKHHEE	DGEDELPIAEKDDNHGDGEVDHEHVDDKRLVVDLRLERVVVDSTRALQALGDVPGPAGEWGCGAGQRHDPGEGNAHDGMVLPEAEVAHGLTHHHVPLDGQDHQGPQGNFTAQGGKEALQIASCGPEDEVSMQRGVEGAGEAT*DHEE	1344	1213
+gi|12583664|dbj|AB043817.1|	ENA|BC112106|BC112106.1	54.17	72	33	0	1030	815	1095	880	2e-98	 103	ENA|BC112106|BC112106.1	220	39	53	0	73.61	-3	-2	GGRSILLLEGVLPTAEGGDHAMAELLVHADVDHGVVEG*TLGKEGWNGHEDGPKSASLVGKDVPGHTGIGHP	GSRGLIVTQWVLAAADGGEHAVPELLVHHDIDDRVVDGGALGKERWDGHEDGTEVGALVGEDVECHAGVGHP	1344	1213
+gi|12583664|dbj|AB043817.1|	ENA|BC112106|BC112106.1	47.83	69	36	0	220	14	285	79	2e-98	83.5	ENA|BC112106|BC112106.1	176	33	47	0	68.12	-3	-2	LMFDGDIQGEEVDGESGDNQEEHVGRQS*ESPWFCKVVLRVFKWTPHYTSGV*HGDVEVRSLSAVHGCG	LVLDGDVEREEVDGEAQHDQQKHVGGQHGELPWLSQVVLWVLEGAAYHTRRIGEGHVEVRAFCAIHGCG	1344	1213
+gi|12583664|dbj|AB043817.1|	ENA|BC112106|BC112106.1	72.00	25	7	0	782	708	847	773	2e-98	45.1	ENA|BC112106|BC112106.1	92	18	19	0	76.00	-2	-1	PCG*PPAQPSGWSRTPVAGQRPL*R	PCG*PPSLPSVWWLTPAAGQRPP*R	1344	1213
+gi|12583664|dbj|AB043817.1|	ENA|BC112106|BC112106.1	56.00	75	33	0	532	756	597	821	5e-65	87.7	ENA|BC112106|BC112106.1	185	42	46	0	61.33	1	3	PSLRLVSVHSGRSAVLVRDRLLYPGAWDQQ*VLCDLHVYLPLLHPTRRHLFLLRPTGVHRQRGRCPATGVRDHPE	PTRRLVQVHPRGPAVLVWNRLLHAQAGGQQRVFCHLHVRGPLHHPHDYHLFLLWAARLHRQGGRCPAAGVSHHTE	1344	1213
+gi|12583664|dbj|AB043817.1|	ENA|BC112106|BC112106.1	39.42	104	63	0	19	330	84	395	5e-65	86.8	ENA|BC112106|BC112106.1	183	41	56	0	53.85	1	3	NHERH*GT*LLHPHVKRHWCSEESI*IPAVLPCRTMGFLSSVCLHVLPDYRRIPHQLPHPVCHHRT*ETEDPTELHSAEPGRGRPLHGVWRIHHHDVHLHARLL	SHEWHRRP*LLRALLQCDGCGTQPLRVPTVLPG*AMAVLHAGRLHVSADRAGLPHQLPHALRHRPAQEAAHASQLHPAQPSRG*PLHGPRWLHQHPLHLSAWIL	1344	1213
+gi|12583664|dbj|AB043817.1|	ENA|BC112106|BC112106.1	56.90	58	25	0	829	1002	894	1067	5e-65	81.3	ENA|BC112106|BC112106.1	171	33	40	0	68.97	1	3	QCGLVHLYPPGKHFWAHLHDHSILLCQEFSPLQPHDLHLHEQAVPPLHDHHPLLWEEP	QRGILHLHPPGLQLRSHLHDHPSVLCQERRHLQPCHLYHDEQAVPELHAHHHLLRQEP	1344	1213
+gi|12583664|dbj|AB043817.1|	ENA|BC112106|BC112106.1	50.00	30	15	0	388	477	453	542	5e-65	33.6	ENA|BC112106|BC112106.1	67	15	19	0	63.33	1	3	DCPLVPRCPGH*KVDGRLQASDQFPLR*EP	NCPVVLGGPGHRAVRGGV*AHEQLPLRGEP	1344	1213
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/three_human_mRNA.fasta	Mon Apr 29 09:19:50 2013 -0400
@@ -0,0 +1,183 @@
+>ENA|AB011145|AB011145.1 Homo sapiens mRNA for KIAA0573 protein, partial cds.
+GAGAGGACGAGGTGCCGCTGCCTGGAGAATCCTCCGCTGCCGTCGGCTCCCGGAGCCCAG
+CCCTTTCCTAACCCAACCCAACCTAGCCCAGTCCCAGCCGCCAGCGCCTGTCCCTGTCAC
+GGACCCCAGCGTTACCATGCATCCTGCCGTCTTCCTATCCTTACCCGACCTCAGATGCTC
+CCTTCTGCTCCTGGTAACTTGGGTTTTTACTCCTGTAACAACTGAAATAACAAGTCTTGA
+TACAGAGAATATAGATGAAATTTTAAACAATGCTGATGTTGCTTTAGTAAATTTTTATGC
+TGACTGGTGTCGTTTCAGTCAGATGTTGCATCCAATTTTTGAGGAAGCTTCCGATGTCAT
+TAAGGAAGAATTTCCAAATGAAAATCAAGTAGTGTTTGCCAGAGTTGATTGTGATCAGCA
+CTCTGACATAGCCCAGAGATACAGGATAAGCAAATACCCAACCCTCAAATTGTTTCGTAA
+TGGGATGATGATGAAGAGAGAATACAGGGGTCAGCGATCAGTGAAAGCATTGGCAGATTA
+CATCAGGCAACAAAAAAGTGACCCCATTCAAGAAATTCGGGACTTAGCAGAAATCACCAC
+TCTTGATCGCAGCAAAAGAAATATCATTGGATATTTTGAGCAAAAGGACTCGGACAACTA
+TAGAGTTTTTGAACGAGTAGCGAATATTTTGCATGATGACTGTGCCTTTCTTTCTGCATT
+TGGGGATGTTTCAAAACCGGAAAGATATAGTGGCGACAACATAATCTACAAACCACCAGG
+GCATTCTGCTCCGGATATGGTGTACTTGGGAGCTATGACAAATTTTGATGTGACTTACAA
+TTGGATTCAAGATAAATGTGTTCCTCTTGTCCGAGAAATAACATTTGAAAATGGAGAGGA
+ATTGACAGAAGAAGGACTGCCTTTTCTCATACTCTTTCACATGAAAGAAGATACAGAAAG
+TTTAGAAATATTCCAGAATGAAGTAGCTCGGCAATTAATAAGTGAAAAAGGTACAATAAA
+CTTTTTACATGCCGATTGTGACAAATTTAGACATCCTCTTCTGCACATACAGAAAACTCC
+AGCAGATTGTCCTGTAATCGCTATTGACAGCTTTAGGCATATGTATGTGTTTGGAGACTT
+CAAAGATGTATTAATTCCTGGAAAACTCAAGCAATTCGTATTTGACTTACATTCTGGAAA
+ACTGCACAGAGAATTCCATCATGGACCTGACCCAACTGATACAGCCCCAGGAGAGCAAGC
+CCAAGATGTAGCAAGCAGTCCACCTGAGAGCTCCTTCCAGAAACTAGCACCCAGTGAATA
+TAGGTATACTCTATTGAGGGATCGAGATGAGCTTTAAAAACTTGAAAAACAGTTTGTAAG
+CCTTTCAACAGCAGCATCAACCTACGTGGTGGAAATAGTAAACCTATATTTTCATAATTC
+TATGTGTATTTTTATTTTGAATAAACAGAAAGAAATTTTGGGTTTTTAATTTTTTTCTCC
+CCGACTCAAAATGCATTGTCATTTAATATAGTAGCCTCTTAAAAAAAAAAAAACCTGCTA
+GGATTTAAAAATAAAAATCAGAGGCCTATCTCCACTTTAAATCTGTCCTGTAAAAGTTTT
+ATAAATCAAATGAAAGGTGACATTGCCAGAAACTTACCATTAACTTGCACTACTAGGGTA
+GGGAGGACTTAGGATGTTTCCTGTGTCGTATGTGCTTTTCTTTCTTTCATATGATCAATT
+CTGTTGGTATTTTCAGTATCTCATTTCTCAAAGCTAAAGAGATATACATTCTGGATACTT
+GGGAGGGGAATAAATTAAAGTTTTCACACTGTGTACTGTGTTTTACTGATTGGTTGGATA
+TTGCTTATGAAAATTCCATAGTGGTATTTTTTTGGATTCTTAATGTGTAACTTAAACATA
+CTTTGAAGTGGAGGAGAGTCATAAGACAGAACATTTGGCAGGAATTGTCCTTATGAAACA
+AGAAAAAGAAAATGAAAAGTATTATTAAGCTTCTGTGTTTGTCTAAAAATGTGGCATATG
+GATGGCATTTAAAACTTTGAATGAATTATACCTAAATCTGGGACAGGGAGGTGACAGTGG
+AACAGGCTACCAATCAGAACTAGATGACTTTTAAGGCTCCTCCTATTATGAGACTTCAAT
+TTCCAAAGAGAAGAACTAGCAGAGAAATTGTATTTCAGTAATTTTAAGCTCCTTCTGTCT
+TGTAGAGTCTTGTTATAGTTGTATAAATCAAAAACACAGAATAAGGAACATATTTAACTT
+TTTTTCATTATAAAATGGTTAGAGGACCCTACCCCCTCTAGATTCCCTGATTTCCCCAGG
+CCTGCAGCATACAGTAAGATGGGTCCCTGTGCCAGGCCTCAATACTGCCAGGGAATAAAA
+CCAGAGGGAGAGGACCCTCAGTGTCATATCAGGAAGCCCAGTGCCAGAGGACAGACAGGT
+TCAAAACTGGCTTTTCCTCTGGGCCTGGGTTGGTGCTATAGGCCAAGGGTCATTTTATAC
+TTGGGTATAAATCAATCCCAGTTTGGGAAAAGATTATTTTTAAGCTTAAAAGGCTGACAT
+GTGCCATTATATGTAGTATGTAATATATGTAACATCTTCCAATTCTTTTAAAATAAAATT
+AATATTTATAATGGATATTTAATGATTGTTATTTTTAAAAACCAGCTTATAATTCCTCGT
+TATGCATGATTTATCCAAAGTTTCCATAGTTTTATTCAAAATAATAAATGTTAATAAGGT
+GATAAGGGGTATATTTAATGTATTGTATCAAATTGTGAATAAGAAAGTAGGATGGAGCTT
+TCTAGAGGTTGGGCCTTAGTTCTGTTATCCTCATTGCTTTTAACCAATAAGTTAAATGAA
+GTTAGAGTTATGGTCTTCAGGTTAGATTATGGACCAGATCTGTGAGGGTCAGCATGGAAA
+TTCACATTCAACAAGGTAGCACACAGGACCAAGAGCAGCACATGCAATCAACTGGAATAA
+TATAGTAATCCTGTAACTGGGTTTGAAAAAATAATCAACAAAAGATACAATTCAAGGGTT
+AGGTTGCAGAGAGCTGGCTTGAGAGTAGTTATTATGAAAAAGGCCTCAAGGAGTACGTGT
+TCAGTATGCTCTAAGATGATAAAGTGGCTGTTAAAAAGGGAGTTGATTTGAGGAAGTATT
+ACTTAGCATTCATGCATATTGGGCTTAGGCTCTAGCCCTGCCACTATCATTGTCTTCTCT
+GGACTGTGAAGTCACTGAGGACAAGGAAACTAAATTTAATGTCTGTATCACTAGTGCCTA
+GAATTTCTGGACACTTAGTAGTCACCATCAGGCGTTTATTTAATGAATGAGAAGCAAAGT
+GACCTTGGTTACTTTTTTACCCTGAGGGGCTCAGCACTCATTAGGACTTGGTGCCTAATT
+TTATAAAAAGTCACTAAGCTCAAGTGCTTGGATGAAAGGACAGCGTGGATAAAAAGGTTT
+TTAAAACATGGATGTTAAGGCTGTTTTGCTTGGAGAAGACTTGGGACTGGGACAGTCTTT
+AGATATTATTTGAAATGCTGGCACTGTCTATCTGGATCCCAGGGCTTGAACTAGGATTTG
+AGGAAGTCACAGGGAAGCAGATTTCAGTCTGACATTTATTCAGTGCAAGTTTTTTGGTGC
+TGTAGTATATGATGAAAGATGTAAAGCTGAATAAAGCATTATTTCTGCCCTAGAGTTGTT
+CACAGCCTAGTCAGGCATATGGATATGTAAACAATGACTGTAACGTGTTATAGATGTAAA
+GACAAAATAAAGGTTAAAGAGGGCATAAAGGAGCACTCAATTGCAGAGATTTGAGGACAT
+TATTTTTATTTTGAGCTTTAAAAAGATGAATAGGTGTTCTCAGGAGGTAGGGATCTGGCT
+GAGAGGGAATAATCTGAGCAAAGGTATGAAACAGCCTAATGCATTAGAGAAAAAAGTTCT
+TTTAGTAAGGCATTTGGGGTTGGGGAAGCTAGAAAAAGAAATGGGAGCTGGTCACACAGG
+GCCTTGTGTGCCAGACTAAGGGGTTTGTAGTATATATTGTAGGCAGAAGAGATCCATCAA
+CAGATTGCAAGCAAGGAAGTATGTTCACTTTAAAGTTTGAGAAAGAATAGTGTGGAAGCA
+CGTCTCAAATTTAGACTTACTTGTTCCCCCTCTGAACCGTGAATCAGACCATTTCAGGTA
+GAAGTCTTCCCCGGTTTATCTGATCTACTCGGGGCCTCAGGCTTCTCAGCTGGGAAGAGA
+GGATGCAAGACCAGACTGAAGAACACGGTTGAGTCCCCAGAACCAAAAGGGGGCCTTTCT
+GCTTCTTAGCCAGCTACCTCTTCGAGTTTTTCAAATTGTGAGGGGGACCATAAAAGGATG
+GAAACTTTTAGATGACATTCTACAAATTATTTTTTTCTTTAAATTAAAAGAACCTAGCCA
+ATAAGATAGAGAATGGGCATCTAAGGCATCTCAGAGCTCTCTGATGAAGCCAGGTTGTCA
+AAGATCATTTGCAAAAGAAGGGAAAACTGGCATGACAAAAGCTACAGAGAGGAGAGTGAA
+ATATAGAAGTGTTTGAAATGTTCAAGCTCACAATAAGCTTAAATTTATAGAAAATGCTAA
+GGTTGTCAAGAAGGCTTTTTTTTTTTTCTTTTTTAAACCTGAGGGCAAAAAGGAATGGAT
+AAAGTAGTGTAATGGATTGACAATCAGGAAGAACAGAATAACTCAGTTTTTTTTTCTCCT
+ACAAGGAGATATGGCTGGACCAAAATAAAATGACATGAAATTGCAAAAATGAAAAT
+>ENA|M10051|M10051.1 Human insulin receptor mRNA, complete cds.
+GGGGGGCTGCGCGGCCGGGTCGGTGCGCACACGAGAAGGACGCGCGGCCCCCAGCGCTCT
+TGGGGGCCGCCTCGGAGCATGACCCCCGCGGGCCAGCGCCGCGCGCCTGATCCGAGGAGA
+CCCCGCGCTCCCGCAGCCATGGGCACCGGGGGCCGGCGGGGGGCGGCGGCCGCGCCGCTG
+CTGGTGGCGGTGGCCGCGCTGCTACTGGGCGCCGCGGGCCACCTGTACCCCGGAGAGGTG
+TGTCCCGGCATGGATATCCGGAACAACCTCACTAGGTTGCATGAGCTGGAGAATTGCTCT
+GTCATCGAAGGACACTTGCAGATACTCTTGATGTTCAAAACGAGGCCCGAAGATTTCCGA
+GACCTCAGTTTCCCCAAACTCATCATGATCACTGATTACTTGCTGCTCTTCCGGGTCTAT
+GGGCTCGAGAGCCTGAAGGACCTGTTCCCCAACCTCACGGTCATCCGGGGATCACGACTG
+TTCTTTAACTACGCGCTGGTCATCTTCGAGATGGTTCACCTCAAGGAACTCGGCCTCTAC
+AACCTGATGAACATCACCCGGGGTTCTGTCCGCATCGAGAAGAACAATGAGCTCTGTTAC
+TTGGCCACTATCGACTGGTCCCGTATCCTGGATTCCGTGGAGGATAATCACATCGTGTTG
+AACAAAGATGACAACGAGGAGTGTGGAGACATCTGTCCGGGTACCGCGAAGGGCAAGACC
+AACTGCCCCGCCACCGTCATCAACGGGCAGTTTGTCGAACGATGTTGGACTCATAGTCAC
+TGCCAGAAAGTTTGCCCGACCATCTGTAAGTCACACGGCTGCACCGCCGAAGGCCTCTGT
+TGCCACAGCGAGTGCCTGGGCAACTGTTCTCAGCCCGACGACCCCACCAAGTGCGTGGCC
+TGCCGCAACTTCTACCTGGACGGCAGGTGTGTGGAGACCTGCCCGCCCCCGTACTACCAC
+TTCCAGGACTGGCGCTGTGTGAACTTCAGCTTCTGCCAGGACCTGCACCACAAATGCAAG
+AACTCGCGGAGGCAGGGCTGCCACCAATACGTCATTCACAACAACAAGTGCATCCCTGAG
+TGTCCCTCCGGGTACACGATGAATTCCAGCAACTTGCTGTGCACCCCATGCCTGGGTCCC
+TGTCCCAAGGTGTGCCACCTCCTAGAAGGCGAGAAGACCATCGACTCGGTGACGTCTGCC
+CAGGAGCTCCGAGGATGCACCGTCATCAACGGGAGTCTGATCATCAACATTCGAGGAGGC
+AACAATCTGGCAGCTGAGCTAGAAGCCAACCTCGGCCTCATTGAAGAAATTTCAGGGTAT
+CTAAAAATCCGCCGATCCTACGCTCTGGTGTCACTTTCCTTCTTCCGGAAGTTACGTCTG
+ATTCGAGGAGAGACCTTGGAAATTGGGAACTACTCCTTCTATGCCTTGGACAACCAGAAC
+CTAAGGCAGCTCTGGGACTGGAGCAAACACAACCTCACCACCACTCAGGGGAAACTCTTC
+TTCCACTATAACCCCAAACTCTGCTTGTCAGAAATCCACAAGATGGAAGAAGTTTCAGGA
+ACCAAGGGGCGCCAGGAGAGAAACGACATTGCCCTGAAGACCAATGGGGACAAGGCATCC
+TGTGAAAATGAGTTACTTAAATTTTCTTACATTCGGACATCTTTTGACAAGATCTTGCTG
+AGATGGGAGCCGTACTGGCCCCCCGACTTCCGAGACCTCTTGGGGTTCATGCTGTTCTAC
+AAAGAGGCCCCTTATCAGAATGTGACGGAGTTCGATGGGCAGGATGCGTGTGGTTCCAAC
+AGTTGGACGGTGGTAGACATTGACCCACCCCTGAGGTCCAACGACCCCAAATCACAGAAC
+CACCCAGGGTGGCTGATGCGGGGTCTCAAGCCCTGGACCCAGTATGCCATCTTTGTGAAG
+ACCCTGGTCACCTTTTCGGATGAACGCCGGACCTATGGGGCCAAGAGTGACATCATTTAT
+GTCCAGACAGATGCCACCAACCCCTCTGTGCCCCTGGATCCAATCTCAGTGTCTAACTCA
+TCATCCCAGATTATTCTGAAGTGGAAACCACCCTCCGACCCCAATGGCAACATCACCCAC
+TACCTGGTTTTCTGGGAGAGGCAGGCGGAAGACAGTGAGCTGTTCGAGCTGGATTATTGC
+CTCAAAGGGCTGAAGCTGCCCTCGAGGACCTGGTCTCCACCATTCGAGTCTGAAGATTCT
+CAGAAGCACAACCAGAGTGAGTATGAGGATTCGGCCGGCGAATGCTGCTCCTGTCCAAAG
+ACAGACTCTCAGATCCTGAAGGAGCTGGAGGAGTCCTCGTTTAGGAAGACGTTTGAGGAT
+TACCTGCACAACGTGGTTTTCGTCCCCAGAAAAACCTCTTCAGGCACTGGTGCCGAGGAC
+CCTAGGCCATCTCGGAAACGCAGGTCCCTTGGCGATGTTGGGAATGTGACGGTGGCCGTG
+CCCACGGTGGCAGCTTTCCCCAACACTTCCTCGACCAGCGTGCCCACGAGTCCGGAGGAG
+CACAGGCCTTTTGAGAAGGTGGTGAACAAGGAGTCGCTGGTCATCTCCGGCTTGCGACAC
+TTCACGGGCTATCGCATCGAGCTGCAGGCTTGCAACCAGGACACCCCTGAGGAACGGTGC
+AGTGTGGCAGCCTACGTCAGTGCGAGGACCATGCCTGAAGCCAAGGCTGATGACATTGTT
+GGCCCTGTGACGCATGAAATCTTTGAGAACAACGTCGTCCACTTGATGTGGCAGGAGCCG
+AAGGAGCCCAATGGTCTGATCGTGCTGTATGAAGTGAGTTATCGGCGATATGGTGATGAG
+GAGCTGCATCTCTGCGTCTCCCGCAAGCACTTCGCTCTGGAACGGGGCTGCAGGCTGCGT
+GGGCTGTCACCGGGGAACTACAGCGTGCGAATCCGGGCCACCTCCCTTGCGGGCAACGGC
+TCTTGGACGGAACCCACCTATTTCTACGTGACAGACTATTTAGACGTCCCGTCAAATATT
+GCAAAAATTATCATCGGCCCCCTCATCTTTGTCTTTCTCTTCAGTGTTGTGATTGGAAGT
+ATTTATCTATTCCTGAGAAAGAGGCAGCCAGATGGGCCGCTGGGACCGCTTTACGCTTCT
+TCAAACCCTGAGTATCTCAGTGCCAGTGATGTGTTTCCATGCTCTGTGTACGTGCCGGAC
+GAGTGGGAGGTGTCTCGAGAGAAGATCACCCTCCTTCGAGAGCTGGGGCAGGGCTCCTTC
+GGCATGGTGTATGAGGGCAATGCCAGGGACATCATCAAGGGTGAGGCAGAGACCCGCGTG
+GCGGTGAAGACGGTCAACGAGTCAGCCAGTCTCCGAGAGCGGATTGAGTTCCTCAATGAG
+GCCTCGGTCATGAAGGGCTTCACCTGCCATCACGTGGTGCGCCTCCTGGGAGTGGTGTCC
+AAGGGCCAGCCCACGCTGGTGGTGATGGAGCTGATGGCTCACGGAGACCTGAAGAGCTAC
+CTCCGTTCTCTGCGGCCAGAGGCTGAGAATAATCCTGGCCGCCCTCCCCCTACCCTTCAA
+GAGATGATTCAGATGGCGGCAGAGATTGCTGACGGGATGGCCTACCTGAACGCCAAGAAG
+TTTGTGCATCGGGACCTGGCAGCGAGAAACTGCATGGTCGCCCATGATTTTACTGTCAAA
+ATTGGAGACTTTGGAATGACCAGAGACATCTATGAAACGGATTACTACCGGAAAGGGGGC
+AAGGGTCTGCTCCCTGTACGGTGGATGGCACCGGAGTCCCTGAAGGATGGGGTCTTCACC
+ACTTCTTCTGACATGTGGTCCTTTGGCGTGGTCCTTTGGGAAATCACCAGCTTGGCAGAA
+CAGCCTTACCAAGGCCTGTCTAATGAACAGGTGTTGAAATTTGTCATGGATGGAGGGTAT
+CTGGATCAACCCGACAACTGTCCAGAGAGAGTCACTGACCTCATGCGCATGTGCTGGCAA
+TTCAACCCCAAGATGAGGCCAACCTTCCTGGAGATTGTCAACCTGCTCAAGGACGACCTG
+CACCCCAGCTTTCCAGAGGTGTCGTTCTTCCACAGCGAGGAGAACAAGGCTCCCGAGAGT
+GAGGAGCTGGAGATGGAGTTTGAGGACATGGAGAATGTGCCCCTGGACCGTTCCTCGCAC
+TGTCAGAGGGAGGAGGCGGGGGGCCGGGATGGAGGGTCCTCGCTGGGTTTCAAGCGGAGC
+TACGAGGAACACATCCCTTACACACACATGAACGGAGGCAAGAAAAACGGGCGGATTCTG
+ACCTTGCCTCGGTCCAATCCTTCCTAACAGTGCCTACCGTGGCGGGGGCGGGCAGGGGTT
+CCCATTTTCGCTTTCCTCTGGTTTGAAAGCCTCTGGAAAACTCAGGATTCTCACGACTCT
+ACCATGTCCAGTGGAGTTCAGAGATCGTTCCTATACATTTCTGTTCATCTTAAGGTGGAC
+TCGTTTGGTTACCAATTTAACTAGTCCTGCAGAGGATTTAACTGTGAACCTGGAGGGCAA
+GGGGTTTCCACAGTTGCTGCTCCTTTGGGGCAACGACGGTTTCAAACCAGGATTTTGTGT
+TTTTTCGTTCCCCCCACCCGCCCCCAGCAGATGGAAAGAAAGCACCTGTTTTTACAAATT
+CTTTTTTTTTTTTTTTTTTTTTTTTTTTTGCTGGTGTCTGAGCTTCAGTATAAAAGACAA
+AACTTCCTGTTTGTGGAACAAAATTTCGAAAGAAAAAACCAAA
+>ENA|BC112106|BC112106.1 Homo sapiens rhodopsin, mRNA (cDNA clone MGC:138311 IMAGE:8327574), complete cds.
+CCAGCTGGAGCCCTGAGTGGCTGAGCTCAGGCCTTCGCAGCATTCTTGGGTGGGAGCAGC
+CACGGGTCAGCCACAAGGGCCACAGCCATGAATGGCACAGAAGGCCCTAACTTCTACGTG
+CCCTTCTCCAATGCGACGGGTGTGGTACGCAGCCCCTTCGAGTACCCACAGTACTACCTG
+GCTGAGCCATGGCAGTTCTCCATGCTGGCCGCCTACATGTTTCTGCTGATCGTGCTGGGC
+TTCCCCATCAACTTCCTCACGCTCTACGTCACCGTCCAGCACAAGAAGCTGCGCACGCCT
+CTCAACTACATCCTGCTCAACCTAGCCGTGGCTGACCTCTTCATGGTCCTAGGTGGCTTC
+ACCAGCACCCTCTACACCTCTCTGCATGGATACTTCGTCTTCGGGCCCACAGGATGCAAT
+TTGGAGGGCTTCTTTGCCACCCTGGGCGGTGAAATTGCCCTGTGGTCCTTGGTGGTCCTG
+GCCATCGAGCGGTACGTGGTGGTGTGTAAGCCCATGAGCAACTTCCGCTTCGGGGAGAAC
+CATGCCATCATGGGCGTTGCCTTCACCTGGGTCATGGCGCTGGCCTGCGCCGCACCCCCA
+CTCGCCGGCTGGTCCAGGTACATCCCCGAGGGCCTGCAGTGCTCGTGTGGAATCGACTAC
+TACACGCTCAAGCCGGAGGTCAACAACGAGTCTTTTGTCATCTACATGTTCGTGGTCCAC
+TTCACCATCCCCATGATTATCATCTTTTTCTGCTATGGGCAGCTCGTCTTCACCGTCAAG
+GAGGCCGCTGCCCAGCAGCAGGAGTCAGCCACCACACAGAAGGCAGAGAAGGAGGTCACC
+CGCATGGTCATCATCATGGTCATCGCTTTCCTGATCTGCTGGGTGCCCTACGCCAGCGTG
+GCATTCTACATCTTCACCCACCAGGGCTCCAACTTCGGTCCCATCTTCATGACCATCCCA
+GCGTTCTTTGCCAAGAGCGCCGCCATCTACAACCCTGTCATCTATATCATGATGAACAAG
+CAGTTCCGGAACTGCATGCTCACCACCATCTGCTGCGGCAAGAACCCACTGGGTGACGAT
+GAGGCCTCTGCTACCGTGTCCAAGACGGAGACGAGCCAGGTGGCCCCGGCCTAAGACCTG
+CCTAGGACTCTGTGGCCGACTATAGGCGTCTCCCATCCCCTACACCTTCCCCCAGCCACA
+GCCATCCCACCAG
--- a/tools/ncbi_blast_plus/ncbi_blast_plus.txt	Thu Apr 25 09:36:38 2013 -0400
+++ b/tools/ncbi_blast_plus/ncbi_blast_plus.txt	Mon Apr 29 09:19:50 2013 -0400
@@ -101,6 +101,7 @@
         - Correct case of exception regular expression (for error handling
           fall-back in case the return code is not set properly).
         - Clearer naming of output files.
+v0.0.20 - Added unit tests for BLASTN and TBLASTX.
 
 
 Developers
--- a/tools/ncbi_blast_plus/ncbi_blastn_wrapper.xml	Thu Apr 25 09:36:38 2013 -0400
+++ b/tools/ncbi_blast_plus/ncbi_blastn_wrapper.xml	Mon Apr 29 09:19:50 2013 -0400
@@ -1,4 +1,4 @@
-<tool id="ncbi_blastn_wrapper" name="NCBI BLAST+ blastn" version="0.0.19">
+<tool id="ncbi_blastn_wrapper" name="NCBI BLAST+ blastn" version="0.0.20">
     <description>Search nucleotide database with nucleotide query sequence(s)</description>
     <!-- If job splitting is enabled, break up the query file into parts -->
     <parallelism method="multi" split_inputs="query" split_mode="to_size" split_size="1000" shared_inputs="subject,histdb" merge_outputs="output1"></parallelism>
@@ -148,6 +148,18 @@
             </change_format>
         </data>
     </outputs>
+    <tests>
+        <test>
+            <param name="query" value="rhodopsin_nucs.fasta" ftype="fasta" />
+            <param name="db_opts_selector" value="file" />
+            <param name="subject" value="three_human_mRNA.fasta" ftype="fasta" />
+            <param name="database" value="" />
+            <param name="evalue_cutoff" value="1e-40" />
+            <param name="out_format" value="6" />
+            <param name="adv_opts_selector" value="basic" />
+            <output name="output1" file="blastn_rhodopsin_vs_three_human.tabular" ftype="tabular" />
+        </test>
+    </tests>
     <help>
     
 .. class:: warningmark
--- a/tools/ncbi_blast_plus/ncbi_tblastx_wrapper.xml	Thu Apr 25 09:36:38 2013 -0400
+++ b/tools/ncbi_blast_plus/ncbi_tblastx_wrapper.xml	Mon Apr 29 09:19:50 2013 -0400
@@ -1,4 +1,4 @@
-<tool id="ncbi_tblastx_wrapper" name="NCBI BLAST+ tblastx" version="0.0.19">
+<tool id="ncbi_tblastx_wrapper" name="NCBI BLAST+ tblastx" version="0.0.20">
     <description>Search translated nucleotide database with translated nucleotide query sequence(s)</description>
     <!-- If job splitting is enabled, break up the query file into parts -->
     <parallelism method="multi" split_inputs="query" split_mode="to_size" split_size="1000" shared_inputs="subject,histdb" merge_outputs="output1"></parallelism>
@@ -190,6 +190,18 @@
             </change_format>
         </data>
     </outputs>
+    <tests>
+        <test>
+            <param name="query" value="rhodopsin_nucs.fasta" ftype="fasta" />
+            <param name="db_opts_selector" value="file" />
+            <param name="subject" value="three_human_mRNA.fasta" ftype="fasta" />
+            <param name="database" value="" />
+            <param name="evalue_cutoff" value="1e-40" />
+            <param name="out_format" value="ext" />
+            <param name="adv_opts_selector" value="basic" />
+            <output name="output1" file="tblastx_rhodopsin_vs_three_human_ext.tabular" ftype="tabular" />
+        </test>
+    </tests>
     <help>
     
 .. class:: warningmark