Questions tagged [bash]

Use this tag to refer to the questions related to the Linux and Unix shell abbreviated as bash (Bourne Again SHell). It has notable programmatic strengths over C shell (csh) and touch C shell (tcsh)

Filter by
Sorted by
Tagged with
0
votes
0answers
32 views

Permutation Function of 4 amino acid sequence on Bash shell

How can I write a Programming Language to construct about 160000 tetrapeptide from the 20Amino acids using permutation on Bash Shell.
0
votes
1answer
474 views

curl: (22) The requested URL returned error: 400 Bad Request when using entrez direct

I've just dowloaded entrez direct on my Mac and I was trying to perform some simple queries like: ...
1
vote
2answers
32 views

Is there a way to calculate the Ka/Ks value in R/unix?

I am investigating the rate of synonymous non synonymous substitutions and would like to find the Ka/Ks value also known as dn and ds between multiple alleles? I want to investigate multiple alleles ...
2
votes
1answer
39 views

Automatize classical debugging in a bash script? (normal debugging by hand)

And thanks in advance. I am debugging a text file which is the one the two inputs on the program RGmatch https://bitbucket.org/pfurio/rgmatch/src/master/. The code works. Is that some records give ...
1
vote
2answers
110 views

Retain only part of a file name and fasta header in fasta directory

I have a MainDirectory/hundreds_of_subdirectories/thousands_of_fastas structure. For each fasta in the bottom level, I am trying to change the file name as well as ...
1
vote
0answers
16 views

Nucleotide Frequency - HOMER?

I'm trying to find the frequency of UA dinucleotide in a bed file. I want to convert the bed file to FASTA, count all the UA dinucleotides up, then divide by the total length of the FASTA file. Also, ...
-1
votes
1answer
174 views

Cellranger gives error

I am trying to run cellranger but I get fastq permission denied error ...
2
votes
2answers
41 views

Calling for one specific SNP from multiple sequencing runs

Just learning the basics of bioinformatics and bash shell scripting. Is there a correct way to search for a specific allele with one or even multiple sequencing runs? I am searching for prevalence of ...
1
vote
1answer
59 views

how to sort and join two files based on First Column Id

I have got two files large files like this in Tab-delimited format, trying to merge in R, do we need to sort the file before merging? ...
5
votes
3answers
714 views

Removing duplicate FASTA sequences based on headers with Bash

I used the following command to remove duplicate FASTA sequences based on the header sequence: ...
-1
votes
2answers
72 views

Filter rows of VCF file for Match=EXACT?

How do I remove rows in a VCF file on 1kGenome column, where Match=EXACT using bash ?
0
votes
2answers
63 views

Print specific columns in a matrix on the basis of sample id's in the header

I have a matrix file (expression.txt) in which the first column is a gene_id and from the second column the sample id's start. This matrix has 20,000 columns with sample id's corresponding to ...
0
votes
1answer
182 views

STAR aligner multiple fastq files

I’m using STAR to align fastq files from SMART-seq2. I have raw data folder containing sub-folders with samples names the sub-folders each contain fastq file. How can I make a bash command in order ...
2
votes
2answers
116 views

How to remove sequences from a fasta file using a sequence ID list which contains a space within the id?

I have a fasta file that contains sequence reads and sequence id file that needed to be removed from the fasta file. I have done this earlier, but since id contains a space my piece of code is not ...
1
vote
3answers
64 views

How to filter out partially similar strings from two lists and make one list?

I have following two lists. List1 ...
2
votes
1answer
67 views

Merging several fq.gz files or R1 and R2 classes into a single one

I have the following files: ...
-2
votes
2answers
79 views

How to remove low abundance and less prevalent data from my dataset?

I have a relative abundance table in .tsv format where samples are in column and rows contain the features (pathways). Something like this reproducible example: ...
1
vote
0answers
45 views

How to calculate GenePy score for vcf files in all subdirectories?

In a parent directory, I have several sub-directories, each of them contain one .vcf.gz file on which analysis should be performed. I want to run following commands on all .vcf.gz files in all ...
1
vote
1answer
48 views

Table error when combining counts columns

I have these files that I want to make table of (https://github.com/learnseq/RNAseqfiles01.git), the table that I want (since all the files have the same ID column),, I want merge the ID column with ...
1
vote
2answers
293 views

Trying to create a .bam file without the need for a .sam file

I'm trying to use the code specified in this link to create a .bam file without the need for a .sam file. Here is the code I'm using: ...
2
votes
2answers
112 views

How to get rows with similar values in two different columns using command line?

I have following example data from a large file. ...
1
vote
2answers
158 views

Mapping statistics from bam file using bbtools and sambamba

Below are the statistics for RNA-seq mapped and unmapped paired-end reads to rice genome using reformat.sh from bbtools on bam files. It gives 77% mapped and 5% unmapped, what about the remaining 18% ...
1
vote
2answers
223 views

How to run Jupyter script on Slurm HPC

Now jupyter installed on the server and I am using below code to plot the rarefaction plot but I am still getting some error. could you please suggest how I can get rid of it? ...
0
votes
1answer
93 views

repophlan script to download bacterial genome

I am trying to donwload the microbial genome using the repophlan_get_microbes.py (https://github.com/SegataLab/repophlan) but now it is running more than 10 days and still not finished on slurm HPC. ...
2
votes
5answers
803 views

How to measure the total size of a fastq file in base pairs?

Or Kbps/Gbps. It feels like it should be conceptually very simple, but I can't seem to figure out the right combination of keywords to find it via my search engine. Help would be appreciated! I have ...
0
votes
0answers
42 views

Issues with AutoDock Vina

I am trying to use AutoDock Vina to do docking however I am getting this error. I would really appreciate any help as I have not been able to resolve this error despite trying numerous times.
0
votes
2answers
91 views

Using variables with fasterq-dump?

I am trying to download multiple fastq files from the SRA NCBI database which, conveniently enough, have their IDs in a range. Building a simple loop doesn't seem to do the trick, as bash appears to ...
1
vote
1answer
313 views

Split fasta file based on groups in header information and output as separate files

I have a fasta file containing the sequence of a gene across different species. In total there are around 900 samples and 12 species. (Each sequences is over multiple lines and longer than 100bp.) My ...
3
votes
2answers
291 views

How to subset genes and its nested features from a GFF file using a gene list

I would like to subset a GFF file (gene and nested features) from a gene list. The GFF file looks like this ...
4
votes
1answer
4k views

Calculating average coverage for .bam files (sequence data)

(Full discolosure that this is my first time working with sequence data, and with the bash scripting.) I need to calculate the average coverage for any .bam file. After some searching I wrote the ...
1
vote
4answers
2k views

Edit FASTA header using sed

I need to rename the following headers from a FASTA file. Something like: ...
1
vote
1answer
728 views

remove sequences from fasta file matching a string in the header

I have a file with 16S sequences. some headers contain species information. For my purposes I would like to exclude a number of species from the file, therefore I would like to do a pattern matching ...
2
votes
0answers
318 views

Cannot blast against specific NCBI databases

I am having issues with some prokaryote reference genome databases (exact names : ref_prok_rep_genomes.*), that I downloaded from the NCBI website : https://ftp.ncbi.nlm.nih.gov/blast/db/. Files in ...
-2
votes
2answers
71 views

Calculating the number of probes for a given genomic range

I have a big data frame like this for copy number (exome seq) ...
-4
votes
1answer
48 views

Joinning columns from different files

I have 48 files in which I have r1 ...
1
vote
2answers
45 views

Trying to divide number in fourth TSV column by a certain number according to its corresponding identifier in the first column

I was wondering if anyone could help with my little issue in manipulating tsv files in order to calculate estimated coverage across many tsv files. I have tsv files that look something like this: <...
1
vote
1answer
40 views

Running tophat dockerfile in background

I need to use tophat2 to align some sequences, and I wish to use the docker container. These are large sequences so it will take a long time, plus I'm working on a university server so the chances of ...
-1
votes
1answer
38 views

Appending some columns from one tabular file to another

I have a list of files in which I have this information ...
2
votes
2answers
59 views

Moving file based on their names

I have a list of vcf files; I also have a list of names in a txt file like LP6005409-DNA_F01 LP2000325-DNA_A01 LP6005409-DNA_E02 LP6005500-DNA_C03 What I have in ...
0
votes
2answers
131 views

Looping over several files in bash

I have a bash script: I am wondering how I can change this script to loop over a bunch of .vcf files and give output .txt with the name of corresponding .vcf I tried changes done in similar script ...
0
votes
2answers
126 views

Changing this code in a way to work for my files

I have a bash script which extracts some information from a .vcf file but one .vcf file at each time. How I can change this script in a way to work on a bunch of .vcf files and the output is a .txt ...
-1
votes
1answer
57 views

bash loop for strelka2 processing for multiple files from one directory

In the below bash loop for strelka2, 3 bam files are being stored in a variable. The two questions are: Is this proper way to loop through multi bam and wait for ...
1
vote
1answer
53 views

Fastq-dump script download X spots or all

Im trying to write a script where an optional input of -X flag can be used, or if that info is not available download all reads. my script as follow: ...
0
votes
2answers
101 views

How I can run this code on my files?

I am annotating some .txt files by Annovar software by this code nnovar]$ module load annovar/2016Feb01 [cyan01 annovar]$ table_annovar.pl But I really got ...
1
vote
2answers
131 views

How can I run a command for multiple files?

I have a bunch of .vcf files in a folder and I want to run command below on all of them but doing that one by one manually is really painful. I am seeking for a way ...
-1
votes
1answer
79 views

How I move these files based on their names in R or terminal

EDITED I have list of my samples in first column and corresponding file name for my .vcf files in second column of a file name clin_name.txt like below. For ...
0
votes
2answers
461 views

How to loop multiple function in shell script?

I need to extract sequences one after another consecutively from a large fasta files (multiple fasta files) and each extracted files to be saved in new fasta file (I mean the first sequence extracted ...
2
votes
2answers
884 views

Error with BWA Mem input having multiple fastq files using cat and process substitution

Running BWA Mem on a number of paired end fastq files using process substitution on the inputs results in this error: ...
2
votes
1answer
68 views

Pattern mining from a genomic sequence

I need to find the following pattern from a genomic sequence ...
0
votes
1answer
154 views

How to edit the headers of multiple fasta files from multiple folders? (recursively)

This is an expansion of my previous question, How to edit the headers of multiple fasta files from multiple folders? My directories are organized as follows: one main directory, in which I have ...