923 days ago
Perl script to rename the fasta file !
#Script #1 #!/usr/bin/perl -w use strict; #USAGE #perl extractPattern.pl kmerfasta > uniref100_result_broad my %kHash; local $/ = '>'; my $infile2 = "$ARGV[0]"; # Kmer fasta open( FH2, '919 days ago
911 days ago
911 days ago
899 days ago
Installing SEVA environment in Conda !
...aml Collecting package metadata (repodata.json): done Solving environment: done Downloading and Extr...action: done Executing transaction: done # # To activate this environment, use # # $ conda activa...899 days ago
Extract all fasta sequences except ids !
...umber the kmers perl storeKmer.pl omi_kmer19_formated.fa > omi_kmer19_formated_numbered.fa # cat all *.rd file cat *.rd > all_reads_hits #count the lines in each file -- go...889 days ago
886 days ago
886 days ago
bash script to extract sequence by ids !
Use a Perl one-liner, grep and seqtk subseq to extract the desired fasta sequences: # Create test input: cat > in.fasta BGI_novel_T016313 Solyc03g025570.2.1 TT...884 days ago