Skip to content

Instantly share code, notes, and snippets.

@maasha

maasha/test.rb Secret

Created September 12, 2014 13:41
Show Gist options
  • Save maasha/0b512afcb3ab512ef08d to your computer and use it in GitHub Desktop.
Save maasha/0b512afcb3ab512ef08d to your computer and use it in GitHub Desktop.
Download ZIP
Raw
test.rb
#!/usr/bin/env ruby
require 'biopieces'
p = BP.new.
read_fastq(input: "Bioconval-C-75-Caecum_S2_L001_R1_001.fastq", input2: "Bioconval-C-75-Caecum_S2_L001_R2_001.fastq").
trim_seq.
mean_scores.
grab(select: "ATCG", keys: :SEQ).
plot_scores(output: "bp2.png", force: true).
write_fastq(output: "bp2.fq", force: true).
run(fork: false, verbose: true)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment