Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
enrichment profile maria_mari Bioinformatics 9 09-24-2013 08:21 PM
DESeq Plot Error coralgirl Bioinformatics 2 11-29-2011 10:35 AM
FASTX quality plot error vinay052003 Bioinformatics 2 11-18-2011 08:26 AM
Citation for illumina error profile jstjohn Illumina/Solexa 1 08-10-2011 04:24 AM
RNA-Seq: Sequence-specific error profile of Illumina sequencers. Newsbot! Literature Watch 0 05-18-2011 03:30 AM

Thread Tools
Old 09-19-2012, 06:55 AM   #1
Location: Wageningen, Netherlands

Join Date: Sep 2012
Posts: 19
Default error profile plot

I'm doing quality control of my alignments from bwa-sw. I would like to make a error profile plot for that, but samstat does not generate one (with comparision on insertions, deletions and so on). Also for long reads all the plots from samstat are not very readable.

Also I managed to get the graph of hits per contig from sam file via samtools, but I'm also interested if estimating to how many contigs particular read matches (from blast alignments it was even up to 100 contigs to which read was mapped) so I'm curious how does it look for bwa-sw alignment. Does any one have any python script for that?
nanto is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 10:11 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO