I have gut human microbiome data, I would like to filter out human sequences. What would be the best practice? I thought of megablasting the reads (454) against the human genome and removing 100% identities. Any ideas to make this better / shorter?
Thanks,
Iddo
•
written
by
Iddo
3 answers
If you have shotgun metagenomic data, this sounds like a reasonable plan - though you could swap out megabalst for a different/faster aligner like BLAT, and allow some degree of mismatch for sequencing error. If you have have targeted metagenomic data (e.g. a 16S amplicon), then you could additionally impose a filter for contamination hits to be off-target (e.g. non-16S).
If BLASTN against the NCBI-NT is possible for the reads, you could also try to use CARMA3 to filter out human sequences...
Log in to answer this question.
-
Different results in finding human sequences from metagenomes (bowtie2 and Kraken2)
written by kmat •
-
Reads that used for re-constructing MAGs can not be classified by Kraken2
written by Yifeng •
-
kraken2 mapping only 10% metagenomic reads from the gut microbiome
written by luthien •
-
Newbie question: how to create BWA index from multiple .fasta files
written by phonybone •
-
comparing two genomes
written by dr_mm_sh •
-
Best way to find phage genes in microbiome?
written by yax •
-
How Does The Metagenome Assembly Process Affect Snps?
written by weslfield
-
How To Compare The 454 Vs Illumina Sequences
written by deepthithomaskannan
-
Most Sensitive Method To Search 100-1000Bp Sequences To Genome At Human-Mouse Genetic Distances
written by 2184687-1231-83-
-
Diversity Analysis From 454 Data (No 16S Amplicons)
written by Iddo
