Hi,
I am trying to use MECAT2 to correct PacBio reads. I used the command mecat.pl correct for this, which according to your documentation applies mecat2pw and mecat2cns. The output file cns_reads.fasta, has 1,639,674 lines, but the input file had 4,964,062 lines. Are reads supposed to be dropped like this? Does it have to do with the coverage cutoff? I used the defaults for the contig file except for genome size and threads.
Thanks for providing this software for us, and I appreciate any help you can provide when you have time for it.
Hi,
I am trying to use MECAT2 to correct PacBio reads. I used the command mecat.pl correct for this, which according to your documentation applies mecat2pw and mecat2cns. The output file cns_reads.fasta, has 1,639,674 lines, but the input file had 4,964,062 lines. Are reads supposed to be dropped like this? Does it have to do with the coverage cutoff? I used the defaults for the contig file except for genome size and threads.
Thanks for providing this software for us, and I appreciate any help you can provide when you have time for it.