Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Something wrong with access Bam file using "SeqAn" library sunredrock Bioinformatics 2 04-14-2013 11:59 PM
get wrong when I do de novo assembly by trinity for two 22.6G reads file taoxiang180 Bioinformatics 4 08-20-2012 09:52 PM
SAM file from gsMapper: wrong positions¿? CPCantalapiedra Bioinformatics 2 05-23-2012 05:39 AM
bed file for each read m_elena_bioinfo Bioinformatics 2 01-25-2011 07:52 AM
bfast localalign looking for the wrong reference file baldeberre Bioinformatics 16 09-17-2010 01:56 PM

Thread Tools
Old 07-20-2014, 08:40 AM   #1
Junior Member
Location: Germany

Join Date: Jul 2013
Posts: 9
Default what is wrong with this phylip-file? (Can't be read by fpars)

I'm trying to create phylogenetic trees based on gene-content data (seqeunces as discrete character seqeunces where "1"=presence of gene or homolog, "0"=absence of gene or homolog).
One of the tools I want to use is fpars. This seems to work just fine with multiple discrete character alignments in phylip format.

My test_input-data is in phylip format like in this example:
 6 50
CYS1_DIC  10100110101110011010101010111010111100101010101001
ALEU_HOR  10100110101110011010111011111010111101111010101001
CATH_HUM  10101110111111111010101010111010111110101110101111
CYS2_DIC  10011010111001101010101011001011110000101010100100
ALEU_HOA  11100100101111011010111011111110110101111010101001
CATH_HUB  10101110110011111010101010111110111110101111101111
However, I want to do bootstrapping with this data, and as stated in this thread, the fseqboot-tool just doesnt seem to handle my input data.

So I wrote my own function in python to do resampling of my alignment data and store the resampled data in a multiple alignmentfile using Bio.AlignIO.write().

The resulting file looks like this (first 14 lines as example):
 6 50
CYS1_DIC  00000111110010011011111011111010110111101010001100
ALEU_HOR  00011111110011011011111011111010110111111010001100
CATH_HUM  00000111111010011111111011111110111111101110101111
CYS2_DIC  01100111110011001100111000000110111011100101001011
ALEU_HOA  00011111110011011011111001111110110111111010010100
CATH_HUB  00000111111110111111111011111110111111101100101111
 6 50
CYS1_DIC  11010111100011010101101010110101010100000001110111
ALEU_HOR  11010111110011011111101010110101110100000001111111
CATH_HUM  11110111100011010101111111110101011100011001110111
CYS2_DIC  10100001100011010101000101110100011000000100010011
ALEU_HOA  11110110110110111111101000110111111100000001111111
CATH_HUB  11110011100011010101111111110101011101011001110111
this looks just like the results i would have expected from fseqboot (if it had handled my data correctly). The next steps would now be to run fpars on this data, and then to run fcondense on the fpars-result to get the bootstrap values for my trees.

However, when i try to run fpars on this data (command:"fpars -infile bootstrap_test4.phy -outfile TREES_BS_test4.tree -auto -seed 12345"), I get the following message:
Error: Bad discrete states file b 'bootstrap_test4.phy': read 61 states for '<null>', expected 50
Error: Unable to read discrete states from 'bootstrap_test4.phy'
Died: fpars terminated: Bad value for '-infile' with -auto defined
Can anybody see what's wron here? What can I do?

Last edited by someperson; 07-20-2014 at 08:43 AM.
someperson is offline   Reply With Quote
Old 07-21-2014, 05:15 AM   #2
Senior Member
Location: Research Triangle Park, NC

Join Date: Aug 2009
Posts: 245

Check your bootstrap file. From the error, the program is expecting 61 characters, yet you only have 50. So you must have un-intended spaces or some formatting error in your input text file - the program is misreading that input file, for some reason, and wrapping lines on you.
Michael Black, Ph.D.
ScitoVation LLC. RTP, N.C.
mbblack is offline   Reply With Quote

alignments, emboss, fpars, phylip, python

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 01:26 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO