Question: Pindel error - malformed record detected
0
gravatar for rooksie
12 months ago by
rooksie0
rooksie0 wrote:

Hello,

I have sam2pindel input files that are giving me a malformed record detected error.

My command is:

pindel -f ${refdir}/${refname} -p ${filename}_recal_output4Pindel.txt \
 -c ALL -o ${filename}_recal -T 10 --report_long_insertions true \
 --report_breakpoints true --report_close_mapped_reads true

The output is:

Pindel
Initializing parameters...
Pindel version 0.2.5b5, 20150726.
Loading reference genome ...
Loading reference genome done.
Initializing parameters done.
SearchRegion::SearchRegion
Processing region: NC_006088.4  1   196202544
Chromosome Size: 196202544
NumBoxes: 60007 BoxSize: 6546

Looking at chromosome NC_006088.4 bases 1 to 5000001 of the bed region: chromosome NC_006088.4:1-196202544 
Scanning and processing reads anchored in NC_006088.4
Something wrong with the read name: 
malformed record detected in Sample32_recal_output4Pindel.txt

The pindel input file was generated using sam2pindel because the bam files were from bowtie2 alignment and then GATK was used to realign and recalibrate them.

Here is that command:

samtools view ../5_SNP_Indel/${filename}_recal_reads.bam \
   | ./sam2pindel - ${filename}_recal_Output4Pindel.txt ${length} ${filename} 0

I previously ran pindel using input files generated from bam files immediately after bowtie2 (before processing with GATK), and these files did not have this error. They ran successfully.

I suspect the issue may have arisen due to the realign and/or recalibrating of GATK.

Is there a way to check or correct whatever error(s) this is causing pindel?

Please let me know what additional information/data I can supply.

Thank you, Rooksie

pindel gatk genome • 376 views
ADD COMMENTlink modified 12 months ago • written 12 months ago by rooksie0

Try:

samtools view -h ../5_SNP_Indel/${filename}_recal_reads.bam \
   | ./sam2pindel - ${filename}_recal_Output4Pindel.txt ${length} ${filename} 0

The -h will ensure that the SAM header, which records contig information, is retained. May help.

Kevin

ADD REPLYlink written 4 months ago by Kevin Blighe32k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 962 users visited in the last hour