Count GC Content in nucleotide sequence with Perl
...= -1) { usage(); exit; } $fasta_file = $ARGV[0]; $out_file = "gc_out.txt"; unless ( open(IN, "$fasta_file") ) { print "Got a bad fasta file: $fasta_file\n\n";...2893 days ago
Installing Porechop on Ubuntu !
...XTRA_MIDDLE_TRIM_GOOD_SIDE] [--extra_middle_trim_bad_side EXTRA_MIDDLE_TRIM_BAD_SIDE] [--min...dapters on their "good" side (default: 10) --extra_middle_trim_bad_side EXTRA_MIDDLE_TRIM_BAD_SI...2094 days ago
Perl script to reverse complement a DNA sequence !
.../g; $revcom =~ s/G/C/g; $revcom =~ s/C/G/g; print "Here is the reverse complement DNA: WRONG:\n\n"; print "$revcom\n"; print "\nThat was a bad algorithm, and the reverse co...2031 days ago
Tadpole is 250x faster than SPADes assembler !
...lower coverage for any kmer after correction. markbadbases=0 (mbb) Any base fu...this will have its quality reduced. markdeltaonly=t (mdo) Only mark bad bases adjacent to good bases....964 days ago