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
___________________________________________________________
The Galaxy User list should be used for the discussion of
Galaxy analysis and other features on the public server
at usegalaxy.org. Please keep all replies on the list by
using "reply all" in your mail client. For discussion of
local Galaxy instances and the Galaxy source code, please
use the Galaxy Development list:
http://lists.bx.psu.edu/listinfo/galaxy-dev
To manage your subscriptions to this and other Galaxy lists,
please use the interface at:
http://lists.bx.psu.edu/
To search Galaxy mailing lists use the unified search at:
http://galaxyproject.org/search/mailinglists/