Question: Depth Of Coverage Error
0
Sandrine Imbeaud • 20 wrote:
Hi,
I am using the "Depth of coverage on BAM files" tool from the NGS:
GATK
Tools but encounter problem.
While part of the BAM files proceed successfully, some BAM files end
proceeding with a systematic error (see below).
The entire BAM file dataset was generated with the same deepseq
analysis
pipeline (MiSeq Reporter, PCR amplicon Workflow) => alignemnent is
done
on Hs. hg19 and manifest is focused on specific regions.
Does someone know how to solve this problem?
Also, the chromosome X and Y appear to be exclude from the
calculation.
Is there any settings to select in order to include both sexual
chromosomes? Is it expected?
Kind Regards
/ Sandrine
--
Sandrine Imbeaud
INSERM, UMR U-674, IUH
Université Paris Descartes
Génomique Fonctionnelle des tumeurs solides
27 rue Juliette Dodu
F75010 Paris, France
TEL: +33 (0)1 53 72 51 98
FAX: +33 (0)1 53 72 51 92
MOBILE: +33 (0)6 12 69 80 29
http://www.inserm-u674.net/
Picked up _JAVA_OPTIONS: -Djava.io.tmpdir=/tmp
[Mon Feb 10 07:08:03 CST 2014]
net.sf.picard.sam.CreateSequenceDictionary
REFERENCE=/tmp/tmp-gatk-hPNp7B/gatk_input.fasta
OUTPUT=/tmp/tmp-gatk-hPNp7B/dict1033574942916796988.tmp
TRUNCATE_NAMES_AT_WHITESPACE=true NUM_SEQUENCES=2147483647
VERBOSITY=INFO QUIET=false VALIDATION_STRINGENCY=STRICT
COMPRESSION_LEVEL=5 MAX_RECORDS_IN_RAM=500000 CREATE_INDEX=false
CREATE_MD5_FILE=false
[Mon Feb 10 07:08:03 CST 2014] Executing as
g2main@roundup49.tacc.utexas.edu on Linux
2.6.32-431.1.2.0.1.el6.x86_64
amd64; Java HotSpot(TM) 64-Bit Server VM 1.7.0_40-b43; Picard version:
1.58(1057)
[Mon Feb 10 07:08:04 CST 2014]
net.sf.picard.sam.CreateSequenceDictionary done. Elapsed time: 0.00
minutes.
Runtime.totalMemory()=2025324544
##### ERROR
##### ERROR A USER ERROR has occurred (version 1.4-18-g80a4ce0):
##### ERROR The invalid arguments or inputs must be corrected before
the
GATK can proceed
##### ERROR Please do not post this error to the GATK forum
##### ERROR
##### ERROR See the documentation (rerun with -h) for this tool to
view
allowable command-line arguments.
##### ERROR Visit our wiki for extensive documentation
http://www.broadinstitute.org/gsa/wiki
##### ERROR Visit our forum to view answers to commonly asked
questions
http://getsatisfaction.com/gsa
##### ERROR
##### ERROR MESSAGE: Input files reads and reference have incompatible
contigs: No overlapping contigs found.
##### ERROR reads contigs = [chrM, chr1, chr2, chr3, chr4, chr5,
chr6,
chr7, chr8, chr9, chr10, chr11, chr12, chr13, chr14, chr15, chr16,
chr17, chr18, chr19, chr20, chr21, chr22, chrX, chrY]
##### ERROR reference contigs = [c1_27022371-27108741,
c1_65310990-65312594, c1_103341958-103574198, c1_115256193-115258963,
c11_533546-534529, c11_108093339-108239879, c12_25379959-25398502,
c12_46123245-46302016, c12_49412540-49449288, c12_70910587-71031382,
c12_121416337-121440487, c13_77618564-77901349, c16_337220-402825,
c16_2097254-2138873, c17_7571510-7591038, c17_40474813-40500706,
c19_10596592-10614570, c2_21224182-21267085, c2_178098515-178099162,
c20_57484200-57484647, c3_41240709-41282027, c3_41265294-41266864,
c3_47057688-47205620, c3_178928002-178948337, c4_74269729-74287296,
c4_185308648-185395907, c5_1253083-1294978, c5_1294819-1295755,
c5_55251637-55260189, c5_71402848-71505569, c6_36644019-36655290,
c7_135242597-135333692, c7_140452856-140453348, c9_21967531-21994664,
c9_135766517-135820181]
##### ERROR
mv: cannot stat
`/galaxy-repl/main/files/007/581/dataset_7581988.dat.sample_summary':
No
such file or directory
ADD COMMENT
• link
•
modified 4.8 years ago
by
Jennifer Hillman Jackson ♦ 25k
•
written
4.8 years ago by
Sandrine Imbeaud • 20
