More posts like this
-
Genome Visualization Tools
written by dlera.lozano •Hi Everyone, I was wondering if anyone was familiar with any whole genome visualization tools, I have been using BRIG to visualize some bacterial genomes …
-
Get read counts without GTF/GFF files
written by weather •Hi, I was running some RNAseq for my non-model species, which only have reference genome (FASTA) available in NCBI. Therefore, I was able to run …
-
Import or read Snapgene .dna files into R
written by eggrandio •Hello, I was wondering if there is any package that allows reading Snapgene .dna files directly into R. Even if it does not conserve the …
-
Is it possible to use pseudoalignments.bam to run rMATs?
written by Lillian •Hi I have an urgent question about running rMATs, I usually running rMATs after running STAR alignment with original fastq files. However, I'm unable to …
-
Getting unspliced counts for scRNA-seq data generated from cellranger aggr pipeline
written by Luna_P •Hi everyone, I have 15 scRNA-seq samples that I have aggregated with cellranger aggr pipeline. My problem is that, I would like to use [scVelo][1] …
-
[WGS] how to correctly visualize long indels (50nt+)
written by memento1984_1 •Hi All, I was wondering if there is any tool that can correctly visualize long indels (eg known deletions that are 50bp+ long). Bam files …
-
Converting Bisulfite-seq BAM files to fastq
written by divyswar01 •TCGA has some DNA Methylation files here https://portal.gdc.cancer.gov/legacy-archive/search/f. But they only have BAM files here. I need fastq files for these BAM files. I have …
-
Salmon and reads order in fastq in quasi-mapping-based mode
written by ZheFrenchReally a naïve question but it's a little while I have this one in mind. From Salmon doc : *"If your reads or alignments do …
-
How do I remove files that only contain a header OR get tabix not to output a file when the input d…
written by kyle •Hi, I am new to bioinformatics so I apologise if this is an obvious problem that I have missed. I could not find a similar …
-
how to calculate the coverage over a specif coordinate
written by Lila MHi everybody, I have bam files for RNAseq, and I know the coordinates for specific exons of interest. I would like to know if there …
You can export fastq data from your BAM's, trim and realign.
If you don't want to do that then this may be an option: Trim bam file reads to remove all nucleotide over a certain length
How to trim alignments in a BAM file to a restricted genomic window?
What are you going to base your trimming on?
In addition to GenoMax 's answer, I'll like to mention another tool called iVar.
I am basing my trimming on base quality. I found a tool that does trimming of bam files directly (BamUtil, trimBAM) but you have to manually enter how many bases you want to trim from each end. Is there anything like fastp that does the trimming automatically based on quality?
Also, iVar seems like it's for viral amplicon-based sequencing, which I don't think would be appropriate here.