View Single Post
Old 12-21-2010, 02:01 PM   #2
boetsie
Senior Member
 
Location: NL, Leiden

Join Date: Feb 2010
Posts: 245
Default

Hi harrb,

have a look at this older thread about the header;
http://seqanswers.com/forums/showthr...abyss+coverage

About the coverage (per position on the contig i think you mean?) you can probably run ABYSS with the option

--coverage-hist=FILE

This option (as stated at http://seqanswers.com/wiki/ABySS) records the k-mer coverage histogram in FILE. Though, I'm not 100% sure, since i can not test it @tm .

Otherwise, you can map your reads to the contigs (Bowtie, BWA), producing a .sam file. Convert the .sam file to .bam file and do;

samtools pileup -f your_ref.fa your.bam

Have a look at this thread;

http://seqanswers.com/forums/archive...hp/t-7679.html

Good luck,
Boetsie
boetsie is offline   Reply With Quote