组装细菌基因组

一:细菌基因组下载

  • 上Genome Announcements网站(https://mra.asm.org/)找一篇细菌基因组文章;找到文章记载的SRA号**
    这里找了一篇文章:
    https://mra.asm.org/content/7/18/e01118-18
  • 点击一篇文章,查看SRR序列
  • 查看一下大小

    二:下载序列
    prefetch SRR7774037

    三:解压
    fastq-dump --gzip --split-files SRR7774037.sra

    四:fastqc 质量评价
    fastqc SRR7774037_1.fastq.gz

    五:Trimmomatic去接头

mkdir trim_out
java -jar ~/Downloads/Biosofts/Trimmomatic038/Trimmomatic-0.38/trimmomatic-0.38.jar SE 4 -phred33 SRR7774037_1.fastq.gz ./trim_out/output.fq.gz ILLUMINACLIP:/home/xiaoming/Downloads/Biosofts/Trimmomatic038/Trimmomatic-0.38/adapters/TruSeq3-SE.fa:2:30:10 SLIDINGWINDOW:4:15 LEADING:3 TRAILING:3 MINLEN:36

六:SPAdes组装基因组
七:quast评价
八:下载report.html
选择了单端测序,不会后面的基因组组装

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容