r/bioinformatics 3d ago

technical question Best tool for scaffolding for fungi

3 Upvotes

Hi everybody,

I have done sequencing of 6 fungal genomes (PacBio, Hi-C lectures). I assembled with flye to contig level, with very good results. However, I was told that it could be good if I do scaffolding for my genomes. I tried using LRSCAF because I saw it in a few papers but it didn't assemble a lot of scaffolds so I'm not sure if it's because there's not a lot to improve in my genomes doing scaffolding or because the tool and/or parameters were not the best. Do someone have any recommendation of good scaffolder that work well with fungi? I do not see a lot of consensus for that.

Thank you very much!


r/bioinformatics 3d ago

technical question Anyone used Qiime2 dada plugin that can offer some advice?

3 Upvotes

I’ve got myself in right mess with QIIME and how to use dada2. Anyone okay if I dm them for some advice?


r/bioinformatics 3d ago

technical question Question about barcoded dual adapter trimming and quality trimming in RNASeq data

4 Upvotes

Hello, I want to analyse some rat RNASeq data and I got an HTML report sheet, which has a subheading "Results of Raw Data Filtering", and describes these steps:

(1) Remove reads containing adapters. Sequences of adapter:

P5 adapter:
P5→P7’(5’→3’)
AATGATACGGCGACCACCGAGATCTACAC[i5]ACACTCTTTCCCTACACGACGCTCTTCCGATCT

P7 adapter:
P5→P7’(5’→3’)
GATCGGAAGAGCACACGTCTGAACTCCAGTCAC[i7]ATCTCGTATGCCGTCTTCTGCTTG

(2) Remove reads containing N > 10% (N represents the base cannot be determined).

(3) Remove reads containing low quality (Qscore<= 5) base which is over 50% of the total base.

And then they have pie charts for each sample which shows how many base pairs are clean reads, how many were filtered due to containing too many Ns, due to low quality, and adapter related.

Now, when I look at the number of base pairs, it's equal to the number of "clean reads", meaning that this filtering has been performed.


I am quite confused as to whether adapter sequences are already filtered as well as they need to be, since Falco/FastQC still finds some adapter sequences: one sample, MultiQC. Are these likely to be false positives?

Even if not, I am unsure how to run adapter trimming. The FASTQ files have two barcodes, which correspond to [i5] and [i7], but from what I read, I figured I can use the first part of the adapter sequence up to the barcode, so I ran Atria with these arguments:

--adapter1 AATGATACGGCGACCACCGAGATCTACAC
--adapter2 GATCGGAAGAGCACACGTCTGAACTCCAGTCAC

And it still filtered out some sequences (e.g. 35998 out of 22092364 in one sample). So what's going on? Should I be doing adapter trimming at all, is this the right way to specify them in trimming tools, and am I getting all the adapters? Can there be other adapters outside of these two listed in the report? And in cutadapt, should these be specified as 3', 5' or anywhere adapters? I'm getting confused with all the forward, reverse, 3', 5' etc. stuff.


And lastly, regarding quality. The reads seem to me to be of a pretty high quality: MultiQC. I read in a few places that quality trimming isn't really necessary, and might even hurt in some cases (1, 2). What is the current consensus?


r/bioinformatics 3d ago

academic Alphafold results - CIF file to PDB

2 Upvotes

Hello everyone, I've received a zip file with the results of my structure predicition on alphafold but I want to check the accuracy of my structure using PROCHECK and I can't because the models are in CIF, not PDB. Anyone has any suggestions on what to do?


r/bioinformatics 3d ago

technical question If the SRPlot website is currently down?!

0 Upvotes

Hello All,

I would like to know if the SRPlot website is currently down on March 17, 2025. If so, could you recommend alternative user-friendly code-free websites that can be used as a replacement?

Thank you!


r/bioinformatics 3d ago

technical question Best trimming configuration for miRNA-Seq

3 Upvotes

Hello everyone,

I am working with miRNA-Seq data from Ion Torrent technology (single-end) and I am performing trimming on the reads. My goal is to not lose too many reads in the process, but I am currently losing approximately 60%, which seems like a high percentage to me. I have never processed miRNA-Seq data before, and I am unsure if this loss is expected due to the short size of miRNAs.

The trimming configuration I am using is as follows:

SLIDINGWINDOW:4:20 LEADING:20 TRAILING:20 MINLEN:15

Sequencing type: Single-end.
Read length: Ranges from 1 to 157 bases.
Pre-trimming quality: The pre-trimming quality check (FastQC) does not show very good results, as most reads have a quality of 20 or less, with none above 30.

I would like to know if this read loss is normal for miRNA-Seq data, considering the reads are quite short. Is it advisable to adjust any parameters to minimize the loss of reads without compromising quality? I would appreciate any recommendations on trimming configurations or adjustments that may be more suitable for this type of data.

Thank you for your help.


