Skip to content
Snippets Groups Projects
Commit 12544f57 authored by Mike Smith's avatar Mike Smith
Browse files

bwa example

parent f1911724
No related branches found
No related tags found
No related merge requests found
#!/bin/bash
#SBATCH --mem 150
#SBATCH --nodes 8
## create an index
bwa index reference.fa
## perform alignment
bwa aln -I reference.fa reads.txt > out.sai
## create SAM file
bwa samse reference.fa out.sai reads.txt > out.sam
This diff is collapsed.
>example_reference_1
GCTGCCCCTGCTTTCCCAGAACCCGGAAGCTCTTCCCCACTTTTCCCAACCCCATGTCCCTGCCTCCCCT
CAGTTGTGGAGTTACAAGGACAGGCTGTGCTCATGCCAGGTTTGAACTGTGCTCTGGTCTCTCCCCAGTG
GCCCCTGTGGGTTACGGTTCAGGCAAAACCCACAGGGTGGTGTCCGCATCGTCGGCGGGAAGGCTGCACA
GCATGGGGCCTGGCCCTGGATGGTCAGCCTCCAGATCTTCACGTACAACAGCCACAGGTACCACACATGT
GGAGGCAGCTTGCTGAATTCACGATGGGTGCTCACTGCTGCTCACTGCTTCGTCGGCAAAAAGTACGTGT
AGGGATGCACTGAGGGAGGTCTTCAGAACGGCTCTTCTCAGAGAGGGGCGTTCCCCGGGGATGCTGTGCA
GCGTCTCCCTGGGGCTCTGGGCCAAGTGGCTGCAAGACTCCGGGGGCTGGTCCAGACCTTTGCTAGGGGA
AGGCCCTGAGGGTCGCTGTCACCAGGCTTTTGTCCAGCCGGTTGTGACCTGGCTTACCTTTGTGCCCACA
GTAATGTGCATGACTGGAGACTGGTTTTCGGAGCAAAGGAAATTACATATGGGAACAATAAACCAGTAAA
GGCGCCTCTGCAAGAGAGATATGTGGAGAAAATCATCATTCATGAAAAATACAACTCTGCGACAGAGGGA
AATGACATTGCCCTCGTGGAGATCACCCCTCCCATTTCGTGTGGGCGCTTCATTGGGCCGGGCTGCCTGC
CCCACTTTAAGGCAGGCCTCCCCAGAGGCTCCCAGAGCTGCTGGGTGGCCGGCTGGGGATATATAGAAGA
GAAAGGTGAGTATGGGAGCGCCTCCAAGGGGGGACGCTGCTGGCCATTCTCCTGGTGGTCTTTGAGGTGC
AGCGGTCACTTGTTGACACCCAGCCAGGCTGCTTTCATCCTCCTCACGGCGCTACACGTAGAGCCATCAC
TGTGGCCTTCCACAGTCCCCTGTGCCAGGTCACGTGATGGGTGACTCGTCTGGCTGTCTACGGGGGGGCT
GACAGCAGGTGCAGGCAGAGCGCAGCGTTGCTTAGAATGGGGTTGAGGCTGTGTCTGTATTTGGCACCTG
GGGCCCAAGCTGCTGGATCCTGGAAACAAGCAGCAACCCTGAAACAGCCTTTGGAGAAGAGCTCCAGATA
AGGTTATGTGGCCGTATGACAGTGCCTTCCACTCTCTGGGCCTTGGTTTCCTCATCTGTAAATGGAGGGG
CTGGATCAGATGGTCTCCAAGCTGACCATTTCCTAACACCTGTGATTCCAGGGCCCAGCCTTCCTGCTTC
CAGGCTGTCTTTGTCACCATCAGCAGCCACCTCCGGTGCAGTCATAGCACCAAATCTACCCTCTCCATCT
GAGTTACTTGACATTTGCTCTGTAGACTGGGAAGGGGACCGGTGGTTGGTGTGGCTTTCTGTGGGTCATA
GGCATCCCTGTATTTCAGGCTTCCTGGAGTGTCCGGCAGGTACCTCTCCATTCAAGACGGGTGTAGCTGT
CGCTCTAGAGCTCTGTCCATAACCAAGCTGTCTCTCTCTCCCCTGGCTCCCAGAAGATCTCTTGCCCATC
CAACACTACCTGATACACTCGACACCAGATGAGTTCTAGAAGAGTAGGGGGAATTTCTTCCGTACCAGAC
CAGGCTCACGGAGGGCCCCAGGTGTGACCTAGAACAAATTGTCATTACTGACTCCGTGGTCAGAGACGCC
CAGACACCTCTCCCTACATACTTAGGCACACTAGAGATACCATTTATCCCAGAAGCAACAGGTGACCCCA
TTAAATGGCAGCAGTGAGCAGCAAGGTTTTGCCAGCAGGACAAGACCCAGTAGGACTTGATCTCTGCATG
TGACCTATACTTCCAGCTTTCCATAAGTGGTGCATCTCAATCTCTTTTTTTCTTGATAATGGATAATGAG
GCTTCTTAATTCAAGGCAATTCAGTCATATGAGTCTAAAACGCAGCCTTTTTTGGGTAACAACAACAAAA
AGTCAATAATTTGAGGATTCTTCCTTTAACTCTTATTTCCTTCTCTGAGGGCGTGGTCTAAGTTTCTGAC
AGGTTTGTCACCTCTCCTGTCCCTCAGGGAAGGGAGTCTCATCCCAGCTGTCCTTGTTCTCACTGGCCTC
ACTTGACACGTGGCCATCCCCTTCTGGCCTGGAGTGACAGCCTCTATTCCAGCTGCCGATACCCATGGCT
CTGCCCTGGCCTCTGTCCAGCCCCCTAGCCTGCCATGCTCTGTCCCCTGTCTGATGATGTCCACAGCTCT
TGGGGGTCCCTATGTTCCTGGGAGGCTGATCTTAGGGTCCCCACTGCCAACCGTTGCTGTCCTATTCTAT
TACTCAGACAGATTGTAGTGTCCTAGCTATATCTTGGATGCCATGCAAAACCTTTTTTGTTCTGGAATTC
CCATAAACCTAATAGGATCACTCTGTTTCTGTTTGGGGTTTTTTTGTTTTTCATTTTTGCAAACACATCC
TGTACCCATAAGCCCAAGTCTTAAAGAAACCCCCAGGCTCCAACTCCTGTGATCCCTTCAAACCTTCCAT
GTCTTACCACCTCCTTCAACTCATCTGCTTTGGGGCTGGAAGGGTATTGCTTGGGTATTGCTTCTTTCCC
ATTGGGTTTTGAGTCCCTGCAGACTTTTAAGCATTTGTTATTTAAAGCCCAGGTGTCGGCTGGGTGCAGT
TGCTCATTCCTGTAATTTACCCAGCACTTTGGGAGGCCGAGGCAGGTGGATCACCTGAGGTCAGGAGTTT
GAGACCAGCCTGGCCAACATGGTGAAACCCAATCTCTACTAAAAATATAAAACATTAGCAGGGTGTGGTG
GTGGATGCCTGTTGTCCCAGCTACTCGGGAGGCTGAGGCAGGAGAATCACTTGAACCCGGGAGGTGAAGG
TTGCAGTGAGCCAAGATCACGCCACTGCACTCCAGCCTGGGTAACAAGAGCGAAACTTCATCTCAAAAAA
AACCCGCCAAGGTGTTGGGATTTAGAAAACCCTTTTCTGTTTGAGGACACCTGGCTCTTCACAATAAAAC
CCATGGCTTTCAAGACTATGGGCTCTGCCCTAGACTTAAAAAGTCATCTTGGCCGGGCGCGGTGGCTCAC
TCCTGTAATCCCAGCACTTTGGGAGGCCAAGGAGGGTGGATCACAAGGTCGTGAGATCGAGACCATCCTG
GCTAACACGGTGAAACCTCATTTCTACTAAAAATACAGAAAATTAGCCGGGCATGGTGGCGGGCGCCTGT
AGTCCCAGCTACTCGGGAGGCTGAGGCAGGAGAATGGCGTGAACCCAGGAGGCGGAGCTTGCAGTGAGCC
GAGATGGCGCCACCGCACTCCAGCCTGGGCGACAGAGCAAGACTCCGTCTCAAAAAAAAAAAAAAAAAAA
ATCGTCTTAATGGCTAATAAGCCAGGATTGTTGTCATTGTTGTTTCAATGTCACCCCTACCCCAAGGGCA
GTGCATGTACAAAACTCTCACTGCTGCTTGGATGGTAAGAAAAGGTGGGGGCTAAGGGGATTAAAAAACG
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment