基因数据处理64之bwamem处理500bp和1000bp的记录
来源:互联网 发布:淘宝网怎么找人工客服 编辑:程序博客网 时间:2024/06/17 00:05
xubo@xubo:~/xubo/data/alignment/cs-bwamem$ bwa mem bwaindex/GRCH38chr1L3556522.fasta g38l500N10000.fq >g38l500N10000.bwamem.sam[M::bwa_idx_load_from_disk] read 0 ALT contigs[M::process] read 10000 sequences (5000000 bp)...[M::mem_process_seqs] Processed 10000 reads in 20.158 CPU sec, 20.183 real sec[main] Version: 0.7.12-r1039[main] CMD: /home/xubo/xubo/tools/bwa-0.7.12/bwa mem bwaindex/GRCH38chr1L3556522.fasta g38l500N10000.fq[main] Real time: 26.217 sec; CPU: 20.800 secxubo@xubo:~/xubo/data/alignment/cs-bwamem$ samtools g38l500N10000.bwamem.sam [main] unrecognized command 'g38l500N10000.bwamem.sam'xubo@xubo:~/xubo/data/alignment/cs-bwamem$ samtools flagstatg38l500N10000.bwamem.sam [main] unrecognized command 'flagstatg38l500N10000.bwamem.sam'xubo@xubo:~/xubo/data/alignment/cs-bwamem$ samtools flagstat g38l500N10000.bwamem.sam 10000 + 0 in total (QC-passed reads + QC-failed reads)0 + 0 secondary0 + 0 supplementary0 + 0 duplicates10000 + 0 mapped (100.00% : N/A)0 + 0 paired in sequencing0 + 0 read10 + 0 read20 + 0 properly paired (N/A : N/A)0 + 0 with itself and mate mapped0 + 0 singletons (N/A : N/A)0 + 0 with mate mapped to a different chr0 + 0 with mate mapped to a different chr (mapQ>=5)xubo@xubo:~/xubo/data/alignment/cs-bwamem$ bwa mem bwaindex/GRCH38chr1L3556522.fasta g38l5000N10000.fq >g38l500N100000.bwamem.sam[M::bwa_idx_load_from_disk] read 0 ALT contigs[E::main_mem] fail to open file `g38l5000N10000.fq'.xubo@xubo:~/xubo/data/alignment/cs-bwamem$ bwa mem bwaindex/GRCH38chr1L3556522.fasta g38l500N100000.fq >g38l500N100000.bwamem.sam[M::bwa_idx_load_from_disk] read 0 ALT contigs[M::process] read 20000 sequences (10000000 bp)...[M::process] read 20000 sequences (10000000 bp)...[M::mem_process_seqs] Processed 20000 reads in 41.470 CPU sec, 41.454 real sec[M::process] read 20000 sequences (10000000 bp)...[M::mem_process_seqs] Processed 20000 reads in 40.513 CPU sec, 40.457 real sec[M::process] read 20000 sequences (10000000 bp)...[M::mem_process_seqs] Processed 20000 reads in 41.056 CPU sec, 40.989 real sec[M::process] read 20000 sequences (10000000 bp)...[M::mem_process_seqs] Processed 20000 reads in 42.112 CPU sec, 42.049 real sec[M::mem_process_seqs] Processed 20000 reads in 41.096 CPU sec, 41.084 real sec[main] Version: 0.7.12-r1039[main] CMD: /home/xubo/xubo/tools/bwa-0.7.12/bwa mem bwaindex/GRCH38chr1L3556522.fasta g38l500N100000.fq[main] Real time: 206.985 sec; CPU: 206.697 secxubo@xubo:~/xubo/data/alignment/cs-bwamem$ samtools flagstat g38l500N100000.bwamem.sam 100000 + 0 in total (QC-passed reads + QC-failed reads)0 + 0 secondary0 + 0 supplementary0 + 0 duplicates100000 + 0 mapped (100.00% : N/A)0 + 0 paired in sequencing0 + 0 read10 + 0 read20 + 0 properly paired (N/A : N/A)0 + 0 with itself and mate mapped0 + 0 singletons (N/A : N/A)0 + 0 with mate mapped to a different chr0 + 0 with mate mapped to a different chr (mapQ>=5)xubo@xubo:~/xubo/data/alignment/cs-bwamem$ bwa mem bwaindex/GRCH38chr1L3556522.fasta g38l1000N10000.fq >g38l1000N10000.bwamem.sam[M::bwa_idx_load_from_disk] read 0 ALT contigs[M::process] read 10000 sequences (10000000 bp)...[M::mem_process_seqs] Processed 10000 reads in 57.290 CPU sec, 57.343 real sec[main] Version: 0.7.12-r1039[main] CMD: /home/xubo/xubo/tools/bwa-0.7.12/bwa mem bwaindex/GRCH38chr1L3556522.fasta g38l1000N10000.fq[main] Real time: 58.228 sec; CPU: 57.736 secxubo@xubo:~/xubo/data/alignment/cs-bwamem$ samtools flagstat g38l1000N10000.bwamem.sam 10000 + 0 in total (QC-passed reads + QC-failed reads)0 + 0 secondary0 + 0 supplementary0 + 0 duplicates10000 + 0 mapped (100.00% : N/A)0 + 0 paired in sequencing0 + 0 read10 + 0 read20 + 0 properly paired (N/A : N/A)0 + 0 with itself and mate mapped0 + 0 singletons (N/A : N/A)0 + 0 with mate mapped to a different chr0 + 0 with mate mapped to a different chr (mapQ>=5)xubo@xubo:~/xubo/data/alignment/cs-bwamem$ bwa mem bwaindex/GRCH38chr1L3556522.fasta g38l1000N100000.fq >g38l1000N100000.bwamem.sam[M::bwa_idx_load_from_disk] read 0 ALT contigs[M::process] read 10000 sequences (10000000 bp)...[M::process] read 10000 sequences (10000000 bp)...[M::mem_process_seqs] Processed 10000 reads in 55.334 CPU sec, 55.333 real sec[M::process] read 10000 sequences (10000000 bp)...[M::mem_process_seqs] Processed 10000 reads in 57.380 CPU sec, 57.341 real sec[M::process] read 10000 sequences (10000000 bp)...[M::mem_process_seqs] Processed 10000 reads in 56.189 CPU sec, 56.150 real sec[M::process] read 10000 sequences (10000000 bp)...[M::mem_process_seqs] Processed 10000 reads in 59.506 CPU sec, 59.470 real sec[M::process] read 10000 sequences (10000000 bp)...[M::mem_process_seqs] Processed 10000 reads in 57.034 CPU sec, 56.994 real sec[M::process] read 10000 sequences (10000000 bp)...[M::mem_process_seqs] Processed 10000 reads in 58.887 CPU sec, 58.849 real sec[M::process] read 10000 sequences (10000000 bp)...[M::mem_process_seqs] Processed 10000 reads in 57.591 CPU sec, 57.553 real sec[M::process] read 10000 sequences (10000000 bp)...[M::mem_process_seqs] Processed 10000 reads in 55.443 CPU sec, 55.412 real sec[M::process] read 10000 sequences (10000000 bp)...[M::mem_process_seqs] Processed 10000 reads in 56.343 CPU sec, 56.304 real sec[M::mem_process_seqs] Processed 10000 reads in 58.834 CPU sec, 58.846 real sec[main] Version: 0.7.12-r1039[main] CMD: /home/xubo/xubo/tools/bwa-0.7.12/bwa mem bwaindex/GRCH38chr1L3556522.fasta g38l1000N100000.fq[main] Real time: 573.158 sec; CPU: 572.982 secxubo@xubo:~/xubo/data/alignment/cs-bwamem$ samtools flagstat g38l1000N100000.bwamem.sam 100000 + 0 in total (QC-passed reads + QC-failed reads)0 + 0 secondary0 + 0 supplementary0 + 0 duplicates100000 + 0 mapped (100.00% : N/A)0 + 0 paired in sequencing0 + 0 read10 + 0 read20 + 0 properly paired (N/A : N/A)0 + 0 with itself and mate mapped0 + 0 singletons (N/A : N/A)0 + 0 with mate mapped to a different chr0 + 0 with mate mapped to a different chr (mapQ>=5)
参考
【1】https://github.com/xubo245/AdamLearning【2】https://github.com/bigdatagenomics/adam/ 【3】https://github.com/xubo245/SparkLearning【4】http://spark.apache.org【5】http://stackoverflow.com/questions/28166667/how-to-pass-d-parameter-or-environment-variable-to-spark-job 【6】http://stackoverflow.com/questions/28840438/how-to-override-sparks-log4j-properties-per-driver
研究成果:
【1】 [BIBM] Bo Xu, Changlong Li, Hang Zhuang, Jiali Wang, Qingfeng Wang, Chao Wang, and Xuehai Zhou, "Distributed Gene Clinical Decision Support System Based on Cloud Computing", in IEEE International Conference on Bioinformatics and Biomedicine. (BIBM 2017, CCF B)【2】 [IEEE CLOUD] Bo Xu, Changlong Li, Hang Zhuang, Jiali Wang, Qingfeng Wang, Xuehai Zhou. Efficient Distributed Smith-Waterman Algorithm Based on Apache Spark (CLOUD 2017, CCF-C).【3】 [CCGrid] Bo Xu, Changlong Li, Hang Zhuang, Jiali Wang, Qingfeng Wang, Jinhong Zhou, Xuehai Zhou. DSA: Scalable Distributed Sequence Alignment System Using SIMD Instructions. (CCGrid 2017, CCF-C).【4】more: https://github.com/xubo245/Publications
Help
If you have any questions or suggestions, please write it in the issue of this project or send an e-mail to me: xubo245@mail.ustc.edu.cnWechat: xu601450868QQ: 601450868
阅读全文
0 0
- 基因数据处理64之bwamem处理500bp和1000bp的记录
- 基因数据处理65之bwa处理500bp和1000bp的记录
- 基因数据处理61之idea运行cs-bwamem处理single-end(1条100bp的reads)
- 基因数据处理63之snap修改默认设置后处理大于400bp的记录
- 基因数据处理52之cs-bwamem集群版运行(1千万条100bp的reads)
- 基因数据处理53之cs-bwamem集群版运行paird-end(1千万条100bp的reads)
- 基因数据处理62之snap默认无法处理大于400bp的reads
- 基因数据处理84之cs-bwamem处理小数据集
- 基因数据处理57之BWA-MEM运行single-end(1千万条100bp的reads)
- 基因数据处理59之snap运行single-end(1千万条100bp的reads)
- 基因数据处理60之bwa运行single-end(1千万条100bp的reads)
- 基因数据处理55之cs-bwamem安装记录(idea maven ,没有通过pl)
- 基因数据处理82之cs-bwamem处理SRR003161(参考基因组为GRCH38chr1)
- 基因数据处理82之cs-bwamem处理SRR003161(参考基因组为GRCH38chr1)
- 基因数据处理54之bwa-mem运行paird-end(1千万条100bp的reads)
- 基因数据处理56之bwa运行paird-end(1千万条100bp的reads).md
- 基因数据处理58之snap运行paired-end(1千万条100bp的reads对)
- 基因数据处理44之cloud-scale-bwamem安装
- [非技术]《设计与生存》读书笔记
- Java学习与技术总结——(二)神奇的排序算法
- SSH整合
- 谈谈单例模式(三)
- LeetCode27. Remove Element
- 基因数据处理64之bwamem处理500bp和1000bp的记录
- 基因数据处理65之bwa处理500bp和1000bp的记录
- 基因数据处理66之avocado集群运行
- 基因数据处理67之bwa建立索引时间
- 基因数据处理68之avocado的配置文件默认无法从hdfs读取
- Scrapy_part1_初识SCRAPY和安装
- 基因数据处理69之bowtie安装与使用
- ll -bash: ls: command not found
- 02爬虫---浏览器的模拟Headers属性