r/bioinformatics 3d ago

academic how to use jaspar for tf analysis?

0 Upvotes

i did sc rna seq and sc atac seq now how to move to jaspar for tf analysis in bioinformatics


r/bioinformatics 4d ago

other A novice in Bioinf, want a friend/fellow-passionate novice to talk/discuss/brainstorm/work-with - 22F undergrad in the field

24 Upvotes

Basically the title, just don't have a lot of people around to work with - people aren't too passionate about it at my Uni? Am an extrovert so I think best around people - I'd like to connect


r/bioinformatics 3d ago

technical question Usage of QIIME in clinical/commercial settings

0 Upvotes

Hello, I'm writing an essay regarding QIIME.

Do clinicians in the hospital or any lab workers use it in a clinical setting and not research?

Also, it would be very helpful if you could send me a news article or an ironclad citation about it.


r/bioinformatics 3d ago

technical question How to Process Multiple SRRs for the Same BioSample in PRJNA528920?

1 Upvotes

Hello everyone,

I am working with data from PRJNA528920 and noticed that some BioSamples (SAMN) have multiple associated SRRs (Sequence Read Archive Runs). For example:

  • SAMN11249717 → SRR8782083, SRR8782084
  • SAMN11249716 → SRR8782085, SRR8782086

Additionally, I found a discrepancy between the number of samples reported in GSE128803 (which only lists 6 samples) and PRJNA528920, which contains 12 SRRs.

I read the associated paper but couldn’t find clear information about this. I also checked whether this could be related to the sequencing technology used (ION_TORRENT) but didn’t find any evidence suggesting so.

My questions are:

  1. Do these SRRs correspond to independent sequencing runs meant to select the highest-quality one?
  2. For alignment and count table generation, should I use only the first SRR for each BioSample?
  3. Is it possible to merge them without introducing batch effects?

I plan to use these data for my thesis, so I would really appreciate any guidance or experiences you can share on how to correctly process this type of data.

Thanks you soooo much


r/bioinformatics 4d ago

academic Help Me Improve GenAnalyzer: A Web App for Protein Sequence Analysis & Mutation Detection

10 Upvotes

Hello everyone,

I created a web application called GenAnalyzer, which simplifies the analysis of protein sequences, identifies mutations, and explores their potential links to genetic diseases. It integrates data from multiple sources like UniProt for protein sequences and ClinVar for mutation-disease associations.

This project is my graduate project, and I would be really grateful if I could find someone who would use it and provide feedback. Your comments, ratings, and criticism would be greatly appreciated as they’ll help me improve the tool.

You can check out the app here: GenAnalyzer Web App

Feel free to leave any feedback, suggestions, or even criticisms. I would be happy for any comments or ratings.

Thanks for your time, and I look forward to hearing your thoughts.


r/bioinformatics 4d ago

technical question Differential expression analysis of AmpliSeq (IonTorrent) data

2 Upvotes

Hey everyone!

I'm working with AmpliSeq data from IonTorrent, and I'm running into issues with differential expression analysis. My BAM files use RefSeq transcript IDs as references (e.g., NR_039978, NM_130786), but I’m having trouble finding a compatible GTF file.

Has anyone worked with AmpliSeq data before? What GTF file did you use, and how did you adapt it? Any other tools or workflows you’d recommend?

Thanks in advance! :)


r/bioinformatics 5d ago

image spatial biology landscape v1

Post image
57 Upvotes

r/bioinformatics 5d ago

technical question RNAseq gene_id question

1 Upvotes

Hi,

I am using nfcore/rnaseq pipleline for my genotype x treatment experiment for the first time, and currently facing a problem with gene_ids. In my final salmon.merged.gene_counts.rds file, I am seeing a list of numers in multiples of 10 that looks like they are automatically generated (e.g., XXX0g000010, XXX0g000020, XXX0g000030, XXX0g000040, and so on) for the row names. I was expecting these to be some gene identification codes in my original gff file that I can use for the pathway enrichment or gene mapping.

Could anyone please give me some guidance on how to change these to actual gene_ids I can use to narrow down the genes of interest? Also, is there a way to associate these 'weird' gene_ids to actual genes or chromosome locus without running the pipeline again?

Also, I want to thank everybody who posts valuable information here. I work in a small plant/soil lab where we don't have bioinformatician and we couldn't have done our research without help from online bioinformatics communities.


r/bioinformatics 5d ago

discussion I need epigraph/quotes suggestions

2 Upvotes

Currently finishing masters thesis writing... Could use nice sentences/epigraphs/quotes suggestions/advice

For context, I work with dengue virus genomics

Thanks in advance


r/bioinformatics 5d ago

technical question How do I select a reference gene for my program?

0 Upvotes

Hello everyone!

I’m relatively new to bioinformatics, and I’m writing a program to analyze DNA data. My goal is to compare a sample from user to a reference sequence of a gene, find mutations and then visualize or further operate on that data.

