Question: [ti_index_core] the chromosome blocks not continuous at line xxxx, is the file sorted? [pos 1]
0
gravatar for shubhra.bhattacharya
12 months ago by
shubhra.bhattacharya120 wrote:

While using tabix I encounter the following error.

[ti_index_core] the chromosome blocks not continuous at line 5891584, is the file sorted? [pos 1]

Command: /data/ngs/programs/tabix-0.2.6/tabix -p vcf SRR12954_1.fastq.gz_20180228112323.GRCH38.vcf.gz

I have zipped the vcf file using bgzip data/ngs/programs/tabix-0.2.6/bgzip -c SRR12954_1.fastq.gz_20180228112323.GRCH38.vcf > SRR12954_1.fastq.gz_20180228112323.GRCH38.vcf.gz

I have also sorted the vcf using vcf-sort but run into the same problem

index variant vcf • 825 views
ADD COMMENTlink written 12 months ago by shubhra.bhattacharya120

what is the context around the line 5891584 ?

gunzip -c SRR12954_1.fastq.gz_20180228112323.GRCH38.vcf.gz | cut -f 1,2,4,5 | nl | awk '(NR> 5891580 && NR< 5891590)'
ADD REPLYlink written 12 months ago by Pierre Lindenbaum118k

5891581 chrMT 16129 G A 5891582 chrMT 16129 G A 5891583 chrMT 16179 CAA C 5891584 chr1 3088286013 AAA CCC 5891585 chrMT 16182 A C 5891586 chrMT 16183 A C 5891587 chrMT 16186 C T 5891588 chrMT 16188 C T 5891589 chrMT 16189 T C

This is the output of the command. The chr1 is just dropping by in b/w the chrMT chromosome block.

ADD REPLYlink written 12 months ago by shubhra.bhattacharya120
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1726 users visited in the last hour