Questions tagged [filtering]

The tag has no usage guidance.

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

Filtering out nonarthropod sequences using Blobtools2

I am using Blobtools2 and am working my way through the sample dataset that they provide before using my own. The filtering tutorial I am following is --> https://blobtoolkit.genomehubs.org/...
0
votes
2answers
38 views

Extract sequences from partial Header

I have a library I would like to split them based on their headers. Some are related to RNA and DNA. The header contains much information, but the most important the presence of DNA or RNA/LTR.., ...
0
votes
0answers
15 views

Select for a specific proportion of a genomic feature

Currently have a bed file with exon coordinates for a list of about 20 genes, and I want to select the exons that fall within the first 3/4 or 2/3 of the coding region for each gene. Does anyone know ...
0
votes
1answer
43 views

Filtering raw sequencing reads

I have to filter the raw sequencing reads based on the following criteria: Remove reads containing adapters Remove reads containing N > 10% (N represents base that could not be determined) Remove ...
1
vote
3answers
141 views

Edit all the fasta headers using awk

kindly your comment highly appreciated have more than 4k header sequences look like: ...
1
vote
1answer
48 views

How do I create a filtered gene list using expression medians

Forgive the simple noob question I have TPM data of ~50k genes (rows) across ~1k cell lines (columns). In R, I would like to output an "intermediate expression" gene list for each cell line, like: <...
1
vote
0answers
76 views

Empty values in FILTER field after running Mutect2

I was expecting something like in this image: But when I ran Mutect2, the FILTER field was empty in the output VCF file. How can I enable PASS filters? ...
3
votes
1answer
180 views

Fast processing of fastq data

I am trying to write python script for customized filtering for fastq file (size >3 GB). My proposed script is as follows: ...
1
vote
0answers
35 views

ExpressionSet with a custom gene list

After an ExpressionSet is established, I want it to contain only the genes of my choice ("selected.txt") What is the best approach?? Thanks
3
votes
1answer
200 views

Output from vcftools missingness

I'm new to data filtering on vcf data and vcftools. I performed variant calling on my dataset, CHR22, homo sapiens. I'd like to remove sites that are missing in more than 5% of individuals. ...
2
votes
2answers
428 views

How can I extract the longest N isoforms per gene from a fasta file?

The question of how one can extract the longest isoform per gene from a multi-fasta file has been answered expertly in this thread: How can longest isoforms (per gene) be extracted from a FASTA file? ...
6
votes
5answers
841 views

Delete all 4 lines of a fastq read from a fastq file using read ID

I have the following error when running bowtie2: Error: Read HWI-D00466:116:CC62WANXX:3:1102:7363:63646 1:N:0:GCACACG has more read characters than quality values. I now want to remove all 4 ...
1
vote
0answers
214 views

Error filtering SNPs with Plink

Im trying to filter SNPs using PLINK with this command. Also, the .map and .ped files were converted to a plink format from a vcf file (...
2
votes
1answer
234 views

Select top 100 genes ranked by variance in read counts

I have a gene expression data (read counts) table with many samples. Showing some of them below: ...
3
votes
1answer
141 views

How to filter Ensembl cDNA and ncRNA FASTA files by primary assembly?

I'm currently performing differential expression analysis using alignment-free quantification using Kallisto. To do this, I need to create a Kallisto index using Ensembl's cDNA and ncRNA annotations ...
5
votes
1answer
404 views

Single-cell RNA sequencing (scRNA-seq): filtering cells by transcript counts, how to choose cutoffs?

I am running a notebook with example for the MAGIC algorithm. In the data preprocessing step, a filtering operation is required to filter out cells with a small count of transcripts. My question is ...
5
votes
1answer
268 views

Filter out outliers of the scRNA-seq (heterogenous cells)

I am new to data science. I have a dataset of single-cell gene expression from multiple cell types in C. Elegans. The dataset is from the paper Comprehensive single-cell transcriptional profiling of a ...
11
votes
1answer
1k views

Why does a very strong BLAST hit get lost when I change num_alignments, num_descriptions or max_target_seqs parameter?

Disclaimer: This is a self answered question for documentation purpose and I adapted this from the following github gist. Especially from users terrycojones and peterjc as well as sujaikumar who ...
9
votes
6answers
4k views

How do I find identical sequences in a FASTA file?

I want to create a database for a proteomics study. Therefore, the mapping from a given sequence to a protein ID has to be unique. I am wondering whether there is already a built-in function in ...
1
vote
1answer
224 views

Filtering out all seqs with mutations of list2 from list1

I have 2 lists, list1 and list2, of protein sequences of the same given gene in different strains. In ...
2
votes
3answers
751 views

Filter BAM file for read pairs where one or both of the reads starts with a given sequence pattern

What's the most straightforward way to filter a BAM file for read pairs where one or both of the reads starts with a given sequence pattern? Would a tool that deals with UMIs work for this? Which one ...
6
votes
4answers
400 views

Any fast options to query large VCF bed intervals?

I'm doing some analysis and I need to subset a large VCF file (~8GB gziped) given a bed interval and identify within a subset of rsid. Unfortunately, both my normal choices to do this analysis (<...
7
votes
2answers
2k views

Filtering step for read counts data

I have around 1200 samples as columns and 60,000 genes with Htseq-Counts data. Before normalization with voom function I want to do filtering step. I want to remove genes whose expression is == 0 in ...
6
votes
2answers
647 views

How to safely and efficiently convert subset of bam to fastq?

Question How can I extract reads from a bam file (produced by bwa-mem) to fastq given a ...
9
votes
5answers
2k views

How can longest isoforms (per gene) be extracted from a FASTA file?

Is there a convenient way to extract the longest isoforms from a transcriptome fasta file? I had found some scripts on biostars but none are functional and I'm having difficulty getting them to work. ...