zoukankan      html  css  js  c++  java
  • 生物信息bowtie比对软件的结果格式说明

           生物信息分析中会用到很多的比对软件,比较常用的有bowtie、bowtie2、bwa等,比对文件的标准格式是sam格式,但是bowtie比对默认输出的格式却不是sam格式,由于bowtie适用于短序列比对,并且看突变碱基比较方便,因此它的默认输出格式还是有一定优势的,下面就来说明一下它的默认输出格式。

    example:

    说明:

    1. query id,也就是你自己用于比对的序列的ID。
    2. + 表示正向比对上这条reads;- 表示对反向比对上这条reads。
    3. 对应比对上的参考序列的ID,ref。
    4. 自己用于比对的reads比对到参考序列的起始位置,注意下标是从0开始计数,例如10代表从第11个碱基开始比对上。第2列为"+"时,表示query第一个碱基map到reference上的位置;第2列为"-"时,表示query的反向互补序列第一个碱基map到reference上的位置。
    5. 如果第2列为"+",则和query序列一致;否则,和query序列反向互补。
    6. 质量文件,如果query文件为fasta格式,则无法获取质量文件,用I代替,I的数量与query序列长度一致。
    7. 当前query能map到GENOME的4个不同位置
    8. 如果存在第8列,表示有错配,bowtie的mismatch最大可设置为3。第8列可以分为三个部分,最左端的数字,中间的碱基为reference碱基,最右端的碱基为query碱基,下面分情况讨论:
    第2列为"+"时:最左端的数字9表示query从5'端数起,第10个碱基为"T",而对应的reference为"G";
    第2列为"-"时:最左端的数字9表示query先作反向互补,然后从3'端数起,第10个碱基为"T",而对应的reference为"G";

  • 相关阅读:
    Linux部署Spingboot项目
    Linux Centos7yum安装Mysql8.0.21
    Linux配置网络yum源,提高下载速度
    Linux安装jdk1.8
    Spring的AOP
    Spring的事务管理
    Maven项目中,使用mybatis,根据数据库自动生成pojo实体类、dao、mapper
    Ubuntu14.04中使用docker容器部署tomcat镜像+java web项目
    mybatis
    spl
  • 原文地址:https://www.cnblogs.com/liyichen/p/10843855.html
Copyright © 2011-2022 走看看