Question: Read length from HT-Seq
0
gravatar for gupta
2.4 years ago by
gupta • 30
gupta • 30 wrote:

Hi,

Sorry for asking this very basic question. I have paired-end RNA-Seq data. After initial adapter and quality trimming, I used Tophat for alignment. Later, I obtained the read count using HT-Seq from the alignment file.

My questions is: Can I also somehow obtain the length of the reads that were counted/filtered by HT-Seq? This would allow me to see the distribution of read length in the data.

Thanks in advance for any help.

Best, Pooja

ht-seq rna-seq • 683 views
ADD COMMENT • link • modified 2.4 years ago by Jennifer Hillman Jackson ♦ 25k • written 2.4 years ago by gupta • 30
2
gravatar for Jennifer Hillman Jackson
2.4 years ago by
United States
Jennifer Hillman Jackson ♦ 25k wrote:

Hello Pooja,

The htseq_count tool form option "Additional BAM Output" can be used to produce a dataset that can be parsed for both feature and read length. Search for the keyword "bam" in the tool panel to review the tools to perform the summary (are several choices). Or, convert the BAM-to-SAM, then parse the SAM file as a tabular dataset.

Hopefully this helps! Jen, Galaxy team

ADD COMMENT • link written 2.4 years ago by Jennifer Hillman Jackson ♦ 25k

Thank you so much Jen. Its really helpful :)

ADD REPLY • link written 2.4 years ago by gupta • 30
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 167 users visited in the last hour