扫二维码与项目经理沟通
我们在微信上24小时期待你的声音
解答本文疑问/技术咨询/运营咨询/技术建议/互联网交流
现代操作系统使用LBA而不是CHS来记录硬盘分区.如果用扇区代替柱面,我们将看到:
专注于为中小企业提供成都网站建设、网站设计服务,电脑端+手机端+微信端的三站合一,更高效的管理,为中小企业临湘免费做网站提供优质的服务。我们立足成都,凝聚了一批互联网行业人才,有力地推动了1000多家企业的稳健成长,帮助中小企业通过网站建设实现规模扩充和转变。
sfdisk -uS -l /dev/sda
fdisk –l 命令是按照柱面来进行分区显示的,所以需要按照垂直的空间进行理解,但现在的计算机在显示时仍然会按照扇区来理解(因为现在磁盘是使用LBA(使用扇区进行寻址)取代以前的CHS(柱面、磁头、扇区)方式记录硬盘分区的),所以会判断分区没有在柱面结束;如果我们使用命令 sudo fdisk –uS -l ,这个命令的意思是:“give size in sectors instead of sylinders”即是按照扇区的来进行分区显示,这样就可以看到以扇区方式显示分区的了。可以看到相邻逻辑分区之间都有 64 扇区的空间,那里存放逻辑分区表。
一般来说,Linux创建分区使用fdisk命令,可以自动对齐磁盘。然而,fdisk无法处理大于2TB的磁盘文件。这种情况下,需要使用parted命令去创建分区。
parted创建分区需要手动指定分区开始与结束位置,可能会面临磁盘不对齐的问题(提示:Warning: The resulting partition is not properly aligned for best performance.)
所以需要手动去计算分区起点来对齐磁盘。
输入0.00T有时候可以免于计算起点
所以分区命令是
检查是否正确对齐
如果还有问题,则需要进一步修正,可以参考
linux的命令行默认采用行缓冲模式,换行符\n起到刷新输出缓冲区的作用。
(一)安装bowtie
Bowtie可以在个人计算机上使用,也可以在CSC服务器上使用终端连接。请参阅以下文档的第一部分,了解如何在笔记本电脑上安装Bowtie。特别是对他们的计算机没有管理员权限的那些应该确保软件的正确安装和功能。Bowtie也可以在服务器计算机上远程使用。我们将提供临时帐户访问CSC,但你将需要一个安全Shell终端程序进行通信。默认情况下,Mac和Linux上都有这样的程序,但需要安装Windows。普遍的实现是PuTTY。即使终端程序不用于读取映射,也将需要其他练习,并且应该可用。Bowtie的安装:从下载页面下载相应的版本(Linux,Mac或Win,小编使用的是在Linux下进行)。将zip文件解压缩到新的目录中,并转到该目录。下载的bowtie包装包含大肠杆菌基因组的预先建立的指数,以及从该基因组模拟的一组1000个35bp的读数。要使用Bowtie对齐这些读取,请键入以下命令。bowtiee_colireads/e_coli_1000.fqmap_result.txt
如果你收到错误消息"commandnotfound",请尝试在"bowtie"(./bowtie)之前添加"./"。
(二)使用Bowtie
(1)Mapping
要使用Bowtie对齐示例读取,请发出以下命令。bowtiee_colireads/e_coli_1000.fqmap_result.txt
如果你收到错误消息"commandnotfound",请尝试在"bowtie"(./bowtie)之前添加"./"。"e_coli"与"indexes/e_coli"相同。你可以在文本编辑器中打开map_result.txt。每行都是一个读取对齐。对齐读取的名称显示在第一列中。对于Mac和Linux,使用"少"会更好。
lessmap_result.txt#extrareading
ReadthemanualinthefolderorwebsitetogetadeeperunderstandinghowBowtieworksandfurtheroptionsinBowtie.
我们来看看Bowtie在1中使用的一些不同的选项,报告所有有效的对齐方式与一些不匹配。
./bowtie-a-v2e_coli--suppress1,5,6,7-cATGCATCATGCGCCAT-a/--all报告每个读取或对的所有有效对齐(默认值:off)
-v
最多不相匹配的报告对齐
-c
查询序列在命令行
--suppress
上以默认输出模式抑制输出列
2限制对齐
$./bowtie-k3-v2e_coli--suppress1,5,6,7-cATGCATCATGCGCCAT-k
每次读取或配对时报告有效对齐(默认值:1)。
3不匹配排名
$./bowtie-a--best-v2e_coli--suppress1,5,6,7-cATGCATCATGCGCCAT
所有相同的对齐方式按最佳到最坏的顺序进行报告
4只有最不匹配
$./bowtie-a--best--strata-v2--suppress1,5,6,7e_coli-cATGCATCATG
(2)配对对齐
当使用-1和-2选项指定正确配对的读取文件时,Bowtie可以对齐配对端读取(对于原始,FASTA或FASTQ读取文件)
./bowtiee_coli-1reads/e_coli_1000_1.fq-2reads/e_coli_1000_2.fqmap_paired.txt
SAMtools()是一套用于存储,操纵和分析对齐方式的工具,例如Bowtie输出的对齐方式。bowtie-Se_colireads/e_coli_1000.fqec.sam
我们可以再次检查sam文件以查看与txt文件的区别(也是在r4,r5中未映射的读取)。接下来,我们将SAM文件转换为BAM以准备排序。
samtoolsview-bS-oec.bamec.sam
接下来,我们对BAM文件进行排序,
samtoolssortec.bamec.sorted
这样我们就简单的对bam文件中的基因组进行配对对齐。
sda3 不是逻辑分区,扩展分区没对齐没事。真正存储数据的是逻辑分区,不是扩展分区。
所以sda5 和 sda6 对齐了就行了。
当然, sda3 能对齐当然更好了。
我们在微信上24小时期待你的声音
解答本文疑问/技术咨询/运营咨询/技术建议/互联网交流