How to get multiple genes from bam files
Entering edit mode
6 days ago
沛煒 • 0

Hello guys, I am working on a pool-seq data with 80 genes sequenced from 400 individuals. I have aligned them to the reference and got 400 bam files. I'd like to get all of the 80 genes and write them into 80 fasta files (each contains 400 orthologous sequences).
I already knew the location of each gene on the reference (e.g., gene1 at C1:12211-13100, gene2 at C5: 56733-56823...). Does anyone have any suggestion to output the fasta files separated with genes? Thanks

bam samtools fasta bioinformatics • 99 views

Login before adding your answer.

Traffic: 1883 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6