Hello,
Is fasta2DNAbin giving a different result?
Cheers
Thibaut
Sony Xperia U on O2
markharriso...@aol.com wrote:
>Hello,
>
>
>I've tried to read a fasta file containing 7 sequences of equal length (5733
>bases) and get the following output:
>
>
>
>> df<-read.dna("AT1G01040.fasta", format="fasta")
>> df
>7 DNA sequences in binary format stored in a list.
>
>
>Mean sequence length: 4224.571
> Shortest sequence: 4223
> Longest sequence: 4226
>
>
>Labels: Aa_0 Abd_0 Ag_0 Ak_1 Altai_5 Amel_1 ...
>
>
>Base composition:
> a c g t
>0.393 0.237 0.000 0.370
>>
>
>
>
>
>
>It appears that Gs are not being read which would also account for the shorter
>and uneven sequence lengths.
>Can anyone imagine what has happened here?
>
>
>Thanks for any help,
>Mark
>
>
> [[alternative HTML version deleted]]
>
>_______________________________________________
>R-sig-phylo mailing list - R-sig-phylo@r-project.org
>https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
>Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/
_______________________________________________
R-sig-phylo mailing list - R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/
