{"id":73258,"date":"2026-11-15T15:33:47","date_gmt":"2026-11-15T07:33:47","guid":{"rendered":"https:\/\/www.biocloudservice.com\/wordpress\/?p=73258"},"modified":"2026-11-15T15:33:47","modified_gmt":"2026-11-15T07:33:47","slug":"bwa-meth-restored-sequence-nm-qc-flag-filter-practice","status":"publish","type":"post","link":"https:\/\/www.biocloudservice.com\/wordpress\/bwa-meth-restored-sequence-nm-qc-flag-filter-practice\/","title":{"rendered":"BWA-meth \u5b9e\u8dd1\u590d\u76d8\uff0c\u6bd4\u5bf9\u6210\u529f\u7684\u516d\u5341\u6761\u8bfb\u6bb5\uff0c\u4e3a\u4ec0\u4e48\u4e0b\u6e38\u53ea\u5269\u4e09\u5341\u6761"},"content":{"rendered":"<div class=\"ysx-method-guide\" style=\"max-width:820px;min-width:0;margin:0 auto;color:#26374a;line-height:1.86;overflow-wrap:anywhere\">\n<p>\u4e00\u4efd\u6bd4\u5bf9\u6587\u4ef6\u91cc\u6709\u516d\u5341\u6761\u5df2\u6bd4\u5bf9\u8bfb\u6bb5\uff0c\u6362\u4e00\u4e2a\u8fc7\u6ee4\u6761\u4ef6\u540e\u53ea\u5269\u4e09\u5341\u6761\uff0c\u5e76\u4e0d\u4e00\u5b9a\u662f\u6bd4\u5bf9\u7a81\u7136\u5931\u8d25\u3002\u672c\u6b21\u5b9e\u9645\u8fd0\u884c BWA-meth \u540e\uff0c\u4fdd\u6301\u8f93\u5165\u548c\u5750\u6807\u4e0d\u53d8\uff0c\u53ea\u589e\u52a0\u4e00\u4e2a\u7279\u5b9a\u65b9\u5411\u7684\u8d28\u91cf\u6807\u8bb0\u8bbe\u7f6e\uff0c\u5c31\u5f97\u5230\u8fd9\u6837\u7684\u7ed3\u679c\u3002\u88ab\u6807\u8bb0\u7684\u8bb0\u5f55\u4ecd\u7136\u5177\u6709\u9ad8\u6bd4\u5bf9\u8d28\u91cf\u548c\u6709\u6548\u5750\u6807\uff0c\u4f46\u4e0b\u6e38\u662f\u5426\u63a5\u7eb3\u5b83\u4eec\u53d6\u51b3\u4e8e\u8fc7\u6ee4\u89c4\u5219\u3002<\/p>\n<p>\u672c\u6587\u5b89\u88c5\u8fd0\u884c BWA-meth 0.2.10\uff0c\u540e\u7aef\u4e3a BWA 0.7.19-r1273\uff0c\u5e76\u4f7f\u7528 Samtools 1.24 \u8f6c\u6362\u3001\u6392\u5e8f\u3001\u7d22\u5f15\u548c\u8ba1\u6570\u3002\u5bb9\u5668\u5185 Python \u4e3a 3.14.7\u3002\u8f93\u5165\u662f\u81ea\u884c\u6784\u9020\u7684\u53c2\u8003\u4e0e\u53cc\u7aef\u8bfb\u6bb5\uff0c\u672c\u6587\u6ca1\u6709\u6267\u884c\u7532\u57fa\u5316\u63d0\u53d6\u6216\u771f\u5b9e\u6837\u672c\u6027\u80fd\u8bc4\u4f30\u3002\u5de5\u5177\u5b9a\u4f4d\u4e0e\u5efa\u8bae\u7684\u4e0b\u6e38\u8854\u63a5\u53ef\u67e5\u9605<a style=\"color:#7155b5;text-decoration:underline;text-underline-offset:2px\" href=\"https:\/\/github.com\/brentp\/bwa-meth\">\u5b98\u65b9\u4ed3\u5e93<\/a>\u3002<\/p>\n<figure class=\"guide-figure\" style=\"margin:24px 0 28px\"><img loading=\"lazy\" decoding=\"async\" src=\"\/wordpress\/wp-content\/uploads\/2026\/11\/bwa-meth-restored-sequence-nm-qc-flag-filter-practice-guide.png\" alt=\"\u516d\u5341\u6761\u5df2\u6bd4\u5bf9\uff0c\u4e3a\u4f55\u53ea\u4fdd\u7559\u4e09\u5341\u6761\u3002BWA-meth 0.2.10 \u5b9e\u8dd1 \u00b7 \u6bd4\u5bf9\u4e0e\u4e0b\u6e38\u8fc7\u6ee4\u300260 \u2192 30\uff0c\u5750\u6807\u672a\u53d8\uff0c\u8d28\u91cf\u6807\u8bb0\u6539\u53d8\uff0c\u9ad8\u6bd4\u5bf9\u8d28\u91cf\u4e0d\u80fd\u8986\u76d6\u5176\u4ed6\u8fc7\u6ee4\u6761\u4ef6\u3002NM = 0\uff0c\u4e0d\u7b49\u4e8e\u6062\u590d\u5e8f\u5217\u4e0e\u539f\u59cb\u53c2\u8003\u76f8\u540c\uff0c\u672c\u6b21\u6807\u7b7e\u4fdd\u7559\u8f6c\u6362\u7a7a\u95f4\u7684\u6bd4\u5bf9\u4fe1\u606f\u3002SAM \u2192 BAM\uff0c\u91cd\u5b9a\u5411\u4e0d\u4f1a\u81ea\u52a8\u8f6c\u6362\u6587\u4ef6\u683c\u5f0f\uff0c\u660e\u786e\u6392\u5e8f\u3001\u5efa\u7acb\u7d22\u5f15\u5e76\u91cd\u65b0\u8bfb\u53d6\u6838\u9a8c\u3002\" width=\"1200\" height=\"840\" style=\"display:block;width:100%;height:auto;border-radius:10px\"><figcaption style=\"font-size:14px;line-height:1.6;margin:9px 2px 0;color:#475569\">BWA-meth 0.2.10 \u5b9e\u8dd1 \u00b7 \u6bd4\u5bf9\u4e0e\u4e0b\u6e38\u8fc7\u6ee4<\/figcaption><\/figure>\n<h2 style=\"font-size:1.38em;line-height:1.5;margin:1.55em 0 .55em;color:#182d43\">\u4e00\u3001\u5148\u660e\u786e\u5b83\u4e3b\u8981\u89e3\u51b3\u7684\u662f\u6bd4\u5bf9\u95ee\u9898<\/h2>\n<p>\u4e9a\u786b\u9178\u6c22\u76d0\u6d4b\u5e8f\u8bfb\u6bb5\u4e0e\u539f\u59cb\u53c2\u8003\u4e4b\u95f4\u53ef\u80fd\u51fa\u73b0\u5927\u91cf\u7531\u5904\u7406\u8fc7\u7a0b\u5bfc\u81f4\u7684\u78b1\u57fa\u5dee\u5f02\u3002\u5982\u679c\u628a\u8fd9\u4e9b\u5dee\u5f02\u90fd\u5f53\u6210\u666e\u901a\u6d4b\u5e8f\u9519\u914d\uff0c\u5b9a\u4f4d\u95ee\u9898\u4f1a\u53d8\u5f97\u4e0d\u540c\u3002BWA-meth \u4f7f\u7528\u76f8\u5e94\u7684\u8f6c\u6362\u8868\u793a\u5b8c\u6210\u6bd4\u5bf9\uff0c\u518d\u6062\u590d\u8f93\u51fa\u4e2d\u9700\u8981\u4fdd\u7559\u7684\u8bfb\u6bb5\u4fe1\u606f\u3002\u8fd9\u4e2a\u8fc7\u7a0b\u7684\u76ee\u6807\u662f\u627e\u5230\u4f4d\u7f6e\uff0c\u4e0d\u7b49\u4e8e\u5df2\u7ecf\u7edf\u8ba1\u51fa\u6bcf\u4e2a\u80de\u5627\u5576\u7684\u7532\u57fa\u5316\u6bd4\u4f8b\u3002<\/p>\n<p>\u56e0\u6b64\uff0c\u5224\u65ad\u6d41\u7a0b\u662f\u5426\u5b8c\u6210\u81f3\u5c11\u8981\u5206\u4e24\u5c42\u3002\u7b2c\u4e00\u5c42\u68c0\u67e5\u8bfb\u6bb5\u662f\u5426\u5b9a\u4f4d\u3001\u65b9\u5411\u4e0e\u914d\u5bf9\u662f\u5426\u5408\u7406\u3001\u539f\u59cb\u6d4b\u91cf\u5e8f\u5217\u662f\u5426\u4fdd\u7559\u3002\u7b2c\u4e8c\u5c42\u624d\u662f\u6839\u636e\u5408\u9002\u7684\u89c4\u5219\u6c47\u603b\u4f4d\u70b9\u72b6\u6001\uff0c\u5e76\u8003\u8651\u8d28\u91cf\u3001\u8986\u76d6\u3001\u91cd\u53e0\u914d\u5bf9\u53ca\u91cd\u590d\u5904\u7406\u3002\u62ff\u5230\u4e00\u4e2a\u53ef\u4ee5\u6253\u5f00\u7684 BAM \u6587\u4ef6\uff0c\u53ea\u80fd\u8bf4\u660e\u83b7\u5f97\u4e86\u6bd4\u5bf9\u8f7d\u4f53\uff0c\u4e0d\u80fd\u81ea\u52a8\u8bf4\u660e\u7b2c\u4e8c\u5c42\u4e5f\u5df2\u7ecf\u5b8c\u6210\u3002<\/p>\n<p>BWA-meth \u7684\u5b98\u65b9\u8bf4\u660e\u5c06\u540e\u7eed\u4f4d\u70b9\u7edf\u8ba1\u4e0e\u5176\u4ed6\u4efb\u52a1\u4ea4\u7ed9\u76f8\u5e94\u5de5\u5177\u3002\u672c\u6587\u4e5f\u9075\u5b88\u8fd9\u4e2a\u8fb9\u754c\uff0c\u4e3b\u8981\u9a8c\u8bc1\u7d22\u5f15\u3001\u6bd4\u5bf9\u3001\u5e8f\u5217\u6062\u590d\u548c\u8d28\u91cf\u6807\u8bb0\u3002\u5bf9\u4e8e\u7532\u57fa\u5316\u6bd4\u4f8b\u3001\u8f6c\u5316\u6548\u7387\u6216\u5dee\u5f02\u533a\u57df\uff0c\u672c\u6587\u6ca1\u6709\u8f93\u51fa\u7ed3\u679c\uff0c\u4e5f\u4e0d\u4f1a\u4ece\u6bd4\u5bf9\u7387\u95f4\u63a5\u63a8\u65ad\u8fd9\u4e9b\u6307\u6807\u3002<\/p>\n<h2 style=\"font-size:1.38em;line-height:1.5;margin:1.55em 0 .55em;color:#182d43\">\u4e8c\u3001\u4e00\u6761\u53c2\u8003\uff0c\u4e3a\u4ec0\u4e48\u751f\u6210\u4e24\u4e2a\u8f6c\u6362\u7248\u672c<\/h2>\n<p>\u672c\u4f8b\u53c2\u8003\u662f\u4e00\u6761\u4e8c\u4e07\u56db\u5343\u78b1\u57fa\u7684\u865a\u6784\u67d3\u8272\u4f53\u3002\u5b9e\u9645\u7d22\u5f15\u547d\u4ee4\u5982\u4e0b\uff0c\u5b83\u4f1a\u6839\u636e\u539f\u59cb\u53c2\u8003\u521b\u5efa\u8f6c\u6362\u6587\u4ef6\u4e0e\u76f8\u5e94\u7d22\u5f15\uff0c\u4e0d\u662f\u53ea\u5bf9\u539f\u59cb\u53c2\u8003\u6267\u884c\u666e\u901a\u6bd4\u5bf9\u7d22\u5f15\u3002<\/p>\n<pre style=\"max-width:100%;box-sizing:border-box;overflow-x:auto;margin:18px 0;padding:14px 16px;background:#f3f6fa;color:#1a293c;border:1px solid #d8e1ea;border-radius:8px;line-height:1.55;white-space:pre;overflow-wrap:normal;word-break:normal\"><code class=\"language-bash\">bwameth.py index reference.fa\r\n<\/code><\/pre>\n<p>\u68c0\u67e5\u751f\u6210\u7684\u8f6c\u6362\u53c2\u8003\uff0c\u53d1\u73b0\u5305\u542b\u4e24\u6761\u5404\u4e8c\u4e07\u56db\u5343\u78b1\u57fa\u7684\u5e8f\u5217\u3002\u5176\u4e2d\u4e00\u6761\u628a\u80de\u5627\u5576\u8f6c\u6362\u4e3a\u80f8\u817a\u5627\u5576\uff0c\u53e6\u4e00\u6761\u628a\u9e1f\u560c\u5464\u8f6c\u6362\u4e3a\u817a\u560c\u5464\u3002\u811a\u672c\u4e0d\u662f\u53ea\u68c0\u67e5\u6587\u4ef6\u540d\u5b58\u5728\uff0c\u800c\u662f\u9010\u5b57\u7b26\u9a8c\u8bc1\u4e24\u6761\u5e8f\u5217\u7b49\u4e8e\u539f\u59cb\u53c2\u8003\u7684\u76f8\u5e94\u8f6c\u6362\u3002\u524d\u4e00\u6761\u80de\u5627\u5576\u6570\u91cf\u4e3a\u96f6\uff0c\u540e\u4e00\u6761\u9e1f\u560c\u5464\u6570\u91cf\u4e3a\u96f6\uff0c\u4e0e\u6784\u9020\u89c4\u5219\u4e00\u81f4\u3002<\/p>\n<div class=\"guide-table\" style=\"max-width:100%;overflow-x:auto;margin:18px 0\">\n<table style=\"display:table;width:100%;max-width:100%;min-width:0;table-layout:fixed;border-collapse:collapse;font-size:.92em;line-height:1.6\">\n<thead>\n<tr>\n<th scope=\"col\" style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;background:#f7f3ff;color:#1e293b;vertical-align:top;text-align:left;padding:6px;border:1px solid #c4cfdb\">\u6587\u4ef6\u4e2d\u7684\u5bf9\u8c61<\/th>\n<th scope=\"col\" style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;background:#f7f3ff;color:#1e293b;vertical-align:top;text-align:left;padding:6px;border:1px solid #c4cfdb\">\u672c\u4f8b\u957f\u5ea6<\/th>\n<th scope=\"col\" style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;background:#f7f3ff;color:#1e293b;vertical-align:top;text-align:left;padding:6px;border:1px solid #c4cfdb\">\u6838\u9a8c\u91cd\u70b9<\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">\u539f\u59cb\u53c2\u8003 chrToy<\/td>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">24000<\/td>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">\u4fdd\u7559\u539f\u59cb\u5e8f\u5217\u4e0e\u6821\u9a8c\u503c<\/td>\n<\/tr>\n<tr>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">\u8f6c\u6362\u53c2\u8003 fchrToy<\/td>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">24000<\/td>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">C \u8f6c\u4e3a T<\/td>\n<\/tr>\n<tr>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">\u8f6c\u6362\u53c2\u8003 rchrToy<\/td>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">24000<\/td>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">G \u8f6c\u4e3a A<\/td>\n<\/tr>\n<tr>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">\u6700\u7ec8\u6bd4\u5bf9\u53c2\u8003\u540d\u79f0<\/td>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">chrToy<\/td>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">\u56de\u5230\u539f\u59cb\u5750\u6807\u547d\u540d<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<\/div>\n<p>\u8fd9\u4e2a\u68c0\u67e5\u53ef\u4ee5\u9632\u6b62\u4e00\u79cd\u5e38\u89c1\u6df7\u6dc6\uff0c\u76ee\u5f55\u4e2d\u6709\u4e00\u5957\u7d22\u5f15\uff0c\u5c31\u5047\u8bbe\u5b83\u4e00\u5b9a\u5bf9\u5e94\u5f53\u524d\u53c2\u8003\u3002\u53c2\u8003\u5e8f\u5217\u66f4\u65b0\u3001\u6587\u4ef6\u88ab\u8986\u76d6\u6216\u8005\u4e0d\u540c\u6784\u5efa\u5171\u7528\u76f8\u4f3c\u540d\u79f0\u65f6\uff0c\u4ec5\u770b\u6269\u5c55\u540d\u65e0\u6cd5\u5224\u65ad\u3002\u5206\u6790\u8bb0\u5f55\u5e94\u7ed1\u5b9a\u539f\u59cb\u53c2\u8003\u5185\u5bb9\u3001\u8f6c\u6362\u6587\u4ef6\u548c\u540e\u7aef\u7248\u672c\u3002\u91cd\u65b0\u751f\u6210\u7d22\u5f15\u4e0d\u662f\u4e3a\u4e86\u5f62\u5f0f\u5b8c\u6574\uff0c\u800c\u662f\u4e3a\u4e86\u4fdd\u8bc1\u6240\u6709\u5750\u6807\u548c\u8f6c\u6362\u89c4\u5219\u6307\u5411\u540c\u4e00\u4efd\u5e8f\u5217\u3002<\/p>\n<p>\u672c\u6587\u6ca1\u6709\u6bd4\u8f83\u4e0d\u540c\u540e\u7aef\u7d22\u5f15\u7684\u6027\u80fd\uff0c\u4e5f\u6ca1\u6709\u6d4b\u8bd5\u66f4\u5927\u53c2\u8003\u7684\u5185\u5b58\u9700\u6c42\u3002\u4e24\u6761\u5c0f\u5e8f\u5217\u80fd\u591f\u5feb\u901f\u5b8c\u6210\u7d22\u5f15\uff0c\u53ea\u8bf4\u660e\u6f14\u793a\u6d41\u7a0b\u6210\u7acb\u3002\u771f\u5b9e\u9879\u76ee\u7684\u6210\u672c\u53d6\u51b3\u4e8e\u53c2\u8003\u5927\u5c0f\u3001\u6570\u636e\u91cf\u3001\u540e\u7aef\u9009\u62e9\u548c\u8fd0\u884c\u73af\u5883\uff0c\u4e0d\u80fd\u4ece\u8fd9\u4e2a\u5c0f\u6848\u4f8b\u5916\u63a8\u5168\u57fa\u56e0\u7ec4\u8017\u65f6\u3002<\/p>\n<h2 style=\"font-size:1.38em;line-height:1.5;margin:1.55em 0 .55em;color:#182d43\">\u4e09\u3001\u53cc\u7aef\u8bfb\u6bb5\u7684\u65b9\u5411\uff0c\u4e0d\u80fd\u53ea\u9760\u6587\u4ef6\u540d\u731c<\/h2>\n<p>\u6211\u4eec\u751f\u6210\u4e09\u5341\u5bf9\u6709\u6548\u8bfb\u6bb5\uff0c\u5176\u4e2d\u5341\u4e94\u5bf9\u6765\u81ea\u539f\u59cb\u9876\u90e8\u94fe\uff0c\u5341\u4e94\u5bf9\u6765\u81ea\u539f\u59cb\u5e95\u90e8\u94fe\u3002\u6bcf\u4e2a\u7247\u6bb5\u957f\u5ea6\u4e8c\u767e\u516d\u5341\u78b1\u57fa\uff0c\u4e24\u7aef\u5404\u6d4b\u4e00\u767e\u4e94\u5341\u78b1\u57fa\uff0c\u56e0\u6b64\u4e24\u7aef\u6709\u56db\u5341\u78b1\u57fa\u91cd\u53e0\u3002\u518d\u52a0\u4e24\u5bf9\u5168\u672a\u77e5\u78b1\u57fa\u5bf9\u7167\uff0c\u603b\u8ba1\u4e09\u5341\u4e8c\u5bf9\u3001\u516d\u5341\u56db\u6761\u8bb0\u5f55\u3002<\/p>\n<p>\u751f\u6210\u65f6\u5148\u6309\u8bbe\u5b9a\u8f6c\u6362\u7247\u6bb5\u4e2d\u7684\u80de\u5627\u5576\uff0c\u518d\u63d0\u53d6\u7b2c\u4e00\u7aef\u4e0e\u7b2c\u4e8c\u7aef\u7684\u76f8\u5e94\u5e8f\u5217\u3002\u7b2c\u4e8c\u7aef\u5305\u542b\u65b9\u5411\u8f6c\u6362\uff0c\u4e0d\u80fd\u7b80\u5355\u590d\u5236\u7b2c\u4e00\u7aef\u6216\u628a\u540c\u4e00\u4e2a\u5b57\u7b26\u4e32\u4f5c\u4e3a\u4e24\u7aef\u3002\u811a\u672c\u4e3a\u6bcf\u4e00\u5bf9\u4fdd\u5b58\u4e86\u4e24\u7aef\u7684\u539f\u59cb\u5e8f\u5217\u3001\u9884\u671f\u8d77\u70b9\u548c\u6765\u6e90\u65b9\u5411\uff0c\u4fbf\u4e8e\u6bd4\u5bf9\u540e\u9010\u6761\u6838\u67e5\u3002<\/p>\n<p>\u5b9e\u9645\u8f93\u51fa\u4e2d\u516d\u5341\u6761\u6709\u6548\u8bfb\u6bb5\u5747\u6bd4\u5bf9\uff0c\u56db\u6761\u5168\u672a\u77e5\u78b1\u57fa\u8bb0\u5f55\u672a\u6bd4\u5bf9\u3002\u516d\u5341\u6761\u90fd\u5177\u6709\u9884\u671f\u5750\u6807\u548c\u5b8c\u6574\u7684\u4e00\u767e\u4e94\u5341\u78b1\u57fa\u5339\u914d\uff0c\u57fa\u7840\u8fd0\u884c\u4e2d\u4e5f\u90fd\u5e26\u6709\u7b26\u5408\u914d\u5bf9\u5173\u7cfb\u7684\u6807\u8bb0\u3002\u8f93\u5165\u8bfb\u540d\u4f7f\u7528\u7b2c\u4e00\u7aef\u3001\u7b2c\u4e8c\u7aef\u540e\u7f00\uff0c\u8f93\u51fa\u5f52\u5e76\u5230\u5171\u540c\u7684\u7247\u6bb5\u540d\u79f0\uff0c\u914d\u5bf9\u5173\u7cfb\u4e0e\u751f\u6210\u8bb0\u5f55\u4e00\u81f4\u3002<\/p>\n<p>\u8fd9\u91cc\u7684\u6210\u529f\u4f9d\u8d56\u4e8e\u672c\u4f8b\u975e\u5e38\u5e72\u51c0\u7684\u8f93\u5165\u3002\u6ca1\u6709\u63a5\u5934\u3001\u63d2\u5165\u7f3a\u5931\u3001\u590d\u6742\u91cd\u590d\u6216\u771f\u5b9e\u6d4b\u5e8f\u8d28\u91cf\u5206\u5e03\uff0c\u4e0d\u80fd\u636e\u6b64\u5ba3\u79f0\u8f6f\u4ef6\u5728\u5b9e\u9645\u6570\u636e\u4e0a\u5177\u6709\u767e\u5206\u4e4b\u767e\u5b9a\u4f4d\u51c6\u786e\u7387\u3002\u5c0f\u6848\u4f8b\u7684\u5408\u7406\u7528\u9014\uff0c\u662f\u4e3a\u81ea\u5df1\u7684\u5b89\u88c5\u548c\u6574\u7406\u6d41\u7a0b\u5efa\u7acb\u4e00\u4e2a\u80fd\u660e\u786e\u68c0\u67e5\u5bf9\u9519\u7684\u8d77\u70b9\u3002<\/p>\n<h2 style=\"font-size:1.38em;line-height:1.5;margin:1.55em 0 .55em;color:#182d43\">\u56db\u3001\u8f93\u51fa\u4e3a\u4ec0\u4e48\u4ecd\u8981\u4fdd\u7559\u6ca1\u6709\u8f6c\u6362\u5b8c\u7684\u5e8f\u5217<\/h2>\n<p>\u6bd4\u5bf9\u65f6\u4f7f\u7528\u8f6c\u6362\u8868\u793a\uff0c\u5e76\u4e0d\u610f\u5473\u7740\u6700\u7ec8\u8f93\u51fa\u5e94\u8be5\u4e22\u6389\u539f\u59cb\u6d4b\u5e8f\u4e2d\u7684\u78b1\u57fa\u4fe1\u606f\u3002\u5982\u679c\u628a\u6240\u6709\u80de\u5627\u5576\u90fd\u6c38\u4e45\u6539\u6210\u80f8\u817a\u5627\u5576\uff0c\u540e\u7eed\u5c31\u65e0\u6cd5\u518d\u533a\u5206\u54ea\u4e9b\u4f4d\u7f6e\u5728\u539f\u59cb\u8bfb\u6bb5\u91cc\u4fdd\u7559\u4e86\u80de\u5627\u5576\u3002\u8f6c\u6362\u4fbf\u4e8e\u5b9a\u4f4d\uff0c\u539f\u59cb\u89c2\u6d4b\u5219\u662f\u540e\u7eed\u72b6\u6001\u5224\u65ad\u7684\u4f9d\u636e\uff0c\u4e24\u8005\u9700\u8981\u540c\u65f6\u6b63\u786e\u5904\u7406\u3002<\/p>\n<p>\u6211\u4eec\u9010\u6761\u6bd4\u8f83\u6700\u7ec8\u6bd4\u5bf9\u4e2d\u7684\u5e8f\u5217\u4e0e\u8f93\u5165 FASTQ\u3002\u5bf9\u4e8e\u4ee5\u6b63\u65b9\u5411\u5b58\u50a8\u7684\u8bb0\u5f55\u76f4\u63a5\u6bd4\u8f83\uff0c\u5bf9\u4e8e\u5e26\u53cd\u5411\u6807\u8bb0\u7684\u8bb0\u5f55\uff0c\u5148\u6309\u683c\u5f0f\u7ea6\u5b9a\u4f7f\u7528\u76f8\u5e94\u53cd\u5411\u4e92\u8865\u5e8f\u5217\u6bd4\u8f83\u3002\u516d\u5341\u6761\u6709\u6548\u8bb0\u5f55\u5168\u90e8\u4e00\u81f4\u3002\u8fd9\u4e00\u6b65\u5f88\u91cd\u8981\uff0c\u82e5\u4e0d\u8003\u8651\u6bd4\u5bf9\u65b9\u5411\uff0c\u76f4\u63a5\u628a\u4e24\u4e2a\u5b57\u7b26\u4e32\u9010\u5b57\u6bd4\u8f83\uff0c\u53cd\u800c\u4f1a\u628a\u6b63\u5e38\u7684\u683c\u5f0f\u8868\u793a\u8bef\u5224\u6210\u5e8f\u5217\u635f\u574f\u3002<\/p>\n<p>\u672c\u6b21\u8f93\u51fa\u8fd8\u542b\u6709\u4e0e\u8f6c\u6362\u65b9\u5411\u76f8\u5173\u7684\u6807\u7b7e\u3002\u5b83\u4eec\u5e2e\u52a9\u89e3\u91ca\u8bb0\u5f55\u5982\u4f55\u6765\u81ea\u8f6c\u6362\u7a7a\u95f4\uff0c\u4f46\u4e0d\u80fd\u66ff\u4ee3\u5b9e\u9645\u5e8f\u5217\u6838\u9a8c\u3002\u5224\u65ad\u6062\u590d\u662f\u5426\u6b63\u786e\uff0c\u6700\u76f4\u63a5\u7684\u4f9d\u636e\u4ecd\u7136\u662f\u539f\u59cb FASTQ\u3001\u65b9\u5411\u6807\u8bb0\u4ee5\u53ca\u6700\u7ec8\u5e8f\u5217\uff0c\u800c\u4e0d\u662f\u770b\u5230\u4e00\u4e2a\u719f\u6089\u7684\u6807\u7b7e\u5c31\u8ba4\u4e3a\u6d41\u7a0b\u5fc5\u7136\u6b63\u786e\u3002<\/p>\n<p>\u5b9e\u9645\u6574\u7406\u65f6\u5e94\u4fdd\u7559\u4e00\u5c0f\u90e8\u5206\u53ef\u56de\u67e5\u7684\u8bfb\u540d\uff0c\u4f7f\u6837\u672c\u5904\u7406\u524d\u540e\u7684\u8bb0\u5f55\u80fd\u591f\u8fde\u63a5\u3002\u82e5\u4e2d\u9014\u91cd\u547d\u540d\u3001\u5220\u9664\u540e\u7f00\u6216\u6df7\u5408\u6587\u4ef6\u800c\u4e0d\u7559\u6620\u5c04\uff0c\u51fa\u73b0\u914d\u5bf9\u5f02\u5e38\u540e\u4f1a\u5f88\u96be\u8ffd\u6eaf\u3002\u8bfb\u540d\u6574\u7406\u770b\u8d77\u6765\u662f\u683c\u5f0f\u5de5\u4f5c\uff0c\u5374\u53ef\u80fd\u76f4\u63a5\u5f71\u54cd\u914d\u5bf9\u8bc6\u522b\u4e0e\u540e\u7eed\u53bb\u91cd\u3002<\/p>\n<h2 style=\"font-size:1.38em;line-height:1.5;margin:1.55em 0 .55em;color:#182d43\">\u4e94\u3001\u4e3a\u4ec0\u4e48 NM \u4e3a\u96f6\uff0c\u539f\u59cb\u53c2\u8003\u5374\u6709\u51e0\u5341\u4e2a\u4e0d\u540c\u78b1\u57fa<\/h2>\n<p>\u672c\u6b21\u4e00\u4e2a\u5f88\u6709\u4ee3\u8868\u6027\u7684\u8bb0\u5f55\uff0c\u539f\u59cb\u8f93\u51fa\u5e8f\u5217\u4e0e\u672a\u8f6c\u6362\u53c2\u8003\u9010\u4f4d\u7f6e\u6bd4\u8f83\u6709\u4e09\u5341\u4e94\u5904\u5dee\u5f02\uff0c\u8bb0\u5f55\u4e2d\u7684 NM \u6807\u7b7e\u5374\u662f\u96f6\uff0cMD \u663e\u793a\u6574\u6bb5\u5339\u914d\u3002\u516d\u5341\u6761\u6709\u6548\u8bb0\u5f55\u7684 NM \u90fd\u4e3a\u96f6\uff0c\u800c\u5b83\u4eec\u4e0e\u539f\u59cb\u53c2\u8003\u7684\u9010\u5b57\u7b26\u5dee\u5f02\u8303\u56f4\u4e3a\u5341\u4e5d\u5230\u4e94\u5341\u5904\u3002<\/p>\n<p>\u8fd9\u8868\u660e\u5728\u672c\u6b21\u7248\u672c\u548c\u8fd0\u884c\u8def\u5f84\u4e2d\uff0c\u8fd9\u4e9b\u6807\u7b7e\u4fdd\u7559\u4e86\u8f6c\u6362\u7a7a\u95f4\u7684\u6bd4\u5bf9\u4fe1\u606f\uff0c\u4e0d\u80fd\u76f4\u63a5\u5f53\u4f5c\u6062\u590d\u540e\u5e8f\u5217\u76f8\u5bf9\u672a\u8f6c\u6362\u53c2\u8003\u7684\u5dee\u5f02\u603b\u6570\u3002\u4e0a\u8ff0\u5dee\u5f02\u662f\u6211\u4eec\u5728\u751f\u6210\u8bfb\u6bb5\u65f6\u660e\u786e\u52a0\u5165\u7684\u8f6c\u6362\u7ed3\u679c\uff0c\u4e0d\u662f\u5b9e\u9a8c\u53d1\u73b0\u4e86\u51e0\u5341\u4e2a\u771f\u5b9e\u53d8\u5f02\uff0c\u4e5f\u4e0d\u662f\u7a0b\u5e8f\u51ed\u7a7a\u6539\u53d8\u4e86\u5e8f\u5217\u3002<\/p>\n<p><a style=\"color:#7155b5;text-decoration:underline;text-underline-offset:2px\" href=\"https:\/\/samtools.github.io\/hts-specs\/SAMtags.pdf\">SAM \u6807\u7b7e\u89c4\u8303<\/a>\u8bf4\u660e\u4e86\u6807\u51c6\u5b57\u6bb5\u7684\u542b\u4e49\uff0c\u4f46\u5177\u4f53\u6587\u4ef6\u7ecf\u8fc7\u8f6c\u6362\u4e0e\u6062\u590d\u4ee5\u540e\uff0c\u4ecd\u9700\u68c0\u67e5\u751f\u6210\u5de5\u5177\u5982\u4f55\u7ef4\u62a4\u76f8\u5173\u6807\u7b7e\u3002BWA-meth \u7684<a style=\"color:#7155b5;text-decoration:underline;text-underline-offset:2px\" href=\"https:\/\/raw.githubusercontent.com\/brentp\/bwa-meth\/master\/bwameth.py\">\u5b9e\u73b0\u4ee3\u7801<\/a>\u4e0e\u672c\u6b21\u9010\u6761\u6838\u9a8c\u7ed3\u679c\u53ef\u5171\u540c\u5e2e\u52a9\u8ffd\u8e2a\u8fd9\u4e00\u8fc7\u7a0b\u3002\u4e0d\u80fd\u62ff\u666e\u901a\u6d4b\u5e8f\u4e2d\u719f\u6089\u7684\u89e3\u91ca\uff0c\u4e0d\u52a0\u68c0\u67e5\u5730\u5957\u5728\u7ecf\u8fc7\u7279\u6b8a\u5904\u7406\u7684\u6587\u4ef6\u4e0a\u3002<\/p>\n<p>\u5982\u679c\u540e\u7eed\u7a0b\u5e8f\u9700\u8981\u76f8\u5bf9\u4e8e\u539f\u59cb\u53c2\u8003\u91cd\u65b0\u8ba1\u7b97\u7684\u5dee\u5f02\u4fe1\u606f\uff0c\u5e94\u5148\u6838\u5bf9\u5b83\u7684\u8f93\u5165\u8981\u6c42\uff0c\u518d\u51b3\u5b9a\u662f\u5426\u6d3e\u751f\u4e00\u4efd\u91cd\u65b0\u8ba1\u7b97\u7684\u6587\u4ef6\uff0c\u5e76\u4fdd\u7559\u539f\u59cb\u4ea7\u7269\u3002\u76f4\u63a5\u8986\u76d6\u539f\u59cb\u6807\u7b7e\u4f1a\u8ba9\u540e\u6765\u8005\u96be\u4ee5\u8fd8\u539f\u4e0a\u6e38\u8bed\u4e49\u3002\u66f4\u4e0d\u80fd\u6839\u636e\u8fd9\u91cc\u7684\u5dee\u5f02\u5b57\u7b26\u4e32\u76f4\u63a5\u8fdb\u884c\u666e\u901a\u53d8\u5f02\u89e3\u91ca\uff0c\u56e0\u4e3a\u5b83\u4eec\u5305\u542b\u6d4b\u5e8f\u5904\u7406\u65b9\u5f0f\u5e26\u6765\u7684\u9884\u671f\u53d8\u5316\u3002<\/p>\n<h2 style=\"font-size:1.38em;line-height:1.5;margin:1.55em 0 .55em;color:#182d43\">\u516d\u3001\u6587\u4ef6\u6269\u5c55\u540d\u4e0d\u4f1a\u66ff\u4f60\u5b8c\u6210\u683c\u5f0f\u8f6c\u6362<\/h2>\n<p>BWA-meth \u7684\u672c\u6b21\u5b9e\u9645\u8f93\u51fa\u662f\u6807\u51c6\u8f93\u51fa\u4e2d\u7684 SAM \u6587\u672c\u3002\u628a\u91cd\u5b9a\u5411\u6587\u4ef6\u547d\u540d\u6210 BAM\uff0c\u5e76\u4e0d\u4f1a\u8ba9\u5b83\u81ea\u52a8\u53d8\u6210\u538b\u7f29\u4e8c\u8fdb\u5236\u6bd4\u5bf9\u6587\u4ef6\u3002\u6211\u4eec\u660e\u786e\u4fdd\u5b58 SAM\uff0c\u518d\u4f7f\u7528 Samtools \u6392\u5e8f\u751f\u6210 BAM\uff0c\u6700\u540e\u5efa\u7acb\u7d22\u5f15\u3002<\/p>\n<pre style=\"max-width:100%;box-sizing:border-box;overflow-x:auto;margin:18px 0;padding:14px 16px;background:#f3f6fa;color:#1a293c;border:1px solid #d8e1ea;border-radius:8px;line-height:1.55;white-space:pre;overflow-wrap:normal;word-break:normal\"><code class=\"language-bash\">bwameth.py --threads 2 --reference reference.fa \\\r\n  R1.fq R2.fq &gt; normal.sam\r\n\r\nsamtools sort -@ 1 -o normal.bam normal.sam\r\nsamtools index normal.bam\r\nsamtools flagstat -O json normal.bam\r\n<\/code><\/pre>\n<p>\u5206\u6b65\u9aa4\u4fdd\u5b58\u7684\u597d\u5904\u662f\u6bcf\u4e00\u6b65\u90fd\u53ef\u4ee5\u68c0\u67e5\u9000\u51fa\u7801\u548c\u4ea7\u7269\u3002\u5c24\u5176\u5728\u6279\u91cf\u5904\u7406\u65f6\uff0c\u4e0d\u80fd\u56e0\u4e3a\u91cd\u5b9a\u5411\u521b\u5efa\u4e86\u4e00\u4e2a\u975e\u7a7a\u6587\u4ef6\uff0c\u5c31\u628a\u8be5\u6837\u672c\u6807\u6210\u5b8c\u6210\u3002\u9700\u8981\u786e\u8ba4\u6587\u4ef6\u5934\u3001\u53c2\u8003\u540d\u79f0\u3001\u8bb0\u5f55\u6570\u4ee5\u53ca\u540e\u7eed\u8bfb\u53d6\u662f\u5426\u6210\u529f\u3002\u547d\u4ee4\u65e5\u5fd7\u4e5f\u5e94\u4e0e\u7ed3\u679c\u5206\u5f00\u4fdd\u5b58\uff0c\u907f\u514d\u9519\u8bef\u4fe1\u606f\u6df7\u8fdb\u6bd4\u5bf9\u6b63\u6587\u3002<\/p>\n<p>\u82e5\u5c06\u8fd9\u4e9b\u6b65\u9aa4\u6539\u6210\u7ba1\u9053\uff0c\u9700\u8981\u540c\u65f6\u68c0\u67e5\u4e0a\u6e38\u9000\u51fa\u72b6\u6001\uff0c\u4e0d\u80fd\u53ea\u770b\u6700\u540e\u4e00\u4e2a\u7a0b\u5e8f\u662f\u5426\u9000\u51fa\u6210\u529f\u3002\u672c\u6b21\u590d\u73b0\u811a\u672c\u91c7\u7528\u5206\u6b65\u6267\u884c\uff0c\u5e76\u68c0\u67e5\u6bcf\u4e2a\u8fd4\u56de\u503c\uff0c\u6240\u4ee5\u672c\u6587\u6570\u5b57\u6765\u81ea\u5b8c\u6210\u6392\u5e8f\u540e\u91cd\u65b0\u8bfb\u53d6\u7684\u7ed3\u679c\uff0c\u800c\u4e0d\u662f\u6839\u636e\u547d\u4ee4\u5f62\u5f0f\u63a8\u6d4b\u5b83\u5e94\u8be5\u751f\u6210\u4ec0\u4e48\u3002<\/p>\n<h2 style=\"font-size:1.38em;line-height:1.5;margin:1.55em 0 .55em;color:#182d43\">\u4e03\u3001\u9ad8\u6bd4\u5bf9\u8d28\u91cf\u4e3a\u4ec0\u4e48\u6ca1\u6709\u963b\u6b62\u8bb0\u5f55\u88ab\u8fc7\u6ee4<\/h2>\n<p>\u7b2c\u4e8c\u8f6e\u8fd0\u884c\u4f7f\u7528\u540c\u4e00\u8f93\u5165\u548c\u7d22\u5f15\uff0c\u589e\u52a0\u5bf9\u6307\u5b9a\u6765\u6e90\u65b9\u5411\u8bbe\u7f6e\u8d28\u91cf\u5931\u8d25\u6807\u8bb0\u7684\u9009\u9879\u3002\u672c\u6587\u8fd9\u6837\u505a\u662f\u4e3a\u4e86\u6f14\u793a\u6807\u8bb0\u7684\u5f71\u54cd\uff0c\u4e0d\u662f\u5efa\u8bae\u4e00\u822c\u6587\u5e93\u542f\u7528\u5b83\u3002<\/p>\n<pre style=\"max-width:100%;box-sizing:border-box;overflow-x:auto;margin:18px 0;padding:14px 16px;background:#f3f6fa;color:#1a293c;border:1px solid #d8e1ea;border-radius:8px;line-height:1.55;white-space:pre;overflow-wrap:normal;word-break:normal\"><code class=\"language-bash\">bwameth.py --threads 2 --reference reference.fa \\\r\n  --set-as-failed f R1.fq R2.fq &gt; qc_forward.sam\r\n<\/code><\/pre>\n<p>\u7ed3\u679c\u4ecd\u6709\u516d\u5341\u6761\u6709\u6548\u5750\u6807\u8bb0\u5f55\uff0c\u6240\u6709\u6bd4\u5bf9\u8d28\u91cf\u4ecd\u4e3a\u516d\u5341\u3002\u4f46\u5176\u4e2d\u4e09\u5341\u6761\u88ab\u8bbe\u7f6e\u8d28\u91cf\u5931\u8d25\u6807\u8bb0\uff0c\u540c\u65f6\u672c\u6b21\u5b9e\u73b0\u6e05\u9664\u4e86\u5b83\u4eec\u7684\u6b63\u5e38\u914d\u5bf9\u6807\u8bb0\u3002\u4ee5\u7b2c\u4e00\u5bf9\u4e3a\u4f8b\uff0c\u4e24\u4e2a\u65d7\u6807\u503c\u7531\u4e5d\u5341\u4e5d\u3001\u4e00\u767e\u56db\u5341\u4e03\uff0c\u53d8\u6210\u516d\u767e\u96f6\u4e5d\u3001\u516d\u767e\u4e94\u5341\u4e03\u3002\u4f4d\u7f6e\u548c\u5e8f\u5217\u672a\u53d8\uff0c\u53d8\u5316\u53d1\u751f\u5728\u8bb0\u5f55\u5c5e\u6027\u3002<\/p>\n<p>\u6211\u4eec\u518d\u7528 Samtools \u5b9e\u9645\u8ba1\u6570\u3002\u53ea\u6392\u9664\u672a\u6bd4\u5bf9\u8bb0\u5f55\uff0c\u4ecd\u7136\u5f97\u5230\u516d\u5341\u6761\u3002\u518d\u6392\u9664\u8d28\u91cf\u5931\u8d25\u8bb0\u5f55\uff0c\u7ed3\u679c\u4e3a\u4e09\u5341\u6761\u3002\u8981\u6c42\u5177\u6709\u6b63\u5e38\u914d\u5bf9\u6807\u8bb0\uff0c\u4e5f\u53ea\u5269\u4e09\u5341\u6761\u3002\u8fd9\u4e09\u4e2a\u6570\u5b57\u90fd\u6b63\u786e\uff0c\u53ea\u662f\u7b5b\u9009\u6761\u4ef6\u4e0d\u540c\u3002<\/p>\n<div class=\"guide-table\" style=\"max-width:100%;overflow-x:auto;margin:18px 0\">\n<table style=\"display:table;width:100%;max-width:100%;min-width:0;table-layout:fixed;border-collapse:collapse;font-size:.92em;line-height:1.6\">\n<thead>\n<tr>\n<th scope=\"col\" style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;background:#f7f3ff;color:#1e293b;vertical-align:top;text-align:left;padding:6px;border:1px solid #c4cfdb\">\u8ba1\u6570\u6761\u4ef6<\/th>\n<th scope=\"col\" style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;background:#f7f3ff;color:#1e293b;vertical-align:top;text-align:left;padding:6px;border:1px solid #c4cfdb\">\u57fa\u7840\u8fd0\u884c<\/th>\n<th scope=\"col\" style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;background:#f7f3ff;color:#1e293b;vertical-align:top;text-align:left;padding:6px;border:1px solid #c4cfdb\">\u6307\u5b9a\u65b9\u5411\u6807\u8bb0\u540e<\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">\u5df2\u6bd4\u5bf9\u7684\u4e3b\u8bb0\u5f55<\/td>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">60<\/td>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">60<\/td>\n<\/tr>\n<tr>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">\u5df2\u6bd4\u5bf9\u4e14\u901a\u8fc7\u8d28\u91cf\u6807\u8bb0<\/td>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">60<\/td>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">30<\/td>\n<\/tr>\n<tr>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">\u5e26\u6b63\u5e38\u914d\u5bf9\u6807\u8bb0\u7684\u8bb0\u5f55<\/td>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">60<\/td>\n<td style=\"white-space:normal;overflow-wrap:anywhere;word-break:normal;padding:6px;vertical-align:top;border:1px solid #c4cfdb\">30<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<\/div>\n<p>\u9ad8\u6bd4\u5bf9\u8d28\u91cf\u4e0d\u80fd\u81ea\u52a8\u8986\u76d6\u5176\u4ed6\u8d28\u91cf\u6807\u8bb0\u3002\u6bd4\u5bf9\u8d28\u91cf\u3001\u672a\u6bd4\u5bf9\u72b6\u6001\u3001\u8d28\u91cf\u5931\u8d25\u72b6\u6001\u548c\u914d\u5bf9\u5173\u7cfb\u662f\u4e0d\u540c\u7ef4\u5ea6\uff0c\u4e0b\u6e38\u7a0b\u5e8f\u53ef\u80fd\u540c\u65f6\u4f7f\u7528\u5176\u4e2d\u51e0\u4e2a\u3002\u9047\u5230\u201c\u4e0a\u6e38\u516d\u5341\u6761\u3001\u4e0b\u6e38\u4e09\u5341\u6761\u201d\u65f6\uff0c\u5148\u6bd4\u8f83\u7b5b\u9009\u89c4\u5219\uff0c\u6bd4\u91cd\u65b0\u8c03\u6574\u6bd4\u5bf9\u53c2\u6570\u66f4\u6709\u9488\u5bf9\u6027\u3002<a style=\"color:#7155b5;text-decoration:underline;text-underline-offset:2px\" href=\"https:\/\/www.htslib.org\/doc\/samtools-flagstat.html\">Samtools \u8ba1\u6570\u6587\u6863<\/a>\u4e5f\u660e\u786e\u533a\u5206\u8d28\u91cf\u901a\u8fc7\u4e0e\u5931\u8d25\u7684\u7c7b\u522b\u3002<\/p>\n<p>\u4e0b\u9762\u4e24\u6761\u662f\u672c\u4f8b\u5b9e\u9645\u6267\u884c\u8fc7\u7684\u8fc7\u6ee4\u8ba1\u6570\u3002\u7b2c\u4e8c\u6761\u7684\u6392\u9664\u503c\u5305\u542b\u672a\u6bd4\u5bf9\u4e0e\u8d28\u91cf\u5931\u8d25\u4e24\u79cd\u6807\u8bb0\u3002\u672c\u4f8b\u6ca1\u6709\u6b21\u8981\u6216\u8865\u5145\u6bd4\u5bf9\uff0c\u56e0\u6b64\u7ed3\u679c\u53ef\u76f4\u63a5\u5bf9\u5e94\u4e0a\u8868\uff0c\u771f\u5b9e\u6570\u636e\u82e5\u5b58\u5728\u8fd9\u4e9b\u8bb0\u5f55\uff0c\u8fd8\u9700\u8981\u53e6\u884c\u660e\u786e\u662f\u5426\u53ea\u7edf\u8ba1\u4e3b\u6bd4\u5bf9\u3002<\/p>\n<pre style=\"max-width:100%;box-sizing:border-box;overflow-x:auto;margin:18px 0;padding:14px 16px;background:#f3f6fa;color:#1a293c;border:1px solid #d8e1ea;border-radius:8px;line-height:1.55;white-space:pre;overflow-wrap:normal;word-break:normal\"><code class=\"language-bash\">samtools view -c -F 4 qc_forward.bam\r\nsamtools view -c -F 516 qc_forward.bam\r\n<\/code><\/pre>\n<h2 style=\"font-size:1.38em;line-height:1.5;margin:1.55em 0 .55em;color:#182d43\">\u516b\u3001\u80fd\u8bfb\u53d6 BAM\uff0c\u4e0d\u7b49\u4e8e\u5177\u5907\u4e0b\u6e38\u6240\u9700\u4fe1\u606f<\/h2>\n<p>\u672c\u6b21 BWA-meth \u8f93\u51fa\u6ca1\u6709 Bismark \u4f7f\u7528\u7684\u7532\u57fa\u5316\u8c03\u7528\u5b57\u7b26\u4e32\u6807\u7b7e\u3002\u56e0\u6b64\uff0c\u4e0d\u80fd\u56e0\u4e3a\u4e24\u4e2a\u5de5\u5177\u90fd\u8f93\u51fa BAM\uff0c\u5c31\u5047\u8bbe\u5b83\u4eec\u7684\u7532\u57fa\u5316\u63d0\u53d6\u7a0b\u5e8f\u53ef\u4ee5\u76f4\u63a5\u4e92\u6362\u3002BAM \u662f\u6570\u636e\u5bb9\u5668\uff0c\u5de5\u5177\u8fd8\u53ef\u80fd\u8981\u6c42\u7279\u5b9a\u6807\u7b7e\u3001\u6392\u5e8f\u65b9\u5f0f\u3001\u53c2\u8003\u4fe1\u606f\u6216\u914d\u5bf9\u7ec4\u7ec7\u3002<\/p>\n<p>\u672c\u4f8b\u7684\u4e24\u7aef\u6709\u56db\u5341\u78b1\u57fa\u91cd\u53e0\uff0c\u4e5f\u63d0\u793a\u540e\u7eed\u4f4d\u70b9\u7edf\u8ba1\u9700\u8981\u5904\u7406\u540c\u4e00\u7247\u6bb5\u91cd\u590d\u8986\u76d6\u7684\u95ee\u9898\u3002\u8bfb\u6bb5\u5c42\u9762\u6709\u516d\u5341\u6761\u8bb0\u5f55\uff0c\u5e76\u4e0d\u610f\u5473\u7740\u6bcf\u4e2a\u8986\u76d6\u4f4d\u7f6e\u90fd\u83b7\u5f97\u516d\u5341\u4e2a\u72ec\u7acb\u5206\u5b50\u7684\u8bc1\u636e\u3002\u672c\u6587\u6ca1\u6709\u5b9e\u9645\u8fdb\u884c\u4f4d\u70b9\u6c47\u603b\uff0c\u6240\u4ee5\u4e0d\u4f1a\u628a\u8fd9\u6bb5\u91cd\u53e0\u8f6c\u6210\u4e00\u4e2a\u672a\u7ecf\u8fd0\u884c\u7684\u7532\u57fa\u5316\u8ba1\u6570\u7ed3\u8bba\u3002<\/p>\n<p>BWA-meth \u7684\u4f18\u52bf\u662f\u80fd\u591f\u5229\u7528\u6210\u719f\u7684\u6bd4\u5bf9\u540e\u7aef\uff0c\u5e76\u8f93\u51fa\u4fbf\u4e8e\u5e38\u89c4\u6bd4\u5bf9\u5de5\u5177\u68c0\u67e5\u7684\u683c\u5f0f\u3002\u5b83\u7684\u5c40\u9650\u662f\u540e\u7eed\u7532\u57fa\u5316\u89e3\u91ca\u4ecd\u9700\u76f8\u5e94\u6d41\u7a0b\uff0c\u7279\u6b8a\u8f6c\u6362\u6807\u7b7e\u4e0e\u8d28\u91cf\u6807\u8bb0\u4e5f\u9700\u8981\u88ab\u4e0b\u6e38\u6b63\u786e\u7406\u89e3\u3002\u5224\u65ad\u4e00\u4e2a\u5206\u6790\u7ec4\u5408\u662f\u5426\u9002\u5408\u81ea\u5df1\uff0c\u5e94\u68c0\u67e5\u8fd9\u4e9b\u63a5\u53e3\u8981\u6c42\uff0c\u800c\u4e0d\u662f\u53ea\u786e\u8ba4\u6587\u4ef6\u6269\u5c55\u540d\u80fd\u591f\u8854\u63a5\u3002<\/p>\n<h2 style=\"font-size:1.38em;line-height:1.5;margin:1.55em 0 .55em;color:#182d43\">\u4e5d\u3001\u8fd9\u6b21\u5b9e\u8dd1\u6700\u503c\u5f97\u4fdd\u7559\u7684\u68c0\u67e5<\/h2>\n<p>\u672c\u6b21\u73af\u5883\u66fe\u51fa\u73b0\u4f9d\u8d56\u5305\u5728\u5f53\u524d Python \u7248\u672c\u4e0b\u7684\u8bed\u6cd5\u8b66\u544a\uff0c\u4f46\u7d22\u5f15\u3001\u6bd4\u5bf9\u3001\u6392\u5e8f\u548c\u9010\u6761\u6838\u9a8c\u5747\u5b8c\u6210\u3002\u6211\u4eec\u4fdd\u7559\u539f\u59cb\u8b66\u544a\uff0c\u6ca1\u6709\u4e3a\u4e86\u8ba9\u65e5\u5fd7\u770b\u8d77\u6765\u5e72\u51c0\u800c\u5c4f\u853d\u5168\u90e8\u6807\u51c6\u9519\u8bef\u3002\u662f\u5426\u6210\u529f\u5e94\u7531\u8fd4\u56de\u503c\u3001\u6587\u4ef6\u53ef\u8bfb\u6027\u548c\u6570\u503c\u6838\u9a8c\u5171\u540c\u5224\u65ad\uff0c\u4e0d\u80fd\u4ec5\u51ed\u65e5\u5fd7\u91cc\u6709\u6ca1\u6709\u63d0\u793a\u6587\u5b57\u3002<\/p>\n<p><a style=\"color:#7155b5;text-decoration:underline;text-underline-offset:2px\" href=\"https:\/\/www.biocloudservice.com\/wordpress\/wp-content\/uploads\/2026\/11\/ysx-bioinformatics-practice-168-170-4bad75ff9d6f.zip\">\u4e0b\u8f7d\u672c\u7ec4\u5b8c\u6574\u811a\u672c\u4e0e\u5b9e\u8dd1\u8bb0\u5f55<\/a>\uff0c\u5728\u672c\u7bc7\u76ee\u5f55\u6267\u884c\u4e0b\u9762\u547d\u4ee4\uff0c\u53ef\u4ee5\u91cd\u65b0\u751f\u6210\u53c2\u8003\u4e0e\u4e09\u5341\u4e8c\u5bf9\u8bfb\u6bb5\u3002\u811a\u672c\u56fa\u5b9a\u955c\u50cf\u6458\u8981\uff0c\u9650\u5236\u6bcf\u4e2a\u5bb9\u5668\u4e24\u6838\u3001\u56db GiB \u5185\u5b58\uff0c\u4fdd\u5b58\u539f\u59cb\u8f93\u5165\u3001\u8f6c\u6362\u53c2\u8003\u68c0\u67e5\u3001\u9010\u6761\u5750\u6807\u548c\u5e8f\u5217\u6bd4\u8f83\u3001\u4e24\u8f6e\u8d28\u91cf\u6807\u8bb0\u7ed3\u679c\u4ee5\u53ca\u8fc7\u6ee4\u8ba1\u6570\uff0c\u5efa\u8bae\u6bcf\u6b21\u4f7f\u7528\u65b0\u76ee\u5f55\u3002<\/p>\n<pre style=\"max-width:100%;box-sizing:border-box;overflow-x:auto;margin:18px 0;padding:14px 16px;background:#f3f6fa;color:#1a293c;border:1px solid #d8e1ea;border-radius:8px;line-height:1.55;white-space:pre;overflow-wrap:normal;word-break:normal\"><code class=\"language-bash\">python3 run.py output\r\n<\/code><\/pre>\n<p>\u628a\u8fd9\u5957\u68c0\u67e5\u8fc1\u79fb\u5230\u771f\u5b9e\u9879\u76ee\u65f6\uff0c\u4f18\u5148\u4fdd\u7559\u53c2\u8003\u5185\u5bb9\u6821\u9a8c\u3001\u8bfb\u540d\u5173\u7cfb\u3001\u8f6c\u6362\u4e0e\u6062\u590d\u7684\u4e00\u81f4\u6027\u3001\u4e0b\u6e38\u8fc7\u6ee4\u6761\u4ef6\u56db\u9879\u3002\u5b83\u4eec\u80fd\u56de\u7b54\u5f88\u591a\u6bd4\u5355\u72ec\u4e00\u4e2a\u6bd4\u5bf9\u7387\u66f4\u5173\u952e\u7684\u95ee\u9898\uff0c\u4e5f\u80fd\u8ba9\u540e\u7eed\u7532\u57fa\u5316\u5206\u6790\u5efa\u7acb\u5728\u6e05\u695a\u7684\u8f93\u5165\u542b\u4e49\u4e0a\u3002<\/p>\n<p class=\"guide-end\" style=\"margin-top:1.6em;padding-top:1.0em;border-top:1px solid #d7e0e9\">\u5f53\u4efb\u52a1\u6269\u5c55\u5230\u591a\u6279\u5168\u57fa\u56e0\u7ec4\u8bfb\u6bb5\uff0c\u9700\u8981\u7edf\u4e00\u4fdd\u5b58\u7d22\u5f15\u3001\u6bd4\u5bf9\u548c\u8fc7\u6ee4\u8bb0\u5f55\u65f6\uff0c\u53ef\u901a\u8fc7<a data-ysx-event=\"contact_click\" style=\"color:#7155b5;text-decoration:underline;text-underline-offset:2px\" href=\"https:\/\/www.biocloudservice.com\/wordpress\/compute-services\/shared-research-workspace\/\">\u4e91\u751f\u4fe1\u751f\u7269\u5206\u6790\u670d\u52a1\u5668\u670d\u52a1<\/a>\u5f00\u5c55\u53ef\u590d\u73b0\u7684\u6279\u91cf\u5904\u7406\u3002<\/p>\n<\/div>\n","protected":false},"excerpt":{"rendered":"<p>\u5b9e\u9645\u8fd0\u884c BWA-meth \u7684\u8f6c\u6362\u7d22\u5f15\u548c\u53cc\u7aef\u6bd4\u5bf9\uff0c\u9a8c\u8bc1\u539f\u59cb\u5e8f\u5217\u6062\u590d\u3001NM \u6807\u7b7e\u8bed\u4e49\u3001SAM \u8f6c BAM \u4e0e\u8d28\u91cf\u5931\u8d25\u6807\u8bb0\u7684\u5f71\u54cd\u3002<\/p>\n","protected":false},"author":1,"featured_media":73257,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_jetpack_feature_clip_id":0,"_jetpack_memberships_contains_paid_content":false,"footnotes":"","jetpack_post_was_ever_published":false},"categories":[49],"tags":[294,1426,1427,958],"class_list":["post-73258","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-software","tag-bam","tag-bwa-meth","tag-1427","tag-958"],"jetpack_sharing_enabled":true,"jetpack_featured_media_url":"https:\/\/i0.wp.com\/www.biocloudservice.com\/wordpress\/wp-content\/uploads\/2026\/11\/bwa-meth-restored-sequence-nm-qc-flag-filter-practice-guide.png?fit=1200%2C840&ssl=1","_links":{"self":[{"href":"https:\/\/www.biocloudservice.com\/wordpress\/wp-json\/wp\/v2\/posts\/73258","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.biocloudservice.com\/wordpress\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.biocloudservice.com\/wordpress\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.biocloudservice.com\/wordpress\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.biocloudservice.com\/wordpress\/wp-json\/wp\/v2\/comments?post=73258"}],"version-history":[{"count":1,"href":"https:\/\/www.biocloudservice.com\/wordpress\/wp-json\/wp\/v2\/posts\/73258\/revisions"}],"predecessor-version":[{"id":73261,"href":"https:\/\/www.biocloudservice.com\/wordpress\/wp-json\/wp\/v2\/posts\/73258\/revisions\/73261"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.biocloudservice.com\/wordpress\/wp-json\/wp\/v2\/media\/73257"}],"wp:attachment":[{"href":"https:\/\/www.biocloudservice.com\/wordpress\/wp-json\/wp\/v2\/media?parent=73258"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.biocloudservice.com\/wordpress\/wp-json\/wp\/v2\/categories?post=73258"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.biocloudservice.com\/wordpress\/wp-json\/wp\/v2\/tags?post=73258"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}