Parse the NCBI taxonomy database with Perl !
use Bio::DB::Taxonomy; use warnings; my $sps="human"; # Get one from a NCBI taxonomy database my $dbh =...tor_taxon = $dbh->ancestor($nam); print "$ancestor_taxon->id ****"; #my $tree = $dbh->get_tree(...1169 days ago
Perl script to get all the descendent in a tree !
use warnings; use Bio::DB::Taxonomy; #Rotifera 10190 # Get one from a NCBI taxonomy database my $dbh = Bio::DB::Taxonomy->new(-source => 'flatfile',...1161 days ago
1051 days ago
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") basena...1042 days ago
Commandline to Extract a list of specific read IDs from a bam file
#Save all the IDs in IDs.txt file #Run it on BAM file samtools view file.bam | fgrep -w -f IDs.txt993 days ago
978 days ago
947 days ago
661 days ago
Genome Scaffolding and gap filling !
scaffolding with ARCS v1.0.3 (−c3, −l,4, −a,0.9, −z500, −m50, −20 000, −e30000, −s90). https://github.com/bcgsc/arcs Next, automated gap filling was performed using Sealer v2.0.1 (−L150, -P10, −k75-115 [step = 10]) https://github.com/bcgsc/abyss/tree/sealer-release617 days ago
Extract the mapped and unmapped reads !
PROCESSORS=20 #Single_End_Layout: samtools view --threads $PROCESSORS -b -F 4 in.bam > mapped.bam samtools view --threads $PROCESSORS -b -f 4 in.bam > unmapped.ba...587 days ago