##选取部分组装结果
##bwa把原始数据比对到组装结果上,samtools比对排序
bwa index -p ref part.genome.fa
cat sample.list |while read a;do echo "bwa mem -t 8 -M ref R1.fastq.gz R2.fastq.gz -R \"@RG\tID:XX\tPL:ILLUMINA\tLB:$a\tSM:XX\" 2>$a.bwamem.log |samtools sort -@ 8 -T /XXX/tmp/$a - >$a.sorted.bam &";done >run.bwa.sh
##比对结果去Duplicates
awk_2 run.bwa.sh|sed 's/^>//;s/.sorted.bam//'|while read a;do echo "java -Xmx100g -jar /XXX/picard.jar MarkDuplicates TMP_DIR=/XXX/tmp MAX_FILE_HANDLES_FOR_READ_ENDS_MAP=800 INPUT=$a.sorted.bam OUTPUT=$a.sorted.dedup.bam METRICS_FILE=$a.sorted.dedup.metrics REMOVE_DUPLICATES=true ASSUME_SORTED=true &>$a.dedup.log && echo "$a.sorted.bam done" && rm $a.sorted.bam";done >run.dedup.sh
##写脚本算出part.genome.fa各条contig下特定bin（例：2000bp）下的深度，gc深度
##写脚本将part.genome.fa序列分成特定bin长度的序列subseq.fa，比对nt数据库，取tophit
##写脚本整合结果，subseq.fa的depth,GC content,及Taxon
##写脚本可视化
评估组装中的污染
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- 下面选了最近十年里，十位名人所做的毕业演讲。那么多的故事与经历，其实只想告诉你一件事： 面对迷茫和不确定的未来，我...
- 主CP:Zayn Malik/Liam Payne 副CP:Louis Tomlinson/Harry Style...
