T
T
Tutorials
Search…
Mapping pipeline
First log into the serve
Copied!
Then go to the correct directory with the data:
1
cd /opt/storage5/monica/NGS_data
Copied!
If you have a new reference you have to index it using:
1
bwa index <reference_file>
Copied!
Perform mapping using the following code
1
bwa mem -t 10 -R "@RG\tID:AGT12\tSM:ATG12\tPL:Illumina" Anopheles_gambiae.AgamP4.dna.toplevel.fa AGT12_1.fastq.gz AGT12_2.fastq.gz | samtools view [email protected] 10 -b - | samtools sort [email protected] 10 -o AGT12.bam -
Copied!
Then index the bam file with
1
samtools index AGT12.bam
Copied!
Create a table of variants
1
bcftools mpileup -f Anopheles_gambiae.AgamP4.dna.toplevel.fa AGT12.bam | bcftools call -mv | bcftools query -f '%CHROM\t%POS\t%REF\t%ALT[\t%TGT]\n'
Copied!
Copy link