2902 days ago
2295 days ago
Perl script to check fastq reads qualities !
#!/usr/bin/env perl use strict; use warnings; sub readfq { my ($fh, $aux) = @_; @$aux = [undef, 0] if (!defined(@$aux)); return if ($aux->[1]); if (!def...2185 days ago
2240 days ago
Perl subroutine to read genome/reads fasta file !
sub readSeqFromFasta{ my $file = $_[0]; my (%map,$id,$seq); open(IN,$file) or die("Cannot open file for reading $file:$!\n"); while(){ chomp; if($_ =~ /^>/...2164 days ago
2158 days ago
2157 days ago
Long reads mapper bash script !
#!/bin/bash #only for LONG READS mapping #USAGE: runMapper.sh minimap2 ref.fa reads.fa/fq 48 ont echo -e "This i...$bwaMemLoc mem $refFasta $longReads -t $thread > $fileName.out.sa...pping with $toolName" if [ $readsType == "ont" ]; then $n...2153 days ago
2154 days ago
Installing Porechop on Ubuntu !
➜ Tools git:(master) ✗ git clone https://github...[-v VERBOSITY] [-t THREADS] [-b BARCODE_DIR]..._THRESHOLD] [--check_reads CHECK_READS] [--scoring_schem...adapters found at the ends of reads (default: 2) --end_thresho...2136 days ago