Hi,
I encounter a error when I am running the bowtie2 for genome mapping.
I use the latest version, which is bowtie2. 2.3
After installing bowtie2, I build the index for it, which run successfully and the index files are generated. Then I began to mapping using the following command:
bowtie2 \
-N 1 \
-x Atbt2index \
-1 /home/zxu/document/GATA_ChIP-seq/CGA1_INPUT_Rep1_R1.fastq \
-2 /home/zxu/document/GATA_ChIP-seq/CGA1_INPUT_Rep1_R2.fastq \
--un /home/zxu/document/GATA_ChIP-seq/Unaligned_CGA1_INPUT_Rep1.fq \
-S /home/zxu/document/GATA_ChIP-seq/CGA1_INPUT_Rep1.sam
At the beginning the program run successfully, and I can see the both the sam file and the file for unaligned read are generated. And I can also see the file size of .sam is increasing as the program running. However, when the file size increased around 1GB, the program died with the following information.
Error, fewer reads in file specified with -1 than in file specified with -2
terminate called after throwing an instance of 'int'
bowtie2-align died with signal 6 (ABRT) (core dumped)
The sam file can still be opened and I can see the data in there, I assume that the data is not complete. And the file for the unaligned reads are still empty with noting in it.
It seems like it is the problem of my input files. Anyone knows what is going on with my file data? Many thanks.
Zhenhua
I believe you have files with different numbers of reads. You should check this and verify the existence of singletons.
What's the output of
wc -l /home/zxu/document/GATA_ChIP-seq/CGA1_INPUT_Rep1_R?.fastq?I am getting the same problem... while using Trinity though...
output of `wc -l`
127349656 /home/kaalindi/3N/3Nleft.fastq
Presumably the other file has more/fewer reads. You should find out why that happened and fix it.