cd /home/andy/data/other_project/xjf/
#fastq convert to fasta:
cat /home/public/sequence_datas/01_zhy_genome/clean_data/cixingtou3DNA1-26_HL3FTCCXY_L1_2.clean.fq | paste - - - - | sed 's/^@/>/g'| cut -f1-2 | tr '\t' '\n' > female1.2.fasta
#make blast database
makeblastdb -dbtype nucl -in female1.2.fasta -out female1.2.fasta_db
#blastn short
blastn -query miRNA -db /home/andy/data/other_project/xjf/database/female1.2.fasta_db -out female_1.2_genome_miR275 -evalue 100 -task blastn-short
#extract target sequence
source activate python2
cd /home/andy/data/other_project/xjf/extract/
filter_fasta.py -f /home/andy/data/other_project/xjf/female1.2.fasta -o /home/andy/data/other_project/xjf/extract/01_miR275_genome -s /home/andy/data/other_project/xjf/extract/list_1
#assemble the candidate sequence using Trinity
source activate RNAseq
cd /home/andy/data/other_project/xjf/trinity_assemble/
Trinity --seqType fa --single /home/andy/data/other_project/xjf/extract/01_miR275_genome --CPU 20 --max_memory 50G --output 01_trinity_miRNA275_genome --min_contig_length 120
Trinity --seqType fa --single /home/andy/data/other_project/xjf/extract/02_miR275_genome --CPU 20 --max_memory 50G --output 02_trinity_assemble/Trinity_miRNA275_genome --min_contig_length 100
2021-11-19 短序列手动组装
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
