Thanks a lot
Hi all,
My ATAC-seq library was generated from a genetic-modified human cell line which carried genes (e.g GFP) that not found in reference human genome. in this case, standard mapping tools such as Bowtie2 and STAR will not be the suitable tools, is it correct? if so, should I map them as a de-novo-genome and use tool like Flye?
Thanks
1 answer
Just because you have a single transgene does not mean you cannot use standard tools. If you want the reads included mapping to it then add the cDNA sequence to the reference genome add make a new index. Or just use the standard index. It in this case does not really matter to align the potentially few reads coming from the transgene (if any, maybe the promoter). I like bowtie2 as it predictably and reliably uses a low memory footprint.
Log in to answer this question.