How to see maximal read length / average read length using STAR
I am trying to understand the length of the short reads in star.png. In the screenshot, the average is 90, but as I know that’s not how this is defined for illumina. If the average is 90 mapped, the actual length is ~100 with the rest going to adapters. Is there a way to see the maximal length of a read (e.g. something like check the max length of the first 5k of the reads)?
• 863 views
•
link
updated
by
Istvan Albert
•
written
by
shinyjj •
•
written
by
shinyjj •
1 answer
You can run a tool like seqkit stats, or run for a quick solution run something like:
cat data.fq | awk 'NR%4==2 { print length($0) }' | datamash mean 1 median 1
• 0 views
•
link
written
by
Istvan Albert
Log in to answer this question.
More posts like this
-
STAR alignment - high % of reads mapped to multiple loci
written by sophie •
-
STAR aligner can't map too short reads
written by Assa Yeroslaviz
-
How to calculate coverage in short-reads and long-reads in RNA-Seq?
written by shinyjj •
-
trouble with 'too short' reads using STAR. each mate maps fine, but the paired end mapping does not
written by james.alan.gregory •
-
Too many unmapped reads - STAR alignment
written by codsorre •
-
Terrible RNA mapping result by STAR
written by 610225668 •
-
cause of unmapped short with STAR
written by lkianmehr
-
STAR outputs interpretation
written by XBria
-
Long Read Length, yet STAR says many reads too short
written by dec986
-
14 contigs recieved from SMRTportal using HGAP.3
written by komwit.s •
