这篇文章主要讲解了“Leucaena trichandra线粒体基因组举例分析”，文中的讲解内容简单清晰，易于学习与理解，下面请大家跟着小编的思路慢慢深入，一起来研究和学习“Leucaena trichandra线粒体基因组举例分析”吧！
将路径改和数据替换为自己的以后运行脚本，遇到报错
[Pomgroup@localhost Pome_Mito_practice]$ bash Iternative_assembly_Pome_Mito.sh
Iternative_assembly_Pome_Mito.sh: line 2: $'\r': command not found
Iternative_assembly_Pome_Mito.sh: line 4: syntax error near unexpected token `$'\r''
'ternative_assembly_Pome_Mito.sh: line 4: `
sed -i 's/\r$//' Iternative_assembly_Pome_Mito.sh
blasr nanopore.fastq reference.fasta --nproc 16 > blasr.out
awk '{a=$8-$7;print $0,a;}' blastr.out
第8列减去第7列赋值给a并且将a添加到文件的最后一列
awk '{a=$8-$7;print $0,a;}' blastr.out | sort -n -r -k14,14
按照第14列倒叙排列
awk '{a=$8-$7;print $0,a;}' blastr.out | sort -n -r -k14,14 | awk '$14>500'
第14列大于500的行
awk '{a=$8-$7;print $0,a;}' blastr.out | sort -n -r -k14,14 | awk '$14>500' | cut -d ' ' -f1,1
以空格作为分隔符分割然后提取第一列 这样就得到了比对长度大于500的fastq的reads的id
grep -F -x -v -f
这行命令是干什么的还不知道
seqtk subseq nanopore.fasta ids.txt > aligned.fastq
canu -p hehuan -d hehuan-oxford genomeSize=2000k -nanopore-raw aligned.fastq感谢各位的阅读，以上就是“Leucaena trichandra线粒体基因组举例分析”的内容了，经过本文的学习后，相信大家对Leucaena trichandra线粒体基因组举例分析这一问题有了更深刻的体会，具体使用情况还需要大家实践验证。这里是亿速云，小编将为大家推送更多相关知识点的文章，欢迎关注！
免责声明：本站发布的内容（图片、视频和文字）以原创、转载和分享为主，文章观点不代表本网站立场，如果涉及侵权请联系站长邮箱：is@yisu.com进行举报，并提供相关证据，一经查实，将立刻删除涉嫌侵权内容。
