[Bioc-sig-seq] Calculating insert size for paired-end /mate pair samples

Tyler Backman tbackman at ucr.edu
Wed Jan 20 19:16:39 CET 2010


I forgot to mention- this is for use with ShortRead.

Sincerely,
Tyler William H Backman
Bioinformatics Analyst
Institute for Integrative Genome Biology (IIGB)
Department of Botany and Plant Sciences
E-mail: tyler.backman at ucr.edu
1207E Genomics Building
University of California
Riverside, CA 92521

On Jan 20, 2010, at 10:13 AM, Tyler Backman wrote:

> Hello Abhi,
> 
> Here is an example for calculating insert size of GA reads via a bowtie alignment:
> 
> # align with bowtie
> system("bowtie -q reference.fasta -1 "pair1filtered.fastq" -2 "pair2filtered.fastq" -n 3 -e 160 -p 5 -I 0 -X 500 --fr > alignment.bowtie")
> # import alignment
> alignedReads <- readAligned("./", pattern="alignment.bowtie", type="Bowtie")
> # calculate insert sizes
> distance <- position(alignedReads[seq(2,length(alignedReads),by=2)]) - position(alignedReads[seq(1,length(alignedReads),by=2)]) + width(alignedReads)[1]
> # look at stats
> median(distance)
> sd(distance)
> 
> Sincerely,
> Tyler William H Backman
> Bioinformatics Analyst
> Institute for Integrative Genome Biology (IIGB)
> Department of Botany and Plant Sciences
> E-mail: tyler.backman at ucr.edu
> 1207E Genomics Building
> University of California
> Riverside, CA 92521
> 
> On Jan 20, 2010, at 9:52 AM, Pratap, Abhishek wrote:
> 
>> Hi All
>> 
>> Recently it has been quiet here. I hope you all had a good start to 2010.
>> 
>> I wanted to know if there is already a method in ShortRead or some other package to calculate the insert size for pair end or mate pair data from GA and plot it on a graph. I guess the latter part is trivial.
>> 
>> Thanks,
>> -Abhi
>> 
>> 
>> 	[[alternative HTML version deleted]]
>> 
>> _______________________________________________
>> Bioc-sig-sequencing mailing list
>> Bioc-sig-sequencing at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing
> 
> _______________________________________________
> Bioc-sig-sequencing mailing list
> Bioc-sig-sequencing at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing



More information about the Bioc-sig-sequencing mailing list