Questions tagged [trimming]

To cut from the end. Trimming is often done on raw Illumina reads to discard low-quality bases or adapter sequences. Questions related to trimming algorithms, when, how to trim, trimming software.

Filter by
Sorted by
Tagged with
0 votes
1 answer
37 views

Paired end reads with different read lengths

I downloaded some raw RNA-seq reads from NCBI using the sra-toolkit (Accession numbers: SRR2053159-64). The reads were generated using ABI SOLiD, and in colour space. I downloaded them in base space ...
Arkajyoti Banerjee's user avatar
1 vote
0 answers
62 views

Can files with different R1 and R2 lengths be trusted?

I received paired end amplicon sequence from a LAB with different lengths for R1 (320) and R2(280). Should I trust this lab to sequence other samples? Also, I had to do a trimming over the R2 at 220(...
abi's user avatar
  • 11
1 vote
1 answer
214 views

Primer trimming-fasta files

My goal: I want to trim off the primers (Forward : CGAGAAGACCCTRTGRAGCT, Reverse : GTTGGGGYGACCNYGG) from a fasta file with a lot of dna sequences allowing for some (e.g. 3) mismatches (identity). ...
Damianos's user avatar
4 votes
1 answer
47 views

Trim reads 1kb upstream of sequence

I need a quick way to trim multiple reads in a FASTA file. I need to trim everything that is 1kbp upstream of this sequence ...
rimo's user avatar
  • 923
3 votes
1 answer
108 views

Does cutadapt trim trailing N's first and then use max_n to filter reads?

Background I want to trim leading and (likely just) trailing N's from my WES (Illumina NextSeq500) reads with cutadapt (--trim-n)...
Dandelion's user avatar
  • 301
1 vote
0 answers
83 views

Steps for trimming aligned peptides using Trimal

I am trimming the aligned peptide sequences and back translating them to codon using Trimal. Here i am using following cmd's: ...
Kishor Kumar Sarker's user avatar
1 vote
1 answer
871 views

Searching for adapter sequences in FASTQ files - metgenomics

I have recently received some metagenomic data from 16S rrna sequencing. The sequencing company claim to have removed primers, however not adapter sequences. Please note that the files have been ...
h3ab74's user avatar
  • 836
2 votes
1 answer
192 views

When using cutadapt, can I specify an R2 adapter as optional when I have a required R1 adapter?

I'm using cutadapt 3.5 to trim adapters and perform some filtering on paired-end data. Both R1 and R2 sequences have 3' adapters that might be found depending on the sequence length, but R1 also has ...
Jesse's user avatar
  • 947
1 vote
0 answers
185 views

Why is cutadapt discarding some paired end reads?

So, I am using cutadapt to remove two primers from some paired end reads. I have 2 files for each sample, i.e. sample_R1.fastq and sample_R2.fastq. I have primer_1 and primer_2 and both of them can be ...
gabt's user avatar
  • 320
1 vote
1 answer
186 views

fastq file format unknown

I have paired-end fastq files some of which seem to be in a weird format (from a collaborator, not a public database). When I cat the file I get what seem to me to ...
iichelhadi's user avatar
4 votes
1 answer
86 views

Telling grep to treat N as [ATCG]

Okay so I'm using grep to try and get a preview of some trimming operations that are not going as expected.. Lets say that my sequence in the FastQ file is: ATNGCNATCG What I want to do is.. ...
RPINerd's user avatar
  • 51
0 votes
0 answers
76 views

FastqCleaner in R, not recognizing paired-end fastq file?

I've recently downloaded the FastqCLeaner to process paired-end fastq files locally using R (4.0.3). I'm having trouble with the shiny app. It is not uploading my second fastq file, and is repeated ...
Najeha Mohamed's user avatar
1 vote
0 answers
184 views

FastQC and Trimmomatic in Galaxy?

I am independently working on data retrieved SRA database, paired-end data as separate inputs. The proceedure I followed is, After running FastQC using Galaxy, the majority of the modules have failed....
Najeha Mohamed's user avatar
1 vote
1 answer
117 views

Trimmomatic trims most of the reverse strand

I have a paired-end data. It is from small part of human genome not a WGS or WES. I use BWA for alignment and do variant calling. Since I had some false positives I wanted to do trimming. I used ...
Hamza Umut Karakurt's user avatar
0 votes
1 answer
524 views

How to run trimmomatic in HPC

Sorry if this has been asked before, I've done a quick search and I don't think there's an easy explanation for me to understand. I'm really new to bioinformatics/RNASEQ analysis, having only taken ...
user7322's user avatar
1 vote
0 answers
291 views

Calculating the trimming parameter when constructing an additive phylogeny

How do you calculate the trimming parameter when constructing an additive phylogeny? I think it has a relationship with the length of the shortest hanging edge, but I can't figure out how this works.
halloiamhere's user avatar
0 votes
1 answer
2k views

What are the right parameters to trim a small RNA transcriptome with trimmomatic?

I'm having some problems with finding the right parameters to trim my small RNA Illumina reads (51 nt long) with Trimmomatic. Before trimming, one of the samples (21M reads) looks like this: So for ...
LinuxBlanket's user avatar
3 votes
2 answers
1k views

Is there a safe catch-all adapter sequence for trimming?

I would like to trim/mark adapters using trimmomatic or picard MarkIlluminaAdapters from a series of Illumina Paired-End read fastqs. The fastq files may have been done using different kits or ...
init_js's user avatar
  • 319
3 votes
3 answers
331 views

Tools for quality trimming at 5 prime?

I'm looking for a tool that is capable to do quality trimming at 5' end and it is configurable. For example, I can choose the quality trheshold, the read length, etc. Any recommendations? I'm ...
aLbAc's user avatar
  • 196
4 votes
2 answers
3k views

"Sequence Duplication Levels" module still fails after pre-processing Illumina data

I want to ask about why the sequence duplication levels are high after I trimmed by using Trimmomatic? I am using the following Trimmomatic operations: ...
yy97's user avatar
  • 43
1 vote
1 answer
242 views

Using sickle for quality trimming

I am using sickle package from bioconda to trim my reads from the Illumina Miseq. When I tried following code, it failed to run ...
L R Joshi's user avatar
  • 719
0 votes
1 answer
234 views

Demultiplexing and preprocessing for custom single-nucleus Drop-seq data

I am trying to reproduce the preprocessing of paired-end sequencing reads in Lake et al. (ref. 1). First, paired-end reads were filtered out if read 1 had more than four non-T bases in the last ten ...
gc5's user avatar
  • 1,783
8 votes
3 answers
1k views

Why does cutadapt remove low quality bases from the ends of reads only?

I use cutadapt to remove low quality bases from my Illumina reads. The algorithm only removes low quality bases from the end until it reaches a good quality base. If there is a bad quality base beyond ...
Biomagician's user avatar
  • 2,459
7 votes
1 answer
1k views

How to trim adapter sequences from GSE65360 in order to map the reads?

I am trying to map single-cell chromatin accessibility data from doi:10.1038/nature14590, obtained using scATAC-seq, to the reference genome. Example paired-end reads are here. The reads contain ...
firefly's user avatar
  • 73
9 votes
1 answer
2k views

FASTQC overrepresented sequences after trimming

I have a set of RNA-seq samples from different experiments (Single and Paired End, depending on the experiment). I ran FASTQC in all the samples and found overrepresented adapter sequences: I removed ...
plat's user avatar
  • 1,022
7 votes
1 answer
655 views

insert size pre and post trimming

I have a problem here with my rna seq data: Sequencing details: rRNA was removed, followed by cDNA preparation and generation of stranded libraries using the TruSeq Stranded Total RNA Sample Prep Kit. ...
novicebioinforesearcher's user avatar
12 votes
4 answers
955 views

How can I systematically detect unknown barcode/adapter sequences within a set of samples?

I have often downloaded datasets from the SRA where the authors failed to mention which adapters were trimmed during the processing. Local alignments tend to overcome this obstacle, but it feels a ...
story's user avatar
  • 1,573