<div dir="ltr">Thanks for your help. I still have problems with parsing. First I run the program in automated mode with $run_blast = 1 which gives me a file file-AA_with_blast_on with 384412 lines. For the same data I used blast on slurm then parsed the record for inparanoid input with $run_blast = 0. The slurm blast/parsing resulted in file-AA_blast_off.fasta file which has only 4894 records.<br><br>I got the Parser Warning during the command<br>junk after document element at line 1, column 6899908, byte 6899908 at <a href="http://blast_parser.pl">blast_parser.pl</a> line 105<br><br>This is the only command I used for parsing blastXML file<br>cat blastall_output | perl <a href="http://blast_parser.pl">blast_parser.pl</a> 40 >> file-AA_blast_off.fasta<br><br>for blastall I used the same parameters given in <a href="http://inparanoid.pl">inparanoid.pl</a>. Sir, please let me know what I am doing wrong.<br><br>Thanks<br>Sammy<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 10, 2014 at 5:43 PM, Erik Sonnhammer <span dir="ltr"><<a href="mailto:erik.sonnhammer@scilifelab.se" target="_blank">erik.sonnhammer@scilifelab.se</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
<div>Looks ok, but I don't know what that
warning means.<span class="HOEnZb"><font color="#888888"><br>
<br>
/Erik</font></span><div><div class="h5"><br>
<br>
On 11/10/2014 05:36 PM, Sammy Desilva wrote:<br>
</div></div></div><div><div class="h5">
<blockquote type="cite">
<div dir="ltr">Thanks, please have a look whether the format of
file-AA is correct for input<br>
<br>
aml_g1t1 aml_g1t1 880.6 423 423 423 423
423 423 q:1-423 h:1-423 <br>
aml_g1t1 aml_g10385t1 150.6 423 301 210 248
210 248 q:1-210 h:39-286 <br>
aml_g1t1 aml_g15430t1 98.6 423 162 139 137
139 137 q:1-139 h:3-139 <br>
aml_g1t1 aml_g10384t1 89.0 423 620 132 128
132 128 q:266-397 h:15-142 <br>
aml_g1t1 aml_g15143t1 85.5 423 246 108 107
108 107 q:301-408 h:16-122 <br>
aml_g1t1 aml_g13567t1 81.3 423 376 112 104
112 104 q:306-417 h:181-284 <br>
aml_g1t1 aml_g15245t1 67.0 423 153 95 94
95 94 q:313-407 h:33-126<br>
<br>
Parser Warning<br>
<br>
junk after document element at line 1, column 6899908, byte
6899908 at <a href="http://blast_parser.pl" target="_blank">blast_parser.pl</a> line 105<br>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Mon, Nov 10, 2014 at 1:49 PM, Erik
Sonnhammer <span dir="ltr"><<a href="mailto:erik.sonnhammer@scilifelab.se" target="_blank">erik.sonnhammer@scilifelab.se</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000"><span>
<div>On 11/10/2014 12:27 PM, Sammy Desilva wrote:<br>
</div>
<blockquote type="cite">
<div dir="ltr">
<div>
<div>As mentioned in <a href="http://inparanoid.pl" target="_blank">inparanoid.pl</a>
$run_blast could be set 0 if we have the blast
files. For the pairwise comparison of a big
dataset I want to run the blastall on cluster.
Could you tell me which parameter-set of
blastall I should use so that there will be no
difference in results (when $run_blast is 1)? </div>
</div>
</div>
</blockquote>
</span> Hi, simply use the same parameters ;)<br>
<br>
You can save compute time by running the "A-A" comparisons
only once. If you run the script for each species
comparisons, A-A and B-B will be rerun every time, so you
almost save half the compute if you just do it once for
the whole batch.<span><br>
<br>
<blockquote type="cite">
<div dir="ltr">
<div>
<div>Next, when setting $run_blast = 0 the
blastall results should already be parsed?<br>
</div>
</div>
</div>
</blockquote>
</span> Yes.<br>
<br>
/Erik<br>
<blockquote type="cite">
<div dir="ltr">
<div>
<div><br>
</div>
Thanks.<br>
</div>
Sammy<br>
</div>
<br>
<fieldset></fieldset>
<br>
<pre>_______________________________________________
InParanoid mailing list
<a href="mailto:InParanoid@lists.su.se" target="_blank">InParanoid@lists.su.se</a>
<a href="https://lists.su.se/mailman/listinfo/inparanoid-at-sbc.su.se" target="_blank">https://lists.su.se/mailman/listinfo/inparanoid-at-sbc.su.se</a>
</pre>
</blockquote>
<br>
</div>
</blockquote>
</div>
<br>
</div>
</blockquote>
<br>
</div></div></div>
</blockquote></div><br></div>