Extract Fasta Sequences Sub Sets by position

【Extract Fasta Sequences Sub Sets by position】的更多相关文章

Extract Fasta Sequences Sub Sets by position

cut -d " " -f 1 sequences.fa | tr -s "\n" "\t"| sed -s 's/>/\n/g' > sequences.tab while read id start end; do \ g=$(grep "$id" sequences.tab | cut -f 2 | cut -c $start-$end);\ echo ">$id";\ echo $…

根据位置信息提取 fasta 文件中的序列 -- extract fasta sequence by their position

#!/usr/bin/env python # usages: python extract_seq_by_pos.py input.fasta id_start_end > result.fasta import sys import re FASTA= sys.argv[1] BED= sys.argv[2] fasta= open(FASTA, 'U') fasta_dict= {} for line in fasta: line= line.strip() if line == '':…

Lua 5.1 参考手册

Lua 5.1 参考手册 by Roberto Ierusalimschy, Luiz Henrique de Figueiredo, Waldemar Celes 云风译 www.codingnow.com Copyright © 2006 Lua.org, PUC-Rio. All rights reserved. 1 - 介绍 Lua 是一个扩展式程序设计语言,它被设计成支持通用的过程式编程,并有相关数据描述的设施. Lua 也能对面向对象编程,函数式编程,数据驱动式编程提供很好的支持.…

Getting started with new I/O (NIO)--reference

The new input/output (NIO) library, introduced with JDK 1.4, provides high-speed, block-oriented I/O in standard Java code. This hands-on tutorial covers the NIO library in great detail, from the high-level concepts to under-the-hood programming deta…

NGS中的一些软件功能介绍

1.bowtie 短序列比对工具,blast也是短序列比对工具,速度快,结果易理解. 输入可以是fastq或者fasta文件. 生成比对结果文件sam格式的吧. 2.bwa 转自:https://www.jianshu.com/p/1552cc6ac3be 将DNA序列比对到参考基因组上的软件,包含三种算法: BWA-backtrack:适合比对长度不超过100bp的序列: BWA-SW:合于长度为70-1M bp的序列: BWA-MEM:合于长度为70-1M bp的序列,高质量的测序数据,其比…

T-SQL Recipes之Separating elements

Separating elements Separating elements is a classic T-SQL challenge. It involves a table called Arrays with strings holding comma-separated lists of values in a column called arr. Run the following code to create the Arrays table, and populate it wi…