• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    迪恩网络公众号

Golang simplestream.SimpleStream类代码示例

原作者: [db:作者] 来自: [db:来源] 收藏 邀请

本文整理汇总了Golang中github.com/abeconnelly/simplestream.SimpleStream的典型用法代码示例。如果您正苦于以下问题:Golang SimpleStream类的具体用法?Golang SimpleStream怎么用?Golang SimpleStream使用的例子?那么恭喜您, 这里精选的类代码示例或许可以为您提供帮助。



在下文中一共展示了SimpleStream类的10个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于我们的系统推荐出更棒的Golang代码示例。

示例1: emit_nocall_ref

func emit_nocall_ref(start_pos int64, n int64, ref_ain *simplestream.SimpleStream, aout *bufio.Writer) (int64, error) {

	end_pos := start_pos + n
	for ; start_pos < end_pos; start_pos++ {

		if ref_ain.Pos >= ref_ain.N {
			if e := ref_ain.Refresh(); e != nil {
				return 0, e
			}
		}

		bp := ref_ain.Buf[ref_ain.Pos]
		ref_ain.Pos++

		switch bp {
		case 'a', 'A':
			aout.WriteByte('A')
		case 'c', 'C':
			aout.WriteByte('C')
		case 'g', 'G':
			aout.WriteByte('G')
		case 't', 'T':
			aout.WriteByte('T')
		default:
			aout.WriteByte(bp)
		}

		gCounter++

	}

	return start_pos, nil

}
开发者ID:abeconnelly,项目名称:pasta,代码行数:34,代码来源:gvcf2pasta.go


示例2: convert

func convert(pa_ain *simplestream.SimpleStream, aout *os.File) error {
	var bp byte

	bufout := bufio.NewWriter(aout)
	defer bufout.Flush()

	for {
		if pa_ain.Pos >= pa_ain.N {
			if e := pa_ain.Refresh(); e != nil {
				return e
			}
		}

		bp = pa_ain.Buf[pa_ain.Pos]
		pa_ain.Pos++

		switch bp {
		case '=', '#', '*', '\'', 'a', 'Q':
			bufout.WriteByte('a')
		case '~', '&', '+', '"', 'c', 'S':
			bufout.WriteByte('c')
		case '?', ':', '-', ',', 'g', 'W':
			bufout.WriteByte('g')
		case '@', ';', '%', '_', 't', 'd':
			bufout.WriteByte('t')
		case 'A', 'C', 'G', 'T', 'n', 'N':
			bufout.WriteByte('n')
		}

	}

	return nil
}
开发者ID:abeconnelly,项目名称:pasta,代码行数:33,代码来源:pasta2fasta.go


示例3: emit_alt

func emit_alt(start_pos int64, n int64, alt_seq string, ref_ain *simplestream.SimpleStream, aout *bufio.Writer) (int64, error) {

	end_pos := start_pos + n
	for ; start_pos < end_pos; start_pos++ {
		if ref_ain.Pos >= ref_ain.N {
			if e := ref_ain.Refresh(); e != nil {
				return 0, e
			}
		}
		ref_ain.Pos++
	}

	sub_len := n
	if n > int64(len(alt_seq)) {
		sub_len = int64(len(alt_seq))
	}

	for i := 0; i < len(alt_seq); i++ {

		if int64(i) < sub_len {
			switch alt_seq[i] {
			case 'a', 'A':
				aout.WriteByte('A')
			case 'c', 'C':
				aout.WriteByte('C')
			case 'g', 'G':
				aout.WriteByte('G')
			case 't', 'T':
				aout.WriteByte('T')
			default:
				aout.WriteByte(alt_seq[i])
			}
		} else {
			switch alt_seq[i] {
			case 'a', 'A':
				aout.WriteByte('b')
			case 'c', 'C':
				aout.WriteByte('d')
			case 'g', 'G':
				aout.WriteByte('h')
			case 't', 'T':
				aout.WriteByte('u')
			default:
				aout.WriteByte(alt_seq[i])
			}
		}

		//aout.WriteByte(alt_seq[i])

	}

	return start_pos, nil

}
开发者ID:abeconnelly,项目名称:pasta,代码行数:54,代码来源:gvcf2pasta.go


