Extract fasta sequences with ids in another file !
#Ids are in test.txt - one ids per line #sequences are in test.fa grep -w -A 2 -f test.txt test.fa --no-group-separator # seqtk seqtk subseq test.fa test.txt #faSomeRecods faSomeRecords in.fa listFile out.fa # seqkit seqkit grep -n -f list.txt sequences.fas > newfile2.fas948 days ago
Run Pango on your multifasta file !
...ova_all_samples_combined.fa (pangolin) [jnarayan@hn1 FASTA]$...n1 FASTA]$ pangolin --update pangolin already latest release...equences detected: 320 Total passing QC: 293 Data files fo...1 1 pangolearn 1...946 days ago
Extract fasta header with ids !
...txt minimap2 -t 36 -k19 -w5 -A1 -B2 -O3,13 -E2,1 -s200 -z200 -N50 --min-occ-floor=100 finaal_output.fasta finaal_output.fasta > finaal_self_align.paf936 days ago
921 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, '917 days ago
909 days ago
897 days ago
Installing SEVA environment in Conda !
...Downloading and Extracting Packages libtiff-4.3.0 |...#################### | 100% pandoc-2.17 | 12.7 MB...############ | 100% r-colorspace-2.0_1 | 2.5 MB | #####...############### | 100% r-rcpparmadillo-0.10 | 1.2 MB | ##...897 days ago
884 days ago
884 days ago