Virus+ Sequence Masked Human Reference Genome (hg19)

<p>A version of the human genome (hg19) originally masked for ribosomal, plant, animal, fungal and&nbsp;low-entropy sequences&nbsp;by Brian Bushnell (<a href="https://zenodo.org/record/1208052#.X5BuTy9h3UI">Bushnell Masked Human Genome</a>) additionally masked for all possible viral sequences.</p> <p>The following commands were used to generate the additional virus sequence masked reference database:</p> <p><strong>1) Download all RefSeq and Neighbor nucleotide records:</strong></p> <p><a href="https://www.ncbi.nlm.nih.gov/nuccore/?term=Viruses[Organism]%20NOT%20cellular%20organisms[ORGN]%20NOT%20wgs[PROP]%20NOT%20gbdiv%20syn[prop]%20AND%20(srcdb_refseq[PROP]%20OR%20nuccore%20genome%20samespecies[Filter])">https://www.ncbi.nlm.nih.gov/nuccore/?term=Viruses[Organism]%20NOT%20cellular%20organisms[ORGN]%20NOT%20wgs[PROP]%20NOT%20gbdiv%20syn[prop]%20AND%20(srcdb_refseq[PROP]%20OR%20nuccore%20genome%20samespecies[Filter])</a></p> <p><strong>2) Shred the downloaded viral genomes using shred.sh from the <a href="https://jgi.doe.gov/data-and-tools/bbtools/">bbtools</a> package</strong></p> <p>shred.sh in=refseq_virus_reformated.fasta out=virus_shred.fasta.gz length=85 minlength=75 overlap=30</p> <p><strong>3) Map shredded virus sequence to the hg19-masked human genome using bbmap.sh&nbsp;from the&nbsp;<a href="https://jgi.doe.gov/data-and-tools/bbtools/">bbtools</a> package</strong></p> <p>bbmap.sh ref=hg19_main_mask_ribo_animal_allplant_allfungus.fa.gz in=virus_shred.fasta.gz outm=map_human_all_viruses.sam minid=0.90</p> <p><strong>4) Mask virus sequenced mapped regions from the hg19-masked human genome using bbmask.sh from the&nbsp;<a href="https://jgi.doe.gov/data-and-tools/bbtools/">bbtools</a> package</strong></p> <p>bbmask.sh in=hg19_main_mask_ribo_animal_allplant_allfungus.fa.gz out=human_virus_masked.fasta.gz sam=map_human_all_viruses<br> .sam</p> <p><strong>5) Remove all N&#39;s to further reduce file size using <a href="https://bioinf.shenwei.me/seqkit/">seqkit</a></strong><br> seqkit -is replace -p &quot;n&quot; -r &quot;&quot; human_virus_masked.fasta.gz &nbsp;&gt; human_virus_masked.fasta_Ns_removed.gz</p> <p><strong>Additional References:</strong></p> <ol> <li><a href="http://seqanswers.com/forums/showthread.php?t=42552">http://seqanswers.com/forums/showthread.php?t=42552</a> for additional information on the original masking of hg19</li> <li><a href="https://jgi.doe.gov/data-and-tools/bbtools/">bbtools</a></li> <li><a href="https://bioinf.shenwei.me/seqkit/">seqkit</a></li> <li><a href="https://www.ncbi.nlm.nih.gov/genome/viruses/">NCBI Virus Genome RefSeq</a></li> </ol>