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

Join Date: Feb 2010
Posts: 245

Hi harrb,

have a look at this older thread about the header;

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


This option (as stated at 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;

Good luck,
boetsie is offline   Reply With Quote