This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Calculate Distribution Of Mapped Reads To Transcriptome Assembly

Hi,

I have a de novo transcriptome assembly done with trinity and afterwards I mapped the reads to it using bowtie2. Does anyone knows of a tool/script to calculate the distribution of the reads in the contigs using the information in the sam file? Basically I would like to know if on average reads map to the beggining or the end of the contigs (~10⁶ contigs)

Thanks

bowtie2 trinity assembly transcriptome
• 3,800 views • link
updated by hbw hbw • written by Biojl Biojl

1 answer

You could use HTSeq to read in the SAM file and then just count the number of reads mapping to each location. You could then accumulate the reads mapped to each percent length of the transcripts. If you don't mind asking, why would you expect to find more reads in the 3' or 5' end?

• 345 views • link
written by hbw hbw

For example low RIN scores are supposed to cause that.

• 0 views • link
written by Biojl Biojl

Log in to answer this question.

More posts like this
  • Sam To Bam
    written by freddy • freddy
    <p>I sequenced a mito genome using paired end method on the mi-seq I mapped all the reads back to reference gnome using <a href='http://bio-bwa.sourceforge.net/'>BWA</a> and …
  • How To Sort Sam File
    written by Chen Sun Chen Sun
    <p>I want to process the <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file created by <a href='http://bio-bwa.sourceforge.net/'>BWA</a>, if it is sorted by the sequence coordinate on reference, it will be …
  • Sam To Bam - Loss Of Data Or Just Great Compression?
    written by Clare Clare
    <p>I use picard tools to take a <a href='http://bio-bwa.sourceforge.net/'>BWA</a> alignment <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file and convert it into a sorted <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file. Normally this works …
  • Extract Raw Counts From Sam File
    written by mad.cichlids mad.cichlids
    <p>Hi, All</p> <p>I was performing a de novo assembly of RNA seq data. I aligned my reads to the assembly for each of my samples …
  • Error In Converting Sam To Bam By Samtools
    written by Raghav Raghav
    <p>hello every one,</p> <p>after mapping,i encountered an error during converting my <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> out file into <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> out put file.</p> <p>[cdac@nbri <a href='http://samtools.sourceforge.net/'>samtools</a>-0.1.19]$ <a …
  • Filtering A Sam File For Quality Scores
    written by Varun Gupta Varun Gupta
    <p>Hi I have a <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file. The 5th column of the <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file tells us the quality score. In my <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> alignment …
  • In Bowtie2'S Sam Output, What Does The Yt Attribute Mean?
    written by Eric Fournier Eric Fournier
    <p>In Bowtie2's <a href='http://samtools.sourceforge.net/SAM1.pdf'>SAM</a> output, what does the YT flag mean? I can't seem to find a mention of it in Bowtie2's documentation, and the …
  • Soap Result Converted To Bam Not Working In Samtool
    written by Sarwar • Sarwar
    <p>I converted the output of soap2 which is map file to <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file using map2sam.pl and then converted this <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> to <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> …
  • Reads Per Contig From Sam File?
    written by Newvin • Newvin
    <p>After recently assembling a transcriptome, I took the reads used to generate the assembly and aligned them to the assembled contigs using <a href='http://bio-bwa.sourceforge.net/'>BWA</a>. The …
  • Using Gatk, Removing Read Contigs On Bam Level
    written by Mdeng Mdeng
    <p>Hey everyone,</p> <p>ATM I am working with GATK for the first time and a get an error, wich everybody knows:</p> <pre><code>##### ERROR MESSAGE: Input files …