959 days ago
959 days ago
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.fas946 days ago
Run Pango on your multifasta file !
...va_all_samples_combined.fa (pangolin) [jnarayan@hn1 FASTA]$ pangolin --update pangolin already latest release (v3.1.16) pangolearn updated to 2021-11-18 constella...944 days ago
919 days ago
917 days ago
Extract all fasta sequences except ids !
awk 'BEGIN{while((getline0)l[">"$1]=1}/^>/{f=!l[$1]}f' genomic.fna > filtered_without_omi.fasta #extract subseq seqtk subseq omi_ids.fa omi_single_id.txt > omi_si...885 days ago
882 days ago
882 days ago
880 days ago