Question: Preprocessing for GATK
0
gravatar for anailg92
3.6 years ago by
anailg920
Spain
anailg920 wrote:

Hi, I'm doing a GWAs study in Arabidopsis thaliana, I'm using for preprocessing the samples de command bwa aln and sampe, then samtools import reference.fa sample.sam sample.bam.

But when the script continous, it has an error and I don't know how to fix. This is my error:

[samopen] SAM header is present: 7 sequences.
[sam_read1] reference 'ID:bwa   PN:bwa  VN:0.7.10-r789  CL:/home/bwa-0.7.10/

bwa sampe ./thaliana.fa SRR519713_1.sai SRR519713_2.sai SRR519713_1.fastq SRR519713_2.fastq
' is recognized as '*'.
[main_samview] truncated file.

If anyone knows where I am wrong?

And another question which reference genome I have to put in this commandline: 

java -Xmx4g -jar $HOME/opt/GATK2/GenomeAnalysisTK.jar -T RealignerTargetCreator -R $2 -I $3.bam -o $3

Thank you!

snp • 1.1k views
ADD COMMENTlink modified 3.5 years ago by Biostar ♦♦ 20 • written 3.6 years ago by anailg920
1

What command are you giving that produces that error message? It looks like the header is completely messed up. BTW, use samtools view instead of samtools import.
 

ADD REPLYlink written 3.6 years ago by Devon Ryan84k
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: 618 users visited in the last hour