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
1
vote
2answers
28 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
37 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
56 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
15 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
69 views

Cellranger gives error

I am trying to run cellranger but I get fastq permission denied error ...
2
votes
2answers
38 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
58 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
394 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
65 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
59 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
0answers
83 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
92 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
53 views

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

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

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

I have the following files: ...
-2
votes
2answers
71 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
39 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
44 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
195 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
104 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
105 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
187 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
69 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
416 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
36 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
64 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
187 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 ...
2
votes
2answers
167 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
3k 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
560 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
229 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
47 views

Joinning columns from different files

I have 48 files in which I have r1 ...
1
vote
2answers
44 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
35 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
56 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
120 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
117 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
54 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
50 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
89 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
129 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
78 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
386 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
641 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
64 views

Pattern mining from a genomic sequence

I need to find the following pattern from a genomic sequence ...
0
votes
1answer
124 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 ...
2
votes
1answer
287 views

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 multiple directories that end with a number ranging from 314 to 727, followed by .3 . For example, 'mgm4761314.3'. Within ...
1
vote
2answers
132 views

MiXCR: only create a single export file for all clonotypes

I am using following command from MiXCR to both align, assemble and export my input files: ...