Entering edit mode
4.2 years ago
A_heath
▴
170
Hi all,
I have paired-end reads for a genome and a specific contig of interest. I would like to prolong my contig as much as possible using PRICE. I tried the following command line:
./PriceTI -fp L_raff_FAM23217_raw_reads_1.fq L_raff_FAM23217_raw_reads_2.fq 150 -icf contig.fasta 1 1 5 -o contig_extended.fasta -nc 5 -mpi 95
However, I obtained a very weird output file with additional contigs. Here is an example of an output file obtained after cycle 1:
>contig_1 (14073nt)
TGTAATAGAGGGGGCTTTATCT...
>contig_2 (340nt)
TCATTATTGTCG...
>contig_3 (307nt)
CTCTTCGGGTTTTC...
I don't understand why there are additional contigs. I just want an extension of my original contig (initial size = 13 035 bp).
Is there a problem with my command line?
Thank you in advance for your help!