soapdenovo2

Reply to soapdenovo2 on Wed, 13 Jul 2022 06:07:54 GMT

anneng — Wed, 13 Jul 2022 06:07:54 GMT

a. *.contig
contig sequences without using mate pair information.
b. *.scafSeq
scaffold sequences (final contig sequences can be extracted by breaking down scaffold sequences at gap regions).

Reply to soapdenovo2 on Wed, 13 Jul 2022 05:42:48 GMT

anneng — Wed, 13 Jul 2022 05:42:48 GMT

对组装的结果进行质控
docker pull staphb/quast

 quast.py test_data/contigs_1.fasta \
           test_data/contigs_2.fasta \
        -r test_data/reference.fasta.gz \
        -g test_data/genes.txt \
        -1 test_data/reads1.fastq.gz -2 test_data/reads2.fastq.gz \
        -o quast_test_output

Reply to soapdenovo2 on Tue, 12 Jul 2022 12:23:29 GMT

anneng — Tue, 12 Jul 2022 12:23:29 GMT

Highly variable chloroplast genome from two endangered Papaveraceae lithophytes Corydalis tomentella and Corydalis saxicola
https://onlinelibrary.wiley.com/doi/full/10.1002/ece3.7312
2.2 Genome assembly and annotation
The cp genome were assembled on a Linux system. First, raw sequencing data were filtered using Trimmomatic (Version 0.36) to get the high-quality clean data (Bolger et al., 2014). In the second step, we used the thirteen chloroplast genome sequences of Papaveraceae species which were downloaded from GenBank to establish a Basic Local Alignment Search Tool (BLASTn) database. Then the clean data were mapped to the BLAST database, and the mapped reads which were considered as reads from chloroplast genome were extracted. Next step, the extracted reads were assembled to contigs using SOAPdenovo2 (Luo et al., 2012). At last, SSPACE was used to construct the scaffold of the chloroplast genome (Boetzer et al., 2011), and GapCloser was used to fill gaps (Luo et al., 2012). The completed genomes were annotated using CPGAVAS2 (Shi et al., 2019), and the results were modified for starter and terminator revisions by Apollo software (Lee et al., 2009). CPGAVAS2 software was used to convert revised GFF3 format annotation results into a sqn format for NCBI submission. Sequin software was used to check and correct unsatisfactory comments in the sqn file, and the corrected results were submitted to the NCBI database. Physical maps of the cp genomes were drawn by GenomeDRAW (Marc et al., 2013) using a GB format file exported from the sqn file by sequin software.

Reply to soapdenovo2 on Tue, 12 Jul 2022 07:41:00 GMT

anneng — Tue, 12 Jul 2022 07:41:00 GMT

cd /home/bioinfo/soapdenovo2/SOAPdenovo2-r242
./SOAPdenovo-63mer all -s tobacco.config -K 63 -R -o tobacco

这个命令的输入文件是在配置文件里面写的建议用shell封装下把 R1 R2 暴露出来其他参数用配置文件也可以