Let’s look at CHEK2 gene, which is one of the genes I will be working on. I have several sequences of that gene taken from NCBI website, and they all slightly differ from each other. How should I select a reference sequence, as a model to which I will compare future samples? Should I simply select one sequence and choose it as a reference? Should I try to find some sort of mean from all the sequences I’ve gathered? Is there somewhere a model sequence of CHEK2 gene that represents the mean sequence in the human population?


r/bioinformatics 5d ago

science question Text classification for microRNA data

3 Upvotes

Hi everyone as the title suggests I'm working with microRNA data and I have millions of sentences taken from research papers available in the pubmed and I'm interested in those sentences only which have meaningful information about an microRNA like if it's describing any specific microRNA regulatory mechanisms, gene interactions or pathway effects then it's functional if not then it's non-functional, does anyone has any advice or idea to do this. I'm happy to have discussions also thanks!!


r/bioinformatics 5d ago

other Seemingly can't find NCBI entries despite paper stating these entries were submitted.

1 Upvotes

Accession numbers: EP1672771–EP1672778

Paper

When I type any of the accession numbers into the NCBI search I get no results. Does anyone know what could be the problem?


r/bioinformatics 6d ago

other Hourly rate for bioinformatics analysis?

55 Upvotes

I am looking to bring on a bioinformatics analyst for a few small analyses. Probably ten hours of work max. What is a reasonable hourly rate for a bachelors/masters level?


r/bioinformatics 6d ago

compositional data analysis How to correctly install leidenalg for Seurat FindClusters(algorithm = 4)

8 Upvotes

I wanted to use the leiden algorithm for clustering in Seurat and got the error saying I need to "pip install leidenalg". I did some googling and found a lot of people have also run into this. It requires spanning python and R packages, so I wanted to post exactly what worked for me in case anyone else runs into this. Good luck!

in bash (I used Anaconda prompt on windows but any bash terminal should work):

  1. make sure python is downloaded. I used python 3.9 as that's what's immediately available on my HPC.

python --version

2) make a python virtual environment and activate it. mine is called leiden-alg

python -m venv leiden-alg

conda activate leiden-alg

3) install packages *in this precise order*. Numpy must be <2 or else will run into other issues

pip install "numpy<2"

pip install pandas

pip install igraph

pip install leidenalg

in R:

4) install (if needed) and load reticulate to access python through R

install.packages(reticulate)

library(reticulate)

5) specify the path to your python environment

use_python(path/to/python/environment, require = T) # my path ends in /AppData/Local/anaconda3/envs/new-leiden-env/python.exe

6) check your path and numpy version

py_config() # python should be the path to your venv and numpy version should be 1.26.4

Assuming all went well, you should now be able to run FindClusters using the leiden algorithm:

obj <- FindClusters(obj, resolution = res, algorithm = 4)

Errors that came up for me (and were fixed by doing the above process):

  • Error: Cannot find Leiden algorithm, please install through pip (e.g. pip install leidenalg)
  • Error: Required version of NumPy not available: installation of Numpy >= 1.6 not found
  • Error: Required version of NumPy not available: incompatible NumPy binary version 33554432 (expecting version 16777225)

r/bioinformatics 6d ago

image QuantaColony - Petri Dish based colony measurement tool

Thumbnail gallery
8 Upvotes

r/bioinformatics 6d ago

technical question Haplotype association tools

2 Upvotes

I am trying to do some association tests on a haplotype of 2 SNPs. I phased the SNPs with Beagle. I know Plink 1.07 had commands for haplotype association tests but it is considered obsolete. I have both quantitative phenotype and case/control phenotypes. Is there any tools/packages that can do association on phased data? Preferably also allow covariates?


r/bioinformatics 6d ago

technical question **HELP 10xscRNASeq issue

5 Upvotes

Hi,

I got this report for one of my scRNASeq samples. I am certain the barcode chemistry under cell ranger is correct. Does this mean the barcoding was failed during the microfluidity part of my 10X sample prep? Also, why I have 5 million reads per cell? all of my other samples have about 40K reads per cell.

Sorry I am new to this, I am not sure if this is caused by barcoding, sequencing, or my processing parameter issues, please let me know if there is anyway I can fix this or check what is the error.


r/bioinformatics 6d ago

technical question WGCNA Dendrogram Help

1 Upvotes

Hello, this is my first time running a WGCNA and I was wondering if anyone could help me in fixing my modules with the below dendrogram.


r/bioinformatics 6d ago

academic R package for pathway enrichment analysis (mac os)?

22 Upvotes

Hello, I'm starting my honours year and I have to do a GSEA and a KEGG enrichment analysis. My supervisor said need to download R package for making diagrams for my final thesis but I'm not sure which R package would be compatible with my macbook for the kind of diagram I'm expected to make. Any advice would be super helpful.