示例4: peel_ref

func peel_ref(start_pos int64, n int64, ref_ain *simplestream.SimpleStream) (string, int64, error) {
	refseq := []byte{}
	end_pos := start_pos + n
	for ; start_pos < end_pos; start_pos++ {

		if ref_ain.Pos >= ref_ain.N {
			if e := ref_ain.Refresh(); e != nil {
				return "", 0, e
			}
		}

		bp := ref_ain.Buf[ref_ain.Pos]
		ref_ain.Pos++

		refseq = append(refseq, bp)
	}

	return string(refseq), start_pos, nil
}
开发者ID:abeconnelly,项目名称:pasta,代码行数:19,代码来源:gvcf2pasta.go


示例5: emit_nocall

func emit_nocall(start_pos int64, n int64, ref_ain *simplestream.SimpleStream, aout *bufio.Writer) error {

	end_pos := start_pos + n
	for ; start_pos < end_pos; start_pos++ {

		if ref_ain.Pos >= ref_ain.N {
			if e := ref_ain.Refresh(); e != nil {
				return e
			}
		}

		bp := ref_ain.Buf[ref_ain.Pos]
		ref_ain.Pos++

		switch bp {
		case 'A', 'a':
			aout.WriteByte('A')
		case 'C', 'c':
			aout.WriteByte('C')
		case 'G', 'g':
			aout.WriteByte('G')
		case 'T', 't':
			aout.WriteByte('T')
		case 'N', 'n':
			aout.WriteByte('n')
		default:
			if bp != 'n' && bp != 'N' {
				fmt.Printf("!!!! %c ... s%d, n%d\n", bp, start_pos, n)
				panic(bp)
			}

			aout.WriteByte(bp)
		}

	}

	return nil
}
开发者ID:abeconnelly,项目名称:pasta,代码行数:38,代码来源:gff2pasta.go


示例6: _main

func _main(c *cli.Context) {

	if c.String("input") == "" {
		fmt.Fprintf(os.Stderr, "Input required, exiting\n")
		cli.ShowAppHelp(c)
		os.Exit(1)
	}

	gff_ain, err := autoio.OpenReadScannerSimple(c.String("input"))
	_ = gff_ain
	if err != nil {
		fmt.Fprintf(os.Stderr, "%v", err)
		os.Exit(1)
	}
	defer gff_ain.Close()

	ref_ain := simplestream.SimpleStream{}
	ref_fp := os.Stdin
	if c.String("ref-input") != "-" {
		var e error
		ref_fp, e = os.Open(c.String("ref-input"))
		if e != nil {
			fmt.Fprintf(os.Stderr, "%v", err)
			os.Exit(1)
		}
		defer ref_fp.Close()
	}
	ref_ain.Init(ref_fp)

	var ref_start int64
	ref_start = 0
	ss := c.Int("ref-start")
	if ss > 0 {
		ref_start = int64(ss)
	}

	var seq_start int64
	seq_start = 0
	_ = seq_start
	ss = c.Int("seq-start")
	if ss > 0 {
		seq_start = int64(ss)
	}

	aout := os.Stdout
	if c.String("output") != "-" {
		aout, err = os.Open(c.String("output"))
		if err != nil {
			fmt.Fprintf(os.Stderr, "%v", err)
			os.Exit(1)
		}
		defer aout.Close()
	}

	if c.Bool("pprof") {
		gProfileFlag = true
		gProfileFile = c.String("pprof-file")
	}

	if c.Bool("mprof") {
		gMemProfileFlag = true
		gMemProfileFile = c.String("mprof-file")
	}

	gVerboseFlag = c.Bool("Verbose")

	if c.Int("max-procs") > 0 {
		runtime.GOMAXPROCS(c.Int("max-procs"))
	}

	allele := c.Int("allele")

	if gProfileFlag {
		prof_f, err := os.Create(gProfileFile)
		if err != nil {
			fmt.Fprintf(os.Stderr, "Could not open profile file %s: %v\n", gProfileFile, err)
			os.Exit(2)
		}

		pprof.StartCPUProfile(prof_f)
		defer pprof.StopCPUProfile()
	}

	e := convert(&gff_ain, &ref_ain, aout, ref_start, allele)
	if e != nil && e != io.EOF {
		panic(e)
	}

	aout.Sync()

}
开发者ID:abeconnelly,项目名称:pasta,代码行数:91,代码来源:gff2pasta.go


