Extract Description
$ grep -f TATATA.uniq.OsID /home/kimth/data/ggbio/databases/IRGSP/RAP2/gff_RAP2/RAP2.desc.txt
$ while read OsID; do grep $OsID /home/kimth/data/ggbio/databases/IRGSP/RAP2/gff_RAP2/RAP2.desc.txt; done < TATATA.uniq.OsID > TATATA.desc.txt
Extract Unique Sequence ID with Sequence
$ grep Forward ATATAT.txt | cut -f1,7 | sort | uniq > TATATA.2k.txt
Extract matching position
$ cut -f1,2,8 TATATA.txt
$ grep -f TATATA.uniq.OsID /home/kimth/data/ggbio/databases/IRGSP/RAP2/gff_RAP2/RAP2.desc.txt
$ while read OsID; do grep $OsID /home/kimth/data/ggbio/databases/IRGSP/RAP2/gff_RAP2/RAP2.desc.txt; done < TATATA.uniq.OsID > TATATA.desc.txt
Extract Unique Sequence ID with Sequence
$ grep Forward ATATAT.txt | cut -f1,7 | sort | uniq > TATATA.2k.txt
Extract matching position
$ cut -f1,2,8 TATATA.txt
728x90
'Bioinformatics' 카테고리의 다른 글
usage bwa using paired-end reads (0) | 2012.04.05 |
---|---|
how to make reverse complement of sequence in shell (0) | 2012.03.05 |
ABySS 사용시 주의점 (0) | 2012.02.09 |
xys 링크파일만들기 (0) | 2011.12.14 |
엑셀 : 표준편차, 상대표준편차, 표준오차 (0) | 2011.11.02 |