Convert file.lst to file.fasta

Hi,

I am using the command:

seqtk subseq Virsvall.fasta out.file.unmapped.lst > file.unmapped.fasta

And it's not working. The error is:

[E::stk_subseq] could not read the list of regions in the file 'out.file.unmapped.lst'.

The lst file looks like the following

insert image description here

How can I fix it?

seqtk .lst subseq

• 831 views • link

updated 3.8 years ago by GenoMax

• written 3.8 years ago by Gabriel •

Did you make the list file on a windows machine? Can you make sure the line endings are appropriate for unix by doing dos2unix out.file.unmapped.lst?

Another option to try would be filterbyname.sh from BBMap suite.

• 0 views • link

written 3.8 years ago by GenoMax

0 answers

No answers yet.

More posts like this

Decision making on choosing the based on tblastn result
written 1.2 years ago by srestha •

https://blast.ncbi.nlm.nih.gov/Blast.cgi![ 1. enter image description here ][1] ![enter image description here][2] [1]: /media/images/4ee8785b-92df-45fa-9eb4-df6e84e7 [2]: /media/images/a916ec1b-02a1-43fb-beb2-57e08c4e I have performed tblastn using protein accession number against a …
GEOquery getGEO is normalized data or not?
written 1.8 years ago by tyasird •

I'm fetching microarray data using GEOquery. datasets <- c("GSE32676", "GSE43288", "GSE19279", "GSE19280") gset <- getGEO(gse_id, GSEMatrix = TRUE, AnnotGPL = TRUE) So, when I check …
constructing gene-based pangenome from orthofinder output file
written 1.8 years ago by Bushra Bibi

Hi all, I am constructing gene-based pangenome for fungus. I have used orthofinder. Please suggest which file should be used to find core and dispensible …
gatk-funcotator data source
written 1.9 years ago by ashaneev07 •

Hello all, I am seeking clarification regarding the data sources provided with the GATK Funcotator tool. The README.txt file does not clearly specify whether the …
Using samtools with GCS (google cloud storage) on a docker container seems to give "Protocol Not Su…
written 2.3 years ago by abhishekghadge •

![enter image description here][1]Having an issue trying to use the samtools view gs://mybucket/myfile.bam > myfile.bam when using this command with docker. Always get the error …
Increasing per base G content in QC sequencing files
written 2.5 years ago by robertsr •

Hello, I conducted some metagenomic sequencing as follows: - Metagenomic sequencing from human stool samples - PCR-free library prep using NEB kit (450bp insert) - …
DESEQ2 analysis - PCA plot
written 3.0 years ago by shome •

![enter image description here][1]![enter image description here][2] I am not getting a lot of genes from my DESEQ2 analysis hence I was checking the PCA …
Subseq bam file
written 3.7 years ago by shinyjj •

Hi Biostars, Does anyone know how to subsample read from a bam file? The below command gives the read number of this bam file. I …
RNASeq data has 3 files R1, R2 and R3, which one is the barcode file?
written 4.1 years ago by leranwangcs •

Hi, I received some RNASeq files to do analysis. But I'm very confused as a newbie in RNASeq when I checked the files. There are …
format the gtf file
written 4.4 years ago by bright602 •

Hi there, I have a gtf file. Could some one let me know how to move the third column with "transcript" to the first line …