示例7: emit_alt

func emit_alt(start_pos int64, n int64, alt_seq string, ref_ain *simplestream.SimpleStream, aout *bufio.Writer) error {
	ref_pos := 0

	var ref_bp byte

	for i := 0; i < len(alt_seq); i++ {

		if int64(i) < n {
			if ref_ain.Pos >= ref_ain.N {
				if e := ref_ain.Refresh(); e != nil {
					return e
				}
			}
			ref_bp = ref_ain.Buf[ref_ain.Pos]

			// REFERENCE CHECK
			//
			lc_ref_bp := ref_bp
			lc_gff_bp := g_GFF_REF[ref_pos]

			if lc_ref_bp == 'A' {
				lc_ref_bp = 'a'
			} else if lc_ref_bp == 'C' {
				lc_ref_bp = 'c'
			} else if lc_ref_bp == 'G' {
				lc_ref_bp = 'g'
			} else if lc_ref_bp == 'T' {
				lc_ref_bp = 't'
			}

			if lc_gff_bp == 'A' {
				lc_gff_bp = 'a'
			} else if lc_gff_bp == 'C' {
				lc_gff_bp = 'c'
			} else if lc_gff_bp == 'G' {
				lc_gff_bp = 'g'
			} else if lc_gff_bp == 'T' {
				lc_gff_bp = 't'
			}

			if lc_ref_bp != lc_gff_bp {
				fmt.Printf("\nREF MISMATCH: GFF reported %c (at %d+%d) but got %c\n", g_GFF_REF[ref_pos], start_pos, ref_pos, ref_bp)
				panic("!!")
			}
			//
			// REFERENCE CHECK

			ref_ain.Pos++
			ref_pos++

			switch ref_bp {
			case 'a', 'A':
				ref_bp = 'a'
			case 'c', 'C':
				ref_bp = 'c'
			case 'g', 'G':
				ref_bp = 'g'
			case 't', 'T':
				ref_bp = 't'
			case 'n', 'N':
				ref_bp = 'n'
			default:
				return fmt.Errorf("invalid character for reference stream ('%c') at %d", ref_bp, ref_pos)
			}

			// It's considered a sub
			//
			if (ref_bp == 'n' || ref_bp == 'N') && (alt_seq[i] != 'n' && alt_seq[i] != 'N') {
				fmt.Printf("\n\n>>>> start_pos %d, n %d, alt_seq %s, ref_bp %c, alt_seq[%d] %c\n\n", start_pos, n, alt_seq, ref_bp, i, alt_seq[i])
				//panic("whoa!!!")
			}

			switch alt_seq[i] {
			case 'a', 'A':
				aout.WriteByte(gSub[ref_bp]['a'])
			case 'c', 'C':
				aout.WriteByte(gSub[ref_bp]['c'])
			case 'g', 'G':
				aout.WriteByte(gSub[ref_bp]['g'])
			case 't', 'T':
				aout.WriteByte(gSub[ref_bp]['t'])
			case 'n', 'N':
				if ref_bp == 'n' || ref_bp == 'N' {
					fmt.Printf("[email protected]! %c (s%d,n%d) [%s]{%d,%c}\n", ref_bp, start_pos, n, alt_seq, i, alt_seq[i])
					panic("-->")
				}
				aout.WriteByte(gSub[ref_bp]['n'])
			default:
				return fmt.Errorf("invalid character for alt sequence ('%c') at pos %d", alt_seq[i], i)
			}

		} else {

			// It's considered an insertion
			//
			switch alt_seq[i] {
			case 'a', 'A':
				aout.WriteByte('Q')
			case 'c', 'C':
				aout.WriteByte('S')
//.........这里部分代码省略.........
开发者ID:abeconnelly,项目名称:pasta,代码行数:101,代码来源:gff2pasta.go


示例8: convert

func convert(fa_ain, ref_ain *simplestream.SimpleStream, aout *os.File, start_pos int64) error {
	var e error
	var fa_bp byte

	allele_num := 0

	bufout := bufio.NewWriter(aout)
	defer bufout.Flush()

	for {
		if fa_ain.Pos >= fa_ain.N {
			if e := fa_ain.Refresh(); e != nil {
				return e
			}
		}

		if ref_ain.Pos >= ref_ain.N {
			if e := ref_ain.Refresh(); e != nil {
				return e
			}
		}

		fa_bp = fa_ain.Buf[fa_ain.Pos]
		fa_ain.Pos++

		ref_bp = ref_ain.Buf[ref_ain.Pos]
		ref_ain.Pos++

		if fa_bp != 'n' && fa_bp != 'N' && fa_bp == ref_bp {
			switch fa_bp {
			case 'A':
				aout.WriteByte('a')
			case 'C':
				aout.WriteByte('c')
			case 'G':
				aout.WriteByte('g')
			case 'T':
				aout.WriteByte('t')
			default:
				aout.WriteByte(fa_bp)
			}
		} else if fa_bp == 'N' || fa_bp == 'n' {
			switch ref_bp {
			case 'A', 'a':
				aout.WriteByte('!')
			case 'C', 'c':
				aout.WriteByte('#')
			case 'G', 'g':
				aout.WriteByte('\'')
			case 'T', 't':
				aout.WriteByte('4')
			default:
				aout.WriteByte(ref_bp)
			}
		} else if fa_bp != ref_bp {
			switch fa_bp {
			case 'a', 'A':
				aout.WriteByte('A')
			case 'c', 'C':
				aout.WriteByte('C')
			case 'g', 'G':
				aout.WriteByte('G')
			case 't', 'T':
				aout.WriteByte('T')
			default:
				aout.WriteByte(fa_bp)
			}
		} else {
			aout.WriteByte('-')
		}

	}

	return nil
}
开发者ID:abeconnelly,项目名称:pasta,代码行数:75,代码来源:fasta2pasta.go


示例9: _main

func _main(c *cli.Context) {
	var err error

	if c.String("input") == "" {
		fmt.Fprintf(os.Stderr, "Input required, exiting\n")
		cli.ShowAppHelp(c)
		os.Exit(1)
	}

	pasta_ain := simplestream.SimpleStream{}
	pasta_fp := os.Stdin
	if len(c.String("input")) > 0 && c.String("input") != "-" {
		var e error
		pasta_fp, e = os.Open(c.String("input"))
		if e != nil {
			fmt.Fprintf(os.Stderr, "%v", e)
			os.Exit(1)
		}
		defer pasta_fp.Close()
	}
	pasta_ain.Init(pasta_fp)

	aout := os.Stdout
	if c.String("output") != "-" {
		aout, err = os.Open(c.String("output"))
		if err != nil {
			fmt.Fprintf(os.Stderr, "%v", err)
			os.Exit(1)
		}
		defer aout.Close()
	}

	/*
	  aout := os.Stdout
	  if err!=nil {
	    fmt.Fprintf(os.Stderr, "%v", err)
	    os.Exit(1)
	  }
	*/

	if c.Bool("pprof") {
		gProfileFlag = true
		gProfileFile = c.String("pprof-file")
	}

	if c.Bool("mprof") {
		gMemProfileFlag = true
		gMemProfileFile = c.String("mprof-file")
	}

	gVerboseFlag = c.Bool("Verbose")

	if c.Int("max-procs") > 0 {
		runtime.GOMAXPROCS(c.Int("max-procs"))
	}

	if gProfileFlag {
		prof_f, err := os.Create(gProfileFile)
		if err != nil {
			fmt.Fprintf(os.Stderr, "Could not open profile file %s: %v\n", gProfileFile, err)
			os.Exit(2)
		}

		pprof.StartCPUProfile(prof_f)
		defer pprof.StopCPUProfile()
	}

	convert(&pasta_ain, aout)

}
开发者ID:abeconnelly,项目名称:pasta,代码行数:70,代码来源:pasta2fasta.go


示例10: convert

func convert(ref_ain *simplestream.SimpleStream, seq_ain *simplestream.SimpleStream, fout *os.File, start_pos int64, allele_num int) error {

	aout := bufio.NewWriter(fout)
	if g_convert_debug {
		defer func() { fmt.Printf("\n"); aout.Flush() }()
	} else {
		defer aout.Flush()
	}

	icount := 0
	ocount := 0

	//DEBUG
	if g_convert_debug {
		fmt.Printf("start\n")
	}

	for true {

		//DEBUG
		if g_convert_debug {
			fmt.Printf("ref_ain.Pos %d, ref_ain.N %d\n", ref_ain.Pos, ref_ain.N)
		}

		if ref_ain.Pos >= ref_ain.N {
			if e := ref_ain.Refresh(); e != nil {
				return e
			}
		}
		bp_ref := ref_ain.Buf[ref_ain.Pos]
		ref_ain.Pos++

		//DEBUG
		if g_convert_debug {
			fmt.Printf("seq_ain.Pos %d, seq_ain.N %d\n", seq_ain.Pos, seq_ain.N)
		}

		if seq_ain.Pos >= seq_ain.N {
			if e := seq_ain.Refresh(); e != nil {
				return e
			}
		}
		bp_seq := seq_ain.Buf[seq_ain.Pos]
		seq_ain.Pos++

		//DEBUG
		if g_convert_debug {
			fmt.Printf("%d,%d bp_ref %c, bp_seq %c\n", icount, ocount, bp_ref, bp_seq)
		}

		if bp_ref == bp_seq && bp_ref == '\n' {
			break
		}

		if bp_ref == bp_seq {

			switch bp_ref {
			case 'a', 'A':
				aout.WriteByte(gSub[bp_ref]['a'])
			case 'c', 'C':
				aout.WriteByte(gSub[bp_ref]['c'])
			case 'g', 'G':
				aout.WriteByte(gSub[bp_ref]['g'])
			case 't', 'T':
				aout.WriteByte(gSub[bp_ref]['t'])
			case 'n', 'N':
				if bp_ref == 'n' || bp_ref == 'N' {
					return fmt.Errorf("no-call to no-call match: bf_ref:%c i:%d o:%d\n", bp_ref, icount, ocount)
				}
				aout.WriteByte(gSub[bp_ref]['n'])
			default:
				return fmt.Errorf("invalid character for alt sequence ('%c') at i:%d o:%d", bp_ref, icount, ocount)
			}

		} else if bp_seq == '-' {

			switch bp_ref {
			case 'a', 'A':
				aout.WriteByte('!')
			case 'c', 'C':
				aout.WriteByte('$')
			case 'g', 'G':
				aout.WriteByte('7')
			case 't', 'T':
				aout.WriteByte('E')
			case 'n', 'N':
				aout.WriteByte('z')
			default:
				return fmt.Errorf("invalid character for alt sequence ('%c') at i:%d o:%d", bp_ref, icount, ocount)
			}

		} else if bp_ref == '-' {

			switch bp_seq {
			case 'a', 'A':
				aout.WriteByte('Q')
			case 'c', 'C':
				aout.WriteByte('S')
			case 'g', 'G':
				aout.WriteByte('W')
//.........这里部分代码省略.........
开发者ID:abeconnelly,项目名称:pasta,代码行数:101,代码来源:align2pasta.go



注:本文中的github.com/abeconnelly/simplestream.SimpleStream类示例整理自Github/MSDocs等源码及文档管理平台,相关代码片段筛选自各路编程大神贡献的开源项目,源码版权归原作者所有,传播和使用请参考对应项目的License;未经允许,请勿转载。


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
Golang sdl.Renderer类代码示例发布时间:2022-05-24
下一篇:
Golang go-http-auth.DigestAuth类代码示例发布时间:2022-05-24
热门推荐
热门话题
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap