Hey guys, I've been struggling with this question for a while, since I don't program very well.
So, I've mapped some 454 contigs to a very fragmented draft genome (a lot of scaffolds), and I would like to access the quality of this draft by:
1-) Computing how many 454 contigs have mapped to different scaffolds (to investigate how fragmented the draft is)
2-) Computing the alignment length and identity of each aligned 454 contig..
I've mapped with BWA and I have sam, bam and index files.
Can I find these infos with samtools? Or should I use something else?
Thank you so much!
So, I've mapped some 454 contigs to a very fragmented draft genome (a lot of scaffolds), and I would like to access the quality of this draft by:
1-) Computing how many 454 contigs have mapped to different scaffolds (to investigate how fragmented the draft is)
2-) Computing the alignment length and identity of each aligned 454 contig..
I've mapped with BWA and I have sam, bam and index files.
Can I find these infos with samtools? Or should I use something else?
Thank you so much!
Comment