2686 days ago
2105 days ago
910 days ago
Bash script to simulate a genome !
...stq.gz in2=illumina2.fastq.gz out=illumina.int.fastq 2>/dev/null # use KmerGenie 1.7051 to get an idea of k-mer with that produces longest N50 # http://kmergenie.bx.psu.ed...900 days ago
Extract fasta header with ids !
...out.fa more out.fa_class.txt | grep "227859" | awk '{print $2}' > all_real_ids.txt minimap2 -t 36 -k19 -w5 -A1 -B2 -O3,13 -E2,1 -s200 -z200 -N50 --min-occ-floor=100 finaal_ou...859 days ago
Perl script to calculate the basic stats of the assembled genome !
...ngths[0]; my $avg_contig_length = $total_length / $num_contigs; my $median_contig_length = calculate_median(\@contig_lengths); # Calculate N50 my $n50 = calculate_n50(\@contig_leng...76 days ago
Python script for basic stats of the assembled genome !
...an_contig_length = statistics.median(contig_lengths) # Calculate N50 def calculate_n50(lengths): total_size = s...") print(f"Median Contig Length: {median_contig_length}") print(f"N50: {calculate_n50(contig_length...76 days ago