Perl script to extract sequence by Ids from multifasta file !
#!/usr/bin/perl -w use strict; my $idsfile = "$ARGV[0]"; my $seqfile = "$ARGV[1]"; my %ids = (); open FILE, $idsfile; while() { chomp; $ids{$_} += 1; } close FILE; local $/ = "\n>";...as first word in FASTA header if (exists($ids{$id...2106 days ago
2086 days ago
1736 days ago
Samtools commands for bioinformatician !
## count mapped reads samtools view -c -F 260 mapping_file.bam ### converting sam file into fasta samto...reads_mapped.sam > reads.fasta ### converting sam file into bam # -b :...######## # try these recipies: ### convert BAM to FASTA: samtools view...1613 days ago
1563 days ago
To convert just one specific read group to fastq
# Stop script on error. set -uex # The SRR BioProject number for the sequencing...roups to paired files. samtools fastq -t -1 all1.fq -2 all2.fq all.bam # To convert just one specific read group....1553 days ago
1413 days ago
Sequence Ids conversion files !
ftp://ftp.ncbi.nlm.nih.gov/gene/DATA/ Name Size Date Modified ARCHIVE/ 02/01/2020, 05:30:00 ASN_BINARY/ 03/07/2020, 07:49:00 GENE_INFO/ 03/07/2020, 07:48:00 0...1413 days ago
Bash script to handle Multifasta files
#Convert all lowercase residues to uppercase in a FASTA sequence file $ awk 'BEGIN{FS=" "}{if(!/>/){print toupper($0)}else{print $1}}' input.fasta > output.fasta #...1368 days ago
1226 days ago