Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Convert SRA to FASTQ with fastq-dump but problem of read length MurielGB Illumina/Solexa 9 10-07-2014 09:10 AM
Chip-seq tags of unequal length - is this a problem for macs? feralBiologist Bioinformatics 1 01-16-2014 08:24 AM
tophat seq length qual length mismatch error epi Bioinformatics 3 09-10-2012 01:43 AM
read length a problem when mapping with MAQ koustavpal Bioinformatics 3 08-28-2012 06:05 AM
beginning of quality values record not found & qual length () differs from seq length papori Bioinformatics 1 07-09-2012 12:25 PM

Thread Tools
Old 09-26-2015, 03:18 AM   #1
Waqasuddin Khan
Junior Member
Location: Karachi

Join Date: Nov 2013
Posts: 7
Default chr20 length problem

I have an Abi SOLiD WGS csfasta reads. After QC and filtering steps, I aligned the reads via SHRiMP2, specific for csfasta reads. I used hg19 as reference downloaded from SAM file is converted to BAM to sorted BAM to markduplicates (remove duplicates false) to reordersam to validatesam by PICARD. Than, I used GATK for realignertargetcreator but I got the followed error:
ERROR MESSAGE: Input files known and reference have incompatible contigs: Found contigs with the same name but different lengths:
##### ERROR contig known = chr20 / 64444167
##### ERROR contig reference = chr20 / 63025520.

I never ever changed or modified or touched the reference file (any sort of reference file whether it is ucsc.hg19.fasta or ucsc.hg19.dict or ucsc.hg19.fasta.fai). I checked files , such as, in BAM header or all reference files, the length of chr20 is reported as 63025520. I also checked --known files, 1000G_phase1.indels.hg19.sites.vcf and Mills_and_1000G_gold_standard.indels.hg19.sites.vcf.gz
(which I downloaded from where the chr20 length is 63025520. I donít know from where the length of chr20 64444167 comes here? I religned again and again and got the same error.
Thanks in advance
Waqasuddin Khan.
Waqasuddin Khan is offline   Reply With Quote
Old 09-26-2015, 06:18 AM   #2
Senior Member
Location: France

Join Date: Apr 2010
Posts: 141

cross posted
lindenb is offline   Reply With Quote
Old 09-26-2015, 08:27 AM   #3
Richard Finney
Senior Member
Location: bethesda

Join Date: Feb 2009
Posts: 689

grep chr20 chrom.sizes.hg*

chrom.sizes.hg18:chr20 62435964
chrom.sizes.hg19:chr20 63025520
chrom.sizes.hg38:chr20 64444167
Richard Finney is offline   Reply With Quote

chr20 contig length, gatk, reference

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 10:14 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO