I am trying to map the raw reads of a genome back to its assembled genome to generate the pileup. I downloaded the raw reads (fasta and qual files) from NCBI Trace Archive. Because they are long reads, I used BWA to align them against the assembled genome. I used samtools to generate the pileup. Everything was done with default values. The problem I am having is that the assembly I got has a 20X coverage, while the reported coverage in the original paper was ~12x. I am really puzzled by this big difference. Anyone can help? Thanks very much!
Question: Mapping Reads Back To Assembly Contigs To Generate Pileup
6.9 years ago by
Sen • 0
Sen • 0 wrote:
ADD COMMENT • link •
Please log in to add an answer.
Powered by Biostar version 2.3.0
Traffic: 730 users visited in the last hour