bwa比对结果怎么统计比对率

BWA（Burrows-Wheeler Alignment）是一种常用的基因组比对工具，用于将测序数据与参考基因组进行比对。在使用BWA对样本进行比对后，可以根据比对结果来计算比对率。下面是一种常用的统计方法：

首先，从BWA比对输出的文件中提取比对信息。BWA比对结果一般以SAM（Sequence Alignment/Map）格式保存，可以使用SAMtools等相关工具对比对结果进行解析和提取。

从SAM文件中，可以获得每个比对结果的相关信息，如比对的状态（比对成功或失败）、比对上的read数量等。

比对率通常被定义为成功比对的read数量占总read数量的比例。因此，可以通过计算成功比对的read数量除以总read数量得到比对率。

具体的计算步骤如下：
1. 从SAM文件中按行读取每个比对结果。
2. 对于每行结果，判断比对的状态是否为成功（比对状态通常在每行的FLAG列中标注）。
3. 如果比对成功，则成功比对的read数量加1。
4. 继续读取下一行并重复上述步骤，直到读取完所有比对结果。
5. 计算成功比对的read数量除以总read数量，并乘以100%得到比对率（以百分比表示）。

需要注意的是，统计比对率时，应该使用有效的read数量。在有些情况下，比如双端测序，一对read中的其中一条无法比对成功，但另一条可以成功比对，这种情况下，只统计比对成功的read数量。

总之，通过解析BWA比对结果的SAM文件，并统计成功比对的read数量除以总read数量，就可以得到比对率。