how to check sequence length distribution of trimmed reads? (fasta). is there any tool to calculate it?
•
written
by
evashan96 •
4 answers
BBMap/BBTools has a readlength.sh tool that you can you to make a readlength histogram of the reads after trimming.
A bit of an old tool but you can also use PRINSEQ. They also have a webversion, but never used it.
website: http://prinseq.sourceforge.net/
Try OSTRFPD: Multifunctional tool for genome-wide short tandem repeat analysis for DNA, transcripts and amino acid sequences with integrated primer designer link: https://github.com/vivekmathema/OSTRFPD
Lawrence et al published FAST: FAST Analysis of Sequences Toolbox (paper | github).
There you have a tool called faslen, which annotates the fasta header with its length. To only get the length, you can run:
faslen my_test.fasta | grep \>
Cheers,
Michael
Log in to answer this question.
-
rna long read length distibution
written by shinyjj •
-
MultiQC output: opposing output
written by S •
-
Tools for read length distribution
written by shinyjj •
-
Calculate distribution of reads over the reference genome
written by pixie@bioinfo
-
Different Sequence Length Distribution between forward and reverse reads
written by Linda •
-
How to generate sequence length distribution from Fasta file
written by 2822462298 •
-
group reads by sequence length
written by genya35 •
-
Partially map reads to a reference genome
written by kspata
-
Extract a specific subsequence from a fasta file
written by PAn •
-
Script To Calculate Length Distribution Of A Collapsed Fasta File
written by liran0921 •
What have you tried so far?
why not run FastQC (again) on your cleaned reads? that report includes a length distribution plot.
EDIT. Oh, ok, you're asking it for fasta format as input, then FastQC might not be the answer, unless you still have the fastq files as well then use those as input