Corona variant calling steps !
#!/bin/bash FILE_PATH="data/trimmed_fastq_small/" find "$FILE_PATH" -name "*.fasta" | while IFS= read -r my_file do filename=$(basename "$my_file") basename=$(basename "$my_file" .fasta)...1102 days ago
1053 days ago
1047 days ago
Commandline for paired end reads simulation with BBMap !
...ref=mixed.fa out=reads_BBMAP250.fq paired interleaved reads=100k le...ref=mixed.fa out=reads_BBMAP250.fq paired interleaved reads=100k le...ef=mixed.fa, out=reads_BBMAP250.fq, paired, interleaved, reads=100k,..._QUALITY=true REPLACE_NOREF=false paired=true read length=250 re...1047 days ago
List of string comparison algorithms !
String comparison: Levenshtein Distance Damerau-Levenshtein Distance Jaro Distance Jaro-Winkler Distance Match Rating Approach Comparison Hamming Distance More at https://jellyfish.readthedocs.io/en/latest/comparison.html1039 days ago
1038 days ago
Python script to read FASTA and FASTQ file !
...ysam import FastxFile def read_fasta_q_file(fasta_q_file): """Parse FASTA/Q file using `pysam.FastxFile`. Args: fasta_q_file (str): Path to FASTA/Q file. ""...1033 days ago
Tadpole is 250x faster than SPADes assembler !
lege@jit-Lenovo-ideapad-320-15ISK:~/Downloads/MyTool...Second input file for paired data. extra= Extr...inprob for the prefilter. prepasses=1 Use this many p...enome size estimation. contigpasses=16 Build contigs with...thbad=f (rbb) Only discard pairs if both reads are low-dept...1033 days ago
1019 days ago
Run multiple bash command in screen !
#login to screen screen -r 123 #bash.sh srun --partition=compute --nodes=1 --ntasks-per-node=40 --pty bash.sh #Run and check the status in screen contol +A +D1019 days ago