1. Merge re-sequenced FastQ files (cat)
2. Read QC (FastQC)
3. Adapter trimming (fastp)
4. Removal of host reads (Kraken 2; optional)
5. Variant calling
   1. Read alignment (Bowtie 2)
   2. Sort and index alignments (SAMtools)
   3. Primer sequence removal (iVar; amplicon data only)
   4. Duplicate read marking (picard; optional)
   5. Alignment-level QC (picard, SAMtools)
   6. Genome-wide and amplicon coverage QC plots (mosdepth)
   7. Choice of multiple variant calling and consensus sequence generation routes (iVar variants and consensus; default for amplicon data || BCFTools, BEDTools; default for metagenomics data)
      • Variant annotation (SnpEff, SnpSift)
      • Consensus assessment report (QUAST)
      • Lineage analysis (Pangolin)
      • Clade assignment, mutation calling and sequence quality checks (Nextclade)
      • Individual variant screenshots with annotation tracks (ASCIIGenome)
   8. Intersect variants across callers (BCFTools)
6. De novo assembly
   1. Primer trimming (Cutadapt; amplicon data only)
   2. Choice of multiple assembly tools (SPAdes || Unicycler || minia)
      • Blast to reference genome (blastn)
      • Contiguate assembly (ABACAS)
      • Assembly report (PlasmidID)
      • Assembly assessment report (QUAST)
7. Present QC and visualisation for raw read, alignment, assembly and variant calling results (MultiQC)

We first operate a rapid inspection of the different BAM files using samtools flagstat. Illumina provided chr21 read mapping obtained with their GA IIx deep sequencing platform <ftp://webdata:webdata@ussd-ftp.illumina.com/Data/SequencingRuns/NA18507_GAIIx_100_chr21.bam>, aligned to the b36/hg18 reference genome)

Address of the bookmark: https://wiki.bits.vib.be/index.php/NGS_Exercise.6#compare_aln_.26_mem_results_with_vcf-compare

What is PhyloHerb: PhyloHerb is a wrapper program to process genome skimming data collected from plant materials. The outcomes include the plastid genome (plastome) assemblies, mitochondrial genome assemblies, nuclear ribosomal DNAs (NTS+ETS+18S+ITS1+5.8S+ITS2+28S), alignments of gene and intergenic regions, and a species tree. It is designed to be a high throughput program dealing with lower quality data. Examples include low-coverage (5x cpDNA) plastome phylogeny, recycling plastid genes from target enrichment data, retrieving low-copy nuclear genes from medium coverage (5x nucDNA) genome skimming.

License: GNU General Public License

Citation:

• Cai, Liming, Hongrui Zhang, and Charles C. Davis. 2022. PhyloHerb: A high‐throughput phylogenomic pipeline for processing genome‐skimming data. Applications in Plant Sciences 10(3): 1–9. https://doi.org/10.1002/aps3.11475

Address of the bookmark: https://github.com/lmcai/PhyloHerb/

More at http://busco.ezlab.org/

Address of the bookmark: http://busco.ezlab.org/

More at https://sourceforge.net/projects/sb2nhri/files/MyPro/

Address of the bookmark: http://www.sciencedirect.com/science/article/pii/S0167701215001207

http://www.bx.psu.edu/~ratan/

Address of the bookmark: http://www.bx.psu.edu/miller_lab/

Address of the bookmark: https://github.com/nextgenusfs/funannotate

Further information about LACHESIS, including source code, documentation and a user's guide are available at: http://shendurelab.github.io/LACHESIS.

Manuscript describing LACHESIS was published as: Burton JN#, Adey A, Patwardhan RP, Qiu R, Kitzman JO, Shendure J#. Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions. Nature Biotechnology 2013 Dec;31(12):1119-25. doi: 10.1038/nbt.272. PubMed PMID: 24185095.

http://shendurelab.github.io/LACHESIS/

Address of the bookmark: http://shendurelab.github.io/LACHESIS/

A rescently publication paper ( http://jb.asm.org/content/169/12/5429.long )shows significance of an unusual repeated DNA sequences next to a gene in a common bacterium, and their scientific significance. The sequences, it turns out, are part of a sophisticated immune system that bacteria use to fight viruses. And that system, whose very existence was unknown until about seven years ago, may provide scientists with unprecedented power to rewrite the code of life. This means a genome can be edited, much as a writer might change words or fix spelling errors. It allows “customizing the genome of any cell or any species at will,”.

Reference:

http://www.prweb.com/releases/2014/03/prweb11636031.htm

http://www.nytimes.com/2014/03/04/health/a-powerful-new-way-to-edit-dna.html?hpw&rref=health

http://jb.asm.org/content/169/12/5429.long

tax_id:
the unique identifier provided by NCBI Taxonomy
for the species or strain/isolate

GeneID:
the unique identifier for a gene
ASN1: geneid

Symbol:
the default symbol for the gene
ASN1: gene->locus

LocusTag:
the LocusTag value
ASN1: gene->locus-tag

Synonyms:
bar-delimited set of unofficial symbols for the gene

dbXrefs:
bar-delimited set of identifiers in other databases
for this gene. The unit of the set is database:value.
Note that HGNC and MGI include 'HGNC' and 'MGI', respectively,
in the value part of their identifier. Consequently,
dbXrefs for these databases will appear like:
HGNC:HGNC:1100
This would be interpreted as database='HGNC', value='HGNC:1100'
Example for MGI:
MGI:MGI:104537
This would be interpreted as database='MGI', value='MGI:104537'

chromosome:
the chromosome on which this gene is placed.
for mitochondrial genomes, the value 'MT' is used.

map location:
the map location for this gene

description:
a descriptive name for this gene

type of gene:
the type assigned to the gene according to the list of options
provided in https://www.ncbi.nlm.nih.gov/IEB/ToolBox/CPP_DOC/lxr/source/src/objects/entrezgene/entrezgene.asn

Symbol from nomenclature authority:
when not '-', indicates that this symbol is from a
a nomenclature authority

Full name from nomenclature authority:
when not '-', indicates that this full name is from a
a nomenclature authority

Nomenclature status:
when not '-', indicates the status of the name from the
nomenclature authority (O for official, I for interim)

Other designations:
pipe-delimited set of some alternate descriptions that
have been assigned to a GeneID
'-' indicates none is being reported.

Modification date:
the last date a gene record was updated, in YYYYMMDD format

Feature type:
pipe-delimited set of annotated features and their classes or
controlled vocabularies, displayed as feature_type:feature_class
or feature_type:controlled_vocabulary, when appropriate; derived
from select feature annotations on RefSeq(s) associated with the
GeneID

Address of the bookmark: ftp://ftp.ncbi.nih.gov/gene/DATA/GENE_INFO/