Typing of S. epidermidis samples (MD P.B.)

Characterization_of_the_virulence_agr_typing_and_antimicrobial_resistance_profile_of_Staphylococcus_aureus_strains.pdf

Goal of data analyses

S epidermidis genomes
    1. Assemble closed genomes from HDRNA 1, 3, 6, 7, 8, 12, 16, 17, 19, 20 (short read + long read)
        TODO: make a table similar to the paper Characterization of the virulence, agr typing and antimicrobial resistance profile of Staphylococcus aureus strains isolated from food handlers in Brazil Table 2, and draw a tree+heatmaps Figure!
        Based on closed genomes:
        - Sequence type,
        - goeBURST analysis (-->performed is geoBURST TLV-Analysis, see point 2),
        - SCCmec type (https://www.sccmec.org/index.php/en/method-to-identify-sccmmcc-smn-en/review-smn-en) https://cge.cbs.dtu.dk/services/SCCmecFinder/
          SCCmec typing: https://www.genomicepidemiology.org/ --> https://cge.food.dtu.dk/services/SCCmecFinder/
            SCCmecFinder 1.2
            SCCmecFinder identifies SCCmec elements in sequenced S. aureus isolates. The SCCmec element is the defining feature of methicillin-resistant S. aureus isolates, and encodes the single determinant for methicillin resistant, the mecA gene.
            IMPORTANT! SCCmec typing is only available for SCCmec type I-XI and subtyping is currently only available for SCCmec type IV and V
            IMPORTANT! mec gene complex C1 and C2 might produce errors.

        - agr typing: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4187671/ (see point 3)

        - presence phage HH1, SPbeta-like phage, phage related island [referring to A´s paper] (see point 5),

        - absence/presence matrix for selected genes [see attached ppt, see results in Gene_List.pptx] (see point 6)
          gyrB,  fumC, , icd, apsS,
          sigB, sarA, , , ,
          ,  , sdrG(-17), sdrH, ebh, ebp (ebpS), ,
          , , dltA, , lipA,
          , , , , , ,
          --> draw a circle heatmaps including all data with a very big figure and table.

    2. according to assembled genomes describe within host diversity per patient (i.e. compare isolates 2 – 10 to isolate 1).
        python3 /home/jhuang/Scripts/gb_to_excel.py ./gbks/HDRNA_01_K01_conservative_23197.current.gb HDRNA_01_K01_CP133676.xlsx
        python3 /home/jhuang/Scripts/gb_to_excel.py ./gbks/HDRNA_03_K01_bold_bandage_26831.current.gb HDRNA_03_K01_CP133677.xlsx
        python3 /home/jhuang/Scripts/gb_to_excel.py ./gbks/HDRNA_06_K01_conservative_27645.current.gb HDRNA_06_K01_CP133678-CP133679.xlsx
        python3 /home/jhuang/Scripts/gb_to_excel.py ./gbks/HDRNA_07_K01_conservative_27169.current.gb HDRNA_07_K01_CP133680-CP133681.xlsx
        python3 /home/jhuang/Scripts/gb_to_excel.py ./gbks/HDRNA_08_K01_conservative_32455.current.gb HDRNA_08_K01_CP133682-CP133683.xlsx
        python3 /home/jhuang/Scripts/gb_to_excel.py ./gbks/HDRNA_12_K01_bold_37467.current.gb HDRNA_12_K01_CP133684-CP133687.xlsx
        python3 /home/jhuang/Scripts/gb_to_excel.py ./gbks/HDRNA_16_K01_conservative_37834.current.gb HDRNA_16_K01_CP133688-CP133692.xlsx
        python3 /home/jhuang/Scripts/gb_to_excel.py ./gbks/HDRNA_17_K01_conservative_37288.current.gb HDRNA_17_K01_CP133693-CP133695.xlsx
        python3 /home/jhuang/Scripts/gb_to_excel.py ./gbks/HDRNA_19_K01_bold_37377.current.gb HDRNA_19_K01_CP133696-CP133699.xlsx
        python3 /home/jhuang/Scripts/gb_to_excel.py ./gbks/HDRNA_20_K01_conservative_43457.current.gb HDRNA_01_K01_CP133700-CP133701.xlsx

        * SNP INDEL: snippy+spandx --> get the complete list of SNP+INDEL for each isolate group!!! (see point 4)

        * gene absence/presence:
        # A: prepare prokka_HDRNA_01 .. prokka_HDRNA_20 from prokka_remaining

        # B:         #https://github.com/jorvis/biocode/blob/master/gff/convert_genbank_to_gff3.py
        sudo apt-get install -y python3 python3-pip zlib1g-dev libblas-dev liblapack-dev libxml2-dev
        pip3 install biocode
        convert_genbank_to_gff3.py -i HDRNA_01_K01_conservative_23197.current.gb -o ~/DATA/Data_PaulBongarts_S.epidermidis_HDRNA/Data_Holger_S.epidermidis_short/prokka_HDRNA_01/HDRNA_01_K01.gff --with_fasta
        convert_genbank_to_gff3.py -i HDRNA_03_K01_bold_bandage_26831.current.gb -o ~/DATA/Data_PaulBongarts_S.epidermidis_HDRNA/Data_Holger_S.epidermidis_short/prokka_HDRNA_03/HDRNA_03_K01.gff --with_fasta
        convert_genbank_to_gff3.py -i HDRNA_06_K01_conservative_27645.current.gb -o ~/DATA/Data_PaulBongarts_S.epidermidis_HDRNA/Data_Holger_S.epidermidis_short/prokka_HDRNA_06/HDRNA_06_K01.gff --with_fasta
        convert_genbank_to_gff3.py -i HDRNA_07_K01_conservative_27169.current.gb -o ~/DATA/Data_PaulBongarts_S.epidermidis_HDRNA/Data_Holger_S.epidermidis_short/prokka_HDRNA_07/HDRNA_07_K01.gff --with_fasta
        convert_genbank_to_gff3.py -i HDRNA_08_K01_conservative_32455.current.gb -o ~/DATA/Data_PaulBongarts_S.epidermidis_HDRNA/Data_Holger_S.epidermidis_short/prokka_HDRNA_08/HDRNA_08_K01.gff --with_fasta
        convert_genbank_to_gff3.py -i HDRNA_12_K01_bold_37467.current.gb -o ~/DATA/Data_PaulBongarts_S.epidermidis_HDRNA/Data_Holger_S.epidermidis_short/prokka_HDRNA_12/HDRNA_12_K01.gff --with_fasta
        convert_genbank_to_gff3.py -i HDRNA_16_K01_conservative_37834.current.gb -o ~/DATA/Data_PaulBongarts_S.epidermidis_HDRNA/Data_Holger_S.epidermidis_short/prokka_HDRNA_16/HDRNA_16_K01.gff --with_fasta
        convert_genbank_to_gff3.py -i HDRNA_17_K01_conservative_37288.current.gb -o ~/DATA/Data_PaulBongarts_S.epidermidis_HDRNA/Data_Holger_S.epidermidis_short/prokka_HDRNA_17/HDRNA_17_K01.gff --with_fasta
        convert_genbank_to_gff3.py -i HDRNA_19_K01_bold_37377.current.gb -o ~/DATA/Data_PaulBongarts_S.epidermidis_HDRNA/Data_Holger_S.epidermidis_short/prokka_HDRNA_19/HDRNA_19_K01.gff --with_fasta
        convert_genbank_to_gff3.py -i HDRNA_20_K01_conservative_43457.current.gb -o ~/DATA/Data_PaulBongarts_S.epidermidis_HDRNA/Data_Holger_S.epidermidis_short/prokka_HDRNA_20/HDRNA_20_K01.gff --with_fasta

        # C: check gff.files if containing the repeated gene names using check_duplicate_cds.py, remove manually, recheck the removed one.
        python3 ~/Scripts/check_duplicate_cds.py HDRNA_01_K01.gff
        python3 ~/Scripts/check_duplicate_cds.py HDRNA_03_K01.gff
        python3 ~/Scripts/check_duplicate_cds.py HDRNA_06_K01.gff
        python3 ~/Scripts/check_duplicate_cds.py HDRNA_07_K01.gff
        python3 ~/Scripts/check_duplicate_cds.py HDRNA_08_K01.gff
        python3 ~/Scripts/check_duplicate_cds.py HDRNA_12_K01.gff
        python3 ~/Scripts/check_duplicate_cds.py HDRNA_16_K01.gff
        python3 ~/Scripts/check_duplicate_cds.py HDRNA_17_K01.gff
        python3 ~/Scripts/check_duplicate_cds.py HDRNA_19_K01.gff
        python3 ~/Scripts/check_duplicate_cds.py HDRNA_20_K01.gff
            Found duplicates for the following CDS IDs:
            RE430_09810.mRNA.0.CDS.1
            RE430_08580.mRNA.0.CDS.1
            RE430_09730.mRNA.0.CDS.1
            RGR13_00575.mRNA.0.CDS.1
            RGR10_09570.mRNA.0.CDS.1
            RGR06_10025.mRNA.0.CDS.1
            RGR06_12305.mRNA.0.CDS.1
            #xxxx
            RGR12_09490.mRNA.0.CDS.1
            RGR12_00425.mRNA.0.CDS.1
            RGR12_11570.mRNA.0.CDS.1
            RGR09_12270.mRNA.0.CDS.1
            RGR09_09635.mRNA.0.CDS.1
            RGR09_01280.mRNA.0.CDS.1
            RGR08_01845.mRNA.0.CDS.1
            RGR08_13135.mRNA.0.CDS.1
            RGR08_10360.mRNA.0.CDS.1
            RGR08_09845.mRNA.0.CDS.1
            RGR14_00635.mRNA.0.CDS.1
            RGR14_00355.mRNA.0.CDS.1
            RGR14_09870.mRNA.0.CDS.1
            RGR14_01550.mRNA.0.CDS.1

            RGR07_09075.mRNA.0.CDS.1
            RGR07_00190.mRNA.0.CDS.1
            RGR07_11485.mRNA.0.CDS.1
            RGR07_00090.mRNA.0.CDS.1
            RGR07_11545.mRNA.0.CDS.1
            RGR07_01915.mRNA.0.CDS.1
            RGR11_01280.mRNA.0.CDS.1
            RGR11_09700.mRNA.0.CDS.1

        rsync -a -P prokka_HDRNA_01/HDRNA_01_K01.gff jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_01/
        rsync -a -P prokka_HDRNA_03/HDRNA_03_K01.gff jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_03/
        rsync -a -P prokka_HDRNA_06/HDRNA_06_K01.gff jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_06/
        rsync -a -P prokka_HDRNA_07/HDRNA_07_K01.gff jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_07/
        rsync -a -P prokka_HDRNA_08/HDRNA_08_K01.gff jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_08/
        rsync -a -P prokka_HDRNA_12/HDRNA_12_K01.gff jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_12/
        rsync -a -P prokka_HDRNA_16/HDRNA_16_K01.gff jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_16/
        rsync -a -P prokka_HDRNA_17/HDRNA_17_K01.gff jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_17/
        rsync -a -P prokka_HDRNA_19/HDRNA_19_K01.gff jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_19/
        rsync -a -P prokka_HDRNA_20/HDRNA_20_K01.gff jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_20/

        cd prokka_HDRNA_01
        roary -p 5 -f ./roary -i 95 -cd 99 -s -e -n -v ./HDRNA_01_K01.gff ./HDRNA_01_K02/HDRNA_01_K02.gff ./HDRNA_01_K03/HDRNA_01_K03.gff ./HDRNA_01_K04/HDRNA_01_K04.gff ./HDRNA_01_K05/HDRNA_01_K05.gff ./HDRNA_01_K06/HDRNA_01_K06.gff ./HDRNA_01_K07/HDRNA_01_K07.gff ./HDRNA_01_K08/HDRNA_01_K08.gff ./HDRNA_01_K09/HDRNA_01_K09.gff ./HDRNA_01_K10/HDRNA_01_K10.gff
        cd ../prokka_HDRNA_03
        roary -p 5 -f ./roary -i 95 -cd 99 -s -e -n -v ./HDRNA_03_K01.gff ./HDRNA_03_K02/HDRNA_03_K02.gff ./HDRNA_03_K03/HDRNA_03_K03.gff ./HDRNA_03_K04/HDRNA_03_K04.gff ./HDRNA_03_K05/HDRNA_03_K05.gff ./HDRNA_03_K06/HDRNA_03_K06.gff ./HDRNA_03_K07/HDRNA_03_K07.gff ./HDRNA_03_K08/HDRNA_03_K08.gff ./HDRNA_03_K09/HDRNA_03_K09.gff ./HDRNA_03_K10/HDRNA_03_K10.gff
        cd ../prokka_HDRNA_06
        roary -p 5 -f ./roary -i 95 -cd 99 -s -e -n -v ./HDRNA_06_K01.gff ./HDRNA_06_K02/HDRNA_06_K02.gff ./HDRNA_06_K03/HDRNA_06_K03.gff ./HDRNA_06_K04/HDRNA_06_K04.gff ./HDRNA_06_K05/HDRNA_06_K05.gff ./HDRNA_06_K06/HDRNA_06_K06.gff ./HDRNA_06_K07/HDRNA_06_K07.gff ./HDRNA_06_K08/HDRNA_06_K08.gff ./HDRNA_06_K09/HDRNA_06_K09.gff ./HDRNA_06_K10/HDRNA_06_K10.gff
        cd ../prokka_HDRNA_07
        roary -p 5 -f ./roary -i 95 -cd 99 -s -e -n -v ./HDRNA_07_K01.gff ./HDRNA_07_K02/HDRNA_07_K02.gff ./HDRNA_07_K03/HDRNA_07_K03.gff ./HDRNA_07_K04/HDRNA_07_K04.gff ./HDRNA_07_K05/HDRNA_07_K05.gff ./HDRNA_07_K06/HDRNA_07_K06.gff ./HDRNA_07_K07/HDRNA_07_K07.gff ./HDRNA_07_K08/HDRNA_07_K08.gff ./HDRNA_07_K09/HDRNA_07_K09.gff ./HDRNA_07_K10/HDRNA_07_K10.gff
        cd ../prokka_HDRNA_08
        roary -p 5 -f ./roary -i 95 -cd 99 -s -e -n -v ./HDRNA_08_K01.gff ./HDRNA_08_K02/HDRNA_08_K02.gff ./HDRNA_08_K03/HDRNA_08_K03.gff ./HDRNA_08_K04/HDRNA_08_K04.gff ./HDRNA_08_K05/HDRNA_08_K05.gff ./HDRNA_08_K06/HDRNA_08_K06.gff ./HDRNA_08_K07/HDRNA_08_K07.gff ./HDRNA_08_K08/HDRNA_08_K08.gff ./HDRNA_08_K09/HDRNA_08_K09.gff ./HDRNA_08_K10/HDRNA_08_K10.gff
        cd ../prokka_HDRNA_12
        roary -p 5 -f ./roary -i 95 -cd 99 -s -e -n -v ./HDRNA_12_K01.gff ./HDRNA_12_K02/HDRNA_12_K02.gff ./HDRNA_12_K03/HDRNA_12_K03.gff ./HDRNA_12_K04/HDRNA_12_K04.gff ./HDRNA_12_K05/HDRNA_12_K05.gff ./HDRNA_12_K06/HDRNA_12_K06.gff ./HDRNA_12_K07/HDRNA_12_K07.gff ./HDRNA_12_K08/HDRNA_12_K08.gff ./HDRNA_12_K09/HDRNA_12_K09.gff ./HDRNA_12_K10/HDRNA_12_K10.gff
        cd ../prokka_HDRNA_16
        roary -p 5 -f ./roary -i 95 -cd 99 -s -e -n -v ./HDRNA_16_K01.gff ./HDRNA_16_K02/HDRNA_16_K02.gff ./HDRNA_16_K03/HDRNA_16_K03.gff ./HDRNA_16_K04/HDRNA_16_K04.gff ./HDRNA_16_K05/HDRNA_16_K05.gff ./HDRNA_16_K06/HDRNA_16_K06.gff ./HDRNA_16_K07/HDRNA_16_K07.gff ./HDRNA_16_K08/HDRNA_16_K08.gff ./HDRNA_16_K09/HDRNA_16_K09.gff ./HDRNA_16_K10/HDRNA_16_K10.gff
        cd ../prokka_HDRNA_17
        roary -p 5 -f ./roary -i 95 -cd 99 -s -e -n -v ./HDRNA_17_K01.gff ./HDRNA_17_K02/HDRNA_17_K02.gff ./HDRNA_17_K03/HDRNA_17_K03.gff ./HDRNA_17_K04/HDRNA_17_K04.gff ./HDRNA_17_K05/HDRNA_17_K05.gff ./HDRNA_17_K06/HDRNA_17_K06.gff ./HDRNA_17_K07/HDRNA_17_K07.gff ./HDRNA_17_K08/HDRNA_17_K08.gff ./HDRNA_17_K09/HDRNA_17_K09.gff ./HDRNA_17_K10/HDRNA_17_K10.gff
        cd ../prokka_HDRNA_19
        roary -p 5 -f ./roary -i 95 -cd 99 -s -e -n -v ./HDRNA_19_K01.gff ./HDRNA_19_K02/HDRNA_19_K02.gff ./HDRNA_19_K03/HDRNA_19_K03.gff ./HDRNA_19_K04/HDRNA_19_K04.gff ./HDRNA_19_K05/HDRNA_19_K05.gff ./HDRNA_19_K06/HDRNA_19_K06.gff ./HDRNA_19_K07/HDRNA_19_K07.gff ./HDRNA_19_K08/HDRNA_19_K08.gff ./HDRNA_19_K09/HDRNA_19_K09.gff ./HDRNA_19_K10/HDRNA_19_K10.gff
        cd ../prokka_HDRNA_20
        roary -p 5 -f ./roary -i 95 -cd 99 -s -e -n -v ./HDRNA_20_K01.gff ./HDRNA_20_K02/HDRNA_20_K02.gff ./HDRNA_20_K03/HDRNA_20_K03.gff ./HDRNA_20_K04/HDRNA_20_K04.gff ./HDRNA_20_K05/HDRNA_20_K05.gff ./HDRNA_20_K06/HDRNA_20_K06.gff ./HDRNA_20_K07/HDRNA_20_K07.gff ./HDRNA_20_K08/HDRNA_20_K08.gff ./HDRNA_20_K09/HDRNA_20_K09.gff
        cd ..

        rsync -a -P jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_01/roary prokka_HDRNA_01
        rsync -a -P jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_03/roary prokka_HDRNA_03
        rsync -a -P jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_06/roary prokka_HDRNA_06
        rsync -a -P jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_07/roary prokka_HDRNA_07
        rsync -a -P jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_08/roary prokka_HDRNA_08
        rsync -a -P jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_12/roary prokka_HDRNA_12
        rsync -a -P jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_16/roary prokka_HDRNA_16
        rsync -a -P jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_17/roary prokka_HDRNA_17
        rsync -a -P jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_19/roary prokka_HDRNA_19
        rsync -a -P jhuang@hamm:~/DATA/Data_Holger_S.epidermidis_short/prokka_HDRNA_20/roary prokka_HDRNA_20

        cp prokka_HDRNA_01/roary/gene_presence_absence.csv gene_presence_absence__HDRNA_01.csv
        cp prokka_HDRNA_03/roary/gene_presence_absence.csv gene_presence_absence_HDRNA_03.csv
        cp prokka_HDRNA_06/roary/gene_presence_absence.csv gene_presence_absence_HDRNA_06.csv
        cp prokka_HDRNA_07/roary/gene_presence_absence.csv gene_presence_absence_HDRNA_07.csv
        cp prokka_HDRNA_08/roary/gene_presence_absence.csv gene_presence_absence_HDRNA_08.csv
        cp prokka_HDRNA_12/roary/gene_presence_absence.csv gene_presence_absence_HDRNA_12.csv
        cp prokka_HDRNA_16/roary/gene_presence_absence.csv gene_presence_absence_HDRNA_16.csv
        cp prokka_HDRNA_17/roary/gene_presence_absence.csv gene_presence_absence_HDRNA_17.csv
        cp prokka_HDRNA_19/roary/gene_presence_absence.csv gene_presence_absence_HDRNA_19.csv
        cp prokka_HDRNA_20/roary/gene_presence_absence.csv gene_presence_absence_HDRNA_20.csv

        #Wenn man will open the files mit libreoffice, needs "," --> "|"; "\n"-->\n; the first and last " in the text; in the kate, then open file with libreoffice with delimiter '|'.

        sed -i 's/\.mRNA\.0\.CDS\.1//g' gene_presence_absence__HDRNA_01.csv

        group_10||Y_phosphoryl: pyrimidine-nucleoside phosphorylase|10|11|1.11|1|416||||362|1301|1165|RE430_03860|HDRNA_01_K02_02111|HDRNA_01_K03_02181|HDRNA_01_K04_02110|HDRNA_01_K05_02181|HDRNA_01_K06_02156|HDRNA_01_K07_02196     HDRNA_01_K07_02197|HDRNA_01_K08_01754|HDRNA_01_K09_02226|HDRNA_01_K10_02107

        group_10||Y_phosphoryl: pyrimidine-nucleoside phosphorylase|10|11|1.11|1|416||||362|1301|1165|RE430_03860.mRNA.0.CDS.1|HDRNA_01_K02_02111|HDRNA_01_K03_02181|HDRNA_01_K04_02110|HDRNA_01_K05_02181|HDRNA_01_K06_02156|HDRNA_01_K07_02196    HDRNA_01_K07_02197|HDRNA_01_K08_01754|HDRNA_01_K09_02226|HDRNA_01_K10_02107

        for file in gene_presence_absence_HDRNA_*; do
            sed -i 's/\.mRNA\.0\.CDS\.1//g' "$file"
        done

        ~/Tools/csv2xls-0.4/csv_to_xls.py gene_presence_absence__HDRNA_01.csv gene_presence_absence_HDRNA_03.csv gene_presence_absence_HDRNA_06.csv gene_presence_absence_HDRNA_07.csv gene_presence_absence_HDRNA_08.csv gene_presence_absence_HDRNA_12.csv gene_presence_absence_HDRNA_16.csv gene_presence_absence_HDRNA_17.csv gene_presence_absence_HDRNA_19.csv gene_presence_absence_HDRNA_20.csv -d'|' -o gene_presence_absence.xls

        * genomic rearrangements (e.g. SCCmec deletions, ACME deletions, agr insertions)
          TODOs using Easyfig!
          #Staphylococcal Cassette Chromosome mec
          #arginine catabolic mobile element (ACME)

        #Prevalence and genetic diversity of arginine catabolic mobile element (ACME) in clinical isolates of coagulase-negative staphylococci: identification of ACME type I variants in Staphylococcus epidermidis.
        Fig. 1. A schematic drawing of genetic structures of ACME (a region from the arc to opp3 cluster, or corresponding genetic components) among the three DI subtypes (DI.1, DI.2, and DI.3: strains CNS266, CNS115, and CNS149, respectively), type I (strain USA300-FPR3757, accession number CP000255), type II (strain ATCC12228, accession number AE015929) and type DII (strain M08/0126, accession number FR753166). Putative ORFs of genes are represented by arrows colored with green (arc cluster), red (opp3 cluster), blue (a region between the arc and opp3 clusters in ACME I), or dark blue (genes in ACME II). The regions in light pink including the arc cluster indicate genetically identical areas to both ATCC12228 and USA300-FPR3757. The regions with light blue are identical to only ATCC12228, while those with light orange to USA300-FPR3757. White space regions between argR and SAUSA300_0072 show no sequence homology either to ATCC12228 or to USA300_FPR3757; however, these regions show 91–99% nucleotide identity among the three ACME subtypes. Regions colored with dark orange in the three ACME DI subtypes show=98% nucleotide sequence identity to each other. Regions colored with grey (type DI.1), purple or cyan (type DI.3) do not show high nucleotide identity (<98%) to cognate genes in other ACME types (Table S2.2). Positions of primers used for PCR profile (Tables 1 and 4) are shown with arrowheads under ACME I sequence. Collapse

          #Smash++: https://academic.oup.com/gigascience/article/9/5/giaa048/5841055
          #https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7597632/

        Artemis Comparison Tool (ACT): Allows for the visual comparison of genomes and can be used to investigate the presence or absence of genomic regions (such as SCCmec or ACME) and other structural variations.
        #https://journals.plos.org/plospathogens/article?id=10.1371/journal.ppat.1009304
        #https://journals.plos.org/plospathogens/article/figure?id=10.1371/journal.ppat.1009304.g010

        The image you uploaded appears to be a schematic representation of genomic rearrangements, typically found in scientific publications or reports. These kinds of images are often created using bioinformatics visualization tools or general graphic design software. While it's not possible to determine the exact software used to create this specific image without more context, I can suggest several tools commonly used for such purposes:

    3. Epidome: (using R?) Create bar plot showing ST distribution in noses from patients HDRNA 1, 3, 6,   --> 7, 8, 12, 16, 17, 19, 20 <--
    (Easy using R) Create table showing presence / absence of STs / per patient
    NOTE: Epidome data have not been processed due to potential missing of the data; At first sending the results without Epidome results, ask again where is the epidome data for the 10 patients?

generate goeBURST

5,1,1,1,2,2,1,1
87,7,1,1,2,2,1,1
35,2,1,2,2,4,1,1
69,1,18,6,2,2,1,1

23,7,1,2,1,3,3,1
224,19,16,19,6,3,19,10
640,28,3,13,5,8,9,11

'-',1,13,2,1,2,1,29

in goeBURST-1.2.1.jar, I have gelesen "Edge level to define group SLV DLV TLV". Wie kann ich generate SLV, DLV and TLV files?

The input you've provided seems to represent allelic profiles or sequence types (STs) used in microbial typing, particularly in methods like Multilocus Sequence Typing (MLST). When you input these profiles into software like goeBURST (implemented in a tool like PHYLOViZ), the software uses these profiles to construct a phylogenetic network. The network shows relationships between different strains or isolates based on their allelic similarity.

In goeBURST, groups are defined based on their allelic differences:

SLV (Single-Locus Variants): Strains or STs differing by only one locus.
DLV (Double-Locus Variants): Strains or STs differing by two loci.
TLV (Triple-Locus Variants): Strains or STs differing by three loci.

agr typing S. epidermidis

> Comparison of the amino acid sequences of a region of the N-terminus of AgrC (A), and AgrB (B) of S. epidermidis (S.e.), S. aureus (S.a.) and S. lugdunensis (S.l.).
> https://academic.oup.com/femsle/article/163/1/1/625220: Cloning and characterization of an accessory gene regulator (agr)-like locus from Staphylococcus epidermidis

https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5064449/
https://brieflands.com/articles/archcid-62833#4.-Results

https://www.ncbi.nlm.nih.gov/Taxonomy/Browser/wwwtax.cgi?id=1282
Staphylococcus epidermidis Agr Operon
https://academic.oup.com/femspd/article/51/1/220/501159

#pAgr (S. epidermidis agr operon promoter)
https://parts.igem.org/Part:BBa_K212003#

- An overall homology of 68% was found between the agr locus from S. epidermidis and S. aureus.
- The agr locus from S. epidermidis was organized similar to those from S. aureus and S. lugdunensis.
- The putative RNAII molecule contains four open reading frames, agrA, B, C and D. AgrA was a response regulator.
- AgrB showed homology with transducer and translocase molecules.
- AgrC is expected to act as a histidine protein kinase in which a leucine zipper is present.
- AgrD is presumably processed into an autoinducer peptide.

- For a long time, Staphylococcus epidermidis, as a member of the coagulase-negative staphylococci, was considered as part of the physiological skin flora of the human being with no pathogenic significance.
- Today, we know that S. epidermidis is one of the most prevalent causes for implant-associated and nosocomial infections.
- We performed pheno- and genotypic analysis (ica, IS256, SCCmec types, agr groups) of biofilm formation in 200 isolates.
- Fifty percent were genetically ica-positive and produced biofilm.
- Among all studied isolates, agr II and III and SCCmec type I were the most prevalent, whereas within the selected multi-resistant isolates (29%), agr I and III and SCCmec type II dominated.
- SCCmec type I and mecA-negative S. epidermidis isolates were associated with agr II.
- The majority of the blood culture and biopsy isolates were assigned to agr III and SCCmec type I, whereas agr II was predominantly detected in mecA-negative S. epidermidis isolated from catheter and implant materials.
- MLST analysis revealed the major clonal lineages of ST2, ST5, ST10, and ST242 (total 13 STs).
- ST2 isolates from blood cultures were icaA/D-positive and harbored SCCmec types II and III and IS256, whereas the icaA/D- and IS256-positive ST23 isolates were assigned to SCCmec types I and IV.

# -- OPTION 1 using AgrVATE gp1-4-operon_ref.fasta (failed since the database contains only agr from Staphylococcus aureus) --

> Species-Wide Phylogenomics of the Staphylococcus aureus Agr Operon Revealed Convergent Evolution of Frameshift Mutations

makeblastdb -in HDRNA_K01.fna -dbtype nucl
-perc_identity 90 -qcov_hsp_perc 90
blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query /home/jhuang/Tools/AgrVATE/agrvate_databases/references/gp1-4-operon_ref.fasta -evalue 0.1 -num_threads 15 -outfmt "6 sseqid qseqid evalue pident sstart send" -strand both -max_target_seqs 1

blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query gp1-4-operon_ref.fasta -evalue 0.1 -num_threads 15 -outfmt "6 sseqid qseqid evalue pident sstart send" -strand both > HDRNA_01.blastn
blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query /home/jhuang/Tools/AgrVATE/agrvate_databases/references/gp1-4-operon_ref.fasta -evalue 0.1 -num_threads 15 -outfmt "6 sseqid qseqid evalue pident sstart send" -strand both > HDRNA_03.blastn

blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query /home/jhuang/Tools/AgrVATE/agrvate_databases/gp1234_all_motifs.fna -evalue 0.1 -num_threads 15 -outfmt "6 sseqid qseqid evalue pident sstart send" -strand both > HDRNA_01_.blastn
blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query /home/jhuang/Tools/AgrVATE/agrvate_databases/gp1234_all_motifs.fna -evalue 0.1 -num_threads 15 -outfmt "6 sseqid qseqid evalue pident sstart send" -strand both > HDRNA_03_.blastn

# -- OPTION 2 using agr-1-3_Se.fasta (failed!) --

> High Genetic Variablity of the agr Locus in Staphylococcus Species

>gi|3320006|emb|Z49220.1| Staphylococcus epidermidis hld and agr[A,B,C,D] genes

>gi|18251022|gb|AF346724.1| Staphylococcus epidermidis strain N910160 AgrB (agrB) gene, partial cds; AgrD (agrD) gene, complete cds; and AgrC (agrC) gene, partial cds

>gi|18251026|gb|AF346725.1| Staphylococcus epidermidis strain N910191 AgrB (agrB) gene, partial cds; AgrD (agrD) gene, complete cds; and AgrC (agrC) gene, partial cds

blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query agr-1-3_Se.fasta -evalue 0.1 -num_threads 15 -outfmt 6 -strand both > HDRNA_01__.blastn

# -- OPTION 3 using agrD_I-III.fasta (successful!) --

> Staphylococcus epidermidis agr Quorum-Sensing System: Signal Identification, Cross Talk, and Importance in Colonization

> https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4187671/

> https://www.researchgate.net/publication/264391883

tblastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna  -query agrD_I-III.fasta -evalue 0.1 -num_threads 15  > HDRNA_01.tblastn
tblastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna  -query agrD_I-III.fasta -evalue 0.1 -num_threads 15  > HDRNA_03.tblastn
tblastn -db HDRNA_06_K01_conservative_27645.current.gb_converted.fna  -query agrD_I-III.fasta -evalue 0.1 -num_threads 15  > HDRNA_06.tblastn
tblastn -db HDRNA_07_K01_conservative_27169.current.gb_converted.fna  -query agrD_I-III.fasta -evalue 0.1 -num_threads 15  > HDRNA_07.tblastn
tblastn -db HDRNA_08_K01_conservative_32455.current.gb_converted.fna  -query agrD_I-III.fasta -evalue 0.1 -num_threads 15  > HDRNA_08.tblastn
tblastn -db HDRNA_12_K01_bold_37467.current.gb_converted.fna  -query agrD_I-III.fasta -evalue 0.1 -num_threads 15  > HDRNA_12.tblastn
tblastn -db HDRNA_16_K01_conservative_37834.current.gb_converted.fna  -query agrD_I-III.fasta -evalue 0.1 -num_threads 15  > HDRNA_16.tblastn
tblastn -db HDRNA_17_K01_conservative_37288.current.gb_converted.fna  -query agrD_I-III.fasta -evalue 0.1 -num_threads 15  > HDRNA_17.tblastn
tblastn -db HDRNA_19_K01_bold_37377.current.gb_converted.fna  -query agrD_I-III.fasta -evalue 0.1 -num_threads 15  > HDRNA_19.tblastn
tblastn -db HDRNA_20_K01_conservative_43457.current.gb_converted.fna  -query agrD_I-III.fasta -evalue 0.1 -num_threads 15  > HDRNA_20.tblastn

II?
II
I
II
I

II
II
II
III
II

>AgrD_I
MENIFNLFIKFFTTILEFIGTVAGDSVCASYFDEPEVPEELTKLYE
>AgrD_II
MNLLGGLLLKIFSNFMAVIGNASKYNPCSNYLDEPQVPEELTKLDE
>AgrD_III
MNLLGGLLLKLFSNFMAVIGNAAKYNPCASYLDEPQVPEELTKLDE

Variant (SNP+INDEL) calling

Input files:
    HDRNA_01_K01_conservative_23197.current.gb

    HDRNA_01_K01_conservative_23197.current.gb:LOCUS       CP133676             2502964 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_03_K01_bold_bandage_26831.current.gb:LOCUS       CP133677             2590275 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_06_K01_conservative_27645.current.gb:LOCUS       CP133678             2465260 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_06_K01_conservative_27645.current.gb:LOCUS       CP133679               19348 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_07_K01_conservative_27169.current.gb:LOCUS       CP133680             2544074 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_07_K01_conservative_27169.current.gb:LOCUS       CP133681                2241 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_08_K01_conservative_32455.current.gb:LOCUS       CP133682             2425353 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_08_K01_conservative_32455.current.gb:LOCUS       CP133683                6358 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_12_K01_bold_37467.current.gb:LOCUS       CP133684             2490139 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_12_K01_bold_37467.current.gb:LOCUS       CP133685               43849 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_12_K01_bold_37467.current.gb:LOCUS       CP133686                6642 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_12_K01_bold_37467.current.gb:LOCUS       CP133687                2241 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_16_K01_conservative_37834.current.gb:LOCUS       CP133688             2626440 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_16_K01_conservative_37834.current.gb:LOCUS       CP133689               24906 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_16_K01_conservative_37834.current.gb:LOCUS       CP133690               20828 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_16_K01_conservative_37834.current.gb:LOCUS       CP133691                4567 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_16_K01_conservative_37834.current.gb:LOCUS       CP133692                2242 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_17_K01_conservative_37288.current.gb:LOCUS       CP133693             2503072 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_17_K01_conservative_37288.current.gb:LOCUS       CP133694               29861 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_17_K01_conservative_37288.current.gb:LOCUS       CP133695                4439 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_19_K01_bold_37377.current.gb:LOCUS       CP133696             2362062 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_19_K01_bold_37377.current.gb:LOCUS       CP133697               55320 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_19_K01_bold_37377.current.gb:LOCUS       CP133698               46464 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_19_K01_bold_37377.current.gb:LOCUS       CP133699               11827 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_20_K01_conservative_43457.current.gb:LOCUS       CP133700             2490778 bp    DNA     circular BCT 05-SEP-2023
    HDRNA_20_K01_conservative_43457.current.gb:LOCUS       CP133701                2241 bp    DNA     circular BCT 05-SEP-2023

ln -s /home/jhuang/Tools/spandx/ spandx

nextflow run spandx/main.nf --fastq "trimmed_HDRNA_01/*_P_{1,2}.fastq.gz" --ref db/CP133676.fasta --annotation --database CP133676 -resume
mv work work_CP133676
mv Outputs Outputs_CP133676

for fasta_file in CP133677 CP133678 CP133679 CP133680 CP133681 CP133682 CP133683 CP133684 CP133685 CP133686 CP133687 CP133688 CP133689 CP133690 CP133691 CP133692 CP133693 CP133694 CP133695 CP133696 CP133697 CP133698 CP133699 CP133700 CP133701; do
  echo "nextflow run spandx/main.nf --fastq "trimmed_HDRNA_01/*_P_{1,2}.fastq.gz" --ref db/${fasta_file}.fasta --annotation --database ${fasta_file} -resume"
  echo "mv work work_${fasta_file}"
  echo "mv Outputs Outputs_${fasta_file}"
done

for file in *.fastq.gz; do mv $file $(echo $file | cut -d'_' -f1)-$(echo $file | cut -d'_' -f1)-$(echo $file | cut -d'_' -f3)_$(echo $file | cut -d'_' -f6); done

for file in *.fastq.gz; do mv $file $(echo $file | cut -d'_' -f3)_$(echo $file | cut -d'_' -f6); done
for file in *.fastq.gz; do mv $file $(echo $file | cut -d'_' -f1)-$(echo $file | cut -d'_' -f2); done
for file in *.fastq.gz; do mv $file $(echo $file | cut -d'-' -f1)_$(echo $file | cut -d'-' -f2); done

Input read files could not be found.
Have you included the read files in the current directory and do they have the correct naming?
With the parameters specified, SPANDx is looking for reads named *_{1,2}.fastq.gz.
To fix this error either rename your reads to match this formatting or specify the desired format
when initializing SPANDx e.g. --fastq "*_{1,2}_sequence.fastq.gz"

cd trimmed_HDRNA_01
nextflow run ../spandx/main.nf --ref ../db/CP133676.fasta --annotation --database CP133676 -resume
mv Outputs Outputs_CP133676
cd ..

cd trimmed_HDRNA_03
nextflow run ../spandx/main.nf --ref ../db/CP133677.fasta --annotation --database CP133677 -resume
mv Outputs Outputs_CP133677
cd ..

cd trimmed_HDRNA_06
nextflow run ../spandx/main.nf --ref ../db/CP133678.fasta --annotation --database CP133678 -resume
mv work work_CP133678
mv Outputs Outputs_CP133678
nextflow run ../spandx/main.nf --ref ../db/CP133679.fasta --annotation --database CP133679 -resume
mv work work_CP133679
mv Outputs Outputs_CP133679
cd ..

cd trimmed_HDRNA_07
nextflow run ../spandx/main.nf --ref ../db/CP133680.fasta --annotation --database CP133680 -resume
mv work work_CP133680
mv Outputs Outputs_CP133680
nextflow run ../spandx/main.nf --ref ../db/CP133681.fasta --annotation --database CP133681 -resume
mv work work_CP133681
mv Outputs Outputs_CP133681
cd ..

cd trimmed_HDRNA_08
nextflow run ../spandx/main.nf --ref ../db/CP133682.fasta --annotation --database CP133682 -resume
mv work work_CP133682
mv Outputs Outputs_CP133682
nextflow run ../spandx/main.nf --ref ../db/CP133683.fasta --annotation --database CP133683 -resume
mv work work_CP133683
mv Outputs Outputs_CP133683
cd ..

cd trimmed_HDRNA_12
nextflow run ../spandx/main.nf --ref ../db/CP133684.fasta --annotation --database CP133684 -resume
mv work work_CP133684
mv Outputs Outputs_CP133684
nextflow run ../spandx/main.nf --ref ../db/CP133685.fasta --annotation --database CP133685 -resume
mv work work_CP133685
mv Outputs Outputs_CP133685
nextflow run ../spandx/main.nf --ref ../db/CP133686.fasta --annotation --database CP133686 -resume
mv work work_CP133686
mv Outputs Outputs_CP133686
nextflow run ../spandx/main.nf --ref ../db/CP133687.fasta --annotation --database CP133687 -resume
mv work work_CP133687
mv Outputs Outputs_CP133687
cd ..

cd trimmed_HDRNA_16
nextflow run ../spandx/main.nf --ref ../db/CP133688.fasta --annotation --database CP133688 -resume
mv work work_CP133688
mv Outputs Outputs_CP133688
nextflow run ../spandx/main.nf --ref ../db/CP133689.fasta --annotation --database CP133689 -resume
mv work work_CP133689
mv Outputs Outputs_CP133689
nextflow run ../spandx/main.nf --ref ../db/CP133690.fasta --annotation --database CP133690 -resume
mv work work_CP133690
mv Outputs Outputs_CP133690
nextflow run ../spandx/main.nf --ref ../db/CP133691.fasta --annotation --database CP133691 -resume
mv work work_CP133691
mv Outputs Outputs_CP133691
nextflow run ../spandx/main.nf --ref ../db/CP133692.fasta --annotation --database CP133692 -resume
mv work work_CP133692
mv Outputs Outputs_CP133692
cd ..

cd trimmed_HDRNA_17
nextflow run ../spandx/main.nf --ref ../db/CP133693.fasta --annotation --database CP133693 -resume
mv work work_CP133693
mv Outputs Outputs_CP133693
nextflow run ../spandx/main.nf --ref ../db/CP133694.fasta --annotation --database CP133694 -resume
mv work work_CP133694
mv Outputs Outputs_CP133694
nextflow run ../spandx/main.nf --ref ../db/CP133695.fasta --annotation --database CP133695 -resume
mv work work_CP133695
mv Outputs Outputs_CP133695
cd ..

cd trimmed_HDRNA_19
nextflow run ../spandx/main.nf --ref ../db/CP133696.fasta --annotation --database CP133696 -resume
mv work work_CP133696
mv Outputs Outputs_CP133696
nextflow run ../spandx/main.nf --ref ../db/CP133697.fasta --annotation --database CP133697 -resume
mv work work_CP133697
mv Outputs Outputs_CP133697
nextflow run ../spandx/main.nf --ref ../db/CP133698.fasta --annotation --database CP133698 -resume
mv work work_CP133698
mv Outputs Outputs_CP133698
nextflow run ../spandx/main.nf --ref ../db/CP133699.fasta --annotation --database CP133699 -resume
mv work work_CP133699
mv Outputs Outputs_CP133699
cd ..

cd trimmed_HDRNA_20
nextflow run ../spandx/main.nf --ref ../db/CP133700.fasta --annotation --database CP133700 -resume
mv work work_CP133700
mv Outputs Outputs_CP133700
nextflow run ../spandx/main.nf --ref ../db/CP133701.fasta --annotation --database CP133701 -resume
mv work work_CP133701
mv Outputs Outputs_CP133701
cd ..

#-------------------------

rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_01/Outputs_CP133676 .

rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_03/Outputs_CP133677 .

rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_06/Outputs_CP133678 .
rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_06/Outputs_CP133679 .

rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_07/Outputs_CP133680 .
rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_07/Outputs_CP133681 .

rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_08/Outputs_CP133682 .
rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_08/Outputs_CP133683 .

rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_12/Outputs_CP133684 .
rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_12/Outputs_CP133685 .
rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_12/Outputs_CP133686 .
rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_12/Outputs_CP133687 .

rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_16/Outputs_CP133688 .
rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_16/Outputs_CP133689 .
rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_16/Outputs_CP133690 .
rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_16/Outputs_CP133691 .
rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_16/Outputs_CP133692 .

rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_17/Outputs_CP133693 .
rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_17/Outputs_CP133694 .
rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_17/Outputs_CP133695 .

rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_19/Outputs_CP133696 .
rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_19/Outputs_CP133697 .
rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_19/Outputs_CP133698 .
rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_19/Outputs_CP133699 .

rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_20/Outputs_CP133700 .
rsync -a -P jhuang@hamm:/home/jhuang/DATA/Data_Holger_S.epidermidis_short/trimmed_HDRNA_20/Outputs_CP133701 .

cut -f2 -d$'\t' snippy.core.tab > f2
cut -f3 -d$'\t' snippy.core.tab > f3
cut -f4 -d$'\t' snippy.core.tab > f4

diff snippy/merged_snp.vcf.id variants/f2
5d5
< 138824
7d6
< 139197
9d7
< 139844
61d58
< 2475573

# -- merging vcf-files using bcftools --

results_HDRNA_01/snippy
bcftools merge HDRNA_01_K01/HDRNA_01_K01.vcf.gz HDRNA_01_K02/HDRNA_01_K02.vcf.gz HDRNA_01_K03/HDRNA_01_K03.vcf.gz HDRNA_01_K04/HDRNA_01_K04.vcf.gz HDRNA_01_K05/HDRNA_01_K05.vcf.gz HDRNA_01_K06/HDRNA_01_K06.vcf.gz HDRNA_01_K07/HDRNA_01_K07.vcf.gz HDRNA_01_K08/HDRNA_01_K08.vcf.gz HDRNA_01_K09/HDRNA_01_K09.vcf.gz HDRNA_01_K10/HDRNA_01_K10.vcf.gz -o merged.vcf
#bcftools index merged.vcf.gz
cp merged.vcf merged_CP133676.vcf

cd results_HDRNA_03/snippy
bcftools merge HDRNA_03_K01/HDRNA_03_K01.vcf.gz HDRNA_03_K02/HDRNA_03_K02.vcf.gz HDRNA_03_K03/HDRNA_03_K03.vcf.gz HDRNA_03_K04/HDRNA_03_K04.vcf.gz HDRNA_03_K05/HDRNA_03_K05.vcf.gz HDRNA_03_K06/HDRNA_03_K06.vcf.gz HDRNA_03_K07/HDRNA_03_K07.vcf.gz HDRNA_03_K08/HDRNA_03_K08.vcf.gz HDRNA_03_K09/HDRNA_03_K09.vcf.gz HDRNA_03_K10/HDRNA_03_K10.vcf.gz -o merged.vcf
cp merged.vcf merged_CP133677.vcf

cd results_HDRNA_06/snippy
bcftools merge HDRNA_06_K01/HDRNA_06_K01.vcf.gz HDRNA_06_K02/HDRNA_06_K02.vcf.gz HDRNA_06_K03/HDRNA_06_K03.vcf.gz HDRNA_06_K04/HDRNA_06_K04.vcf.gz HDRNA_06_K05/HDRNA_06_K05.vcf.gz HDRNA_06_K06/HDRNA_06_K06.vcf.gz HDRNA_06_K07/HDRNA_06_K07.vcf.gz HDRNA_06_K08/HDRNA_06_K08.vcf.gz HDRNA_06_K09/HDRNA_06_K09.vcf.gz HDRNA_06_K10/HDRNA_06_K10.vcf.gz -o merged.vcf
#split merged.vcf to merged_CP133678.vcf and merged_CP133679.vcf

cd results_HDRNA_07/snippy
bcftools merge HDRNA_07_K01/HDRNA_07_K01.vcf.gz HDRNA_07_K02/HDRNA_07_K02.vcf.gz HDRNA_07_K03/HDRNA_07_K03.vcf.gz HDRNA_07_K04/HDRNA_07_K04.vcf.gz HDRNA_07_K05/HDRNA_07_K05.vcf.gz HDRNA_07_K06/HDRNA_07_K06.vcf.gz HDRNA_07_K07/HDRNA_07_K07.vcf.gz HDRNA_07_K08/HDRNA_07_K08.vcf.gz HDRNA_07_K09/HDRNA_07_K09.vcf.gz HDRNA_07_K10/HDRNA_07_K10.vcf.gz -o merged.vcf
cp merged.vcf merged_CP133680.vcf
#Note that merged_CP133681.vcf is empty.

cd results_HDRNA_08/snippy
bcftools merge HDRNA_08_K01/HDRNA_08_K01.vcf.gz HDRNA_08_K02/HDRNA_08_K02.vcf.gz HDRNA_08_K03/HDRNA_08_K03.vcf.gz HDRNA_08_K04/HDRNA_08_K04.vcf.gz HDRNA_08_K05/HDRNA_08_K05.vcf.gz HDRNA_08_K06/HDRNA_08_K06.vcf.gz HDRNA_08_K07/HDRNA_08_K07.vcf.gz HDRNA_08_K08/HDRNA_08_K08.vcf.gz HDRNA_08_K09/HDRNA_08_K09.vcf.gz HDRNA_08_K10/HDRNA_08_K10.vcf.gz -o merged.vcf
#split merged.vcf to merged_CP133682.vcf and merged_CP133683.vcf.
#----ERROR: IGNORING the record---->
#CP133683        1718    .       G       A       155736  .       QR=0;RO=0;ANN=A|intron_variant|MODIFIER|RGR12_11570|GENE_RGR12_11570|transcript|TRANSCRIPT_RGR12_11570|protein_coding|1/1|c.358-506C>T||||||WARNING_TRANSCRIPT_NO_START_CODON;DP=12890;AB=0;AO=3768;QA=121701;TYPE=snp  GT:DP:RO:QR:AO:QA:GL    ./.:.:.:.:.:.:. 1/1:5361:0:0:5347:173706:-15627.9,-1609.61,0    1/1:3750:0:0:3747:121191:-10903.3,-1127.96,0    ./.:.:.:.:.:.:. 1/1:3779:0:0:3768:121701:-10949.3,-1134.28,0    ./.:.:.:.:.:.:. ./.:.:.:.:.:.:. ./.:.:.:.:.:.:. ./.:.:.:.:.:.:. ./.:.:.:.:.:.:.

cd results_HDRNA_12/snippy
bcftools merge HDRNA_12_K01/HDRNA_12_K01.vcf.gz HDRNA_12_K02/HDRNA_12_K02.vcf.gz HDRNA_12_K03/HDRNA_12_K03.vcf.gz HDRNA_12_K04/HDRNA_12_K04.vcf.gz HDRNA_12_K05/HDRNA_12_K05.vcf.gz HDRNA_12_K06/HDRNA_12_K06.vcf.gz HDRNA_12_K07/HDRNA_12_K07.vcf.gz HDRNA_12_K08/HDRNA_12_K08.vcf.gz HDRNA_12_K09/HDRNA_12_K09.vcf.gz HDRNA_12_K10/HDRNA_12_K10.vcf.gz -o merged.vcf
#split merged.vcf to merged_CP133684.vcf and merged_CP133685.vcf.
#Note that merged_CP133686.vcf and merged_CP133687.vcf are empty.

cd results_HDRNA_16/snippy
bcftools merge HDRNA_16_K01/HDRNA_16_K01.vcf.gz HDRNA_16_K02/HDRNA_16_K02.vcf.gz HDRNA_16_K03/HDRNA_16_K03.vcf.gz HDRNA_16_K04/HDRNA_16_K04.vcf.gz HDRNA_16_K05/HDRNA_16_K05.vcf.gz HDRNA_16_K06/HDRNA_16_K06.vcf.gz HDRNA_16_K07/HDRNA_16_K07.vcf.gz HDRNA_16_K08/HDRNA_16_K08.vcf.gz HDRNA_16_K09/HDRNA_16_K09.vcf.gz HDRNA_16_K10/HDRNA_16_K10.vcf.gz -o merged.vcf
cp merged.vcf merged_CP133688.vcf
#Note that merged_CP133689.vcf - merged_CP133692.vcf are empty.

cd results_HDRNA_17/snippy
bcftools merge HDRNA_17_K01/HDRNA_17_K01.vcf.gz HDRNA_17_K02/HDRNA_17_K02.vcf.gz HDRNA_17_K03/HDRNA_17_K03.vcf.gz HDRNA_17_K04/HDRNA_17_K04.vcf.gz HDRNA_17_K05/HDRNA_17_K05.vcf.gz HDRNA_17_K06/HDRNA_17_K06.vcf.gz HDRNA_17_K07/HDRNA_17_K07.vcf.gz HDRNA_17_K08/HDRNA_17_K08.vcf.gz HDRNA_17_K09/HDRNA_17_K09.vcf.gz HDRNA_17_K10/HDRNA_17_K10.vcf.gz -o merged.vcf
cp merged.vcf merged_CP133693.vcf
#Note that merged_CP133694.vcf - merged_CP133695.vcf are empty.

cd results_HDRNA_19/snippy
bcftools merge HDRNA_19_K01/HDRNA_19_K01.vcf.gz HDRNA_19_K02/HDRNA_19_K02.vcf.gz HDRNA_19_K03/HDRNA_19_K03.vcf.gz HDRNA_19_K04/HDRNA_19_K04.vcf.gz HDRNA_19_K05/HDRNA_19_K05.vcf.gz HDRNA_19_K06/HDRNA_19_K06.vcf.gz HDRNA_19_K07/HDRNA_19_K07.vcf.gz HDRNA_19_K08/HDRNA_19_K08.vcf.gz HDRNA_19_K09/HDRNA_19_K09.vcf.gz HDRNA_19_K10/HDRNA_19_K10.vcf.gz -o merged.vcf
cp merged.vcf merged_CP133696.vcf
#Note that merged_CP133697.vcf - merged_CP133699.vcf are empty.

cd results_HDRNA_20/snippy
bcftools merge HDRNA_20_K01/HDRNA_20_K01.vcf.gz HDRNA_20_K02/HDRNA_20_K02.vcf.gz HDRNA_20_K03/HDRNA_20_K03.vcf.gz HDRNA_20_K04/HDRNA_20_K04.vcf.gz HDRNA_20_K05/HDRNA_20_K05.vcf.gz HDRNA_20_K06/HDRNA_20_K06.vcf.gz HDRNA_20_K07/HDRNA_20_K07.vcf.gz HDRNA_20_K08/HDRNA_20_K08.vcf.gz HDRNA_20_K09/HDRNA_20_K09.vcf.gz -o merged.vcf
cp merged.vcf merged_CP133700.vcf
#Note that merged_CP133701.vcf is empty.
#Note not enough reads exist in 'HDRNA_20 S K10 365_BB84_S87_R1_001.fastq.gz' and 'HDRNA_20 S K10 365_BB84_S87_R2_001.fastq.gz'.

python3 ~/Scripts/merge_snippy_and_spandx_results.py results_HDRNA_01/snippy/merged_CP133676.vcf Outputs_CP133676/Phylogeny_and_annotation/All_SNPs_indels_annotated.txt HDRNA_01_not_in_vcf_file_output.txt HDRNA_01_not_in_txt_file_output.txt HDRNA_01_common_records_output.txt
python3 ~/Scripts/merge_snippy_and_spandx_results.py results_HDRNA_03/snippy/merged_CP133677.vcf Outputs_CP133677/Phylogeny_and_annotation/All_SNPs_indels_annotated.txt HDRNA_03_not_in_vcf_file_output.txt HDRNA_03_not_in_txt_file_output.txt HDRNA_03_common_records_output.txt

python3 ~/Scripts/merge_snippy_and_spandx_results.py results_HDRNA_06/snippy/merged_CP133678.vcf Outputs_CP133678/Phylogeny_and_annotation/All_SNPs_indels_annotated.txt HDRNA_06_1_not_in_vcf_file_output.txt HDRNA_06_1_not_in_txt_file_output.txt HDRNA_06_1_common_records_output.txt
python3 ~/Scripts/merge_snippy_and_spandx_results.py results_HDRNA_06/snippy/merged_CP133679.vcf Outputs_CP133679/Phylogeny_and_annotation/All_SNPs_indels_annotated.txt HDRNA_06_2_not_in_vcf_file_output.txt HDRNA_06_2_not_in_txt_file_output.txt HDRNA_06_2_common_records_output.txt

python3 ~/Scripts/merge_snippy_and_spandx_results.py results_HDRNA_07/snippy/merged_CP133680.vcf Outputs_CP133680/Phylogeny_and_annotation/All_SNPs_indels_annotated.txt HDRNA_07_not_in_vcf_file_output.txt HDRNA_07_not_in_txt_file_output.txt HDRNA_07_common_records_output.txt

python3 ~/Scripts/merge_snippy_and_spandx_results.py results_HDRNA_08/snippy/merged_CP133682.vcf Outputs_CP133682/Phylogeny_and_annotation/All_SNPs_indels_annotated.txt HDRNA_08_not_in_vcf_file_output.txt HDRNA_08_not_in_txt_file_output.txt HDRNA_08_common_records_output.txt
#ERROR: python3 ~/Scripts/merge_snippy_and_spandx_results.py results_HDRNA_08/snippy/merged_CP133683.vcf Outputs_CP133683/Phylogeny_and_annotation/All_SNPs_indels_annotated.txt HDRNA_08_2_not_in_vcf_file_output.txt HDRNA_08_2_not_in_txt_file_output.txt HDRNA_08_2_common_records_output.txt

python3 ~/Scripts/merge_snippy_and_spandx_results.py results_HDRNA_12/snippy/merged_CP133684.vcf Outputs_CP133684/Phylogeny_and_annotation/All_SNPs_indels_annotated.txt HDRNA_12_1_not_in_vcf_file_output.txt HDRNA_12_1_not_in_txt_file_output.txt HDRNA_12_1_common_records_output.txt
python3 ~/Scripts/merge_snippy_and_spandx_results.py results_HDRNA_12/snippy/merged_CP133685.vcf Outputs_CP133685/Phylogeny_and_annotation/All_SNPs_indels_annotated.txt HDRNA_12_2_not_in_vcf_file_output.txt HDRNA_12_2_not_in_txt_file_output.txt HDRNA_12_2_common_records_output.txt

python3 ~/Scripts/merge_snippy_and_spandx_results.py results_HDRNA_16/snippy/merged_CP133688.vcf Outputs_CP133688/Phylogeny_and_annotation/All_SNPs_indels_annotated.txt HDRNA_16_not_in_vcf_file_output.txt HDRNA_16_not_in_txt_file_output.txt HDRNA_16_common_records_output.txt

python3 ~/Scripts/merge_snippy_and_spandx_results.py results_HDRNA_17/snippy/merged_CP133693.vcf Outputs_CP133693/Phylogeny_and_annotation/All_SNPs_indels_annotated.txt HDRNA_17_not_in_vcf_file_output.txt HDRNA_17_not_in_txt_file_output.txt HDRNA_17_common_records_output.txt

python3 ~/Scripts/merge_snippy_and_spandx_results.py results_HDRNA_19/snippy/merged_CP133696.vcf Outputs_CP133696/Phylogeny_and_annotation/All_SNPs_indels_annotated.txt HDRNA_19_not_in_vcf_file_output.txt HDRNA_19_not_in_txt_file_output.txt HDRNA_19_common_records_output.txt

python3 ~/Scripts/merge_snippy_and_spandx_results.py results_HDRNA_20/snippy/merged_CP133700.vcf Outputs_CP133700/Phylogeny_and_annotation/All_SNPs_indels_annotated.txt HDRNA_20_not_in_vcf_file_output.txt HDRNA_20_not_in_txt_file_output.txt HDRNA_20_common_records_output.txt

mv HDRNA_01_common_records_output.txt _HDRNA_01.txt
mv HDRNA_03_common_records_output.txt  HDRNA_03.txt
mv HDRNA_06_common_records_output.txt  HDRNA_06.txt
mv HDRNA_07_common_records_output.txt  HDRNA_07.txt
mv HDRNA_08_common_records_output.txt  HDRNA_08.txt
mv HDRNA_12_common_records_output.txt  HDRNA_12.txt
mv HDRNA_16_common_records_output.txt  HDRNA_16.txt
mv HDRNA_17_common_records_output.txt  HDRNA_17.txt
mv HDRNA_19_common_records_output.txt  HDRNA_19.txt
mv HDRNA_20_common_records_output.txt  HDRNA_20.txt

sed -i '1s/_trimmed_P//g' _HDRNA_01.txt
sed -i '1s/_trimmed_P//g' HDRNA_03.txt HDRNA_06.txt HDRNA_07.txt HDRNA_08.txt HDRNA_12.txt HDRNA_16.txt HDRNA_17.txt HDRNA_19.txt HDRNA_20.txt

# -- check if f3==f6 --
cut -f3 -d$'\t' HDRNA_17.txt > f3
cut -f6 -d$'\t' HDRNA_17.txt > f6
diff f3 f6

# -- check if f6==f7 in HDRNA_7.txt since they have the sample names --
cut HDRNA_07.txt -d$'\t' -f6  > f6
cut HDRNA_07.txt -d$'\t' -f7  > f7
diff d6 f7
--> delete the column HDRNA_07_K01-BB28 in variant_calling.xls.

~/Tools/csv2xls-0.4/csv_to_xls.py _HDRNA_01.txt HDRNA_03.txt HDRNA_06.txt HDRNA_07.txt HDRNA_08.txt HDRNA_12.txt HDRNA_16.txt HDRNA_17.txt HDRNA_19.txt HDRNA_20.txt -d$'\t' -o variant_calling.xls

processing commands for presence phage HH1, SPbeta-like phage, phage related island

#makeblastdb -in HDRNA_K01.fna -dbtype nucl
blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query MT880870.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880870_on_01.blastn
blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query MT880871.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880871_on_01.blastn
blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query MT880872.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880872_on_01.blastn

blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query MT880870.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880870_on_03.blastn
blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query MT880871.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880871_on_03.blastn
blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query MT880872.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880872_on_03.blastn

blastn -db HDRNA_06_K01_conservative_27645.current.gb_converted.fna -query MT880870.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880870_on_06.blastn
blastn -db HDRNA_06_K01_conservative_27645.current.gb_converted.fna -query MT880871.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880871_on_06.blastn
blastn -db HDRNA_06_K01_conservative_27645.current.gb_converted.fna -query MT880872.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880872_on_06.blastn

blastn -db HDRNA_07_K01_conservative_27169.current.gb_converted.fna -query MT880870.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880870_on_07.blastn
blastn -db HDRNA_07_K01_conservative_27169.current.gb_converted.fna -query MT880871.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880871_on_07.blastn
blastn -db HDRNA_07_K01_conservative_27169.current.gb_converted.fna -query MT880872.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880872_on_07.blastn

blastn -db HDRNA_08_K01_conservative_32455.current.gb_converted.fna -query MT880870.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880870_on_08.blastn
blastn -db HDRNA_08_K01_conservative_32455.current.gb_converted.fna -query MT880871.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880871_on_08.blastn
blastn -db HDRNA_08_K01_conservative_32455.current.gb_converted.fna -query MT880872.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880872_on_08.blastn

blastn -db HDRNA_12_K01_bold_37467.current.gb_converted.fna -query MT880870.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880870_on_12.blastn
blastn -db HDRNA_12_K01_bold_37467.current.gb_converted.fna -query MT880871.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880871_on_12.blastn
blastn -db HDRNA_12_K01_bold_37467.current.gb_converted.fna -query MT880872.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880872_on_12.blastn

blastn -db HDRNA_16_K01_conservative_37834.current.gb_converted.fna -query MT880870.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880870_on_16.blastn
blastn -db HDRNA_16_K01_conservative_37834.current.gb_converted.fna -query MT880871.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880871_on_16.blastn
blastn -db HDRNA_16_K01_conservative_37834.current.gb_converted.fna -query MT880872.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880872_on_16.blastn

blastn -db HDRNA_17_K01_conservative_37288.current.gb_converted.fna -query MT880870.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880870_on_17.blastn
blastn -db HDRNA_17_K01_conservative_37288.current.gb_converted.fna -query MT880871.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880871_on_17.blastn
blastn -db HDRNA_17_K01_conservative_37288.current.gb_converted.fna -query MT880872.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880872_on_17.blastn

blastn -db HDRNA_19_K01_bold_37377.current.gb_converted.fna -query MT880870.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880870_on_19.blastn
blastn -db HDRNA_19_K01_bold_37377.current.gb_converted.fna -query MT880871.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880871_on_19.blastn
blastn -db HDRNA_19_K01_bold_37377.current.gb_converted.fna -query MT880872.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880872_on_19.blastn

blastn -db HDRNA_20_K01_conservative_43457.current.gb_converted.fna -query MT880870.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880870_on_20.blastn
blastn -db HDRNA_20_K01_conservative_43457.current.gb_converted.fna -query MT880871.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880871_on_20.blastn
blastn -db HDRNA_20_K01_conservative_43457.current.gb_converted.fna -query MT880872.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > MT880872_on_20.blastn

ΦSepi-HH1(MT880870): 34053 bp in ST2, ST83 (34053 bt)
PI-Sepi-HH2(MT880871): not in ST2, but in ST290, ST297 and ST487 (36164 bt)
ΦSPbeta-like(MT880872): in ST2 and ST22 (147057 bt)
# In the new isolates, we have the MLST

-
ST130
ST224
ST23
ST35
ST487 -->
ST5
ST640
ST69
ST86
ST87
--> HDRNA_17_K01 (ST69) contains MT880871!
It agrees with the description in the Anna's paper, ST487 has PI-Sepi-HH2!
#shovill/HDRNA_11_K01/contigs.fa sepidermidis    487     arcC(1) aroE(1) gtr(1)  mutS(5) pyrR(2) tpiA(1) yqiL(1)

# -- 01 --
33342   34053 --> am Grenzen
33634   34132 --> am Grenzen
131829  147057 + 1       3810 --> am Grenzen
30857   34053 --> NA
33634   34132 --> NA
131829  147057 + 1       3807 --> NA
NA
NA
NA
07: NA
08: NA
12: NA
16: NA
17: MT880871_on_17.blastn YES
19: NA
20: NA

shovill/HDRNA_11_K01/contigs.fa sepidermidis    487     arcC(1) aroE(1) gtr(1)  mutS(5) pyrR(2) tpiA(1) yqiL(1)

a Biofilm
formation compared to S. epidermidis 1457: –, ⬍30%; ⫹, 30% to 59%; ⫹⫹, 60% to 84%; ⫹⫹⫹, ⱖ85%.
b ST,
sequence type determined by MLST.
c GC, genetic cluster based on Bayesian analysis of population structure (BAPS) of MLST.
d ND, not determined. In columns 2 to 5, ⫹ indicates presence and ⫺ indicates absence.

ND, not determined or defined. '+' indicates presence and '-' indicates absence.

SPbeta-like staphylococcus phage (NC_029119.1, 86% sequence identity)[47] and will here be referred to as
ΦSPbeta-like, (GenBank accession number MT880872) the remaining two regions did not correspond to any previously described Staphylococcus phage[48]
ΦSepi-HH1 (MT880870): 34053 bp
phage-related island PI-Sepi-HH2 (MT880871).
MT880870.1: Staphylococcus phage PhiSepi-HH1, complete genome
MT880871.1: Staphylococcus phage PI-Sepi-HH2, complete genome: 36164 bp
MT880872.1: Staphylococcus phage PhiSepi-HH3, complete genome

processing commands for the other genes from Gene_List.pptx

        # AND "Staphylococcus epidermidis"[porgn:__txid1282]

mv HDRNA_16_K01_conservative_37834.current.gb ../gbk
mv ./HDRNA_17_K01_conservative_37288.current.gb ../gbk
mv ./HDRNA_03_K01_bold_bandage_26831.current.gb ../gbk
mv ./HDRNA_06_K01_conservative_27645.current.gb ../gbk
mv ./HDRNA_12_K01_bold_37467.current.gb ../gbk
mv ./HDRNA_01_K01_conservative_23197.current.gb ../gbk
mv ./HDRNA_07_K01_conservative_27169.current.gb ../gbk
mv ./HDRNA_08_K01_conservative_32455.current.gb ../gbk
mv ./HDRNA_19_K01_bold_37377.current.gb ../gbk
mv ./HDRNA_20_K01_conservative_43457.current.gb ../gbk

(agrC) AND "Staphylococcus"
https://www.ncbi.nlm.nih.gov/nuccore/OQ828637.1

gltA
agrC

gene            2184799..2186091
                    /gene="agrC"
                    /locus_tag="SAR2125"
    CDS             2184799..2186091
                    /gene="agrC"
                    /locus_tag="SAR2125"
                    /note="Signal dectecting component of the agr autoinducer
                    peptide-quorum sensing system. Two-component regulatory
                    system family, sensor kinase protein. Similar to
                    Staphylococcus aureus accessory gene regulator C AgrC
                    TR:Q53644 (EMBL:X52543) (423 aa) fasta scores: E():
                    1.7e-101, 75.23% id in 424 aa, and to Staphylococcus
                    epidermidis histidine kinase AgrC TR:O68159
                    (EMBL:AF012132) (429 aa) fasta scores: E(): 6.3e-76,
                    54.93% id in 426 aa"
                    /codon_start=1
                    /transl_table=11
                    /product="autoinducer sensor protein"
                    /protein_id="CAG41106.1"
                    /db_xref="EnsemblGenomes-Gn:SAR2125"
                    /db_xref="EnsemblGenomes-Tr:CAG41106"
                    /translation="MEALNDYNYVLFVIVQVSLMFFISAFISGIRYKKSDYIYIIGIV
                    LSSVYFFDKIRSISLVVITIFIIIFLYFKIRLYSVFLVMVTQIILYCANFVYIIIFSY
                    IITISHSVFIVLPIFLVVYVSISYALAYILNRILKRINGTYLSLNKKFLTVITIVIVI
                    TFSLLFAYSQIDASDASTIKQYSLLFLGIIILLSILIFIYSQFTLKEMKYKRNQEEIE
                    TYYEYTLKIEAINNEMRKFRHDYVNILTTLSEYIREDDMTGLRDYFNKNIVPMKDNLQ
                    MNALKLNGIENLKVREIKGLLTAKILRAQEMNIPISIEIPDEVTRINLNMIDLSRSIG
                    IILDNAIEASSEIDDPIIRVAFIESENSVTFIVMNKCADDIPRIHELFQESFSTKGEG
                    RGLGLSTLKEIADNADNVLLDTIIENGFFIQKVEIINN"
MEALNDYNYVLFVIVQVSLMFFISAFISGIRYKKSDYIYIIGIVLSSVYFFDKIRSISLVVITIFIIIFLYFKIRLYSVFLVMVTQIILYCANFVYIIIFSYIITISHSVFIVLPIFLVVYVSISYALAYILNRILKRINGTYLSLNKKFLTVITIVIVITFSLLFAYSQIDASDASTIKQYSLLFLGIIILLSILIFIYSQFTLKEMKYKRNQEEIETYYEYTLKIEAINNEMRKFRHDYVNILTTLSEYIREDDMTGLRDYFNKNIVPMKDNLQMNALKLNGIENLKVREIKGLLTAKILRAQEMNIPISIEIPDEVTRINLNMIDLSRSIGIILDNAIEASSEIDDPIIRVAFIESENSVTFIVMNKCADDIPRIHELFQESFSTKGEGRGLGLSTLKEIADNADNVLLDTIIENGFFIQKVEIINN
Staphylococcus_aureus_MRSA252

samtools faidx Staphylococcus_aureus_MRSA252.fasta "gi|49240382|emb|BX571856.1|":2184799-2186091

(yycG) AND "Staphylococcus"
    gene            25617..27443
                    /gene="yycG"
                    /locus_tag="SAR0019"
                    /gene_synonym="vicK"
    CDS             25617..27443
                    /gene="yycG"
                    /locus_tag="SAR0019"
                    /gene_synonym="vicK"
                    /note="Two-component regulatory system family, sensor
                    kinase protein. Previously sequenced as Staphylococcus
                    aureus two-component sensor histidine kinase YycG
                    TR:Q9XCM6 (EMBL:AF136709) (608 aa) fasta scores: E():
                    5.5e-214, 99.836% id in 608 aa. Similar to Bacillus
                    subtilis probable two-component sensor histidine kinase
                    YycG TR:Q45614 (EMBL:D78193) (611 aa) fasta scores: E():
                    2e-98, 46.217% id in 608 aa"
                    /codon_start=1
                    /transl_table=11
                    /product="sensor kinase protein"
                    /protein_id="CAG39047.1"
                    /db_xref="EnsemblGenomes-Gn:SAR0019"
                    /db_xref="EnsemblGenomes-Tr:CAG39047"
                    /db_xref="GOA:Q6GKS6"
                    /db_xref="InterPro:IPR000014"
                    /db_xref="InterPro:IPR000700"
                    /db_xref="InterPro:IPR003594"
                    /db_xref="InterPro:IPR003660"
                    /db_xref="InterPro:IPR003661"
                    /db_xref="InterPro:IPR004358"
                    /db_xref="InterPro:IPR005467"
                    /db_xref="InterPro:IPR029151"
                    /db_xref="UniProtKB/Swiss-Prot:Q6GKS6"
                    /translation="MKWLKQLQSLHTKLVIVYVLLIIIGMQIIGLYFTNNLEKELLDN
                    FKKNITQYAKQLEISIEKVYDEKGSVNAQKDIQNLLSEYANRQEIGEIRFIDKDQIII
                    ATTKQSNRSLINQKANDSSVQKALSLGQSNDHLILKDYGGGKDRVWVYNIPVKVDKKV
                    IGNIYIESKINDVYNQLNNINQIFIVGTAISLLITVILGFFIARTITKPITDMRNQTV
                    EMSRGNYTQRVKIYGNDEIGELALAFNNLSKRVQEAQANTESEKRRLDSVITHMSDGI
                    IATDRRGRIRIVNDMALKMLGMAKEDIIGYYMLSVLSLEDEFKLEEIQENNDSFLLDL
                    NEEEGLIARVNFSTIVQETGFVTGYIAVLHDVTEQQQVERERREFVANVSHELRTPLT
                    SMNSYIEALEEGAWKDEELAPQFLSVTREETERMIRLVNDLLQLSKMDNESDQINKEI
                    IDFNMFINKIINRHEMSTKDTTFIRDIPKKTIFTEFDPDKMTQVFDNVITNAMKYSRG
                    DKRVEFHVKQNPLYNRMTIRIKDNGIGIPINKVDKIFDRFYRVDKARTRKMGGTGLGL
                    AISKEIVEAHNGRIWANSVEGQGTSIFITLPCEVIEDGDWDE"

samtools faidx Staphylococcus_aureus_MRSA252.fasta "gi|49240382|emb|BX571856.1|":25617-27443

psmß1: https://www.ncbi.nlm.nih.gov/nuccore/380448412
hlb: https://www.ncbi.nlm.nih.gov/nuccore/KC242859.1

atlE

tagB

    gene            694999..696102
                    /gene="tagB"
                    /locus_tag="SAR0649"
    CDS             694999..696102
                    /gene="tagB"
                    /locus_tag="SAR0649"
                    /note="Similar to Bacillus subtilis teichoic acid
                    biosynthesis protein B TagB SW:TAGB_BACSU (P27621) (381
                    aa) fasta scores: E(): 4.3e-26, 31.302% id in 361 aa, and
                    to Lactococcus lactis teichoic acid biosynthesis protein B
                    TagB TR:Q9CH14 (EMBL:AE006327) (371 aa) fasta scores: E():
                    2.8e-06, 24.834% id in 302 aa. Lack of similarity at the
                    N-terminus in comparison to other orthologues"
                    /codon_start=1
                    /transl_table=11
                    /product="teichoic acid biosynthesis protein"
                    /protein_id="CAG39666.1"
                    /db_xref="EnsemblGenomes-Gn:SAR0649"
                    /db_xref="EnsemblGenomes-Tr:CAG39666"
                    /translation="MNVLIKKFYHLVVRILSKMITPQVIDKPHIVFMMTFPEDIKPII
                    KALNNSLYQKTVLTTPKQAPYLSELSDDVNVIEMTNRTLVKQIKALKSAQMIIIDNYY
                    LLLGGYNKTSNQHIVQTWHASGALKNFGLTDHQVDVSDKAMVQQYRKVYQATDFYLVG
                    CEQMSQCFKQSLGATEEQMLYFGLPRINKYYTADRATVKAELKDKYGITNKLALYVPT
                    YREDKADNRAIDKAYFEKCLPGYTLINKLHPSIEHSDIDDVSSIDTSILMLMSDIIIS
                    DYSSLPIEASLLDIPTIFYVYDEGTYDKVRGLNQFYKAIPDSYKVYTEEDLIMTIQEK
                    EHLLSPLFKDWHKYNTDKSLHQLTEYIDKMVTK"

samtools faidx Staphylococcus_aureus_MRSA252.fasta "gi|49240382|emb|BX571856.1|":694999-696102 > tagB.fasta

    gene            2137..2586
                    /gene="capC"
    CDS             2137..2586
                    /gene="capC"
                    /codon_start=1
                    /transl_table=11
                    /product="CapC"
                    /protein_id="BAB13485.1"
                    /translation="MFGSDLYIALILGVLLSLIFAEKTGIVPAGLVVPGYLGLVFNQP
                    VFILLVLLVSLLTYVIVKYGLSKFMILYGRRKFAAMLITGIVLKIAFDFLYPIVPFEI
                    AEFRGIGIIVPGLIANTIQKQGLTITFGSTLLLSGATFAIMFVYYLI"

samtools faidx capBCA_ywtC.fasta "gi|10119860|dbj|AB039950.1|":2137-2586 > capC.fasta

    gene            2783..3256
                    /gene="sepA"
    CDS             2783..3256
                    /gene="sepA"
                    /function="multidrug resistance"
                    /experiment="experimental evidence, no additional details
                    recorded"
                    /note="similar to BAB43260.1 of S. aureus N315"
                    /codon_start=1
                    /transl_table=11
                    /product="SepA"
                    /protein_id="BAB83937.1"
                    /translation="MIVNYLKHKFYNLLTTMIVLFIFVLSGAIFLTFLGFGLYGLSRI
                    LIYFRLGDFTYNRSMYDNLLYYGSYIIFGYFIIFAVEHLMDYFRKMLPENAYFRGATF
                    HLISYTVATTLFYFIIHLNYVYINIDFWVIMVIIGFLYVCKLQFYPESKNLNNRK"

samtools faidx ORF123_sepA_ORF5.fasta "gi|18250967|dbj|AB078343.1|":2783-3256 > sepA.fasta

    gene            1427047..1429569
                    /gene="mprF"
                    /locus_tag="SAR1372"
                    /gene_synonym="fmtC"
    CDS             1427047..1429569
                    /gene="mprF"
                    /locus_tag="SAR1372"
                    /gene_synonym="fmtC"
                    /note="Similar to Staphylococcus aureus putative membrane
                    protein MprF TR:AAK58115 (EMBL:AF145699) (840 aa) fasta
                    scores: E(): 0, 96.190% id in 840 aa. Similar to
                    Staphylococcus xylosus putative membrane protein MprF
                    TR:AAK58113 (EMBL:AF145698) (841 aa) fasta scores: E():
                    5.4e-208, 62.530% id in 838 aa. Mutations in the CDS have
                    reduced resistance to human defensins and evasion of
                    neutrophil killing"
                    /codon_start=1
                    /transl_table=11
                    /product="putative membrane protein"
                    /protein_id="CAG40370.1"
                    /db_xref="EnsemblGenomes-Gn:SAR1372"
                    /db_xref="EnsemblGenomes-Tr:CAG40370"
                    /db_xref="GOA:Q6GH45"
                    /db_xref="InterPro:IPR016181"
                    /db_xref="InterPro:IPR022791"
                    /db_xref="InterPro:IPR024320"
                    /db_xref="UniProtKB/Swiss-Prot:Q6GH45"
                    /translation="MNQEVKNKIFSILKITFATALFIFVVITLYRELSGINFKDTLVE
                    FSKINRMSLVLLFIGGGASLVILSMYDVILSRALKMDISLGKVLRVSYIINALNAIVG
                    FGGFIGAGVRAMVYKNYTHDKKKLVHFISLILISMLTGLSLLSLLIVFHVFDASLILN
                    KITWVRWVLYAVSLFLPLFIIYSMVRPPDKNNRYVGLYCTLVSCVEWLAAAVVLYFCG
                    VIVDVHVSFMSFIAIFIIAALSGLVSFIPGGFGAFDLVVLLGFKTLGVPEEKVLLMLL
                    LYRFAYYFVPVIIALILSSFEFGTSAKKYIEGSKYFIPAKDVTSFLMSYQKDIIAKIP
                    SLSLAILVFFTSMIFFVNNLTIVYDALYDGNHLTYYLLLAIHTSACLLLLLNVVGIYK
                    QSRRAIIYAMISIILIIVATLFTYASYILITWLVIIFALLIVAFRRARRLKRPIRMRN
                    LVAMLLFSIFILYINHIFIAGTFYALDVYTIEMHTSVLKYYFWITILIIAIIVGAIAW
                    LFDYQFSKVRISSNIEECEEIIDQYGGNYLSHLIYSGDKQFFTNEDKNAFLMYRYKAS
                    SLVVLGDPIGDENAFDELLEAFYNYAEYLGYDVIFYQVTDQHMPLYHNFGNQFFKLGE
                    EAIIDLTQFSTSGKKRRGFRATLNKFDELNISFEIIEPPFSTEFINELQHVSDLWLDN
                    RQEMHFSVGQFNETYLSKAPIGVMRNENNEVIAFCSLMPTYFNDAISVDLIRWLPELD
                    LPLMDGLYLHMLLWSKEQGYTKFNMGMATLSNVGQLHYSYLRERLAGRVFEHFNGLYR
                    FQGLRRYKSKYNPNWEPRFLVYRKDNSLWESLSKVMRVIRHK"

samtools faidx Staphylococcus_aureus_MRSA252.fasta "gi|49240382|emb|BX571856.1|":1427047-1429569 > fmtC.fasta

    gene            1825..2523
                    /gene="sceD"
    CDS             1825..2523
                    /gene="sceD"
                    /codon_start=1
                    /transl_table=11
                    /product="SceD precursor"
                    /protein_id="AAB94657.1"
                    /translation="MKKLLVASSASAALFAVGVGANAHAAEDNNVNQDQLAQTALNNT
                    QQLNDAPVQEGAYNIAFDNSGYNFNFNSDGTNWSWSYNADSSAQQAPAQSTTQEQAPA
                    AQQAPAQSTTQEQAPAAQQAPAQEQTQQPAQQPAQQQTQQPAQQSADSGSNVQVNDHL
                    KAIAQRESGGDIHAINSSSGAAGKYQFLQTTWDSVAPAEYQGKPASEAPEAVQDAAAQ
                    KLYDTAGPSQWVTA"
    sig_peptide     1825..1899
                    /gene="sceD"
    mat_peptide     1900..2520
                    /gene="sceD"
                    /product="SceD"

samtools faidx sceDAE.fasta "gi|6707000|gb|AF109218.1|":1825-2523 > sceD.fasta

    gene            759097..759894
                    /locus_tag="SERP_RS03845"
                    /old_locus_tag="SE0760"
                    /old_locus_tag="SERP0760"
    CDS             759097..759894
                    /locus_tag="SERP_RS03845"
                    /old_locus_tag="SE0760"
                    /old_locus_tag="SERP0760"
                    /inference="COORDINATES: similar to AA
                    sequence:RefSeq:WP_001830106.1"
                    /note="Derived by automated computational analysis using
                    gene prediction method: Protein Homology."
                    /codon_start=1
                    /transl_table=11
                    /product="VOC family protein"
                    /protein_id="WP_002446237.1"
                    /db_xref="GI:488376852"
                    /translation="MFHNKNANFVNGVTINIKDKETIKTFYENVLGFNLINESKTAVQ
                    FEVGDSNQFITFIEIQNGREPLMSEAGLFHIGILLPTLTTLADLLVHLSDFEVPVNGG
                    QQSVATCLFIEDPEGNAIKFYVDRETESWIDEKEGRIRMDIAPINVPRLLQNVSHTQW
                    QGIPDETKLGSLHIKSIRISDVKSYYLHYFGLEESAYMDDYSLFLSSNDYYNHLAVNQ
                    WLSATKRVDNEHTYGLAMIDFHYPKTTHKNLKGPDGIYFRFNRIKEV"

samtools faidx Staphylococcus_epidermidis_RP62A.fasta "gi|57865352|ref|NC_002976.3|":759097-759894 > SE0760.fasta

    gene            3091..8793
                    /gene="esp"
                    /allele="1"
    CDS             3091..8793
                    /gene="esp"
                    /allele="1"
                    /note="enterococcal surface protein; ORF3"
                    /codon_start=1
                    /transl_table=11
                    /product="Esp"
                    /protein_id="AAQ84025.1"
                    /translation="MVSKNNKRVFLEKTKKRVLKYSIKKLSVGVASVLVGVGLVLGTT
                    ELVQAQDEISPSTPLETAISSVQVGDKVASGNTFQENPGYTKNYNFSDLQFSPQELTG
                    DTLKGNTIGFEVYGKHNIAASTKNWEIRLQLDERLAKYVEKIQVDPKKGIGSSRRTFV
                    RINDSLGRPTNIWKVNYIRASDGLFAGAETTDTQTAPNGVITFEKSLDEIFKEIGIDN
                    LKTDRLMYRIYLVSHQDDDKIVPGIDSTGYFLTDSDDFYNSLDVSENNPDQFKHGSVS
                    AKYEEPNTQTKDGSGSTGANGAIILDHKLTKNYNFSYSASAKGTPWYANYKIDERLVP
                    YVAGIQMHMVQADKVTYDVSFESGKKVADLAIERRKDHENYGVGSITDNDLTKLIDFA
                    NASPRPVVIRYVLQLTKPLDEILEDMKATAQVEENKPFGEDFIFDSWLSDTNKKLIQN
                    TYGTGYYYLQDIDGDGNPDDKEESGDTNPYIGKPELEEVYDVDTTVKGKVFIHELAGT
                    GHKAQLVDKEGTVLAEKTIAPNEKDGAPISDTVEFEFTGVDSSKLIAKDELKIQIVSP
                    GFDKPEEGSTVIKESPKAVDKQTVVVGFKPDAKESIRNNKNLPEDAEYSWKTEPDTSN
                    VTDSTKGIVTVKIGNRTFDVDVEFAVKASQAMENDATYVPITTTPETTVQSGKPTFDK
                    PDVPLAKDAFSILDVYNKDFGNASVDANTGIVTFTPAKGVGESEPITGTIPIKIVYQD
                    GSVGTTDLAVTVSKDIYENPGENIPAGYHKVTFTAGEGTSIESGTTVFAVKDGVSLPE
                    DKLPVLKAKDGYTDAKWPEEATQPITADDTEFVSSATKLDDIIENPGENIPAGYHKVT
                    FTAGEGTSIESGTTVFAVKDGVSLPEDKLPVLKAKDGYTDAKWPEEATQPITADDTEF
                    VSSATKLDDIIENPGENIPAGYHKVTFTAGEGTSIESGTTVFAVKDGVSLPEDKLPVL
                    KAKDGYTDAKWPEEATQPITADDTEFVSSATKLDDIIENPGENIPAGYHKVIFTAGEG
                    TSIESGTTVFAVKDGVSLPEDKLPVLKAKDGYTDAKWPEEATQPITADDTEFVSSATK
                    LDDIIENPGENIPAGYHKVIFTAGEGTSIESGTTVFAVKDGVSLPEDKLPVLKAKDGY
                    TDAKWPEEATQPITADDTEFVSSATKLDDKSDADKYNPEGQKVTTELNKEPDASEGIK
                    NKEDLPKDTKYTWKEKVDVSAAGNKKGTVVVTYSDGSSDEVEVDVTVTDNRSDADKYE
                    PTVEGEKVEVGGTVDLTDNVTNLPTLPEGTTVTDVTPDGTIDTNTPGNYEGVIEVTYP
                    DGTKDTVKVPVEVTDNRSDADKYTPMVEGEKVEVGGTVDLTDNVTNLPTLPEGTTVTD
                    VTPDGTIDTNTPGNYEGVIEVTYPDGTKDTVKVPVEVTDNRSDADKYTPMVEGEKVEV
                    GGTVDLTDNVTNLPTLPEGTTVTDVTPGGTIDTNTPGNYEGVIEVTYPDGTKDTVKVP
                    VEVTDNRSDADKYEPTVEGEKVEVGGTVDLTDNVTNLPTLPEGTTVTDVTPGGTIDTN
                    TPGNYEGVIEVTYPDGTKDTVKVPVEVTDNRSDADKYTPMVEGEKVEVGGTVDLTDNV
                    TNLPTLPEGTTVTDVTPDGTIDTNTPGNYEGVIEVTYPDGTKDTVKVPVEVTDNRSDA
                    DKYNPEGQKVTTDLNKEPDASEGIKNKEDLPKGTTYTWKEKVDVSTAGNKKGIVVVTY
                    PDGSKEEVEVTISVEDKKAPNKPQVDPITEGDQIVTGKTEPNAEVTVTLPNGSQYHGT
                    ADKNGNFTVKVPKLEAGTKVIVTATDESGNTSEPTNVVVSSNEKDSEKAVSKDNKTDN
                    QGSKQNTNRGKSSPQKQSSKAYPKTGEIDSNIFTISGGLILLGTLGLLGYKNRKKENE
                    "
samtools faidx Enterococcus_faecium_isolate_E300_pathogenicity_island.fasta "gi|34980227|gb|AY322150.1|":3091-8793 > esp.fasta

mv ~/Downloads/sequence\ \(1\).fasta ecpA.fasta

mv ~/Downloads/sequence\ \(1\).fasta sdrG.fasta

HH1-HP1
HH1-HP2
HH1-TreR

sdrG.fasta

gltA.fasta
agrC.fasta
yycG.fasta
atlE.fasta
psm-beta1.fasta
    gene            366..500
                    /gene="psm beta1"
    CDS             366..500
                    /gene="psm beta1"
                    /note="hemolytic protein; pfam05480"
                    /codon_start=1
                    /transl_table=11
                    /product="Psm beta1"
                    /protein_id="AFD54319.1"
                    /translation="MEGLFNAIKDTVTAAINNDGAKLGTSIVSIVENGVGLLGKLFGF
                    "
    gene            557..691
                    /gene="psm beta2"
    CDS             557..691
                    /gene="psm beta2"
                    /note="hemolytic protein; pfam05480"
                    /codon_start=1
                    /transl_table=11
                    /product="Psm beta2"
                    /protein_id="AFD54320.1"
                    /translation="MTGLAEAIANTVQAAQQHDSVKLGTSIVDIVANGVGLLGKLFGF
                    "
samtools faidx psm-beta.fasta "gi|380448412|gb|JQ066320.1|":366-500 > psm-beta1.fasta

hlb.fasta
samtools faidx hlb_.fasta "gi|441494790|gb|KC242859.1|":16-1008 > hlb.fasta
tagB.fasta
capC.fasta
sepA.fasta
fmtC.fasta
sceD.fasta
SE0760.fasta
esp.fasta
samtools faidx ecpA_.fasta "gi|354620065|gb|JN051494.1|":248-835 > ecpA.fasta
ecpA.fasta

#TODO TOMORROW: add the gene names in the fasta and merge them in a file using the following command line, make once blastn, oder make gene by gene. It is maybe clearer! find the last 3 gene sequences. At first, send the results except for gene arrangements!
#cat sdrG.fasta gltA.fasta agrC.fasta yycG.fasta atlE.fasta psm-beta1.fasta hlb.fasta tagB.fasta capC.fasta sepA.fasta fmtC.fasta sceD.fasta SE0760.fasta esp.fasta ecpA.fasta > 15genes.fasta

#makeblastdb -in HDRNA_K01.fna -dbtype nucl

# -- sdrG (repeated) --
blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query sdrG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sdrG_on_01.blastn
blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query sdrG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sdrG_on_03.blastn
blastn -db HDRNA_06_K01_conservative_27645.current.gb_converted.fna -query sdrG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sdrG_on_06.blastn
blastn -db HDRNA_07_K01_conservative_27169.current.gb_converted.fna -query sdrG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sdrG_on_07.blastn
blastn -db HDRNA_08_K01_conservative_32455.current.gb_converted.fna -query sdrG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sdrG_on_08.blastn
blastn -db HDRNA_12_K01_bold_37467.current.gb_converted.fna -query sdrG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sdrG_on_12.blastn
blastn -db HDRNA_16_K01_conservative_37834.current.gb_converted.fna -query sdrG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sdrG_on_16.blastn
blastn -db HDRNA_17_K01_conservative_37288.current.gb_converted.fna -query sdrG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sdrG_on_17.blastn
blastn -db HDRNA_19_K01_bold_37377.current.gb_converted.fna -query sdrG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sdrG_on_19.blastn
blastn -db HDRNA_20_K01_conservative_43457.current.gb_converted.fna -query sdrG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sdrG_on_20.blastn

# -- gltA --
blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query gltA.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > gltA_on_01.blastn
blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query gltA.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > gltA_on_03.blastn
blastn -db HDRNA_06_K01_conservative_27645.current.gb_converted.fna -query gltA.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > gltA_on_06.blastn
blastn -db HDRNA_07_K01_conservative_27169.current.gb_converted.fna -query gltA.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > gltA_on_07.blastn
blastn -db HDRNA_08_K01_conservative_32455.current.gb_converted.fna -query gltA.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > gltA_on_08.blastn
blastn -db HDRNA_12_K01_bold_37467.current.gb_converted.fna -query gltA.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > gltA_on_12.blastn
blastn -db HDRNA_16_K01_conservative_37834.current.gb_converted.fna -query gltA.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > gltA_on_16.blastn
blastn -db HDRNA_17_K01_conservative_37288.current.gb_converted.fna -query gltA.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > gltA_on_17.blastn
blastn -db HDRNA_19_K01_bold_37377.current.gb_converted.fna -query gltA.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > gltA_on_19.blastn
blastn -db HDRNA_20_K01_conservative_43457.current.gb_converted.fna -query gltA.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > gltA_on_20.blastn
#-->None

# -- agrC --
blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query agrC.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > agrC_on_01.blastn
blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query agrC.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > agrC_on_03.blastn
blastn -db HDRNA_06_K01_conservative_27645.current.gb_converted.fna -query agrC.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > agrC_on_06.blastn
blastn -db HDRNA_07_K01_conservative_27169.current.gb_converted.fna -query agrC.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > agrC_on_07.blastn
blastn -db HDRNA_08_K01_conservative_32455.current.gb_converted.fna -query agrC.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > agrC_on_08.blastn
blastn -db HDRNA_12_K01_bold_37467.current.gb_converted.fna -query agrC.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > agrC_on_12.blastn
blastn -db HDRNA_16_K01_conservative_37834.current.gb_converted.fna -query agrC.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > agrC_on_16.blastn
blastn -db HDRNA_17_K01_conservative_37288.current.gb_converted.fna -query agrC.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > agrC_on_17.blastn
blastn -db HDRNA_19_K01_bold_37377.current.gb_converted.fna -query agrC.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > agrC_on_19.blastn
blastn -db HDRNA_20_K01_conservative_43457.current.gb_converted.fna -query agrC.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > agrC_on_20.blastn
#-->None

# -- yycG (1827 nt) --
blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query yycG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > yycG_on_01.blastn
blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query yycG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > yycG_on_03.blastn
blastn -db HDRNA_06_K01_conservative_27645.current.gb_converted.fna -query yycG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > yycG_on_06.blastn
blastn -db HDRNA_07_K01_conservative_27169.current.gb_converted.fna -query yycG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > yycG_on_07.blastn
blastn -db HDRNA_08_K01_conservative_32455.current.gb_converted.fna -query yycG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > yycG_on_08.blastn
blastn -db HDRNA_12_K01_bold_37467.current.gb_converted.fna -query yycG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > yycG_on_12.blastn
blastn -db HDRNA_16_K01_conservative_37834.current.gb_converted.fna -query yycG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > yycG_on_16.blastn
blastn -db HDRNA_17_K01_conservative_37288.current.gb_converted.fna -query yycG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > yycG_on_17.blastn
blastn -db HDRNA_19_K01_bold_37377.current.gb_converted.fna -query yycG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > yycG_on_19.blastn
blastn -db HDRNA_20_K01_conservative_43457.current.gb_converted.fna -query yycG.fasta -evalue 1e-50 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > yycG_on_20.blastn

#-- atlE.fasta --
blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query atlE.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > atlE_on_01.blastn
blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query atlE.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > atlE_on_03.blastn
blastn -db HDRNA_06_K01_conservative_27645.current.gb_converted.fna -query atlE.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > atlE_on_06.blastn
blastn -db HDRNA_07_K01_conservative_27169.current.gb_converted.fna -query atlE.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > atlE_on_07.blastn
blastn -db HDRNA_08_K01_conservative_32455.current.gb_converted.fna -query atlE.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > atlE_on_08.blastn
blastn -db HDRNA_12_K01_bold_37467.current.gb_converted.fna -query atlE.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > atlE_on_12.blastn
blastn -db HDRNA_16_K01_conservative_37834.current.gb_converted.fna -query atlE.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > atlE_on_16.blastn
blastn -db HDRNA_17_K01_conservative_37288.current.gb_converted.fna -query atlE.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > atlE_on_17.blastn
blastn -db HDRNA_19_K01_bold_37377.current.gb_converted.fna -query atlE.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > atlE_on_19.blastn
blastn -db HDRNA_20_K01_conservative_43457.current.gb_converted.fna -query atlE.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > atlE_on_20.blastn

#-- psm-beta1.fasta --
blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query psm-beta1.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > psm-beta1_on_01.blastn
blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query psm-beta1.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > psm-beta1_on_03.blastn
blastn -db HDRNA_06_K01_conservative_27645.current.gb_converted.fna -query psm-beta1.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > psm-beta1_on_06.blastn
blastn -db HDRNA_07_K01_conservative_27169.current.gb_converted.fna -query psm-beta1.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > psm-beta1_on_07.blastn
blastn -db HDRNA_08_K01_conservative_32455.current.gb_converted.fna -query psm-beta1.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > psm-beta1_on_08.blastn
blastn -db HDRNA_12_K01_bold_37467.current.gb_converted.fna -query psm-beta1.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > psm-beta1_on_12.blastn
blastn -db HDRNA_16_K01_conservative_37834.current.gb_converted.fna -query psm-beta1.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > psm-beta1_on_16.blastn
blastn -db HDRNA_17_K01_conservative_37288.current.gb_converted.fna -query psm-beta1.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > psm-beta1_on_17.blastn
blastn -db HDRNA_19_K01_bold_37377.current.gb_converted.fna -query psm-beta1.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > psm-beta1_on_19.blastn
blastn -db HDRNA_20_K01_conservative_43457.current.gb_converted.fna -query psm-beta1.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > psm-beta1_on_20.blastn

#-- hlb.fasta --
blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query hlb.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > hlb_on_01.blastn
blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query hlb.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > hlb_on_03.blastn
blastn -db HDRNA_06_K01_conservative_27645.current.gb_converted.fna -query hlb.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > hlb_on_06.blastn
blastn -db HDRNA_07_K01_conservative_27169.current.gb_converted.fna -query hlb.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > hlb_on_07.blastn
blastn -db HDRNA_08_K01_conservative_32455.current.gb_converted.fna -query hlb.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > hlb_on_08.blastn
blastn -db HDRNA_12_K01_bold_37467.current.gb_converted.fna -query hlb.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > hlb_on_12.blastn
blastn -db HDRNA_16_K01_conservative_37834.current.gb_converted.fna -query hlb.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > hlb_on_16.blastn
blastn -db HDRNA_17_K01_conservative_37288.current.gb_converted.fna -query hlb.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > hlb_on_17.blastn
blastn -db HDRNA_19_K01_bold_37377.current.gb_converted.fna -query hlb.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > hlb_on_19.blastn
blastn -db HDRNA_20_K01_conservative_43457.current.gb_converted.fna -query hlb.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > hlb_on_20.blastn

#-- tagB.fasta --
blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query tagB.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > tagB_on_01.blastn
blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query tagB.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > tagB_on_03.blastn
blastn -db HDRNA_06_K01_conservative_27645.current.gb_converted.fna -query tagB.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > tagB_on_06.blastn
blastn -db HDRNA_07_K01_conservative_27169.current.gb_converted.fna -query tagB.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > tagB_on_07.blastn
blastn -db HDRNA_08_K01_conservative_32455.current.gb_converted.fna -query tagB.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > tagB_on_08.blastn
blastn -db HDRNA_12_K01_bold_37467.current.gb_converted.fna -query tagB.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > tagB_on_12.blastn
blastn -db HDRNA_16_K01_conservative_37834.current.gb_converted.fna -query tagB.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > tagB_on_16.blastn
blastn -db HDRNA_17_K01_conservative_37288.current.gb_converted.fna -query tagB.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > tagB_on_17.blastn
blastn -db HDRNA_19_K01_bold_37377.current.gb_converted.fna -query tagB.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > tagB_on_19.blastn
blastn -db HDRNA_20_K01_conservative_43457.current.gb_converted.fna -query tagB.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > tagB_on_20.blastn

#-- capC.fasta --
blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query capC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > capC_on_01.blastn
blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query capC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > capC_on_03.blastn
blastn -db HDRNA_06_K01_conservative_27645.current.gb_converted.fna -query capC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > capC_on_06.blastn
blastn -db HDRNA_07_K01_conservative_27169.current.gb_converted.fna -query capC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > capC_on_07.blastn
blastn -db HDRNA_08_K01_conservative_32455.current.gb_converted.fna -query capC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > capC_on_08.blastn
blastn -db HDRNA_12_K01_bold_37467.current.gb_converted.fna -query capC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > capC_on_12.blastn
blastn -db HDRNA_16_K01_conservative_37834.current.gb_converted.fna -query capC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > capC_on_16.blastn
blastn -db HDRNA_17_K01_conservative_37288.current.gb_converted.fna -query capC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > capC_on_17.blastn
blastn -db HDRNA_19_K01_bold_37377.current.gb_converted.fna -query capC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > capC_on_19.blastn
blastn -db HDRNA_20_K01_conservative_43457.current.gb_converted.fna -query capC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > capC_on_20.blastn

#-- sepA.fasta --
blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query sepA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sepA_on_01.blastn
blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query sepA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sepA_on_03.blastn
blastn -db HDRNA_06_K01_conservative_27645.current.gb_converted.fna -query sepA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sepA_on_06.blastn
blastn -db HDRNA_07_K01_conservative_27169.current.gb_converted.fna -query sepA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sepA_on_07.blastn
blastn -db HDRNA_08_K01_conservative_32455.current.gb_converted.fna -query sepA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sepA_on_08.blastn
blastn -db HDRNA_12_K01_bold_37467.current.gb_converted.fna -query sepA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sepA_on_12.blastn
blastn -db HDRNA_16_K01_conservative_37834.current.gb_converted.fna -query sepA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sepA_on_16.blastn
blastn -db HDRNA_17_K01_conservative_37288.current.gb_converted.fna -query sepA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sepA_on_17.blastn
blastn -db HDRNA_19_K01_bold_37377.current.gb_converted.fna -query sepA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sepA_on_19.blastn
blastn -db HDRNA_20_K01_conservative_43457.current.gb_converted.fna -query sepA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sepA_on_20.blastn

fmtC.fasta
blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query fmtC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > fmtC_on_01.blastn
blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query fmtC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > fmtC_on_03.blastn
blastn -db HDRNA_06_K01_conservative_27645.current.gb_converted.fna -query fmtC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > fmtC_on_06.blastn
blastn -db HDRNA_07_K01_conservative_27169.current.gb_converted.fna -query fmtC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > fmtC_on_07.blastn
blastn -db HDRNA_08_K01_conservative_32455.current.gb_converted.fna -query fmtC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > fmtC_on_08.blastn
blastn -db HDRNA_12_K01_bold_37467.current.gb_converted.fna -query fmtC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > fmtC_on_12.blastn
blastn -db HDRNA_16_K01_conservative_37834.current.gb_converted.fna -query fmtC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > fmtC_on_16.blastn
blastn -db HDRNA_17_K01_conservative_37288.current.gb_converted.fna -query fmtC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > fmtC_on_17.blastn
blastn -db HDRNA_19_K01_bold_37377.current.gb_converted.fna -query fmtC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > fmtC_on_19.blastn
blastn -db HDRNA_20_K01_conservative_43457.current.gb_converted.fna -query fmtC.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > fmtC_on_20.blastn

sceD.fasta
blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query sceD.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sceD_on_01.blastn
blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query sceD.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sceD_on_03.blastn
blastn -db HDRNA_06_K01_conservative_27645.current.gb_converted.fna -query sceD.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sceD_on_06.blastn
blastn -db HDRNA_07_K01_conservative_27169.current.gb_converted.fna -query sceD.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sceD_on_07.blastn
blastn -db HDRNA_08_K01_conservative_32455.current.gb_converted.fna -query sceD.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sceD_on_08.blastn
blastn -db HDRNA_12_K01_bold_37467.current.gb_converted.fna -query sceD.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sceD_on_12.blastn
blastn -db HDRNA_16_K01_conservative_37834.current.gb_converted.fna -query sceD.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sceD_on_16.blastn
blastn -db HDRNA_17_K01_conservative_37288.current.gb_converted.fna -query sceD.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sceD_on_17.blastn
blastn -db HDRNA_19_K01_bold_37377.current.gb_converted.fna -query sceD.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sceD_on_19.blastn
blastn -db HDRNA_20_K01_conservative_43457.current.gb_converted.fna -query sceD.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > sceD_on_20.blastn

SE0760.fasta
blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query SE0760.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > SE0760_on_01.blastn
blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query SE0760.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > SE0760_on_03.blastn
blastn -db HDRNA_06_K01_conservative_27645.current.gb_converted.fna -query SE0760.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > SE0760_on_06.blastn
blastn -db HDRNA_07_K01_conservative_27169.current.gb_converted.fna -query SE0760.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > SE0760_on_07.blastn
blastn -db HDRNA_08_K01_conservative_32455.current.gb_converted.fna -query SE0760.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > SE0760_on_08.blastn
blastn -db HDRNA_12_K01_bold_37467.current.gb_converted.fna -query SE0760.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > SE0760_on_12.blastn
blastn -db HDRNA_16_K01_conservative_37834.current.gb_converted.fna -query SE0760.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > SE0760_on_16.blastn
blastn -db HDRNA_17_K01_conservative_37288.current.gb_converted.fna -query SE0760.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > SE0760_on_17.blastn
blastn -db HDRNA_19_K01_bold_37377.current.gb_converted.fna -query SE0760.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > SE0760_on_19.blastn
blastn -db HDRNA_20_K01_conservative_43457.current.gb_converted.fna -query SE0760.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > SE0760_on_20.blastn

esp.fasta
blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query esp.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > esp_on_01.blastn
blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query esp.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > esp_on_03.blastn
blastn -db HDRNA_06_K01_conservative_27645.current.gb_converted.fna -query esp.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > esp_on_06.blastn
blastn -db HDRNA_07_K01_conservative_27169.current.gb_converted.fna -query esp.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > esp_on_07.blastn
blastn -db HDRNA_08_K01_conservative_32455.current.gb_converted.fna -query esp.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > esp_on_08.blastn
blastn -db HDRNA_12_K01_bold_37467.current.gb_converted.fna -query esp.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > esp_on_12.blastn
blastn -db HDRNA_16_K01_conservative_37834.current.gb_converted.fna -query esp.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > esp_on_16.blastn
blastn -db HDRNA_17_K01_conservative_37288.current.gb_converted.fna -query esp.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > esp_on_17.blastn
blastn -db HDRNA_19_K01_bold_37377.current.gb_converted.fna -query esp.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > esp_on_19.blastn
blastn -db HDRNA_20_K01_conservative_43457.current.gb_converted.fna -query esp.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > esp_on_20.blastn

ecpA.fasta
blastn -db HDRNA_01_K01_conservative_23197.current.gb_converted.fna -query ecpA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > ecpA_on_01.blastn
blastn -db HDRNA_03_K01_bold_bandage_26831.current.gb_converted.fna -query ecpA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > ecpA_on_03.blastn
blastn -db HDRNA_06_K01_conservative_27645.current.gb_converted.fna -query ecpA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > ecpA_on_06.blastn
blastn -db HDRNA_07_K01_conservative_27169.current.gb_converted.fna -query ecpA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > ecpA_on_07.blastn
blastn -db HDRNA_08_K01_conservative_32455.current.gb_converted.fna -query ecpA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > ecpA_on_08.blastn
blastn -db HDRNA_12_K01_bold_37467.current.gb_converted.fna -query ecpA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > ecpA_on_12.blastn
blastn -db HDRNA_16_K01_conservative_37834.current.gb_converted.fna -query ecpA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > ecpA_on_16.blastn
blastn -db HDRNA_17_K01_conservative_37288.current.gb_converted.fna -query ecpA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > ecpA_on_17.blastn
blastn -db HDRNA_19_K01_bold_37377.current.gb_converted.fna -query ecpA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > ecpA_on_19.blastn
blastn -db HDRNA_20_K01_conservative_43457.current.gb_converted.fna -query ecpA.fasta -evalue 1e-10 -num_threads 15 -outfmt 6 -strand both -max_target_seqs 1 > ecpA_on_20.blastn

12+15+3-1=29

I only found ebpS instead of ebp.

        #(fumC) AND "complete cds"
        gyrB: https://www.ncbi.nlm.nih.gov/nuccore/MG995415.1 Mycobacterium tuberculosis strain UKR100 GyrB (gyrB) gene, complete cds

        fumC: https://www.ncbi.nlm.nih.gov/nuccore/3153897
        gltA:
        icd:

Bioinformatics Tools:

Artemis: A genome browser and annotation tool that can generate similar kinds of schematic representations.
* GenomeDiagram from Biopython: A toolkit within Biopython that allows for the creation of high-quality genomic graphics.
#--> https://biopython-tutorial.readthedocs.io/en/latest/notebooks/17%20-%20Graphics%20including%20GenomeDiagram.html
IGV (Integrative Genomics Viewer): While primarily a genome browser, IGV can be used to generate snapshots that display genomic regions and variations.
General Graphic Design Software:

Adobe Illustrator: A popular graphic design tool used to create precise and detailed scientific figures.
Inkscape: A free and open-source vector graphics editor that can be used to create diagrams like the one you've shown.
Specialized Genomic Visualization Tools:

SnapGene or Geneious: While mainly used for plasmid mapping and editing, these tools also offer features to create genomic maps and diagrams.
MEGA (Molecular Evolutionary Genetics Analysis): Primarily used for evolutionary analyses, MEGA can also produce graphical representations of genes and genomes.
Online Tools:

PlasMapper: Used primarily for plasmid maps but can be adapted for smaller genomic regions.
ApE (A plasmid Editor): Another tool for visualizing plasmid maps that might be used to create simplified genomic diagrams.
Without specific information, it's hard to pinpoint exactly which one was used for your image, but it's likely that the original creator used one of the bioinformatics tools or a combination of general graphic design software to illustrate the genomic rearrangements. If you're looking to create similar images, you might want to explore some of these options to see which one best suits your needs.

Microbial bioinformatics

Microbial bioinformatics uses computational tools to analyze genomes, track evolution, and study functions in microorganisms, including bacteria and viruses.

Typing of S. epidermidis samples (MD P.B.)

Leave a Reply Cancel reply