I have assembled a transcriptome using Ray. How can I have the numbers of reads by contig ?
I took the reads used to generate the assembly and aligned them to the assembled contigs using BWA to get the unmapped reads, wich are the non-assembled reads (singletons). So I have a SAM format output. Maybe can I use samtools to get the reads number by contigs ?
Thank you for your help.