DavidsonGroup / flexiplex

The Flexible Demultiplexer
https://davidsongroup.github.io/flexiplex/
MIT License
23 stars 2 forks source link

Try trim remainder of the polyT #35

Open ChangqingW opened 4 months ago

ChangqingW commented 4 months ago

Try trimming the remainder of the polyT with seq.find_first_not_of('T', result.endLocations[0]); Replaced pointers with const reference.

curl -sL https://raw.githubusercontent.com/mritchielab/FLAMES/devel/inst/extdata/fastq/musc_rps24.fastq.gz | zcat > musc_rps24.fastq
curl -sL https://raw.githubusercontent.com/mritchielab/FLAMES/devel/inst/extdata/bc_allow.tsv.gz | zcat > bc_allow.tsv
./flexiplex -x CTACACGACGCTCTTCCGATCT -b '????????????????' -u '????????????' -x TTTTTTTTT -e 2 -f 8 -k bc_allow.tsv musc_rps24.fastq > new.fq
vim -d <(sed -n '2~4p' old.fq) <(sed -n '2~4p' new.fq)
Screenshot 2024-03-15 at 3 33 15 pm