Usando o STRUCTURE

 structure_threader run -i Ugibba_filtered1_ref.vcf -o ./results2 -als ~/anaconda3/envs/structure/bin/alstructure_wrapper.R -K 15 -t 3 --ind Ugibba.indfile

Fonte: https://radcamp.github.io/Lisbon2020/05_POPULATION_STRUCTURE.html


Sendo que o arquivo ind:

The “indfile”

Open the file with your favorite text editor. If you don’t have one and want to stay on the CLI, you can try the command nano ~/str_analyses/oyster.indfile in a new terminal window. The first few lines of this file should look like this:

BC1_10_C6	Victoria_BC	4
BC1_20_C6	Victoria_BC	4
BC1_22_C7	Victoria_BC	4
BC1_4_C3	Victoria_BC	4
BC1_7_C5	Victoria_BC	4
BC1_8_C4	Victoria_BC	4
BC1_9_C5	Victoria_BC	4
BC2_10_C5	Klaskino_BC	1

Comentários

Postar um comentário

Postagens mais visitadas deste blog

Instalar o VMPlayer no Linux Ubuntu

Quando fui instalar o VMPlayer no Linux encontrei o seguinte erro: Extracting VMware Installer...done. Installing VMware Player Application 15.5.0 Copying files... Configuring... bora/lib/string/str.c:284 Buffer too small VMware Workstation Error: VMware Workstation unrecoverable error: (host-9802) bora/lib/string/str.c:284 Buffer too small You can request support.  To collect data to submit to VMware support, choose "Collect Support Data" from the Help menu. You can also run the "vm-support" script in the Workstation folder directly. We will respond on the basis of your support entitlement. A solução foi: sudo LC_ALL=C ./VMware-Player-15.5.0-14665864.x86_64.bundle
Leia mais
 Retirado de: https://playingwithgenomics.wordpress.com/2015/07/15/producing-a-bam-file-and-extracting-unique-reads-from-bowtie2-results/ I am always looking for ways to keep my disk usage down. Especially since I’ve been mapping close to 100 Chip-Seq files. I find that 1) piping Bowtie2 output into samtools to create a bam file and 2) keeping only the uniquely mapped reads help a lot. Here is how I do those: I’m dealing with single-end data here. You can modify the command if you are dealing with paired-end data. Text written in red should be substituted according to your data. Converting Bowtie2 output to bam file bowtie2 -x BOWTIE2_INDEX -\ -p NUM_THREADS \ -U INPUT_FILE.fastq.gz | samtools view -bS -t \ YOUR_GENOME_INDEX.fa.fai - > output.bam Counting unique reads You basically can exploit the XS tag that are set by Bowtie2  for reads that can be mapped in multiple places. Therefore, uniquely mapped reads lack the XS tag. The following command exclude u...
Leia mais