2012-06-20 1 views
0

bam 파일에서 일부 FastQ 파일을 추출하려고합니다. Picard는 SamToFastq로이 작업을 수행 할 수 있습니다.이 도구의 설명서에는 bam 또는 sam 파일 중 하나가 허용됩니다.Picard SamToFastq는 읽기 하나만 추출한 다음 오류가 발생합니다.

하지만 실행하면 읽기 하나만 추출한 다음 종료됩니다. 다음은 오류 메시지입니다. 어떤 도움을 주셔서 감사합니다. 그것이 나오는 것에 따라

[[email protected] picard-tools-1.70]$ java -jar SamToFastq.jar I=/home/davy/xxx_trio_data/xxxx-1.bam F=/home/davy/xxx_trio_data/1005-1.fastq 
[Wed Jun 20 14:14:21 BST 2012] net.sf.picard.sam.SamToFastq INPUT=/home/davy/xxx_trio_data/xxxx-1.bam FASTQ=/home/davy/xxxx_trio_data/xxxx-1.fastq OUTPUT_PER_RG=false RE_REVERSE=true INCLUDE_NON_PF_READS=false READ1_TRIM=0 READ2_TRIM=0 INCLUDE_NON_PRIMARY_ALIGNMENTS=false VERBOSITY=INFO QUIET=false VALIDATION_STRINGENCY=STRICT COMPRESSION_LEVEL=5 MAX_RECORDS_IN_RAM=500000 CREATE_INDEX=false CREATE_MD5_FILE=false 
[Wed Jun 20 14:14:21 BST 2012] Executing as [email protected] on Linux 2.6.34.9-69.fc13.x86_64 amd64; OpenJDK 64-Bit Server VM 1.6.0_18-b18; Picard version: 1.70(1215) 
[Wed Jun 20 14:14:21 BST 2012] net.sf.picard.sam.SamToFastq done. Elapsed time: 0.00 minutes. 
Runtime.totalMemory()=2029715456 
FAQ: http://sourceforge.net/apps/mediawiki/picard/index.php?title=Main_Page 
Exception in thread "main" java.lang.IndexOutOfBoundsException: Index: 1, Size: 1 
     at java.util.ArrayList.rangeCheck(ArrayList.java:571) 
     at java.util.ArrayList.get(ArrayList.java:349) 
     at net.sf.picard.sam.SamToFastq.doWork(SamToFastq.java:156) 
     at net.sf.picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:177) 
     at net.sf.picard.sam.SamToFastq.main(SamToFastq.java:118) 

답변

1

, 데이터가 페어링 끝이 단일하지 내가 생각했던 처음에했던대로 읽기 및 피카드는 SECOND_END_FASTQ 옵션으로 지정한이 경우에 두 번째 OUTFILE을 필요로한다.

관련 문제