Thảo luận Thành viên:Phạm Thạch Thảo/Note: Cách đưa thêm annotation vào tree

Từ VLOS
Bước tới: chuyển hướng, tìm kiếm

Hic, em thấy phải nghĩ cách edit cái file annotation thế nào đó đơn giản hơn, edit từng trình tự em thấy cực quá.

Phạm Thạch Thảo, 11:16, 14/5/2011 (UTC)

Em nghĩ có thể dùng:

allSEQ = read.DNAStringSet("sequence.fasta") # Đọc file vào DNAStringSet object.
lenallSEQ = nchar(allSEQ); #Chiều dài tất cả các sequence có mặt.
nSEQlen = sum(lenallSEQ > 1000); #Tổng số sequence dài hơn 1000 bp.
which(lenallSEQ <= 1000); #Định vị vị trí các sequence ngắn hơn 1000.
Phạm Thạch Thảo, 10:59, 14/5/2011 (UTC)

Em tìm giúp anh cách đếm số trình tự trong 1 file fasta có chiều dài lớn hơn 1 số nhất định (ví dụ 1000bp)

Cao Xuân Hiếu, 10:40, 14/5/2011 (UTC)

Em cũng chưa edit xong hoàn toàn, em upload file ví dụ trước vậy (Với Tri thì ok!).

Phạm Thạch Thảo, 08:31, 14/5/2011 (UTC)

em upload file 16S.annot.0.txt của em lên VLOS nhé

Cao Xuân Hiếu, 08:17, 14/5/2011 (UTC)