I used the SplitNCigarReads with the following command.
java -jar GenomeAnalysisTK.jar -T SplitNCigarReads -R test.fasta -I sorted2.bam -o sorted2_split.bam -rf ReassignOneMappingQuality -RMQF 255 -RMQT 60 -U ALLOW_N_CIGAR_READS
But I obtained following log and error messages.
I made index (bai file) using samtools index.
Previously I succeeded in the same way for a different reference sequence.
I wanted to ask you how to solve this problem.
15:08:19,279 HelpFormatter - --------------------------------------------------------------------------------
INFO 15:08:19,284 HelpFormatter - The Genome Analysis Toolkit (GATK) v3.5-0-g36282e4, Compiled 2015/11/25 04:03:56
INFO 15:08:19,284 HelpFormatter - Copyright (c) 2010 The Broad Institute
INFO 15:08:19,284 HelpFormatter - For support and documentation go to http://www.broadinstitute.org/gatk
INFO 15:08:19,288 HelpFormatter - Program Args: -T SplitNCigarReads -R test.fasta -I sorted2.bam -o sorted2_split.bam -rf ReassignOneMappingQuality -RMQF 255 -RMQT 60 -U ALLOW_N_CIGAR_READS
INFO 15:08:19,294 HelpFormatter - Executing as mizunobu@nm5 on Linux 2.6.32-358.el6.x86_64 amd64; Java HotSpot(TM) 64-Bit Server VM 1.8.0_45-b14.
INFO 15:08:19,295 HelpFormatter - Date/Time: 2017/01/23 15:08:19
INFO 15:08:19,295 HelpFormatter - --------------------------------------------------------------------------------
INFO 15:08:19,295 HelpFormatter - --------------------------------------------------------------------------------
INFO 15:08:19,972 GenomeAnalysisEngine - Strictness is SILENT
INFO 15:08:20,056 GenomeAnalysisEngine - Downsampling Settings: No downsampling
INFO 15:08:20,065 SAMDataSource$SAMReaders - Initializing SAMRecords in serial
INFO 15:08:20,098 SAMDataSource$SAMReaders - Done initializing BAM readers: total time 0.03
INFO 15:08:20,211 GenomeAnalysisEngine - Preparing for traversal over 1 BAM files
INFO 15:08:20,217 GenomeAnalysisEngine - Done preparing for traversal
INFO 15:08:20,217 ProgressMeter - [INITIALIZATION COMPLETE; STARTING PROCESSING]
INFO 15:08:20,218 ProgressMeter - | processed | time | per 1M | | total | remaining
INFO 15:08:20,218 ProgressMeter - Location | reads | elapsed | reads | completed | runtime | runtime
INFO 15:08:20,252 ReadShardBalancer$1 - Loading BAM index data
INFO 15:08:20,254 ReadShardBalancer$1 - Done loading BAM index data
INFO 15:08:51,252 ProgressMeter - chr1:88350198 300086.0 31.0 s 103.0 s 0.6% 85.1 m 84.6 m
INFO 17:15:37,400 ProgressMeter - chr14:480923759 2.6647719E7 2.1 h 4.8 m 100.0% 2.1 h 0.0 s
INFO 17:15:57,978 GATKRunReport - Uploaded run statistics report to AWS S3