Entering edit mode
4.7 years ago
dpc
▴
250
Hi there!!!
I want to find gut microbiome composition from V6 region of 16S rRNA sequence. Unidirectional sequencing was done with 917F & 1391R primers. I am using MOTHUR software. Can anybody please tell me how can I extract the V6 region from the SILVA database using pcr.seqs command?
For V3-V4 region we can use the command: pcr.seqs(fasta=silva.bacteria.fasta, start=6388, end=25319, keepdots=F,processors=8)
Similarly, what start and end position can we use for V6 region?
Thanks and Regards
Read this: https://mothur.org/blog/2016/Customization-for-your-region/
Follow same for getting start and end position for V6 region
Thanks a lot gayachit.