BASH script for SelfBLAST a genome
...on settings FASTAFILE=MergedContigs.fasta MYDB=myDB OUTFILE=seeRES THREAD=20 SEQ="" echo "User $USE...Sting"; blastn -task megablast -query $SEQ -db $MYDB -evalue 1e-5 -num_threads $THREAD -max_target_seqs 1 -...2649 days ago
2148 days ago
Perl script to find the distance beetween all the contigs and scaffolds
...Seq->{@$SeqIds[$cnt]}{seq}\n"; close $ofh or die "Could not close file '$tmpFile' $!"; $cnt++; system ("mafft --retree 1 --maxiterate 0 --thread 40 $tmpFile > MAFFT_out");...2137 days ago
Install Parrot Virtual Machine !
...has timespec.......................yes. auto::infnan - Is standard C Inf/NaN handling present..............yes. auto::thread - Does your system have threads.......................yes....1534 days ago
Long reads mapper bash script !
...epted toolName=$1 refFasta=$2 longReads=$3 thread=$4 readsType=$5 if [ $# -...ta $bwaMemLoc mem $refFasta $longReads -t $thread > $fileName.out.sam elif [ $...nimap2Loc -ax map-ont $refFasta $longReads -t $thread > $fileName.out.sam elif [...2119 days ago
Rules to run fastp / Snakemake !
...ple}.html", json=RESULTS + "/fastq_trimmed/{sample}.json" threads: 5 log:...--in2 {input.rev} " "--out1 {output.fwd} --out2 {output.rev} --thread {threads} --cut_tail --html {...898 days ago