understanding bxtools stats output about bxtools HOT 4 CLOSED

Hi Danshu,
I think the best I could advise would be to make a histogram of the BX counts and find a reasonable cutoff for your own data. bxtools was a sort of weekend-type project to help with some 10X analyses led by some collaborators, so I wouldn't presume at this point to give much advice on interpreting the outputs. e.g. the value of "100" was completely ad hoc. I'm not really sure what is meant by their filter either, so you might reach out to the authors.

That said, I'd like to learn more sometime about 10X data. In the meantime, if you have any suggestions for the types of low-level data manipulation (like stats or split) that would be useful from bxtools, do let me know. I'd be interested in expanding this!

Best,
Jeremiah

from bxtools.

Hi Jeremiah,

Thanks for your kind advice. Actually I have tried to plot the BX counts to find a reasonable threshold. But what matters should be the characteristics of the technology, e.g. low frequency of BX count indicates poor quality of that pool, which I'm not sure. Anyway I will try to ask the company for more details.

Actually I'm trying to use my 10x genomics data for scaffolding my assembly. One common requirement of those scaffolders are attaching BX to the read name, while different tools have different styles of attaching BX. For example, ARCS (https://github.com/bcgsc/arcs) will be expecting the barcode at the end of the read name, in this format:
READNAME_TAGCATAGACATCAGA
Great if bxtools may combine filtering with attaching BXs to read names.

I managed to do this job by extracting filtered reads using "seqtk subseq" and then rename them according to the format required: "awk '{ if (NR%4==1) {split($2,a,":"); split(a[3],b,"-"); print $1"_"b[1] } else { print } }' test.fastq > rename_test.fastq".

Best,
Danshu

from bxtools.

I'm in a mood to procrastinate on other things at the moment, decided to give this a shot... Just added relabel. Does this do what you need? You can pipe the output from relabel to samtools view and AWK to get to a fastq from a BAM if you need.

from bxtools.

Thanks and I will try it later!

from bxtools.