python全基因组读取，python基因差异分析

全基因组测序数据获取后应该怎么分析?

如果想用你现有的数据找疾病相关的基因／通路／网络，出门直走差异表达分析。

首先进行基因分类，比如说编码性基因占多大比例，非编码性基因又占多少比例；转录因子占多少比例，蛋白激酶类基因又占多少比例等等。然后将该物种基因组与其它已测序基因组进行比较，包括大小、同源度等等。

如果想用你现有的数据找疾病相关的基因／通路／网络，出门直走差异表达分析。首先，看测序峰图，如果结果的彩图显示峰型是尖锐单一的，碱基所对应的编码是一致的，那么可以判断序列是可以使用的。

．数据量产出总碱基数量、Totally mapped reads、Uniquely mapped reads统计，测序深度分析。

测序得到基因序列后一般都需要进行序列比对，看和目的序列的差异情况，常用的软件有DNAman，还有invitrogen的vectorVI也不错。此外还可以直接在网上进行比对，推荐NCBI网站的BLAST可以直接网上进行。

相关问答

Q1: DNA甲基化数据分析(二)

DML是甲基化差异位点，DMR为甲基化差异区域。

第六列为未甲基化数目 在这里，需要对我们的数据利用R，linux或者python整理成DSS包所需输入文件格式。这里我使用python整理的输入文件格式。

甲基化特异性的PCR（Methylation-specific PCR，MSP）用亚硫酸氢盐处理基因组DNA，所有未发生甲基化的胞嘧啶被转化为尿嘧啶，而甲基化的胞嘧啶不变；随后设计针对甲基化和非甲基化序列的引物进行PCR。

MSAP是利用对DNA甲基化敏感的两种同裂酶Hpa II和Msp I来对基因组DNA进行酶切和连接。由于两种酶能够识别相同的限制性酶切位点，即CCGG位点。

和RNA-seq前期流程类似 -- 质控、去接头、比对参考基因组、排序 后期就是要提取甲基化位点，包括CpG、CHG、CHH三种context，H代表非G位点（A、C、T）。

Q2: 怎么读取千人基因组中的vcf数据格式文件

可以下载一个TextForever软件，进行转码。需要转换编码，下载一个TextForever的转化软件，大小只有550k，在PC上直接运行编码转换，将原VCF文件的Unicode转换为GB即可在电脑上查看了。

vcfR 可以直接读取vcf格式的数据。如果同时读取参照序列fasta格式的序列文件和gff格式文件的注释文件还可以获取更完整的信息(此步骤并非必须，可以只读取vcf数据)。在此处便于重复用到了 pinfsc50 包。

单击菜单栏的“文件”，打开的下拉菜单单击“导入和导出”命令。打开“导入和导出向导”对话框，在请选择要执行的操作列表框中选中“导入 VCARD 文件（.vcf）”，单击“下一步”按钮。

可以使用vcf通讯录编辑器来打开。vcf通讯录编辑器(VCFEditor)能够在电脑上查看编辑手机里通讯录的备份文件。方便实用。

微信如何打开vcf格式的文件 vcf文件也叫vcard是联系人名片文件，用微信是不可以打开的，因为微信无法解析，一般可以用记事本或者outlook打开。

一般的VCF文件都很大，用手动提取里面的信息肯定不大现实。用 vcfR 就可以轻松实现。vcfR 自带测试文件 vcfR_test 。就用这个文件来操作一下吧。

Q3: 用BEDtools/Python序列截取

1、bedtools最基本的用法不用提供正负链信息。其实默认的就是截取基因组的区间序列，既然参考基因组是正链，那么默认截取的就是正链序列。如果bed文件有正负链信息，负链的序列就是软件默认情况下截取的序列的反向互补序列。

2、在Python中，获取序列中某个元素的过程称为索引。

3、序列 序列中的索引操作 序列中的元素都是有序的，每一个元素都带有序号，这个序号叫 索引。索引有正值索引和负值索引之分。加乘操作 切片操作 序列的切片（Slicing）就是从序列中切分出小的子序列。

4、以下运行得到的结果仍然是blast的tabular格式(之后可以经过一些简单的shell命令处理，可转成bed格式，结合bedtools批量提取序列)。

Q4: Day59-用Perl和Python脚本提取FASTA中最长转录本氨基酸序列

我给你perl的解决思路：首先，你要确定你的fasta文件的内容的规律性。比如每段序列的开始是不是都会有一些特殊的标志。那么可以用next函数，将这些不是序列的内容跳过。

cat *.fasta single_all_fasta.fasta 将所有fasta序列整合到一个fasta格式中。

fasta格式是一种基于文本用于表示核酸序列或 多肽 序列的格式。其中 核酸 或 氨基酸 均以单个字母来表示，且允许在序列前添加序列名及注释。该格式已成为 生物信息学 领域的一项标准。

Q5: ab1文件怎么打开

一般可以用Chromas或Bioedit打开。可打开AB1文件的软件：GSLBiotechSnapGene，AppliedBiosystemsSequencingAnalysisSoftware，BioEdit，GeospizaFinchTV，TechnelysiumChromasLiteorChromasPro，CubicDesignDNABaser。

方法1：使用默认的图片编辑软件打开图片文件。当你需要将一个图片文件转换成另一种格式时，最简单的解决方式是使用操作系统默认的图片编辑软件。在Windows系统中，使用画图程序。在Mac电脑上，使用预览程序。

使用Chromas可以打开ab1文件。操作如下：首先找到文件后缀名以“ab1”格式存在的文件。

python全基因组读取的介绍就聊到这里吧，感谢你花时间阅读本站内容，更多关于python基因差异分析、python全基因组读取的信息别忘了在本站进行查找喔。