20 Learning the BAM format
20.1 Introduction
SAMtools provides various (sub)tools for manipulating alignments in the SAM/BAM format. The SAM (Sequence Alignment/Map) format (BAM is just the binary form of SAM) is currently the de facto standard for storing large nucleotide sequence alignments. If you are working with high-throughput sequencing data, at some point you will probably have to deal with SAM/BAM files, so familiarise yourself with them! For the latest information on SAMtools, please refer to the release notes.
20.2 Installing SAMtools
Note: for this code to run with conda samtools, you have to activate conda in every chunk. So the better way to run it would be to copy and paste all of the chunks into your terminal
env /usr/bin/arch -x86_64 /bin/zsh --login
source $HOME/miniconda3/bin/activate
conda activate variantIf you have have installed samtools with homebrew (like I have) that is not necessary.
20.3 Basic usage
If you run samtools on the terminal without any parameters or with --help, all the available utilities are listed:Introduction to RNA-seq
Visualizing and counting of the alignments
20.4 The SAM/BAM format
SAM files, are tab-delimited text files that contain information for each individual read and its alignment to the genome.
The compressed binary version of SAM is called a BAM file. We use this version to reduce size and to allow for indexing, which enables efficient random access of the data contained within the file.
What’s in a SAM/BAM file
A BAM file is just a compressed SAM file. The file begins with a header, which is optional.The header is used to describe source of data, reference sequence, method of alignment, etc., this will change depending on the aligner being used. Following the header is the alignment section. Each line that follows corresponds to alignment information for a single read. Each alignment line has 11 mandatory fields for essential mapping information and a variable number of other fields for aligner specific information. An example entry from a SAM file is displayed below with the different fields highlighted.
Additionally tags (or attribute) can be added to each of of the lines. These tags give some additional information on the alignment. The number and type of tags varies between different alignment tools and the settings within these tools. The ‘CIGAR’ (Compact Idiosyncratic Gapped Alignment Report) string is how the SAM/BAM format represents spliced alignments. Understanding the CIGAR string will help you understand how your query sequence aligns to the reference genome. For example, the position stored is the left most coordinate of the alignment. To get to the right coordinate, you have to parse the CIGAR string.
To start of we’ll have a look at how to use samtools to have a peak at the the contents of the bam files.
SAMtools commands
SAMtools provides the following commands, each invoked as “samtools some_command”.
view The view command filters SAM or BAM formatted data. Using options and arguments it understands what data to select (possibly all of it) and passes only that data through. Input is usually a sam or bam file specified as an argument, but could be sam or bam data piped from any other command. Possible uses include extracting a subset of data into a new file, converting between BAM and SAM formats, and just looking at the raw file contents. The order of extracted reads is preserved.
As these file are binary you can not simply use:
head aln.bam
This will give an unreadable result. SAMtools can help us to make the content readable.
Let’s retrieve a bamfile to work with:
cd RNA-genomics/
source $HOME/miniconda3/bin/activate
conda activate variant
chmod 777 BAM/ERR458493.bam
cat BAM/ERR458493.bam > aln.bam## zsh:cd:2: no such file or directory: RNA-genomics/
source $HOME/miniconda3/bin/activate
conda activate variant
samtools view aln.bam | head20.5 Sorting
The sort command sorts a BAM file based on its position in the reference, as determined by its alignment.
Many downstream tools require sorted BAM files and since they are slightly more compact than unsorted BAM files, you should always sorted BAM files.
You can directly generate a sorted BAM file from a SAM file where -o defines the name of the output file (also a BAM file).
The default for samtools sort is sorting by position.
There are more sorting possibilities to be found with samtools sort –help.
source $HOME/miniconda3/bin/activate
conda activate variant
samtools sort aln.bam -o sorted.bamOnly the input file name needs to be specified, based on this name a .bai (BAM index) is produced.
source $HOME/miniconda3/bin/activate
conda activate variant
ls -lh aln.bam
ls -lh sorted.bam
samtools view aln.bam | head -n 1
samtools view sorted.bam | head -n 1## -rw-r--r-- 2 ggiaever staff 52M Dec 4 21:47 aln.bam
## -rw-r--r-- 2 ggiaever staff 45M Dec 4 21:47 sorted.bam
## ERR458493.748218 16 I 1944 0 51M * 0 0 AGTCACCGTAGTTGAAAACGGCTTCAGCAACTTCGACTGGGTAGGTTTAAG HJFHFJJJJIJJIIJIGGJJJIJJJIJJIIJJJJJHJJHHHHHFFDD=1:@ HI:i:1 NH:i:1 NM:i:1
## ERR458493.748218 16 I 1944 0 51M * 0 0 AGTCACCGTAGTTGAAAACGGCTTCAGCAACTTCGACTGGGTAGGTTTAAG HJFHFJJJJIJJIIJIGGJJJIJJJIJJIIJJJJJHJJHHHHHFFDD=1:@ HI:i:1 NH:i:1 NM:i:1You should use use additional threads (if they are available) to speed up sorting; to use four threads, use -@ 4.
Time taken using one thread (default).
source $HOME/miniconda3/bin/activate
conda activate variant
time samtools sort aln.bam -o sorted.bam## samtools sort aln.bam -o sorted.bam 2.73s user 0.11s system 99% cpu 2.870 totalTime taken using four threads.
source $HOME/miniconda3/bin/activate
conda activate variant
time samtools sort -@ 4 aln.bam -o sorted.bam## [bam_sort_core] merging from 0 files and 4 in-memory blocks...
## samtools sort -@ 4 aln.bam -o sorted.bam 2.89s user 0.21s system 377% cpu 0.822 totalMany of the SAMtools subtools can use additional threads, so make use of them if you have the resources!
20.6 Extracting entries mapping to a specific loci
Use samtools view and the ref:start-end syntax to extract reads mapping within a specific genomic loci; this requires a BAM index file.
20.7 Indexing
The index command creates a new index file that allows fast look-up of data in a (sorted) SAM or BAM. Like an index on a database, the generated .sam.sai or .bam.bai file allows programs that can read it to more efficiently work with the data in the associated files.
An index file is needed to get access rapidly to different alignment regions in the BAM alignment file.
Creating a BAM index file
Various tools require BAM index files, such as IGV, which is a tool that can be used for visualising BAM files.
source $HOME/miniconda3/bin/activate
conda activate variant
samtools index sorted.bamindex outputs two files, the index.bam file and an index.bam.bai file
source $HOME/miniconda3/bin/activate
conda activate variant
samtools index sorted.bam sorted.bai
source $HOME/miniconda3/bin/activate
conda activate variant
samtools view sorted.bam IX:20000-30000 ## ERR458493.829354 0 IX 19964 40 51M * 0 0 GTTGTAGGCATCCTAATCCAAATAACCTCCATAAACAAGTGGTACCAATAC ?@++=BDDHHHHHIIIIIIIIIIIIIIIIIIIIIIIIIIHIIHIIIIIIII HI:i:1 NH:i:1 NM:i:0
## ERR458493.592673 16 IX 20224 40 51M * 0 0 TCTTTGCTTTGCCTGGGCTACATTCATGGTTAGTGGGATGATGTTTGTACC #IJJJJJJIGJJJJJJIJJIJJJJJJJJJJJJJJJJJJHHHHHFFDB=1@B HI:i:1 NH:i:1 NM:i:0
## ERR458493.285385 16 IX 20310 40 51M * 0 0 ATGAGGAAGCGAAACATTCACTATCGAAATCCAACAAAGTCTCAGTTGACG HIJIJJJJJJGIIGJJJIJIIGJHJJJHEHIIJJJIIJHHHHHFFDA+BBB HI:i:1 NH:i:1 NM:i:0
## ERR458493.43268 16 IX 20433 13 51M * 0 0 CATGGTCTGAACTACTCTCCACTAAAACAAAGGTCTTTCAGCGTGTTCTTA IJJJJJJJJIIGFBGIIGGGDIIIIIIJIJJIIJIJIIHGHHHFFFFF@CC HI:i:1 NH:i:1 NM:i:0
## ERR458493.888408 0 IX 20578 20 51M * 0 0 CCCCTTTCAAACTTCGATCATCATTGGTGTGGTCAATTTTTTCTCTTCATT @+=;DDDDHAHFFEHI3AEGII>FGIIFGFEH?FHGHGIIIIGIHGGIIDH HI:i:1 NH:i:1 NM:i:1
## ERR458493.322479 16 IX 20643 0 1S50M * 0 0 CCCATTGAGAGGTTTGGACGCCGTACGTGTCTATTGTGGGGTGCTGCTTCT #A=.7CEEIIFFDC;<(@6:>FFE?G@IIEHEC3@CDCBFC>DDA3DA=?; HI:i:1 NH:i:1 NM:i:0
## ERR458493.284879 16 IX 20906 0 51M * 0 0 GCAATCGCAACAGCTGCAAACTGGATGTGGGGTTTCCTGATTAGTTTCTTT IJJJJJJJIGJJJIIIJIGEJIIEIJJJIIIHHGIJJJHHHHHEDFFD@B@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.472218 0 IX 21069 40 51M * 0 0 CGCTGGAAGAGGTAAATACTATGTGGCTGGAAGGTGTGCCAGCATGGAAAT @B1DDFFFHHHGFGFIIIIJJJJIJJJIIJEIGIDGDBBF?BGBFHFGIIA HI:i:1 NH:i:1 NM:i:0
## ERR458493.623067 0 IX 21639 0 51M * 0 0 TGGGTTTTTCATGATGGAAATGTGATACTACCAGTTCCAATATATTTATCT 111442ADHHFGDGHGBGGHHEGCHIIGHG@GGEEFGCIGGGIADCFDFFG HI:i:1 NH:i:1 NM:i:0
## ERR458493.1087136 0 IX 21646 0 51M * 0 0 TTCATGATGGAAATGTGATACTACCAGTTCCAATATATTTATCTTCCTTAT @@@?DDDDH<+<CD@<EHBH>EFEGGE3+<CEGFGFGEHGIBEH?DGIHFE HI:i:1 NH:i:1 NM:i:0
## ERR458493.98492 0 IX 21660 0 51M * 0 0 GTGATACTACCAGTTCCAATATATTTATCTTCCTTATCTATATGACACGCT ?1:BDFFFHHHHHGHIJJJIIJJJJJIJJJJJJJJIJJJJJJJIJJJJJJI HI:i:1 NH:i:1 NM:i:0
## ERR458493.408716 0 IX 21701 40 51M * 0 0 ATGACACGCTGTTATTTTAGTTCAAGTCAGTGTCCAATTGAAGTGAGTGTC @@@FFFADHFHFHGHIIIJJFFDEIJDFIGCFEE@FHGHIGGGBDBHFEFI HI:i:1 NH:i:1 NM:i:0
## ERR458493.898831 16 IX 21716 40 51M * 0 0 TTTAGTTCAAGTCAGTGTCCAATTGAAGTGAGTGTCACGTGCACAGCGACG EDD8FB9CIIFGG<@IFCFD@E>C<G<FFEEFFA;?GFFC:FF<DDDB:<? HI:i:1 NH:i:1 NM:i:0
## ERR458493.659428 0 IX 22697 0 51M * 0 0 GAACGATTTTACTGTGGACAATAAGTGAAATAAGTGTCCTTAAAATGTGTA @@@?D?DD7D?D<AFE9<C<?+CCBAE+AFH>F9:4CF<<C@D?@B**::? HI:i:1 NH:i:1 NM:i:0
## ERR458493.701162 0 IX 22724 20 51M * 0 0 AAATAAGTGTCCTTAAAATGTGTACGATGTGTACACATCAACCTACTCTCC @@CFFFFBFFHHHIJJGJIEIIIIGHHGIHGIEHIIHGEGIIJIEIGDHHI HI:i:1 NH:i:1 NM:i:0
## ERR458493.620175 0 IX 22961 0 51M * 0 0 GTAAGAACGAAAATACCGACCAATGCATTCAAAGTATTGGCATTCAAACCT @?1:BDDDHHHHHBFE<??CFGCHGIHGBFFHEGFDFEGGDDFF@FGGFGF HI:i:1 NH:i:1 NM:i:1
## ERR458493.617523 0 IX 22991 0 51M * 0 0 AAAGTATTGGCATTCAAACCTGCGGCAATACCCTCGGACTGGACAATGATA 8=+A4BDDFDDFFIIIIFEFFEFIF8FGFFBFFIFFFFEFIICFFIFFFII HI:i:1 NH:i:1 NM:i:0
## ERR458493.987667 16 IX 23004 0 51M * 0 0 TCAAACCTGCGGCAATACCCTCGGACTGGACAATGATACCAGTAGAGGTCT JJIGFIGGGGIGIHGBFHHFDGIIHJIJIHIGGGHEFCGHHFHEDDBD@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.39663 16 IX 23038 20 47M4S * 0 0 GATACCAGTAGAGGTCTTTATTGGAGCTGCAGATTTGGAGTATGTGGAGAT DHHDEIIJIHDEHGEHFIGJIFIIIIJJJIIIJJIIIGHHHHFEDDB=1B? HI:i:1 NH:i:1 NM:i:0
## ERR458493.1037765 0 IX 23084 20 17S34M * 0 0 GCAGATTTGGAGCATGTGGAAATTGCACTAGAAGCCTTGGAAATGGCAGTA 8:?D72B?CA3A+AEECEBDA@>FFA9FEEECAAEECDDIBDC>DDDDD9B HI:i:1 NH:i:1 NM:i:0
## ERR458493.861153 16 IX 23136 0 51M * 0 0 GAACCCTCAGAAGTGACAGTAGTGTATGACTTTGGAGAAGCACTAGTTGCG IDIDDDADEIIIIIFEIIEEDEBFADBA3EAA+4DE?EBDC<D>DADD??? HI:i:1 NH:i:1 NM:i:0
## ERR458493.346191 16 IX 23154 13 10S33M8S * 0 0 TTTTATGGATGAAGTGTATGACTTTGGAGAAGCACTAGTTGCGGTGACGAT GJJIIJJJJFIIGJJJJIIHDIJJJJJJJJJIHGJJJJHHHHHFFFFFC@@ HI:i:1 NH:i:1 NM:i:1
## ERR458493.432364 16 IX 23286 0 51M * 0 0 TAGGGGCAGCACTGGTTTCTGAAGTTTCTGAAGTTTCTTTAGGAGCTTCAG CCD;F?@EEFFG9IGF?B>HHFH?9AC<>F:CBFBABFF??D?=B;DD<@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.972815 0 IX 23297 0 51M * 0 0 CTGGTTTCTGAAGTTTCTGAAGTTTCTTTAGGAGCTTCAGAAGTGACAGTC CBCFFFFFHHHHHIJJJJJJJJJJJJJJJIJJIJIJJJJJIJJGHIIJJHG HI:i:1 NH:i:1 NM:i:0
## ERR458493.380048 16 IX 23298 0 51M * 0 0 TGGTTTCTGAAGTTTCTGAAGTTTCTTTAGGAGCTTCAGAAGTGACAGTCT HIIIIGHFGAEEIIIIIIIIGFEIHGIIHGEAGGGIIIGHHHGFFDDB1C@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.122587 0 IX 23299 0 51M * 0 0 GGTTTCTGAAGTTTCTGAAGTTTCTTTAGGAGCTTCAGAAGTGACAGTCTT B@B4ADEFHHHHHJJJJJJJJJJJJJJJJJIJJIJJJJJJJGHIIIJFHII HI:i:1 NH:i:1 NM:i:0
## ERR458493.1026234 0 IX 23299 0 51M * 0 0 GGTTTCTGAAGTTTCTGAAGTTTCTTTAGGAGCTTCAGAAGTGACAGTCTT =@:44ADEHHHHHJJJJJJJJJJJGJJJJJJJJIGJIJJJIHIIIIJIIIJ HI:i:1 NH:i:1 NM:i:0
## ERR458493.6253 16 IX 23345 20 51M * 0 0 GTCTTGGTGCTACATCCATTGTCATCGCACTGAGTAACAGTAGCGGTGGTG IJJJHGHHAIIIIJJIHHF@HBIGIGEGHJHGFGIJJIHFGGHFFDFF@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.623523 0 IX 23353 40 51M * 0 0 GCTACATCCATTGTCATCGCACTGAGTAACAGTAGCGGTGGTGTATGTCTT ?@@DDFFFGHGHDIIJJJJGGHHJGJJJJJJIHIIIJJ@FHGHFHIJIIJI HI:i:1 NH:i:1 NM:i:0
## ERR458493.814078 0 IX 23353 40 51M * 0 0 GCTACATCCATTGTCATCGCACTGAGTAACAGTAGCGGTGGTGTATGTCTT ?+1:BDFFHFFFHGEGGGEHIIIIHICHIIIIDG@GII?DB@G?@BG?FFF HI:i:1 NH:i:1 NM:i:0
## ERR458493.813613 16 IX 23355 40 51M * 0 0 TACATCCATTGTCATCGCACTGAGTAACAGTAGCGGTGGTGTATGTCTTTG JJIHIGGGIJIGIGIGDDCGIIJIJIIIGIIJIGIHHGGHHGHFFFDD@CC HI:i:1 NH:i:1 NM:i:0
## ERR458493.407288 16 IX 23359 40 51M * 0 0 TCCATTGTCATCGCACTGAGTAACAGTAGCGGTGGTGTATGTCTTTGGAGA 8BGHFBCEFBAGAFBHEIHGEHHGBGIGHAEGGFEHDGHFHHDDDDAB=<@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.701205 0 IX 23369 20 36M15S * 0 0 TCGCACTGAGTAACAGTAGCGGTGGTGTATGTCTTTTTTGAAACAGTAGTG @@@DFFFFGFHFHIEBF@FGGI@EFEF@FIC?FGHEHIJDFHGJGEHHGFF HI:i:1 NH:i:1 NM:i:0
## ERR458493.159467 0 IX 23371 20 51M * 0 0 GCACTGAGTAACAGTAGCGGTGGTGTATGTCTTTGGAGAAACAGTAGTGGT CCCFFFFFHHHHHJIJJIJJCGICGCGIIHIIIJIJGIIJJJJJHIGFHII HI:i:1 NH:i:1 NM:i:0
## ERR458493.934957 16 IX 23390 13 51M * 0 0 GTGGTGTATGTCTTTGGAGAAACAGTAGTGGTTGTGGCTTCAGGGGCCTCA JIJIJJJJJJIHJJIJJJJJIGJJJJJJJJJJJJJIJIHHGGHFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.601349 0 IX 23392 13 51M * 0 0 GGTGTATGTCTTTGGAGAAACAGTAGTGGTTGTGGCTTCAGGGGCCTCAGA =+1=:BDDCFHHHI>BE@HGGIIEGGCBG<CEEGI9DGBB<:BGBG6?### HI:i:1 NH:i:1 NM:i:0
## ERR458493.401885 0 IX 23396 0 51M * 0 0 TTTGTCTTTGGAGAAACAGTAGTGGTTGTGGCTTCAGGGGCCTCAGAGGTG @B@FFFFFHHHHHIJJIJJHJIHIJFHIGIJIIJC>CGGI0?BBBEEGG## HI:i:1 NH:i:1 NM:i:1
## ERR458493.606981 16 IX 23398 0 49M2S * 0 0 TGTCTTTGGAGAAACAGTAGTGGTTGTGGCTTCAGGGGCCTCAGAGGAGCC B>@B/?<BGB?BB??B?GCF;D?:?3FC<?GF;BAF?F=FFAFA=24411: HI:i:1 NH:i:1 NM:i:1
## ERR458493.156699 0 IX 23404 0 51M * 0 0 TGGAGAAACAGTAGTGGTTGTGGCTTCAGGGGCCTCAGAGGTGACAGGCTT ?<+=ADD??FDDBAGBG?CA@CDGFE;*1*:?0@)0****89B4?<8(78F HI:i:1 NH:i:1 NM:i:1
## ERR458493.907482 0 IX 23444 0 51M * 0 0 GTGACAGTCTTAGTGTTACAGCCGTTGTCGTCACAGTGAGTAACGGTAACA @BCFFFFDFHHHHFIJHHHJJJIJHIJJGIIIGIGGDDGI<GHJJJGGIIG HI:i:1 NH:i:1 NM:i:0
## ERR458493.263907 0 IX 23445 0 51M * 0 0 TGACAGTCTTAGTGTTACAGCCGTTGTCGTCACAGTGAGTAACGGTAACAG ?@@DDDDBHDHHAE<CGGICGCEFHHBFHCGIHBE??FHBFHGGGBBCDAC HI:i:1 NH:i:1 NM:i:0
## ERR458493.916041 0 IX 23445 0 51M * 0 0 TGACAGTCTTAGTGTTACAGCCGTTGTCGTCACAGTGAGTAACGGTAACAG CCCFFFFFHHHHHIIIJJJJJJJIJJJJJGIJJJJGHIJFHIIJJHGIJII HI:i:1 NH:i:1 NM:i:0
## ERR458493.661806 0 IX 23451 8 51M * 0 0 TCGTAGTGTTACAGCCGTTGTCGTCACAGGGAGTAACGGTAACAGTAGTGT :@+4AABDDD?FFBFA@C?EGHBGE:CF@):DABFHIIFFEGG=C.=FB@8 HI:i:1 NH:i:1 NM:i:2
## ERR458493.340509 0 IX 23452 13 51M * 0 0 CTTAGTGTTACAGCCGTTGTCGTCACAGTGAGTAACGGTAACAGTAGTGTA BCCFFFFFHHHHHJJJJJGIJJIIIJJJGIGJGHIJJJFHIJIIFIIIIJI HI:i:1 NH:i:1 NM:i:0
## ERR458493.1047071 0 IX 23452 0 51M * 0 0 CTTAGTGTTACAGCCGTTGTCGTCACAGTGAGTAACGGTAACAGTAGTGAA ?@@D;,ABDBF<B@EF@E:EEF+A?;**:C<?1???:0?)?D9*/??#### HI:i:1 NH:i:1 NM:i:1
## ERR458493.338841 0 IX 23453 13 51M * 0 0 TTAGTGTTACAGCCGTTGTCGTCACAGTGAGTAACGGTAACAGTAGTGTAT @@@FDFDFHHDHHEHHGGHGGHJJIGGCGAHFHGGGJDDGGBFFGGHHBFH HI:i:1 NH:i:1 NM:i:0
## ERR458493.445564 16 IX 23460 0 51M * 0 0 TACAGCCGTTGTCGTCACAGTGAGTAACGGTAACAGTAGTGTATGATTTTG C;IGGDIEGEHHCHG>EBIIHEHBCHFDBGHIEIHGIIHHHGFFFFFD@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.267132 16 IX 23460 0 51M * 0 0 TACAGCCGTTGTCGTCACAGTGAGTAACGGTAACAGTAGTGTATGATTTTG FDJIIGIGDJIIIHGIIJIIIIIIHEFIGJHHHIGIIJHHGGHFFFDDCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.213255 16 IX 23460 0 51M * 0 0 TACAGCCGTTGTCGTCACAGTGAGTAACGGTAACAGTAGTGTATGATTTTG HEJIIGIIIIIIJIIHFJJJJJJJJJHJJJJJJJIJJJHHHHFFFFFFCC@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.168583 16 IX 23461 0 51M * 0 0 ACAGCCGTTGTCGTCACAGTGAGTAACGGTAACAGTAGTGTATGATTTTGG AIIGGIIJIIHFHGFDGHIHGHHIIHIIJIHHGIIJIJFGHAHFFFFF@@B HI:i:1 NH:i:1 NM:i:0
## ERR458493.759829 16 IX 23463 0 51M * 0 0 AGCCGTTGTCGTCACAGTGAGTAACGGTAACAGTAGTGTATGATTTTGAAG @@0??0)0GDCF:19AEG?EC<82@:ACF<?@FA:<?A<?FDD?B==A1?: HI:i:1 NH:i:1 NM:i:1
## ERR458493.403662 16 IX 23463 0 51M * 0 0 AGCCGTTGTCGTCACAGTGAGTAACGGTAACAGTAGTGTATGATTTTGGAG BJJIJJIJIIIGIFGJJJJJJJIHJJIGJJJJJJIJJJHHHHHFDDDD:?B HI:i:1 NH:i:1 NM:i:0
## ERR458493.316932 16 IX 23518 13 51M * 0 0 AGTAGTTGCTGAGGTTTCTTCAGGACATTCAGAAGTGACAGTCTTGGTGCT CF@9B?<DCGF;C?11*11<3FA2A:A::?HI<C?CDF?<BFFB?A=4+@= HI:i:1 NH:i:1 NM:i:1
## ERR458493.703696 16 IX 23519 20 51M * 0 0 GTAGTTGCTGAAGTTTCTTCAGGACATTCAGAAGTGACAGTCTTGGTGCTA HIIIHGIHGGGHCIIGIIIHAE@IIDHDIHGGHDEHGBAF?HDED>DF@?@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.377568 16 IX 23531 13 51M * 0 0 GTTTCTTCAGGACATTCAGAAGTGACAGNCTTGGTGCTACATCCATTGACA 9GE@F@EHDBDF@@HB>HCHF@B@AA<4#HCHHC9BF<HFHF<DB:=A1@? HI:i:1 NH:i:1 NM:i:2
## ERR458493.204363 16 IX 23555 20 51M * 0 0 ACAGTCTTGGTGCTACATCCATTGTCATCACACTGAGTAACAGTAGCGGCG ###D88?):4>EDBD3AECIECIIEFEA22<3C22:EA<D<D:@<=:A1<? HI:i:1 NH:i:1 NM:i:1
## ERR458493.579443 16 IX 23568 40 51M * 0 0 TACATCCATTGTCATCACACTGAGTAACAGTAGCGGTGGTGTATGTCTTTG AFIGIGIHFCHJIIEGDGGIHHJJGIGJJJIIIIIJIIHHHHHFFFFF@<B HI:i:1 NH:i:1 NM:i:0
## ERR458493.569900 16 IX 23568 40 51M * 0 0 TACATCCATTGTCATCACACTGAGTAACAGTAGCGGTGGTGTATGTCTTTG IHIIIIIIIIIIIIHGFGIIGIIIIIIIIIIIIIIIIGHHHHHDDDB=;@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.739357 16 IX 23572 40 51M * 0 0 TCCATTGTCATCACACTGAGTAACAGTAGCGGTGGTGTATGTCTTTGGAGA IJJJJJJJJIHH>JIIJJJJJJJJJJJJJJJJJJJJJJHHHHHFFFDDBCB HI:i:1 NH:i:1 NM:i:0
## ERR458493.3109 0 IX 23587 0 51M * 0 0 CTGAGTAACAGTAGCGGTGGTGTATGTCTTTGGAGAAACAGTAGTGGTTGT @@@D?;BDHFHABA>G@@EF@@8ABAEHIIIIGAGEB?BDBDH@?BF0BFG HI:i:1 NH:i:1 NM:i:0
## ERR458493.713771 0 IX 23592 0 51M * 0 0 TAACAGTAGCGGTGGTGTATGTCTTTGGAGAAACAGTAGTGGTTGTGGCTT CCCFFFDFHHHHDGI@FFIIIIIIJJJJFIDHIJJJHIJIIJDGIFGIIGE HI:i:1 NH:i:1 NM:i:0
## ERR458493.740834 16 IX 23592 0 51M * 0 0 TAACAGTAGCGGTGGTGTATGTCTTTGGAGAAACAGTAGTGGTTGTGGCTT JJIJJJJJJJJJJGJJHIIHJJJJJJJJJJJJIJJJJJHHHHHFFFFF@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.414524 16 IX 23595 0 51M * 0 0 CAGTAGCGGTGGTGTATGTCTTTGGAGAAACAGTAGTGGTTGTGGCTTCGG HEJJIHIGHIGHJJJIGGGGGJIIIHGJIHIIIGGGGGHHHGHDDDBA:1@ HI:i:1 NH:i:1 NM:i:1
## ERR458493.486896 16 IX 23596 0 51M * 0 0 AGTAGCGGTGGTGTATGTCTTTGGAGAAACAGTAGTGGTTGTGGCTTCAGG IIIIJJJJJJIJIIJIJIJJIGIJIHIJIJIJJJJJJJHHHGHFFFFF@B@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.154668 0 IX 23601 13 51M * 0 0 CGGTGGTGTATGTCTTTGGAGAAACAGTAGTGGTTGTGGCTTCAGGAGCTT ?+:BDDADCDCD?AEEEEE3CFEEIEI?AACEE:CDDDIC;D30?99;(9D HI:i:1 NH:i:1 NM:i:0
## ERR458493.383391 16 IX 23604 13 51M * 0 0 TGGTGTATGTCTTTGGAGAAACAGTAGTGGTTGTGGCTTCAGGAGCTTCAG IJJJJJJJJJJIIIJJJIIHHJJJJJIJJJIJJJJJJJHHHFHDDFFD@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.872994 16 IX 23611 40 51M * 0 0 TGTCTTTGGAGAAACAGTAGTGGTTGTGGCTTCAGGAGCTTCAGAAGTGAC AHFGEGIIGIGEEHFDHGGHCFFGHHDEGGHEF>GGGIHDHFHFEAFD@?@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.163292 0 IX 23617 40 51M * 0 0 TGGAGAAACAGTAGTGGTTGTGGCTTCAGGAGCTTCAGAAGTGACAGTCTT @@@DDDBDFDD:CG?FDFEF@E=@8F3*11*)))1???DD;:?499*99B# HI:i:1 NH:i:1 NM:i:0
## ERR458493.938810 0 IX 23622 13 51M * 0 0 AAACAGTAGTGGTTGTGGCTTCAGGAGCTTCAGAAGGGACAGTCTTGGTCT @CCFDFEFHHHHFHIHIJJIJ;AFH<?3?CBDF9F@)::6BD?FB??F*9? HI:i:1 NH:i:1 NM:i:1
## ERR458493.50335 16 IX 23625 40 51M * 0 0 CAGTAGTGGTTGTGGCTTCAGGAGCTTCAGAAGTGACAGTCTTGGTCTTAC JJJJJJJJJJIJJJJIJJJJJJJIIJJJJJJJJJIHJJHHHHHFFFDD=:B HI:i:1 NH:i:1 NM:i:0
## ERR458493.385900 0 IX 23629 40 51M * 0 0 AGTGGTTGTGGCTTCAGGAGCTTCAGAAGTGACAGTCTTGGTCTTACATCC =@?DDDDDCFHDHIBEHI@HHIII@GGGG*CDCB?9:4?*:*?@GGECFHF HI:i:1 NH:i:1 NM:i:0
## ERR458493.528857 0 IX 23650 40 51M * 0 0 TTCAGAAGTGACAGTCTTGGTCTTACATCCATTGTCATCACAGTGGGTGAC <@@DFFFDDAFDAE+<CGIJG@GHICH?DGGCGGHIIBFBD*:?DCH:DBC HI:i:1 NH:i:1 NM:i:0
## ERR458493.42328 16 IX 23657 40 51M * 0 0 GTGACAGTCTTGGTCTTACATCCATTGTCATCACAGTGGGTGACGGTGGCA JJIGIGJIGJIJIGGIJIJJGJIJJJJIJIHHHJIJJJHHHHHFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.672512 16 IX 23658 20 51M * 0 0 TGACAGTCTTGGTCTTACATCCATTGTCATCACAGTGGGTGACCGTGGCAG ##FFBFB?2B;3F?FFC9D3FE?*3B@EA+?E>FF;DAC?1C@)DB17;<? HI:i:1 NH:i:1 NM:i:1
## ERR458493.42694 16 IX 23670 20 51M * 0 0 TCTTACATCCATTGTCATCACAGTGGGTGACGGTGGCAGTGGTGTACAATT IJJIHGIIJGIIJJJIIIIIJJJJJJJJJIJJJIJJJJHHHHHFFDA+B@B HI:i:1 NH:i:1 NM:i:1
## ERR458493.481295 16 IX 23673 20 51M * 0 0 TACATCCATTGTCATCACAGTGGGTGACGGTGGCAGTGGTGTACGATTGGG BBFFB@99BD>B949GBGEFFG?=G:C13CEF<DIFF>DFDFFDDDDA1+= HI:i:1 NH:i:1 NM:i:1
## ERR458493.1092013 16 IX 23673 20 51M * 0 0 TACATCCATTGTCATCACAGTGGGTGACGGTGGCAGTGGTGTACGATTAAG DDDCEEDD9IICEEDBD@DDDCEEC>FFF9EFCEB@IIDDBDA;DA==+?? HI:i:1 NH:i:1 NM:i:1
## ERR458493.508969 16 IX 23687 40 51M * 0 0 TCACAGTGGGTGACGGTGGCAGTGGTGTACGATTGAGAAGAAACAGTAGTT HHFGJJJIHEIGIIJJJJJJIJJJJJJJIJIIHGJJJIHHGHHFFEDD@B@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.47894 16 IX 23690 40 51M * 0 0 CAGTGGGTGACGGTGGCAGTGGTGTACGATTGAGAAGAAACAGTAGTTGTG GGHGIGGEFDGIJIGGGJIJIIHFEDHBGEJIIJJJGHHFHDHFFDBD@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.719618 0 IX 23696 40 51M * 0 0 GTGACGGTGGCAGTGGTGTACGATTGAGAAGAAACAGTAGTTGTGGTGGTA @@@FFFFDHHGHHFHIAGFHHJCHIIHIIJJJJIJIJDGHBFHHIJ;BG@F HI:i:1 NH:i:1 NM:i:0
## ERR458493.893751 16 IX 23696 40 51M * 0 0 GTGACGGTGGCAGTGGTGTACGATTGAGAAGAAACAGTAGTTGTGGTGGTA IFD8@BCHG@F>F?9FBFF?;GEFF>HHDF<HEA:FHDD7;<FADD=D<@; HI:i:1 NH:i:1 NM:i:0
## ERR458493.561609 0 IX 23702 40 51M * 0 0 GTGGCAGTGGTGTACGATTGAGAAGAAACAGTAGTTGTGGTGGTAGCTTCA =1=?DDDBFF?DAEHGBBFG>EDHCHIIID>FGGGGEEEG:BF?B?39D?F HI:i:1 NH:i:1 NM:i:0
## ERR458493.491404 16 IX 23704 40 51M * 0 0 GGCAGTGGTGTACGATTGAGAAGAAACAGTAGTTGTGGTGGTAGCTTCAGG 33CGD?GFAHF?F>>GEGG@CGDF<9HGIIIIDIGEEEC3DBD8DD7D@@? HI:i:1 NH:i:1 NM:i:0
## ERR458493.229736 16 IX 23704 40 51M * 0 0 GGCAGTGGTGTACGATTGAGAAGAAACAGTAGTTGTGGTGGTAGCTTCAGG @BHGADAGGHDCGGHCE>F@HFBGGHCGHB>EEBGGGE?>C<DBBADA11? HI:i:1 NH:i:1 NM:i:0
## ERR458493.814112 0 IX 23713 40 51M * 0 0 GTACGATTGAGAAGAAACAGTAGTTGTGGTGGTAGCTTCAGGAGCTTCAGA @1=DDADFFDDHHJJIJJGIBHGIIDFHIFGHGCIEEHCHJG<BG@GIJGC HI:i:1 NH:i:1 NM:i:0
## ERR458493.874578 0 IX 23713 40 51M * 0 0 GTACGATTGAGAAGAAACAGTAGTTGTGGTGGTAGCTTCAGGAGCTTCAGA @@CDFFFFHAFFHIGIIGGG:<CHHIHEG?EGDGGHIIIIII>GIIIIGHI HI:i:1 NH:i:1 NM:i:0
## ERR458493.26779 0 IX 23713 40 51M * 0 0 GTACGATTGAGAAGAAACAGTAGTTGTGGTGGTAGCTTCAGGAGCTTCAGA @:?:DDDFHFHFHIJJJJGIFHIIIJHGICGIEHIFAFHIJJGIIIGIGI< HI:i:1 NH:i:1 NM:i:0
## ERR458493.186740 0 IX 23713 40 51M * 0 0 GTACGATTGAGAAGAAACAGTAGTTGTGGTGGTAGCTTCAGGAGCTTCAGA @B@FFDDFHFHDFIJJIIIIHGHHHHIJIIIG<GH@FHGIHJGIGIJIIHF HI:i:1 NH:i:1 NM:i:0
## ERR458493.174939 16 IX 23735 20 51M * 0 0 GTTGTGGTGGTAGCTTCAGGAGCTTCAGAGGTGACAGTCTTGGTGGTACAG JJJJJIJJJJJJIJJJJIJIJJJJJJJJJJJJJJJJJJHGHHHFFFFFC@C HI:i:1 NH:i:1 NM:i:1
## ERR458493.361354 16 IX 23740 40 51M * 0 0 GGTGGTAGCTTCAGGAGCTTCAGAGGTGACAGTCTTGGTGTTACAGCCATT IIGEGGGGGIIGGHGGHEGBGIIIIIHF?HGHHEIIIIHHHHHDDDB4+@? HI:i:1 NH:i:1 NM:i:0
## ERR458493.1058757 16 IX 23744 40 51M * 0 0 GTAGCTTCAGGAGCTTCAGAGGTGACAGTCTTGGTGTTACAGCCATTGTCG JJJIJJIIIJIGF@IFJGGGIJGGGJIGHCEGJIJJIIHHGHHFFDD7@<@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.1000155 16 IX 23746 40 51M * 0 0 AGCTTCAGGAGCTTCAGAGGTGACAGTCTTGGTGTTACAGCCATTGTCGTC IIIGGJJJJIHEGGJIJJJJJIIJJIGJJIIJJJIHGJFHHGHFFDDFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.992368 16 IX 23746 20 51M * 0 0 AGCTTCAGGAGCTTCAGAGGTGACAGTCTTGGTGTTACAGCCATTGTCGAC IJJJJIJJIJIGIHIIIIJJIHEJJIIIJJJJIJIJJIHHHHHFFFDA:1@ HI:i:1 NH:i:1 NM:i:1
## ERR458493.331766 0 IX 23751 40 51M * 0 0 CAGGAGCTTCAGAGGTGACAGTCTTGGTGTTACAGCCATTGTCGTCACAGT ?@<D1BDDFF<?BGEE=2AA@FFEFHFCFFF>BC99CG>DEG*??0?DDD9 HI:i:1 NH:i:1 NM:i:0
## ERR458493.947474 0 IX 23754 40 51M * 0 0 GAGCTTCAGAGGTGACAGTCTTGGTGTTACAGCCATTGTCGTCACAGTGGG @?<DDFDDDHFFFAHDGGIGIIIJHHGGDHIEHGGIIIHEHFGGIIIBDDB HI:i:1 NH:i:1 NM:i:0
## ERR458493.169174 16 IX 23754 40 51M * 0 0 GAGCTTCAGAGGTGACAGTCTTGGTGTTACAGCCATTGTCGTCACAGTGGG B=>DEED<D<DB4:0**3D;DD??111*+C9F@?3<AED:<DC;D>DD<<? HI:i:1 NH:i:1 NM:i:0
## ERR458493.363103 0 IX 23757 40 51M * 0 0 CTTCAGAGGTGACAGTCTTGGTGTTACAGCCATTGTCGTCACAGTGGGTTA @@=DDDDEHDDHFHIHIIJIIDFHGGHGGFIJIBHIEGGHJJHGBGFG?DH HI:i:1 NH:i:1 NM:i:0
## ERR458493.206013 16 IX 23758 40 51M * 0 0 TTCAGAGGTGACAGTCTTGGTGTTACAGCCATTGTCGTCACAGTGGGTTAC GIEGIHJJHIGJJIGGIEIGIJIIHJIIJJJJJJIJJIHFHHHFFFFD<@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.546100 0 IX 23775 13 51M * 0 0 TGTTGTTACAGCCATTGTAGTCACAGTGGGTTACGGTGACAGTGGTGTATG ?<:4=BDDFFAFFFFFFG,<CCHHHEHHIEGIIIII3CBFFFFGFFGFIFF HI:i:1 NH:i:1 NM:i:2
## ERR458493.8635 0 IX 23775 40 51M * 0 0 TGGTGTTACAGCCATTGTCGTCACAGTGGGTTACGGTGACAGTGGTGTATG @?@?DDDEFDHHHIIG:AFHIG?HGH<ED=C:?DDG?D)99??<B9D;BFG HI:i:1 NH:i:1 NM:i:0
## ERR458493.680060 0 IX 23776 40 51M * 0 0 GGTGTTACAGCCATTGTCGTCACAGTGGGTTACGGTGACAGTGGTGTATGA @@?B;DDFGHHGDGIIGIIHGIIGGIIII:EFHGIIGDH>DBGG?B?FHHE HI:i:1 NH:i:1 NM:i:0
## ERR458493.672858 0 IX 23776 40 51M * 0 0 GGTGTTACAGCCATTGTCGTCACAGTGGGTTACGGTGACAGTGGTGTATGA =;1AABDBBFDDF:CBCGGFHHDEHEFHE):CGHADF=BBFDGIFDD<FG9 HI:i:1 NH:i:1 NM:i:0
## ERR458493.143432 16 IX 23779 40 51M * 0 0 GTTACAGCCATTGTCGTCACAGTGGGTTACGGTGACAGTGGTGTATGATTT FFBCFDDDFFFBFIFF<GCFE??FC8CA3EHHBBGE3IFAFFFDD8DD@<@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.996020 16 IX 23781 20 51M * 0 0 TACAGCCATTGTCGTCACAGTGGGTTACGGTGACAGTGGTGTATGATTTAG =BFB2DB<76@?:?*@G<IHHBGGGFG=IHFABEDBDE?<C>DBD=A=1+? HI:i:1 NH:i:1 NM:i:1
## ERR458493.364566 16 IX 23781 20 51M * 0 0 TACAGCCATTATCGTCACAGTGGGTTACGGTGACAGTGGTGTATGATTTTG GIJJJGIJIHCIGJIEGIIJJJJJJJIJJJJJIJIJJJHHHHHFFFFFC@@ HI:i:1 NH:i:1 NM:i:1
## ERR458493.463468 0 IX 23783 20 51M * 0 0 CTGCCATTGTCGTCACAGTGGGTTACGGTGACAGTGGTGTATGATTTTGGA =1+44ADD>DDDDEIEEIAECFCEEEEAEDDEIECDD?D8B@DDIBDDDIC HI:i:1 NH:i:1 NM:i:1
## ERR458493.795794 16 IX 23792 40 50M1S * 0 0 TCGTCACAGTGGGTTACGGTGACAGTGGTGTATGATTTTGGAGAAGCAGTG ?JJJGHGFIGIHDDDCEIGHHFGIGGHDIIIFEIIHDJFHHHHEDBB=?@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.466977 16 IX 23797 20 51M * 0 0 ACAGTGGGTTACGGTGACAGTGGTGTATGATTTTGGAGAAGCAGTTGTAGC HJJJJJJJJJIJJJIJIJJJJJJJJJJHIJIJJIJJJJHHHHHFFDB=+@@ HI:i:1 NH:i:1 NM:i:1
## ERR458493.752200 16 IX 23797 40 51M * 0 0 ACAGTGGGTTACGGTGACAGTGGTGTATGATTTTGGAGAAGCAGTTGTTGC CGHHDCGJJJJIIIJJHGIHHCJIJJJIGHEIIJJJIIHHFHHFFFFA@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.715681 16 IX 23797 40 51M * 0 0 ACAGTGGGTTACGGTGACAGTGGTGTATGATTTTGGAGAAGCAGTTGTTGC EJJJJJJJIFGEJJJJIJJJJJJJJJJIJHJIJJJJJJHHHHHFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.605693 16 IX 23799 40 51M * 0 0 AGTGGGTTACGGTGACAGTGGTGTATGATTTTGGAGAAGCAGTTGTTGCTG IJJIIGIJIIJIGJJJJJJJJIJJJJJJJJJJJJJIJJHHHHHFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.182821 0 IX 23806 40 51M * 0 0 TACGGTGACAGTGGTGTATGATTTTGGAGAAGCAGTTGTTGCTGATGTTTC ;@@DDBD1=FHD<E2A<CFFHHIIGHG@FBGHB@;??CHG;<???BG9?D4 HI:i:1 NH:i:1 NM:i:0
## ERR458493.244736 16 IX 23806 40 51M * 0 0 TACGGTGACAGTGGTGTATGATTTTGGAGAAGCAGTTGTTGCTGATGTTTC HHHGHIIHIIIIIIIIEIGIIIIIIIIIIIIIIGHIIIHHHHHFFDDD=;@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.261659 16 IX 23808 40 51M * 0 0 CGGTGACAGTGGTGTATGATTTTGGAGAAGCAGTTGTTGCTGATGTTTCTT IJIJJJJIJJIJJJJJJJIJIIJJJJJIJIJJJJJJJJHHHHHFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.35343 0 IX 23810 40 51M * 0 0 GTGACAGTGGTGTATGATTTTGGAGAAGCAGTTGTTGCTGATGTTTCTTTG @@@DFFFDFDHHFHJJIJJJJJJJJIIJJIIGHIGIJIJIHCHFHGIJJJI HI:i:1 NH:i:1 NM:i:0
## ERR458493.380035 0 IX 23816 40 51M * 0 0 GTGGTGTATGATTTTGGAGAAGCAGTTGTTGCTGATGTTTCTTTGGAACAT ?@7BBDFFAADFFIIGEHGEHGIDIIHE@FCEHICHHDGFHIIII@@EHIH HI:i:1 NH:i:1 NM:i:0
## ERR458493.514306 16 IX 23821 20 51M * 0 0 GTATGATTTTGGAGAAGCAGTTGTTGCTGATGTTTCTTTGGAACATTCCGA IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIGE9IIIFFFFFDDBA4+?= HI:i:1 NH:i:1 NM:i:1
## ERR458493.889740 16 IX 23828 40 51M * 0 0 TTTGGAGAAGCAGTTGTTGCTGATGTTTCTTTGGAACATTCAGAAGTGACA JJJJJJJJJJJJJIJJJJJIIJJJJIJJIJJJJJJJJJHHHHHFFFDB4<B HI:i:1 NH:i:1 NM:i:0
## ERR458493.1034697 16 IX 23832 20 51M * 0 0 GAGAAGCAGTTGTTGCTGATGTTTCTTTGGAACATTCAGAAGTGACAGTAT JJJJJIIGJIJJJJJJJJJJJJJIJJJIJJIHJJJJIJHHHHHFFFDB=+B HI:i:1 NH:i:1 NM:i:1
## ERR458493.468248 16 IX 23834 40 51M * 0 0 GAAGCAGTTGTTGCTGATGTTTCTTTGGAACATTCAGAAGTGACAGTCTTG IIIIIIIIIIIIIIIIIGCGIIIIIIIIIIGIIGBIIIHHBHFDDA44+?= HI:i:1 NH:i:1 NM:i:0
## ERR458493.431732 16 IX 23841 40 51M * 0 0 TTGTTGCTGATGTTTCTTTGGAACATTCAGAAGTGACAGTCTTGGTGTTAC BGCIIIIICGHGGIGFAGFGDHCHDIGHHEGDIGGHEHFAFHFBEDDD::? HI:i:1 NH:i:1 NM:i:0
## ERR458493.153413 16 IX 23852 40 51M * 0 0 GTTTCTTTGGAACATTCAGAAGTGACAGTCTTGGTGTTACAGCCATTATCG FGFEFEDBFFG?G>@EGEEEG?C<+AF<+9FFEFBIIIFFFD:BD?DD;@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.1017674 16 IX 23852 10 51M * 0 0 GTTTCTTGGGAACATTCAGAAGTGACAGTCTTGGTGTTACAGCCATTATCG ?;@DB0**CCD:?::>IEE:*19<DEFFC3:+3<2,<<222++4DBDD;?? HI:i:1 NH:i:1 NM:i:1
## ERR458493.120543 16 IX 23852 40 51M * 0 0 GTTTCTTTGGAACATTCAGAAGTGACAGTCTTGGTGTTACAGCCATTATCG EIHBGHEIGHIHGGHEFGBIGIHCHJIIGCGCGGGGHIHDDDBDDD=4+@< HI:i:1 NH:i:1 NM:i:0
## ERR458493.563115 0 IX 23856 40 51M * 0 0 CTTTGGAACATTCAGAAGTGACAGTCTTGGTGTTACAGCCATTATCGTCAC =1?==ADD?BBD?EEDEA33<A<A33<CFF<A??9:??DB3C:?DBDBBD; HI:i:1 NH:i:1 NM:i:0
## ERR458493.845052 16 IX 23868 40 51M * 0 0 CAGAAGTGACAGTCTTGGTGTTACAGCCATTATCGTCACAGTGAGTTACAG GGBGEIHGHCGHIIJJJIGHGFBIHIGGGIJIJJIHHFDAFGGFFFDFCC@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.409922 16 IX 23868 40 51M * 0 0 CAGAAGTGACAGTCTTGGTGTTACAGCCATTATCGTCACAGTGAGTTACAG HGEIGGIHHFGDAGHEDEGCIHEJIGCIHGGHFIHEHEHHHHFDFDDA@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.1015285 0 IX 23874 40 51M * 0 0 TGACAGTCTTGGTGTTACAGCCATTATCGTCACAGTGAGTTACAGTGGCAG ?@@DFFDDHHBHFHIJIJIIJEIHIG@<EGHIIIIHIAH@GGIGGIGEHID HI:i:1 NH:i:1 NM:i:0
## ERR458493.585070 16 IX 23874 40 51M * 0 0 TGACAGTCTTGGTGTTACAGCCATTATCGTCACAGTGAGTTACAGTGGCAG EJHJJJJJIGIGJIIIGIJIIJJJJJIJJJJIJJJJJJHHHHHFFEDB:1B HI:i:1 NH:i:1 NM:i:0
## ERR458493.959661 16 IX 23874 40 51M * 0 0 TGACAGTCTTGGTGTTACAGCCATTATCGTCACAGTGAGTTACAGTGGCAG JJIJIJJJIJJIIGFGGJJJJIJIJIIIJJJJJJJJJJHHHHHFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.103220 16 IX 23874 40 51M * 0 0 TGACAGTCTTGGTGTTACAGCCATTATCGTCACAGTGAGTTACAGTGGCAG 9?9GICGGIGDCGGC?18@EAGGGECIHGHFAEIIHECBDBFDFFFDD<@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.589217 16 IX 23879 40 51M * 0 0 GTCTTGGTGTTACAGCCATTATCGTCACAGTGAGTTACAGTGGCAGTTGTG HFAIIHGHDBDFGIIIGGIGGIGIHGHIIIIGGIHHFHHDDDDDFFFD;@? HI:i:1 NH:i:1 NM:i:0
## ERR458493.994519 0 IX 23882 40 51M * 0 0 TTGGTGTTACAGCCATTATCGTCACAGTGAGTTACAGTGGCAGTTGTGTAT @C1BDDDFGHHHHJJJJJJJJJJJJJJHHIJHIJJJIHGJIJJGHIHIFGH HI:i:1 NH:i:1 NM:i:0
## ERR458493.834276 0 IX 23887 40 51M * 0 0 GTTACAGCCATTATCGTCACAGTGAGTTACAGTGGCAGTTGTGTATGATTT ?@@BDDBDFHDFDGGH@EGEEH<ACFFHG@?@FCGEDHGHGEE3BGD@GE@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.103159 16 IX 23887 20 51M * 0 0 GTTACAGCCATTATCGTCACAGTGAGTTACAGTGGCAGTTGTGTATGAATT GGIGIEIIHIIIHIJGIJHHIIIJIJIIHJJJIJIIIIGGGHHFFEDA+@B HI:i:1 NH:i:1 NM:i:1
## ERR458493.543849 16 IX 23887 40 51M * 0 0 GTTACAGCCATTATCGTCACAGTGAGTTACAGTGGCAGTTGTGTATGATTT CGCHJIGIIJJJGGGHGHHIJIIHGJIIGJJIIGIJJJGHGHHFFFFFCC@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.686780 0 IX 23891 20 51M * 0 0 CAGTCATTATCGTCACAGTGAGTTACAGTGGCAGTTGTGTATGATTTGGAA @@@4=ADDFHHFFH;F?FFH9::A?EHGEEHI@FIGG???*:?@@FGEE9D HI:i:1 NH:i:1 NM:i:1
## ERR458493.90264 0 IX 23898 40 51M * 0 0 TATCGTCACAGTGAGTTACAGTGGCAGTTGTGTATGATTTGGAAGAAACGG ?;:BDDDDHFHHFEHIIIIII?FGGIG9C>?EFCHGIEHHEIIIGGIIIID HI:i:1 NH:i:1 NM:i:0
## ERR458493.597541 16 IX 23899 40 48M3S * 0 0 ATCGTCACAGTGAGTTACAGTGGCAGTTGTGTATGATTTGGAAGAAACAAC H?HGHF?BFHGEFDGEIHFCFBGHFAGBADHIIIHFBFADFHAADB=A11: HI:i:1 NH:i:1 NM:i:0
## ERR458493.894681 0 IX 23902 40 51M * 0 0 GTCACAGTGAGTTACAGTGGCAGTTGTGTATGATTTGGAAGAAACGGTTTC @@;DDFDDDBFFFGGHFCE?7FGCBH+<+<AGCHIIBG9FEHHGIIIHCHI HI:i:1 NH:i:1 NM:i:0
## ERR458493.810936 0 IX 23912 40 51M * 0 0 GTTACAGTGGCAGTTGTGTATGATTTGGAAGAAACGGTTTCAGTTACAACA ?=<DBBDBFBCA;F<C<AA<FH@BAECEGHI9CG;:?1?FC4??D4D9B## HI:i:1 NH:i:1 NM:i:0
## ERR458493.764986 0 IX 23922 40 51M * 0 0 CAGTTGTGTATGATTTGGAAGAAACGGTTTCAGTTACAACAGCACCGGTTG B@@ADDDFHHHHHIJJIJIJJJJJIJJFHHIGIGIGJJJJGIJJJJJJHII HI:i:1 NH:i:1 NM:i:0
## ERR458493.758968 0 IX 23930 40 51M * 0 0 TATGATTTGGAAGAAACGGTTTCAGTTACAACAGCACCGGTTGATGTACTA C@@DFFFFHHHHHIIJIJIFHIJIJDHIIJJJJJFHIJII?GIJJJIGHIJ HI:i:1 NH:i:1 NM:i:0
## ERR458493.202720 0 IX 23942 40 51M * 0 0 GAAACGGTTTCAGTTACAACAGCACCGGTTGATGTACTAACGTGACATCCA ???7DA?A;DFHBHEH@?@F?<BGFFBG:?DHIDFBFG@9D;?B?F<CFHG HI:i:1 NH:i:1 NM:i:0
## ERR458493.523854 0 IX 23942 40 51M * 0 0 GAAACGGTTTCAGTTACAACAGCACCGGTTGATGTACTAACGTGACATCCA ?@@=DD=:CCADHBGEDGHHAE@HGG<F::@B49BGGCBC<;FDGIGEEB3 HI:i:1 NH:i:1 NM:i:0
## ERR458493.180217 0 IX 23944 20 51M * 0 0 AACGGTTTCAGTTACAACAGCACCGGTTGATGTACTAACGGGACATCCAGT 8?@D;=D@CA<CF??A<:CBF+ABC????31?DDD9DBF@).076BFC7=@ HI:i:1 NH:i:1 NM:i:1
## ERR458493.325243 16 IX 23955 40 51M * 0 0 TTACAACAGCACCGGTTGATGTACTAACGTGACATCCAGTTTCATCACAAG 9D?DIGB??00:ADBGCGFDFDC?@HHBHIHFGAEBBEDC?HHD4DDB@@; HI:i:1 NH:i:1 NM:i:0
## ERR458493.811919 16 IX 23955 20 51M * 0 0 TTACAACAGCACCGGTTGATGTACTAACGGGACATCCAGTTTCATCACAAG >FDDFDFEGD:D??;E??AF?BC7GC<@3,FCA9FIIIFDDFC?AADD@@@ HI:i:1 NH:i:1 NM:i:1
## ERR458493.384109 0 IX 23958 40 51M * 0 0 CAACAGCACCGGTTGATGTACTAACGTGACATCCAGTTTCATCACAAGAAG B@?DFFFFHHHHFIIGIEHHHIJJJJJIIJJGIHJJEGIJIIIIIFHIJJI HI:i:1 NH:i:1 NM:i:0
## ERR458493.1048810 0 IX 23959 40 51M * 0 0 AACAGCACCGGTTGATGTACTAACGTGACATCCAGTTTCATCACAAGAAGT @@@DDDD;DFHCFHGH>C?@FEGHIE=F7DCFHFEFHE<?GEHIE<8DFHB HI:i:1 NH:i:1 NM:i:0
## ERR458493.772453 16 IX 23962 40 51M * 0 0 AGCACCGGTTGATGTACTAACGTGACATCCAGTTTCATCACAAGAAGTAAT HEAFEIHHHFD?CIHHCAHFGHEF@DHHFFEIIGIHDIHHDBCBFDDFCC@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.11395 0 IX 23968 20 51M * 0 0 GGGTGATGTACTAACGTGACATCCAGTTTCATCACAAGAAGTAATAGTGGC B@1ADDDFHHHHHJJJEEHIJIJJJICHIHIIJIIJJJIIC:DHIIIHHED HI:i:1 NH:i:1 NM:i:1
## ERR458493.607581 0 IX 23968 20 51M * 0 0 GGGTGATGTACTAACGTGACATCCAGTTTCATCACAAGAAGTAATAGTGGC B@+4ADDFHHHHHJJJHIJJJIJJJJHIJJIJIJJJJJHIG?FGIJJIIIG HI:i:1 NH:i:1 NM:i:1
## ERR458493.309929 0 IX 23968 40 51M * 0 0 GGTTGATGTACTAACGTGACATCCAGTTTCATCACAAGAAGTAATAGTGGC B@1ADDFFHHHGHJJIBHGGIJJJJJIGJIJJJJJJJJGGHCFHIJJGIJG HI:i:1 NH:i:1 NM:i:0
## ERR458493.757281 16 IX 23969 40 51M * 0 0 GTTGATGTACTAACGTGACATCCAGTTTCATCACAAGAAGTAATAGTGGCG JJJJJJJJIJJJJJJJJIIJJJJIJJJIIIJIIJJJJJHHHHHFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.853010 16 IX 23971 40 50M1S * 0 0 TGATGTACTAACGTGACATCCAGTTTCATCACAAGAAGTAATAGTGGCGGA IJJJJIHJJIGJJJJJJJJJJJJJJJJJJJIJJJJJJJHHHHHFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.191725 0 IX 23978 20 51M * 0 0 CTAACGTGACATCCAGTTTCATCACAAGAAGTAATAGTGGCGGTAGGGGAT @CCFDDDDHHGHHIIIGHIEIEIGIEHECCH?FHGII?FHIJG0?@(;A55 HI:i:1 NH:i:1 NM:i:1
## ERR458493.791411 0 IX 23983 40 51M * 0 0 GTGACATCCAGTTTCATCACAAGAAGTAATAGTGGCGGTAGTGGATGAGCA CB@FFFFFHHHHHJJJJJJJJJJJJJHIJJJJIJJJJJ88@9DF@FF>FCH HI:i:1 NH:i:1 NM:i:0
## ERR458493.641045 0 IX 23983 20 51M * 0 0 GAGACATCCAGTTTCATCACAAGAAGTAATAGTGGCGGTAGTGGATGAGCA ;1=;?D+BCBDFDDAFHA>BDC@;CBFECHFHDE@DFD607?/?<C##### HI:i:1 NH:i:1 NM:i:1
## ERR458493.429763 16 IX 23983 40 51M * 0 0 GTGACATCCAGTTTCATCACAAGAAGTAATAGTGGCGGTAGTGGATGAGCA F@FDG>IEBJIHDICGHFCGGF@>GIGGHGGIJJJIIIDHHHHFFDBA@C@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.643682 16 IX 23983 20 51M * 0 0 GTGACATCCAGTTTCATCACAAGAAGTAATAGTGGCGGTAGTGGATGCGCA =<HFB8B@EEGIGGIIIIGHBC:GC?@EFIGIIGAEAAFFAFF@DDDD=<? HI:i:1 NH:i:1 NM:i:1
## ERR458493.622069 0 IX 23985 40 51M * 0 0 GACATCCAGTTTCATCACAAGAAGTAATAGTGGCGGTAGTGGATGAGCATG C<@BFFFFHGHHHJJJJJJJIHIJFHIIJJIIIJJJ:8D9BF<BFBF>DHF HI:i:1 NH:i:1 NM:i:0
## ERR458493.1081963 0 IX 23989 40 51M * 0 0 TCCAGTTTCATCACAAGAAGTAATAGTGGCGGTAGTGGATGAGCATGGAAC ???4ADDDDDDBCAFE;:<C,<?BH?3ABGHGG6?99B>@?########## HI:i:1 NH:i:1 NM:i:0
## ERR458493.153572 0 IX 23993 40 51M * 0 0 GTTTCATCACAAGAAGTAATAGTGGCGGTAGTGGATGAGCATGGAACGGTT CC@FFFFFHHHHHGIIEGHJIJEIJIJJ:?D?DG?<D?DD@DFHEHGGH;; HI:i:1 NH:i:1 NM:i:0
## ERR458493.857364 16 IX 23994 40 51M * 0 0 TTTCATCACAAGAAGTAATAGTGGCGGTAGTGGATGAGCATGGAACGGTTG ###B9*8FAE@EHCDDD9>>HDC6@F@F<IIHFE<<HF>DFD:BBDDD@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.600511 0 IX 23996 40 3S48M * 0 0 GGGTCATCACAAGAAGTAATAGTGGCGGTAGTGGATGAGCATGGAACGGTT =1+4ADDDFHDHFGIIAFHHEHCHHIEE)0?0?D>DE<FC9D@G>G>AG@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.811293 16 IX 24023 40 51M * 0 0 GTGGATGAGCATGGAACGGTTGTGGTAATGGTATAGACGTTACCATTGGAG JJJJJJJJJJJJJJJIJJJJJJJJJJHFA,JJJJJIIJHHHHHFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.23975 0 IX 24035 40 51M * 0 0 GGAACGGTTGTGGTAATGGTATAGACGTTACCATTGGAGTCAGTTGTCTCT :=?D;D?<:@D1C+A::C@C<CC9CEEIADDC9?:9*???DD@DBDE49?@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.87334 0 IX 24035 40 3S48M * 0 0 CTGGGAACGGTTGTGGTAATGGTATAGACGTTACCATTGGAGTCAGTTGTC ?+1=BD>?FF?FHIB?<CFFHE?FFBE><@?FE@<?FB9*9?GHGIB??## HI:i:1 NH:i:1 NM:i:0
## ERR458493.427691 16 IX 24037 40 51M * 0 0 AACGGTTGTGGTAATGGTATAGACGTTACCATTGGAGTCAGTTGTCTCTGA JIJJJJJJJJJJJJJJIJJJJJIJJJIHJJJJJJJJJJHHHHHFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.811877 16 IX 24042 13 51M * 0 0 TTGTGGTAATGGTATAGACGTTACCATTGCAGTCAGTTGTCTCTGAACAAG GJJJJJJJJJJJJJJJIGIIJJIJHHHC<+JJJJJJJJHHHHHFFFFFBCC HI:i:1 NH:i:1 NM:i:1
## ERR458493.395837 0 IX 24054 20 51M * 0 0 TTTAGACGTTACCATTGGAGTCAGTTGTCTCTGAACAAGAAACTACGACTG =+:+=BB1ADFFF@EGBF@><CAFCF+2AC?FE>GCFGFFGII@GFF<GFF HI:i:1 NH:i:1 NM:i:1
## ERR458493.728734 16 IX 24056 40 51M * 0 0 TAGACGTTACCATTGGAGTCAGTTGTCTCTGAACAAGAAACTACGACTGTG EHHFDBB:*EDEGBJJIIHGIIJJIFJIHJIJIIJIHJHHHHHFFFFFCCB HI:i:1 NH:i:1 NM:i:0
## ERR458493.429136 16 IX 24057 40 51M * 0 0 AGACGTTACCATTGGAGTCAGTTGTCTCTGAACAAGAAACTACGACTGTGC 3D@IIIGGGGHGCHEDGBECBIGFHGHF<EGGIIGEFIHHHFHDADD;@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.99357 16 IX 24060 20 51M * 0 0 CGTTCCCATTGGAGTCAGTTGTCTCTGAACAAGAAACTACGACTGTGCTTG @B9B2IHGDHHCECD@>IIHEIGHGIGHABDGHGA;FDFFCADDDDA1@?= HI:i:1 NH:i:1 NM:i:1
## ERR458493.129315 0 IX 24067 20 51M * 0 0 ATCGGAGTCAGTTGTCTCTGAACAAGAAACTACGACTGTGCTTGTGGCATT =?1+4ADDDFHHAFEFGEGIGBCHHHGGGGIDHEBFG?BGIIGGFGGGGGI HI:i:1 NH:i:1 NM:i:1
## ERR458493.252230 16 IX 24068 40 51M * 0 0 TTGGAGTCAGTTGTCTCTGAACAAGAAACTACGACTGTGCTTGTGGCATTT IJJJJJIJJJIIIGIGIIIJIJJJIIIGIIIJJIIJIJGHHHHFFFDA<@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.356552 0 IX 24089 40 51M * 0 0 CAAGAAACTACGACTGTGCTTGTGGCATTTACATAGTCTAGTGTTGTAGTG @<@FFFFFFHBH>FHGIDHGHGGGGGIJGFGIGIII<GHEF?DGGG?BFHG HI:i:1 NH:i:1 NM:i:0
## ERR458493.173448 16 IX 24109 40 50M1S * 0 0 TGTGGCATTTACATAGTCTAGTGTTGTAGTGTAAACAGATGTTGAACCGGG GIIHF<F?BDDHGGGFEE>EC?@HCF;C<?:F?CCGF>FDHFBDA@?7@?= HI:i:1 NH:i:1 NM:i:0
## ERR458493.1011576 16 IX 24109 40 50M1S * 0 0 TGTGGCATTTACATAGTCTAGTGTTGTAGTGTAAACAGATGTTGAACCGGA JJIJIIIJIJIJIHGEIJGHJJJJIJJIJIJJJJJJIJHHHHGFFDFFCC@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.576333 16 IX 24109 40 50M1S * 0 0 TGTGGCATTTACATAGTCTAGTGTTGTAGTGTAAACAGATGTTGAACCGGA JJJJJJJJJJJJJJJJIJIJJIJJJJJJJJJJJJJJJJHHHHHFFDDD=+B HI:i:1 NH:i:1 NM:i:0
## ERR458493.787792 0 IX 24112 20 51M * 0 0 GTCATTTACATAGTCTAGTGTTGTAGTGTAAACAGATGTTGAACCGGTAAC @:=DDFFFHHHHGHIJIJHIIJJIJJJIIJJJJJJIJJJGIJJJJJJGIGG HI:i:1 NH:i:1 NM:i:1
## ERR458493.639898 16 IX 24122 20 51M * 0 0 TAGTCTAGTGTTGTAGTGTAAACAGATGTTGAACCGGTAACGGAAGCAACA IGHIGIIIIIIGIIIIIIIHGIIIIGHH@GED@IGGGFFFFDHEDB=41?? HI:i:1 NH:i:1 NM:i:1
## ERR458493.225857 16 IX 24134 20 50M1S * 0 0 GTAGTGTAAACAGATGTTGAACCGGTAACGGAAGCACCAAATGAAGTGAAG IJIGIHEHG@HGEJJJJHGDEIJJIIHGIGGIIGGHJJHGFGHFEDD=+BB HI:i:1 NH:i:1 NM:i:1
## ERR458493.755592 16 IX 24149 40 51M * 0 0 GTTGAACCGGTAACGGAAGCACCAAATGAAGTGGAAGCAGTTGAACCGGTA CJJHF?JJIH?HGIIIHF<F@GHHEIGIDGHCJJJIIIHHHHHFDFFFCC@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.400063 0 IX 24158 20 51M * 0 0 GAAACGGAAGCACCAAATGAAGTGGAAGCAGTTGAACCGGTAACGGAAGCA 811B?DDD?FAAFEC?GFACDF9CFFHG;EGDG?9D9BGH0?D87;C6@## HI:i:1 NH:i:1 NM:i:1
## ERR458493.559958 16 IX 24173 10 51M * 0 0 AATGAAGTCGAAGCAGTTGAACCGGTAACGGAAGCACCGAATGAAGTGGGA IJIIIGF?0IJJIIGGIGGJIJIJJIIGJJIIJJJIIIHHHHHFFDDA1+? HI:i:1 NH:i:1 NM:i:2
## ERR458493.593573 0 IX 24181 40 50M1S * 0 0 GGAAGCAGTTGAACCGGTAACGGAAGCACCGAATGAAGTGGAAGCAGTTGG @@@FFFFFDFHBH@FADGEHIIIIIGEIII@@FHIIIIGIIIIGHFHFGA# HI:i:1 NH:i:1 NM:i:0
## ERR458493.115927 0 IX 24182 40 51M * 0 0 GAAGCAGTTGAACCGGTAACGGAAGCACCGAATGAAGTGGAAGCAGTTGAA ;??D?DDBDF<B<CGG8CEBA:G=:;?GHGBFHCCC?FHG>FH2BE##### HI:i:1 NH:i:1 NM:i:0
## ERR458493.608840 16 IX 24183 40 51M * 0 0 AAGCAGTTGAACCGGTAACGGAAGCACCGAATGAAGTGGAAGCAGTTGAAC GIIIIIEF@D@IIIIHGEEEIHHHFCGGIFIIIGIIIIHFFDDDDA4+11? HI:i:1 NH:i:1 NM:i:0
## ERR458493.433869 0 IX 24187 20 4S47M * 0 0 ATGAAGTAGAACCGGTAACGGAAGCACCGAATGAAGTGGAAGCAGTTGAAC BCCFFFEEHHHHHJJCGGGIIIIJJIJJJGIGJEIIDFIHIJGIJIJJIII HI:i:1 NH:i:1 NM:i:1
## ERR458493.16479 0 IX 24197 40 51M * 0 0 GTAACGGAAGCACCGAATGAAGTGGAAGCAGTTGAACCGGTAACGGAAGCA @@@FFDDBHFHBHHIEHIGGHICEBAFABG@GH9<@BFFA6AHAGGG@BAE HI:i:1 NH:i:1 NM:i:0
## ERR458493.515745 0 IX 24206 20 51M * 0 0 GCACCGAATGAAGTGGAAGCAGTTGAACCGGTAACGGAAGCACCGAATGAA @@@FFFFFHHHHFEGGIJGIJJEHGGHGIGH@@DHGIGIIIIJIJIJJJII HI:i:1 NH:i:1 NM:i:0
## ERR458493.1043551 0 IX 24207 20 51M * 0 0 CACCGAATGAAGTGGAAGCAGTTGAACCGGTAACGGAAGCACCGAATGAAG @@@DDD<DHHFFHGIEHHDIHHH@HHDCHAGIGIGBGGGGHGIGBHIIIEH HI:i:1 NH:i:1 NM:i:0
## ERR458493.254177 16 IX 24215 20 51M * 0 0 GAAGTGGAAGCAGTTGAACCGGTAACGGAAGCACCGAATGAAGTAGAACCG IIFB?CGFFEEFDD9FBGDEABCGCF<A<3)A+BGC:BFDFFFADDBA=11 HI:i:1 NH:i:1 NM:i:0
## ERR458493.6718 0 IX 24227 40 51M * 0 0 GTTGAACCGGTAACGGAAGCACCGAATGAAGTAGAACCGGTAACGGAAGCG @@;DDEFFHHDFFEHGFIIDAGHIHH=FGIIHFCFG<FDE;@FCGGGHDFD HI:i:1 NH:i:1 NM:i:0
## ERR458493.71907 0 IX 24227 10 31M20S * 0 0 GTTGAACCGGTAACGGAAGCACCGACTGAAGGGGAAGCAGTTGAACCGGTA ???;BDDADA0<C9CEEE@@?EEIA:)8B9?*.?ADD############## HI:i:1 NH:i:1 NM:i:1
## ERR458493.553608 0 IX 24227 40 51M * 0 0 GTTGAACCGGTAACGGAAGCACCGAATGAAGTAGAACCGGTAACGGAAGCG ?@@AAAADHHHAFCHJIJJDHGGGIIGEEEHBFABDABEFAFHGGI@?EHF HI:i:1 NH:i:1 NM:i:0
## ERR458493.568274 0 IX 24237 20 51M * 0 0 TAATGGAAGCACCGAATGAAGTAGAACCGGTAACGGAAGCGGAAGCACCAG @@@:ADEFHGHHHJHIJJJJJEFHIIJJJJFHHIJJGIIGIJDGGIJGHFF HI:i:1 NH:i:1 NM:i:1
## ERR458493.61990 16 IX 24246 40 51M * 0 0 CACCGAATGAAGTAGAACCGGTAACGGAAGCGGAAGCACCAGATGCAGTGG GIIGIIIGFGFB@GGHEGDHGGDFAHBIIIIIIHBGCADHHHHFFFFFC@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.515533 0 IX 24247 40 51M * 0 0 ACCGAATGAAGTAGAACCGGTAACGGAAGCGGAAGCACCAGATGCAGTGGA @@@FDDFFDDFHFGCHJGHGEGHGHHIIIGG@GHIIFHHG?FHFHA=@DAC HI:i:1 NH:i:1 NM:i:0
## ERR458493.748767 0 IX 24247 20 51M * 0 0 ACTGAATGAAGTAGAACCGGTAACGGAAGCGGAAGCACCAGATGCAGTGGA @@+4BDBAFHHFFHGGGIGG@FHHGGGHAHIGIGHGHIGCEHEGIGICGH@ HI:i:1 NH:i:1 NM:i:1
## ERR458493.185392 16 IX 24251 40 51M * 0 0 AATGAAGTAGAACCGGTAACGGAAGCGGAAGCACCAGATGCAGTGGATGAA BBHAIGIHDB@:GBDBB?HEDHC86<3FHF9GFGG>GF??8BHDDD?D?=8 HI:i:1 NH:i:1 NM:i:0
## ERR458493.352994 16 IX 24271 20 12S39M * 0 0 GTAGTTTAGGTAGTAAGCGGAAGCACCAGATGCAGTGGATGAAGAAGCGGC JJIGJJJIJIIHGIHGGIIIIGHIJJJJJJJIIJJJJJHGHHHFFFFFCCC HI:i:1 NH:i:1 NM:i:1
## ERR458493.541216 0 IX 24274 40 51M * 0 0 AGCGGAAGCACCAGATGCAGTGGATGAAGAAGCGGCACTAGTAGCGACAAC CCCFFFFFHHHHHIJJJJJGFGH*:EFGHIIJIJJIBHIGGGGIJJJJHHG HI:i:1 NH:i:1 NM:i:0
## ERR458493.429107 0 IX 24299 20 51M * 0 0 GTAGAAGCGGCACTAGTAGCGACAACACCATTGTAGTTTAGGTAGTAAGCG @:<D;B:?FF<:CF@CCEHFEFBFG9F6?8B@9B4/0=C)=<=BF8@==AE HI:i:1 NH:i:1 NM:i:1
## ERR458493.867298 16 IX 24332 40 49M2S * 0 0 TAGTTTAGGTAGTAAGCGTAAGCACCTGCAAAGATTCCTTCAGAGACGCGC IJJJJJJIIIJJJIIEJJJJJIGEGGGIGJJJIIFIJJHHGHHFFDFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.420671 16 IX 24340 20 51M * 0 0 GTAGTAAGCGTAAGCACCTGCAAAGATTCCTTCAGAGACGCTGAAGTCTGA IIHEIIGIGFIIGIHFIHEGHGIIHFF>GGIIIIIGGFFHFFHEDDA4+@= HI:i:1 NH:i:1 NM:i:1
## ERR458493.535973 16 IX 24403 40 51M * 0 0 GGTACCAATAGCAAAAGAGAATTGTCTGTTACCTTGGGTTGCAGTTAGAAC DIHGHFGIIGIIIGIIIIGIIIIIIGIHIHCEIIHHDIHDFFFBDBBA+:? HI:i:1 NH:i:1 NM:i:0
## ERR458493.474127 16 IX 24407 20 51M * 0 0 CCAATAGCAAAAGAGAATTGTCTGTTACCTTGGGTTGCAGTTAGAACGCCG DBDDDD@B<D?>DB9DA@IDIEECEEC1CIEEEEIIEED>DDC:<0;A+;? HI:i:1 NH:i:1 NM:i:1
## ERR458493.728971 0 IX 24423 40 51M * 0 0 ATTGTCTGTTACCTTGGGTTGCAGTTAGAACACCGGTAGTACCGTCATAAG CCCFFFFFHHHHHJJJJJGHIJJJIJJIJJJJJJJI@GIGGIJJIJJJJJI HI:i:1 NH:i:1 NM:i:0
## ERR458493.208313 0 IX 24446 40 51M * 0 0 GTTAGAACACCGGTAGTACCGTCATAAGAAAGAGCAGTGATGTCAGCGGTA @CCDFFFFHHHHHAGHFGIJJDHIJJIJJIIJGIJIG@GHIIIIIIJJIAA HI:i:1 NH:i:1 NM:i:0
## ERR458493.1037850 0 IX 24446 40 51M * 0 0 GTTAGAACACCGGTAGTACCGTCATAAGAAAGAGCAGTGATGTCAGCGGTA CCCFFFFFHHHHHFHJIGIJIGIJJJJJJJJJHIJJJHIJJJJJJJJJIFH HI:i:1 NH:i:1 NM:i:0
## ERR458493.942872 0 IX 24449 40 51M * 0 0 AGAACACCGGTAGTACCGTCATAAGAAAGAGCAGTGATGTCAGCGGTAATG @@@FFFFFHHHFHFHIJJGIIJJJJJJJJIIFGIGHGIJGGJJIJJAFGGI HI:i:1 NH:i:1 NM:i:0
## ERR458493.794908 0 IX 24455 40 51M * 0 0 CCGGTAGTACCGTCATAAGAAAGAGCAGTGATGTCAGCGGTAATGGCAATT @@<DDDDDDDHDFGGI>GGIII<E93BD?EEEEHHIGIID@DHHGB?D=;D HI:i:1 NH:i:1 NM:i:0
## ERR458493.243029 16 IX 24467 20 51M * 0 0 TCATAAGAAAGAGCAGTGATGTCAGCGGTAATGGCAATTTGGTTCTTACCA HDCIIIHCD@CHGHEHEDIHHIBIIHGIIIGIGHHEFEHABHDDFFFF@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.514807 16 IX 24470 20 51M * 0 0 TAAGAAAGAGCAGTGATGTCAGCGGTAATGGCAATTTGGTTCTTACCAGTG EIIIIIIIIFIIIIIFIIIIIIGIIH>IIIIIIIIIIIFFFFFDDDDD1@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.339061 16 IX 24497 40 51M * 0 0 ATGGCAATTTGGTTCTTACCAGTGTAGCCAACAACAGGAATAGGAGTGGCA JIIIHIIGIIJJIIGHGCIGIJJJJIIIIGGIIGIJIJHHHHHFFFFFC@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.1041852 16 IX 24507 40 51M * 0 0 GGTTCTTACCAGTGTAGCCAACAACAGGAATAGGAGTGGCATTGGTTTCGG 2F<@<9D??9HIGF@<C?FGFGFEHE>B<4GGG@GHEBFHFDBDDD?A<@? HI:i:1 NH:i:1 NM:i:0
## ERR458493.659691 16 IX 24517 40 51M * 0 0 AGTGTAGCCAACAACAGGAATAGGAGTGGCATTGGTTTCGGTTGGGTCAAC GIIJHCGCCHBHGDGIIIJHGDEGGHHGGBJIIIJJJIHGHGHFFFFF@C@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.518919 0 IX 24531 40 51M * 0 0 CAGGAATAGGAGTGGCATTGGTTTCGGTTGGGTCAACAGCAAGGACACCTT @@@DDDDDHHHHHHGGD+A<:<AFFG6:EGGH;GHIIIIIIIEFC(=CCHC HI:i:1 NH:i:1 NM:i:0
## ERR458493.622413 16 IX 24533 20 50M1S * 0 0 GGAATAGGAGTGGCATTGGTTTCGGTTGGGTCAACAGCAAGGACACCTACT IIJIJJJJJIGJJJGJJIIIIEJJJIIIJIIJJJIJJJHHHHHDBDA41B= HI:i:1 NH:i:1 NM:i:1
## ERR458493.992877 0 IX 24542 40 51M * 0 0 GTGGCATTGGTTTCGGTTGGGTCAACAGCAAGGACACCTTCACCCTTGAAT B:?ADFFFHHHHHIJJJJJJJCGHIJJJJJJJJJJJJJJJJIJJJJJJJJJ HI:i:1 NH:i:1 NM:i:0
## ERR458493.435714 16 IX 24547 40 51M * 0 0 ATTGGTTTCGGTTGGGTCAACAGCAAGGACACCTTCACCCTTGAATACAAC GCJJJGIIIJJJJGJIIHHHGJJIJIGFEEFIHGHGAIHHHHHFFFDF@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.873183 0 IX 24552 40 51M * 0 0 TTTCGGTTGGGTCAACAGCAAGGACACCTTCACCCTTGAATACAACAGTTT @CC?DFFFHGGADGHIIEHIIFIIIIIIIIGGGHIIIGHHIIIIGIIICHH HI:i:1 NH:i:1 NM:i:0
## ERR458493.592826 16 IX 24574 40 1S49M1S * 0 0 CGACACCTTCACCCTTGAATACAACAGTTTGGCCAGTAAATGTGTCAGGGA #FCGHGGHF@CIIGIIGEIHFHHFIIIIIIIIIHDEGGHFHDHDDDB+=<@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.410667 0 IX 24616 40 51M * 0 0 GTCAGGGTAGTGTAAGTATAGGTTACCGGAAATGATGTTGACTGTACCTTT ?:1ADAD=DD<C2AEE<CEFEIAEFEEIA@:CCEIII?BBBDD>9DDC>DD HI:i:1 NH:i:1 NM:i:0
## ERR458493.459413 0 IX 24631 40 51M * 0 0 GTATAGGTTACCGGAAATGATGTTGACTGTACCTTTTCCAGAAACTGGTTC @B+BDFFDHHHHHJJJJJJJIJJJIJJJJIIJJJ9FHIJIJIJJJJJJIIJ HI:i:1 NH:i:1 NM:i:0
## ERR458493.699466 16 IX 24634 40 51M * 0 0 TAGGTTACCGGAAATGATGTTGACTGTACCTTTTCCAGAAACTGGTTCGAC EHDIIIHEGIIIIIIGGGIHDF@GHHHFIGIJJIIGGGFHFHFDDFFF@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.856435 16 IX 24644 40 51M * 0 0 GAAATGATGTTGACTGTACCTTTTCCAGAAACTGGTTCGACAATAACATAG IIIGIIIHHIGHFGHFCCJIIJIIIGHEHGIIJIJJJIHHHGHFFFFF@C@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.478349 16 IX 24662 40 51M * 0 0 CCTTTTCCAGAAACTGGTTCGACAATAACATAGGTACTTCCATTGTCTAGG DIHFGECGHFD@@JIIIIIHHEHCGGJHFIGIIHDGIIHFHDBEFDDD@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.851551 16 IX 24662 40 51M * 0 0 CCTTTTCCAGAAACTGGTTCGACAATAACATAGGTACTTCCATTGTCTAGG HEGGD?3@AFDFDDC8BF=FD@HCHHC<?IGCGEDCAEHFHFHDDDDD?@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.789314 0 IX 24674 40 51M * 0 0 ACTGGTTCGACAATAACATAGGTACTTCCATTGTCTAGGTTAATTTCACCG @8=DDADDFFFFFBGIEGHFIF<AEFFFFEEIEEFDFGGCGGIIIIEGGGD HI:i:1 NH:i:1 NM:i:0
## ERR458493.385133 16 IX 24677 40 51M * 0 0 GGTTCGACAATAACATAGGTACTTCCATTGTCTAGGTTAATTTCACCGTTG IIIJIIIBJJIGGJJIIHGFFAD?AGGHFHGHHEHGJGHHDDFB:DFF@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.111125 16 IX 24730 40 51M * 0 0 TACTGAACCTTCAGTGTCGTCTCTACGTTGTAGACCGGAAACAGAACCACC HGGGGGFIDGGGHGFAIIIIHGCGHEFEHHGEAF8GGGC?AFFDDAA?:+? HI:i:1 NH:i:1 NM:i:0
## ERR458493.218886 0 IX 24834 40 51M * 0 0 AAGTGACTTCACCCTTCTTGGACTTTGATAGACTTAAGGAAATGTCACCAC ??;DBBA;DDDD8EEE@ED:?AE3+C4CFE@EC9:*1:9D9??D?4DDEC; HI:i:1 NH:i:1 NM:i:0
## ERR458493.316832 0 IX 24856 40 51M * 0 0 CTTTGATAGACTTAAGGAAATGTCACCACTGTTGTCAAAAGAACTTGGAGT @@@DFFFDGHFFHGIIGIJJGGFHGGGGJJJGIGHHEGGGIGGHCEHIBC9 HI:i:1 NH:i:1 NM:i:0
## ERR458493.937519 16 IX 24872 40 51M * 0 0 GAAATGTCACCACTGTTGTCAAAAGAACTTGGAGTGAAGGAGTAGATGGAT JJJJJJJIHJJIJJJJJJJJJJJJHHHJJJJJIIIJJJHHHHHFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.294190 16 IX 24878 40 51M * 0 0 TCACCACTGTTGTCAAAAGAACTTGGAGTGAAGGAGTAGATGGATGCAGAG HFFJJIIIJJJJJJIIHGGIHJJJIJJJIIJJIIGJIJHHHHGFFFDB:1B HI:i:1 NH:i:1 NM:i:0
## ERR458493.468964 0 IX 24879 20 51M * 0 0 CACCACTGTTGTCAAAAGAACTTGGAGCGAAGGAGTAGATGGATGCAGAGG =;?DD?D3A,2CADEEEABEE4<B>+A+81?DDBDBD<<<DE######### HI:i:1 NH:i:1 NM:i:1
## ERR458493.676636 0 IX 24924 40 51M * 0 0 CAGAGGTGGCAGCGGATTCTTCAGCGTTGAAGTTACCGGTCACGTCAAAGT ?@<D1A=2<BB;<C<@@?3?CEFEC::B::*99BB<DFDFB;7;@###### HI:i:1 NH:i:1 NM:i:0
## ERR458493.748403 0 IX 24926 40 51M * 0 0 GAGGTGGCAGCGGATTCTTCAGCGTTGAAGTTACCGGTCACGTCAAAGTTT CCCFDFFFHHHHHJJJJJJIJJIJGIIJJJGIJJJJJHIIJJHIJJJJEHI HI:i:1 NH:i:1 NM:i:0
## ERR458493.1052357 0 IX 24935 40 51M * 0 0 GCGGATTCTTCAGCGTTGAAGTTACCGGTCACGTCAAAGTTTTCACCGGAG @@@DDD8ADFFAFGEI@83AFHIFFIIIE:@FFIIIIIIFIIIIBFF?E<; HI:i:1 NH:i:1 NM:i:0
## ERR458493.263624 0 IX 24981 40 51M * 0 0 CGGAGATGTCGAATTCGGTACCTGGATAGGTTTCTTTGGACTTTTTCACAA BBCDFFFFHHHHHJJJJJHHIJJJJFHIJJCFHHJJJIIHHIJJJJJIIJJ HI:i:1 NH:i:1 NM:i:0
## ERR458493.113550 0 IX 24988 40 51M * 0 0 GTCGAATTCGGTACCTGGATAGGTTTCTTTGGACTTTTTCACAAGCAAGTC @@<DABDDH?DHHIIDGG9?CFECHHIIIIEGBGIGIIIFGIIIDHIHIFB HI:i:1 NH:i:1 NM:i:0
## ERR458493.37813 0 IX 24995 40 51M * 0 0 TCGGTACCTGGATAGGTTTCTTTGGACTTTTTCACAAGCAAGTCACCGTAG @@@FDFAFHHHAFFGGAFIDHHI<GHIIIIHHIIIIGEHIGIIGGHHGDFD HI:i:1 NH:i:1 NM:i:0
## ERR458493.914000 16 IX 25014 40 51M * 0 0 CTTTGGACTTTTTCACAAGCAAGTCACCGTAGACTTTTAGGTCAGCTGAGA BG>F?*FHGGEGCEGIGHGFEBEEGEA?CEDECDGHBGBFHHHDAEFFCC@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.276156 16 IX 25027 40 51M * 0 0 CACAAGCAAGTCACCGTAGACTTTTAGGTCAGCTGAGAAAGAATATTTGTC ##DB?4>?4DDDD:B?@BC9F@E@EA3:A+C3:A4EEA<2?C<A48BD??8 HI:i:1 NH:i:1 NM:i:0
## ERR458493.670304 16 IX 25077 20 51M * 0 0 CACCGTCAAGTAGAGTCAATTTGGCACCAGCTGCGACCACAACACCTACTG F@IIIIIIIIIIIIIIIIIIIIIIIHIIIIIIIIHFIIHHHHFBDDA4+@@ HI:i:1 NH:i:1 NM:i:1
## ERR458493.661722 16 IX 25081 40 51M * 0 0 GTCAAGTAGAGTCAATTTGGCACCAGCTGCGACCACAACACCTTCTGTTGC <GIIIIIHF<GGEIIHGIGEIIGGIGGGGIHEIGHEC<<<DDADFDFD<@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.799460 0 IX 25084 40 51M * 0 0 AAGTAGAGTCAATTTGGCACCAGCTGCGACCACAACACCTTCTGTTGCAAT B@CFFFFFFHFHHJJJJIGIJJJIJJJJGHJIJGJJJJJJJJJJJJIJJII HI:i:1 NH:i:1 NM:i:0
## ERR458493.514498 0 IX 25096 20 51M * 0 0 TTTGGCACCAGCTGCGACCACAAGACCTTCTGTTGCAATTACTTTACCAGT :=+4ADDDHFHFHGGGGIIIIIG)?CGGBFHHIGIIGGIIF<DHHIIEFG4 HI:i:1 NH:i:1 NM:i:1
## ERR458493.511872 16 IX 25134 40 51M * 0 0 TTACTTTACCAGTAAGAGTTAGATCACCACTGATTGTAGTAGTTGGTAGAT HHHIIIIIIIIIIIIIIIGIHHIHGHIHHIIIIIIIGIGHHHFADDA=:1? HI:i:1 NH:i:1 NM:i:0
## ERR458493.962082 0 IX 25183 40 51M * 0 0 ATAGATAGTACCATCGCTGGAAGAAAGGGAGTTGGAAGCAGTTGCGCTTGC @B@FFFFFCFHHHIJJJJJJEGIIJIJJJFGGHHIEHIHJGHIJJJJJJGI HI:i:1 NH:i:1 NM:i:0
## ERR458493.464504 16 IX 25197 20 51M * 0 0 CGCTGGAAGAAAGGGAGTTGGAAGCAGTTGCGCTTGCACTAGATGCTGGAG GDABIIIG@IIGGEDFGDGHGGGIGGHEGEIGIG@GGCFBHFBDDD=4+=? HI:i:1 NH:i:1 NM:i:1
## ERR458493.44290 0 IX 25199 40 51M * 0 0 CTGGAAGAAAGGGAGTTGGAAGCAGTTGCGCTTGCACTAGATGCTGAAGAA @@@FFFF?AFBBAEFHGHIHIIIEGHIIIDEBGBDD??DHGCEDEH>BBHG HI:i:1 NH:i:1 NM:i:0
## ERR458493.483014 0 IX 25209 40 51M * 0 0 GGGAGTTGGAAGCAGTTGCGCTTGCACTAGATGCTGAAGAAGTAATGGAAC @@@BDDDDHHBFBBFEFEDFEGGGGGIIIIGBGEHGGHIGII9BFHIGHIC HI:i:1 NH:i:1 NM:i:0
## ERR458493.858474 0 IX 25213 40 51M * 0 0 GTTGGAAGCAGTTGCGCTTGCACTAGATGCTGAAGAAGTAATGGAACTCAA CCCFFFFFHHHHHJJJJIJJJIJJJJJJJJJJJJJJJJGHIEIIIJJJJJI HI:i:1 NH:i:1 NM:i:0
## ERR458493.1050228 0 IX 25214 13 51M * 0 0 TTGGAAGAAGTTGCGCTTGCACTAGATGATGAAGAAGTAATGGAACTCAAG :<?D+2A22AA<DE>1@EE8?EDFIC?>*1:**?DEDBBB44B3B*09?## HI:i:1 NH:i:1 NM:i:2
## ERR458493.399420 16 IX 25215 40 51M * 0 0 TGGAAGCAGTTGCGCTTGCACTAGATGCTGAAGAAGTAATGGAACTCAAGG IJJIJJIJJIJJIJJJIIIIJJJJIJJJJJJJJJJJJJHGHHHDDB=4+@? HI:i:1 NH:i:1 NM:i:0
## ERR458493.290736 16 IX 25235 40 51M * 0 0 CTAGATGCTGAAGAAGTAATGGAACTCAAGGAATTGGAAGTGGCAGAACCG FIIIIFEFFFIIIIIIIIIIIHIHIFIIIIIIEIIIIIFFFFFDDB4+;1= HI:i:1 NH:i:1 NM:i:0
## ERR458493.649436 0 IX 25244 40 51M * 0 0 GAAGAAGTAATGGAACTCAAGGAATTGGAAGTGGCAGAACCGGAGGCAGTA B@@FFFFDHHHHHJJJJJJIJJGJJJJJJJJGIIIJJIJJJIIHIJIIJFG HI:i:1 NH:i:1 NM:i:0
## ERR458493.27160 16 IX 25267 40 51M * 0 0 ATTGGAAGTGGCAGAACCGGAGGCAGTAGAAGCGGATGAGGTTGATTGAGG CFC@GIIGEHG>HDD:IIIGGIIIHFH@IGGDHGDEIGGHFFFFFDFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.767969 16 IX 25267 40 51M * 0 0 ATTGGAAGTGGCAGAACCGGAGGCAGTAGAAGCGGATGAGGTTGATTGAGG HF>JJIJHHHGIJIGGIEJJIIIIIIFJJJJJJJJJGJHHHHHFFDBA=@? HI:i:1 NH:i:1 NM:i:0
## ERR458493.18419 0 IX 25276 40 51M * 0 0 GGCAGAACCGGAGGCAGTAGAAGCGGATGAGGTTGATTGAGGGAATGAGCT @@@DFFFFHHGGHIBGIEHGIJIIJJGAHEGHBFH8=FEHIHIAGGIIHHH HI:i:1 NH:i:1 NM:i:0
## ERR458493.199216 0 IX 25297 40 51M * 0 0 AGCGGATGAGGTTGATTGAGGGAATGAGCTTGAGCTACCAGAAGCGGATGA :=@DD@DDFABDFBFB>C+9;?1:EE>F9?D??FFG@F>BGIIDFEB@@EC HI:i:1 NH:i:1 NM:i:0
## ERR458493.410726 0 IX 25298 40 51M * 0 0 GCGGATGAGGTTGATTGAGGGAATGAGCTTGAGCTACCAGAAGCGGATGAG @@@=?D@DHDFAHHDAEHIIIFEEHIIIGIIGIIIJIIIFGGIGH>GBG@D HI:i:1 NH:i:1 NM:i:0
## ERR458493.108369 0 IX 25306 40 51M * 0 0 GGTTGATTGAGGGAATGAGCTTGAGCTACCAGAAGCGGATGAGACTGATTG ?@@DDDDBHADD?FBFEH?BHI@;CFCHFHIG9FCGGBH>FHBBF>EGGF< HI:i:1 NH:i:1 NM:i:0
## ERR458493.1047433 16 IX 25309 40 50M1S * 0 0 TGATTGAGGGAATGAGCTTGAGCTACCAGAAGCGGATGAGACTGATTGGCG JIJJJIIIGGCHEGGGHFDBGHEHGGIIGFEJIJJJJIHDHGHFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.1005912 0 IX 25310 40 51M * 0 0 GATTGAGGGAATGAGCTTGAGCTACCAGAAGCGGATGAGACTGATTGGCTG C@CFFFFFHHGHHIJJJIJJJJJJJJJJIJJJJJIJJJJIIJJIJJJJIJG HI:i:1 NH:i:1 NM:i:0
## ERR458493.121845 0 IX 25312 40 51M * 0 0 TTGAGGGAATGAGCTTGAGCTACCAGAAGCGGATGAGACTGATTGGCTGCC @@+4BDDDHDHHHGIGGBFHJEHGHGB>BGGI@GHHE9DFH3DCEBCGEFE HI:i:1 NH:i:1 NM:i:0
## ERR458493.124644 0 IX 25314 40 51M * 0 0 GAGGGAATGAGCTTGAGCTACCAGAAGCGGATGAGACTGATTGGCTGCCAG @+:BDADDDHDDFDEHGIIIIGIGIFEHEH<GDAEGBGHGGHIIEHI@F>F HI:i:1 NH:i:1 NM:i:0
## ERR458493.502454 16 IX 25322 20 51M * 0 0 GAGCTTGAGCTACCAGAAGCGGATGAGACTGATTGGCTGCCAGAGCTCGAG EIEFFCIIIGFFIIFIGEIIIIIIIIFHEGFIIIIIIIFFDBFDD@4+;@@ HI:i:1 NH:i:1 NM:i:1
## ERR458493.750363 16 IX 25329 40 51M * 0 0 AGCTACCAGAAGCGGATGAGACTGATTGGCTGCCAGAGCTTGAGACACCAG HGBGDJJJJJIGIEGIIHHHFIIIIJIJJIJIJIHIGJHHFGHDDFDFCCB HI:i:1 NH:i:1 NM:i:0
## ERR458493.685910 16 IX 25332 40 47M4S * 0 0 TACCAGAAGCGGATGAGACTGATTGGCTGCCAGAGCTTGAGACACCAACAG CDIIIIIIIHIGGIGEGGIIIIHHIGGIIGGIFIIGEIHBFFC=BB=4+@= HI:i:1 NH:i:1 NM:i:0
## ERR458493.790718 0 IX 25340 40 51M * 0 0 GCGGATGAGACTGATTGGCTGCCAGAGCTTGAGACACCAGAAGTAGAGGAA ?<?1ADDBHAHFDGHEEGHIG@GBFC@A?CDHHG>A9D?GGEC9B?FEDGE HI:i:1 NH:i:1 NM:i:0
## ERR458493.776594 16 IX 25351 20 15S36M * 0 0 AGAAGCTGAGGAAGATGATTGGCTGCCAGAGCTTGAGACACCAGAAGTCGA JJJJJJJJJJGJJJIIJIIJJIIJJJJIJJJJJIJJIHHFHHHFFDB=1@B HI:i:1 NH:i:1 NM:i:1
## ERR458493.878716 0 IX 25352 20 6S45M * 0 0 CTAGGAGATTGGGTGCCAGAGCTTGAGACACCAGAAGTAGAGGAAGCTGAT ???+=D>D3CDDDE1CEEFEFCFEA>EB@EED???DDBBD*?D@DDADDD9 HI:i:1 NH:i:1 NM:i:1
## ERR458493.541162 0 IX 25356 40 51M * 0 0 GGCTGCCAGAGCTTGAGACACCAGAAGTAGAGGAAGCTGATTGACTGACAG CCCFFFFFHHHHHJJJJJJJJJJEGHIFGIIIJGJJIJJIIJIJJJJJJJI HI:i:1 NH:i:1 NM:i:0
## ERR458493.86966 0 IX 25361 40 51M * 0 0 CCAGAGCTTGAGACACCAGAAGTAGAGGAAGCTGATTGACTGACAGAGCTT ?+1AADFFHHGHHJJJJJIIIJEHIIIJHIIGJJGHIJIJJJIGIJJJJII HI:i:1 NH:i:1 NM:i:0
## ERR458493.892577 0 IX 25362 10 51M * 0 0 CAGAGCTTGGGACACCAGAAGTAGAGGAAGCTGATTGACTGACAGAGCTTG ?;?BDDDBF+<3AFGGBC;EGBDEAACHG@FBCGGGH3BF@9?GF;FD@>D HI:i:1 NH:i:1 NM:i:1
## ERR458493.960273 0 IX 25362 40 51M * 0 0 CAGAGCTTGAGACACCAGAAGTAGAGGAAGCTGATTGACTGACAGAGCTTG C@@?DEFFHGHHHJJJJJJJJFHJJIJJJJJJJJJJGIIIIFIJJJJJJJJ HI:i:1 NH:i:1 NM:i:0
## ERR458493.812659 16 IX 25363 40 51M * 0 0 AGAGCTTGAGACACCAGAAGTAGAGGAAGCTGATTGACTGACAGAGCTTGA HHIIIHHGHEAIIIHGGGGGHHIIIGGHGIGGGEIIHIHFHHFDBADD;@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.541541 16 IX 25364 40 51M * 0 0 GAGCTTGAGACACCAGAAGTAGAGGAAGCTGATTGACTGACAGAGCTTGAG CFB0@GFFDCF?B>HE@@E<<FEGHH>D>>@GIFFA>GDHHFBD?DDD@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.797580 16 IX 25364 40 51M * 0 0 GAGCTTGAGACACCAGAAGTAGAGGAAGCTGATTGACTGACAGAGCTTGAG DJJIIIJIHGF?IGIJIGIIIIJIGJIJJJJIGIJIJJHHHHHFFDDA1C@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.669865 16 IX 25371 13 51M * 0 0 AGACACCAGAAGTAGAGGAAGCTGATTGACTGACAGAGCTTGAGACATCAG BAF?:E>ABC91?BEFFA<++HC4<FDA:A:4?GGIIIDFA4FBB4?D<=; HI:i:1 NH:i:1 NM:i:0
## ERR458493.307363 16 IX 25388 13 51M * 0 0 GAAGCTGATTGACTGACAGAGCTTGAGACATCAGAAGTAGAGGAAGCTGAT F>IIGHHD<>EAGHE?@HG>DBIJIHDBGGJIGGHHFBHHDHFDDDDD@C@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.636136 16 IX 25389 13 51M * 0 0 AAGCTGATTGACTGACAGAGCTTGAGACATCAGAAGTAGAGGAAGCTGATT IIJJIJJIHHEJIIIJJJJIJJIIIHJIGGIIIJJIJJGHHHHFFFDD:@B HI:i:1 NH:i:1 NM:i:0
## ERR458493.242207 16 IX 25430 0 51M * 0 0 GAAGCTGATTGACTGACAGAGCTTGAGACATCAGAAGCTGAGGAAGAAGAT IIIIGGIIIHGAIIIIHFEIIIIGHBBDGFFGHGEGGEAHDFDDDA=+=== HI:i:1 NH:i:1 NM:i:0
## ERR458493.937759 16 IX 25434 0 51M * 0 0 CTGATTGACTGACAGAGCTTGAGACATCAGAAGCTGAGGAAGAAGATTGAC IJJIJJIIJIIJJJJJJJJJJJJIJJJJJIJJJJJJJJHHHHHFFFFD?CC HI:i:1 NH:i:1 NM:i:0
## ERR458493.399327 0 IX 25446 0 51M * 0 0 CAGAGCTTGAGACATCAGAAGCTGAGGAAGAAGATTGACTGACAGAGCTTG @@@DDEFFHHGHGGGIIIJJJJIGHHI@GHEHIIIJIJGHCHEHGHGGGGI HI:i:1 NH:i:1 NM:i:0
## ERR458493.148179 0 IX 25455 0 51M * 0 0 AGACATCAGAAGCTGAGGAAGAAGATTGACTGACAGAGCTTGATACATCAG @1;DDEFDHFFHHIIIIICEFFEGGHGDHHIGGIEHHIGHIDHICHEIIIG HI:i:1 NH:i:1 NM:i:0
## ERR458493.167046 0 IX 25474 0 51M * 0 0 AGAAGATTGACTGACAGAGCTTGATACATCAGAAGCTGAAGAAGAAGATTG C@CFFFFFHHHHHJJJJJJJJJJJIJJJJJJJJJJJJJJIIHJJFIIIJJG HI:i:1 NH:i:1 NM:i:0
## ERR458493.363833 0 IX 25474 0 51M * 0 0 AGAAGATTGACTGACAGAGCTTGATACATCAGAAGCTGAAGAAGAAGATTG @@1?BDEFHDHHHIIIIIIIIIIIIIIIIIIIEHIIIIBGIDFHGIIIEIG HI:i:1 NH:i:1 NM:i:0
## ERR458493.883809 16 IX 25474 0 51M * 0 0 AGAAGATTGACTGACAGAGCTTGATACATCAGAAGCTGAAGAAGAAGATTG IJJIJIJJIGIHHGGHJIGGGHJIDDHGIJJJJIIJIJHHGHHFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.905582 16 IX 25474 0 51M * 0 0 AGAAGATTGACTGACAGAGCTTGATACATCAGAAGCTGAAGAAGAAGATTG JJJJJJJIIHJIJIJJJJIJJJJJJIJJIJJIJJJJJJHHHHHFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.889560 16 IX 25477 0 48M3S * 0 0 AGATTGACTGACAGAGCTTGATACATCAGAAGCTGAAGAAGAAGATTGGAT IIIIIGFGFE<IFEIIIIIIIIIFIIHHCEHEHFFEIIFFFFFDDA=41+? HI:i:1 NH:i:1 NM:i:0
## ERR458493.1050903 16 IX 25480 0 51M * 0 0 TTGACTGACAGAGCTTGATACATCAGAAGCTGAAGAAGAAGATTGACTGAC CJHHD99IIGDIGGBCIIGHHGGGJJJJIGGJJIJIG@HBHHDFEDFF@C@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.765941 16 IX 25480 0 51M * 0 0 TTGACTGACAGAGCTTGATACATCAGAAGCTGAAGAAGAAGATTGACTGAC IJIHJJJIJJJJIJJJJJJIJJJJJJJJIJJJJJJJJJHHHHHFFEDD;B@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.308266 16 IX 25484 0 51M * 0 0 CTGACAGAGCTTGATACATCAGAAGCTGAAGAAGAAGATTGACTGACCGAG CGIIHGIHCGHGDFE?EIGIJIIGIGGJIGHGIJIIIHDGHHD@AA=4@@= HI:i:1 NH:i:1 NM:i:1
## ERR458493.174424 16 IX 25512 0 51M * 0 0 AAGAAGAAGATTGACTGACAGAGCTTGAGACATCAGAAGTAGAGGAAGCTG GJIJJJIJJIIIIGJIGFIJJJIIJIJJJIJJJJJJIJHHHHHFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.337570 0 IX 25516 0 51M * 0 0 AGAAGATTGACTGACAGAGCTTGAGACATCAGAAGTAGAGGAAGCTGATTG ;@<1=3B+=ADBD>:?CG>FG@;FEAF@?<<BEGG1::?D@8DD)*?D*?9 HI:i:1 NH:i:1 NM:i:0
## ERR458493.860660 0 IX 25517 0 51M * 0 0 GAAGATTGACTGACAGAGCTTGAGACATCAGAAGTAGAGGAAGCTGATTGA @@@DDBD?CDBFAGIF;CBHGE9FFCGHEE<9?E*9?CDF;DBFG@GIIII HI:i:1 NH:i:1 NM:i:0
## ERR458493.7474 0 IX 25520 0 3S48M * 0 0 GCGGATTGACTGACAGAGCTTGAGACATCAGAAGTAGAGGAAGCTGATTGA =;1+=DDDFH?HHIIIIIIGIIIIIEHGIGGGIICGEGIIIIIIGIIIIGI HI:i:1 NH:i:1 NM:i:0
## ERR458493.130133 0 IX 25538 0 51M * 0 0 GAGACATCAGAAGTAGAGGAAGCTGATTGACTGACAGAGCTTGAGACATCA @+=DDDDDHDF<FG@F?DADFHGIIII@HIIIGIIIIDGIIIGGHGEGIIE HI:i:1 NH:i:1 NM:i:0
## ERR458493.338704 16 IX 25541 0 51M * 0 0 ACATCAGAAGTAGAGGAAGCTGATTGACTGACAGAGCTTGAGACACCAGAA IJIGJJJJJJIJJJJIEIJJJJJIIHHJHHEGJIIJJIHHHHFDDFFF@B@ HI:i:1 NH:i:1 NM:i:1
## ERR458493.417106 0 IX 25545 0 51M * 0 0 CAGAAGTAGAGGAAGCTGATTGACTGACAGAGCTTGAGACATCAGAAGCTG B<@DFFBDHHHHHJIJJJJJJJJJJJIJJJIJJJJJJJJJJJJJJGIJIJJ HI:i:1 NH:i:1 NM:i:0
## ERR458493.687814 16 IX 25550 0 51M * 0 0 GTAGAGGAAGCTGATTGACTGACAGAGCTTGAGACATCAGAAGCTGAGGAA @GGHIIHEHGEC??CC9ABCF<?BFC<?BBC<AGE<EIDHHDHDB=4+;@? HI:i:1 NH:i:1 NM:i:0
## ERR458493.514366 0 IX 25560 13 42M9S * 0 0 CTGATTGACTGACAGAGCTTGAGACATCAGAAGCTGAGGAAGAGGATTGAC @@<:DDDDFHHHFIGGCGIIGHFC>GHIIIIIIIIIHIIEHIE)8?GHHI@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.101742 0 IX 25566 40 51M * 0 0 GACTGACAGAGCTTGAGACATCAGAAGCTGAGGAAGCTGATTGACTGACAG B@@FFFFFHHHHGJJJIJJIJJJJJIJJIJIJJJJIJGIHIJIGHIJGIJI HI:i:1 NH:i:1 NM:i:0
## ERR458493.529915 0 IX 25583 10 51M * 0 0 ACATCAGAAGCTGAGGAAGCTGAGTGACTGACAGAGCTTGAGACATCAGAA CCCFFFFFHHHHHJJJJJJJJJJ+AFGIIIIJJJIJJJIII<DHJJJJJJG HI:i:1 NH:i:1 NM:i:1
## ERR458493.870966 16 IX 25589 13 51M * 0 0 GAAGCTGAGGAAGCTGATTGACTGACAGAGCTTGAGACATCAGAAGCTGAG IJJJJJJJJJJJJJJJJJJIHJJJIJJJIJJJJJJIHJHHHHHFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.740173 16 IX 25687 0 51M * 0 0 AGATTGACTGACAGAGCTTGAGACATCAGAAGTAGAGGAAGCTGATTGACT IJJHHGCGIHGGJIHGGJIIIIHIIGJIIJIIIGGIJJHHFHGFFFFD@CC HI:i:1 NH:i:1 NM:i:0
## ERR458493.683894 0 IX 25689 0 51M * 0 0 ATTGACTGACAGAGCTTGAGACATCAGAAGTAGAGGAAGCTGATTGACTGA @@:ADDDDHFHHBGHIJCEEEDGGIIJHHGFGHIIIIIJJJJ>GHC@HG@F HI:i:1 NH:i:1 NM:i:0
## ERR458493.481488 16 IX 25691 0 51M * 0 0 TGACTGACAGAGCTTGAGACACCAGAAGTAGAGGAAGCTGATTGACTGACA FFFEFFCDEG@FFECGFGAA9IHEEIFIIFIDFFIIGF>DFFDDDBA+@7@ HI:i:1 NH:i:1 NM:i:1
## ERR458493.1009276 16 IX 25691 0 51M * 0 0 TGACTGACAGAGCTTGAGACATCAGAAGTAGAGGAAGCTGATTGACTGAAA B::CDC?DDD<<3?*>DCFCC9>FF@EEEC<+<A?A:?:A:DDDDD=+:1; HI:i:1 NH:i:1 NM:i:1
## ERR458493.324848 16 IX 25691 0 51M * 0 0 TGACTGACAGAGCTTGAGACATCAGAAGTAGAGGAAGCTGATTGACTGACA DF?FEGF>FBFFFD?4EHHBEH<H@@CHBHHIGE<CA+FFFFDDDDB+@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.640111 0 IX 25692 0 51M * 0 0 GACTGACAGAGCTTGAGACATCAGAAGTAGAGGAAGCTGATTGACTGACAG @?@FFDFDDHFFHGHBFHGJIJGIEHGCHHHGHGHIDHG@GHGEHGIIIGG HI:i:1 NH:i:1 NM:i:0
## ERR458493.59532 0 IX 25692 0 51M * 0 0 GACTGACAGAGCTTGAGACATCAGAAGTAGAGGAAGCTGATTGACTGACAG ?@BFFFFFHHHGGJIGGGIIJIGJIII<EEEGHFIIGBHGIGGIEIAFAAB HI:i:1 NH:i:1 NM:i:0
## ERR458493.52354 0 IX 25694 0 51M * 0 0 CTGACAGAGCTTGAGACATCAGAAGTAGAGGAAGCTGATTGACTGACAGAG ?11ABDDFFHHHHJJDEHHJJGHAFAHHDEGEHFGGIFGGFHIIIIHIGHC HI:i:1 NH:i:1 NM:i:0
## ERR458493.688704 16 IX 25697 0 51M * 0 0 ACAGAGCTTGAGACATCAGAAGTAGAGGAAGCTGATTGACTGACAGAGCTT @IGG@IGIIEFF?FIFHEIIHHFIFFFEIIIIDFGFFFFFD:FDDDA+@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.461315 0 IX 25700 0 51M * 0 0 GAGCTTGAGACATCAGAAGTAGAGGAAGCTGATTGACTGACAGAGCTTGAG ?@+BB;;:=CDFA9CB<:?ICHF9?33+A91*?F?*?:CDDGD399*?DCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.422016 0 IX 25702 0 51M * 0 0 GCTTGAGACATCAGAAGTAGAGGAAGCTGATTGACTGACAGAGCTTGAGAC ?@+4=BDFHHHHHJJIJGHIJHIIJJIJJIJJJIJJFHHJJFIIJIJJJIJ HI:i:1 NH:i:1 NM:i:0
## ERR458493.1000255 16 IX 25702 0 51M * 0 0 GCTTGAGACATCAGAAGTAGAGGAAGCTGATTGACTGACAGAGCTTGAAAC EEDCBDDD?1*<FFEEFFDE<FCEIEFDCFEFFADEEEDDDDDD?A4+1?? HI:i:1 NH:i:1 NM:i:1
## ERR458493.702307 0 IX 25716 13 51M * 0 0 AAGTAGAGGAAGCTGATTGACTGACAGAGCTTGAGACATCAGAAGCTGAAG @@@DDDDDHHH?BAC?@FBHHIBGFGGHEEEDGGBGHH99D@?DBB9DFFB HI:i:1 NH:i:1 NM:i:0
## ERR458493.99019 16 IX 25721 13 4S43M4S * 0 0 AGCTGAGGAAGCTGATTGACTGACAGAGCTTGAGACATCAGAAGCTGCGGA JGIIJJJJJJIJJJIJJJIJJJJJJJJIJJJJIIIJJIHHHHHFDDB4;@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.370131 0 IX 25726 40 51M * 0 0 AGCTGATTGACTGACAGAGCTTGAGACATCAGAAGCTGAAGAAGAAGATTG @@@DDBDFBFHHHEHIG@F>HHG@GEDGDGEHEGCHIICGIIGGIGHFHEH HI:i:1 NH:i:1 NM:i:0
## ERR458493.506658 0 IX 25730 40 51M * 0 0 GATTGACTGACAGAGCTTGAGACATCAGAAGCTGAAGAAGAAGATTGACTG ??@+4ADDHFHHFGFGIIIGEGCEHHGIBHIIIIGIJFIIIIHGGGCGGIE HI:i:1 NH:i:1 NM:i:0
## ERR458493.24312 0 IX 25738 40 51M * 0 0 GACAGAGCTTGAGACATCAGAAGCTGAAGAAGAAGATTGACTGACAGAGCT CCCFFFFFHHHHHJJJJJJJJJJJJJJJJEIJJJJJJIIJJJJJJJJJJJJ HI:i:1 NH:i:1 NM:i:0
## ERR458493.738092 16 IX 25749 20 51M * 0 0 AGACATCAGAAGCTGAAGAAGAAGATTGACTGACAGAGCTTGATACATAAG HIGJJJIJIIJIHDJJJIJJJJJJIJJGJIFJJHGIJIHHFHFDDDDD7CB HI:i:1 NH:i:1 NM:i:1
## ERR458493.507408 0 IX 25762 40 51M * 0 0 TGAAGAAGAAGATTGACTGACAGAGCTTGATACATCAGAAACGGATGAAGA @;?DDDDDFHHHHIIBEHHII@GHIICHIIIIIIIIIIIIIIIIFEHCGD> HI:i:1 NH:i:1 NM:i:0
## ERR458493.114246 16 IX 25768 20 51M * 0 0 AGAAGATTGACTGACAGAGCTTGATACATCAGAAACGGATGAAGAGGATCG FIIIGIIIHHHGC:HFIIIIIHHHADIHBIGEDE?GIIHDHF?DDD?B<?= HI:i:1 NH:i:1 NM:i:1
## ERR458493.452870 16 IX 25771 40 51M * 0 0 AGATTGACTGACAGAGCTTGATACATCAGAAACGGATGAAGAGGATTGGCT IIJJJIJIHIGIIIJIJIGJIHFJJIIJJJJIJJJJIJHHHFGFFFFD@CC HI:i:1 NH:i:1 NM:i:0
## ERR458493.826739 16 IX 25774 20 51M * 0 0 TTGACTGACAGAGCTTGATACATCAGAAACGGATGAAGAGGATTGGCTGCC FB<D>AFFAGGGFBFHGC4?FHC9>AF<FEBHCIGGD@FB>HDEFDBA:1; HI:i:1 NH:i:1 NM:i:1
## ERR458493.666318 0 IX 25776 40 51M * 0 0 GACTGACAGAGCTTGATACATCAGAAACGGATGAAGAGGATTGGCTGACAG B@BFFFFFHFDHHGJIIJJJJJFJGGHIJIGIIEGGHIGDIIGGFGIIIIG HI:i:1 NH:i:1 NM:i:0
## ERR458493.63649 0 IX 25776 40 51M * 0 0 GACTGACAGAGCTTGATACATCAGAAACGGATGAAGAGGATTGGCTGACAG @@=DDFFDFHHHGJIIJIJJJJJIGGIJIICHJJIIJJIIIJGIJGIIIIF HI:i:1 NH:i:1 NM:i:0
## ERR458493.72713 0 IX 25797 40 51M * 0 0 CAGAAACGGATGAAGAGGATTGGCTGACAGAGCTTGAGACAACAGAGGAAG @@@DDDDDF?FHABG9FH>F@EHIIGBGCE3DFA@B8D;<;?BFC>FGCH; HI:i:1 NH:i:1 NM:i:0
## ERR458493.240821 16 IX 25817 40 51M * 0 0 TGGCTGACAGAGCTTGAGACAACAGAGGAAGAGGATGGACCGACAGAGCTC EFBB99<HB<F<9D<HGD@HFBGEAGH>HFHBB<<DGEHDFA2BDBD?@?8 HI:i:1 NH:i:1 NM:i:0
## ERR458493.588059 0 IX 25822 40 51M * 0 0 GACAGAGCTTGAGACAACAGAGGAAGAGGATGGACCGACAGAGCTCGAGAC ?<?+4+=BDDDB<EEEEEDI<D>EE8AEE?DE<;<BDDDD@B3=).;8@C# HI:i:1 NH:i:1 NM:i:0
## ERR458493.792673 0 IX 25829 20 51M * 0 0 CGTGAGACAACAGAGGAAGAGGATGGACCGACAGAGCTCGAGACACCAGAA @+:BDDDDFHFHFGII>EGAHIIIII;FHEADHGGIIGIDEGGGGIIIIFD HI:i:1 NH:i:1 NM:i:1
## ERR458493.964993 0 IX 25829 40 51M * 0 0 CTTGAGACAACAGAGGAAGAGGATGGACCGACAGAGCTCGAGACACCAGAA =<+=ADBDFHHFDBFECFGAEGBGGEHIIIIDIGIBBBGAF<@BACEEH=A HI:i:1 NH:i:1 NM:i:0
## ERR458493.988291 0 IX 25831 40 51M * 0 0 TGAGACAACAGAGGAAGAGGATGGACCGACAGAGCTCGAGACACCAGAAGT @@@FDFFDHHH>BGGIIEHGBCGIGGHGEGIHCBF<FHE<;AGCHAHEHH= HI:i:1 NH:i:1 NM:i:0
## ERR458493.1012705 0 IX 25854 40 51M * 0 0 GACCGACAGAGCTCGAGACACCAGAAGTGGATGAAGCTGATTGTGTGATAG @CCDFFFFHHHHHJIJJHIJJJIJHGHDFEHIGIJJIJJ9FHIHGFHGCGE HI:i:1 NH:i:1 NM:i:0
## ERR458493.1080218 0 IX 25862 40 51M * 0 0 GAGCTCGAGACACCAGAAGTGGATGAAGCTGATTGTGTGATAGAACTTAAG ;=++=BD<DDDBDIDEDEF<AE>?ECEEEDE39?D?D:BDDD@DIDDEBEI HI:i:1 NH:i:1 NM:i:0
## ERR458493.741229 0 IX 25868 40 51M * 0 0 GAGACACCAGAAGTGGATGAAGCTGATTGTGTGATAGAACTTAAGATATCG @C@FFFFFHFDHHAEHEIIIIIIIJDHHHDHBGAGEIDFHGIBHIIIJIJJ HI:i:1 NH:i:1 NM:i:0
## ERR458493.870592 16 IX 25868 40 51M * 0 0 GAGACACCAGAAGTGGATGAAGCTGATTGTGTGATAGAACTTAAGATATCG GB90:*3*<GDD@DB@1:*9<B@C:F;?@BB<HB>GFBFC44BD>ADD@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.446988 0 IX 25869 20 41M10S * 0 0 AGACACCAGAAGTGGATGAAGCTGATTGTGTGATAGAGCTTGAGACAACAG C<@FFDDFHGHHFHIHIIJIIJJIJJJJDHFHIJIGGG9FHIIIHIJIIJI HI:i:1 NH:i:1 NM:i:1
## ERR458493.975019 16 IX 25877 40 51M * 0 0 GAAGTGGATGAAGCTGATTGTGTGATAGAACTTAAGATATCGGTAGCGGAC IJJJIIJJJJJIIIGHJJJJJIJJGJJJIIJJJJJJIJHHHHHFFFFF@BC HI:i:1 NH:i:1 NM:i:0
## ERR458493.646693 16 IX 25877 40 51M * 0 0 GAAGTGGATGAAGCTGATTGTGTGATAGAACTTAAGATATCGGTAGCGGAC IJJJJJJJJJJJJJJJJJJJJJJJJIIJIIJJIJJJIIHDHHHFFFFD@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.291361 16 IX 25883 40 51M * 0 0 GATGAAGCTGATTGTGTGATAGAACTTAAGATATCGGTAGCGGACGAGGAT GGGCEGGHG@EEFB*G@BFDDDD?@@HCIGG<FC@FE1FDFC;DDBA+?@= HI:i:1 NH:i:1 NM:i:0
## ERR458493.364483 16 IX 25890 40 51M * 0 0 CTGATTGTGTGATAGAACTTAAGATATCGGTAGCGGACGAGGATGTCTCAG GJJJJIJJJJJIJJJIIIIJJIIGJJJJJJJJIJIHFJHDHHHFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.993454 0 IX 25891 40 51M * 0 0 TGATTGTGTGATAGAACTTAAGATATCGGTAGCGGACGAGGATGTCTCAGC @@@DDFDFHHHHHJJJJJJJJJJJJJJJJHIJIJJJJJIJJIGIGIJJJJI HI:i:1 NH:i:1 NM:i:0
## ERR458493.667893 16 IX 25913 40 51M * 0 0 ATATCGGTAGCGGACGAGGATGTCTCAGCAATAGAACTACTGATCGAAACG EBF@??(???6C?:8C@EEHBE?EH>BB?FBHFA<4EGCFFBBAB;AA1<? HI:i:1 NH:i:1 NM:i:0
## ERR458493.73254 16 IX 25921 40 51M * 0 0 AGCGGACGAGGATGTCTCAGCAATAGAACTACTGATCGAAACGGATCCAGA AIGHEGIIIHF9BIHFC9GD9F9IIHHEGFHEFE<;GBHHHHFDDD7;1<@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.799127 16 IX 25921 40 51M * 0 0 AGCGGACGAGGATGTCTCAGCAATAGAACTACTGATCGAAACGGATCCAGA HGGGEIIEIIGGGGEIGIIIIIGGHBBCHCHGIHGIIGGGHHFFFEAD<@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.5939 0 IX 25927 40 51M * 0 0 CGAGGATGTCTCAGCAATAGAACTACTGATCGAAACGGATCCAGAGGAACT CCCFFFFFHFHHHGIIIGJJIIIJGIIJIIJIJJJJIJJJIIIIGIJGHIH HI:i:1 NH:i:1 NM:i:0
## ERR458493.58277 16 IX 25945 40 51M * 0 0 AGAACTACTGATCGAAACGGATCCAGAGGAACTTGATACCACAGAGGTGGA HHHHHF<GGDBFGGHFCGIHGGGIGJJIIGFGGJJHGAHFFHFFFFFD@@? HI:i:1 NH:i:1 NM:i:0
## ERR458493.205493 16 IX 25946 40 51M * 0 0 GAACTACTGATCGAAACGGATCCAGAGGAACTTGATACCACAGAGGTGGAT #@?HGGGDB:)EGGC8CGHEB?<@HH@BA3B<EDB<<9C<BDFDDDB=@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.703751 16 IX 25949 40 51M * 0 0 CTACTGATCGAAACGGATCCAGAGGAACTTGATACCACAGAGGTGGATGAG <B?@?0*0?DHFCDFBE?>HEGHHFBFBHAH@CGE<AIHFD?ADDDDD@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.268623 0 IX 25955 40 51M * 0 0 ATCGAAACGGATCCAGAGGAACTTGATACCACAGAGGTGGATGAGCTGTTG CCCFFFFFHHHHHJJJEHJHJJJJJJJJJJJJJJIIJ?DGGIJIJJJJJJI HI:i:1 NH:i:1 NM:i:0
## ERR458493.1073989 0 IX 25958 40 51M * 0 0 GAAACGGATCCAGAGGAACTTGATACCACAGAGGTGGATGAGCTGTTGCTT @@@:DDFFHFFHHGIJEHGIIIGGHJJJJIIHGI:BGDHIJJJJJGIIEHH HI:i:1 NH:i:1 NM:i:0
## ERR458493.132828 16 IX 25977 40 51M * 0 0 TTGATACCACAGAGGTGGATGAGCTGTTGCTTGAAATTGAGGTGCTATTGC JJIHGGHGIGEGEHHFGIIIIIIHD?DBDHIIJJJIIJHHFHDDDDD?C@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.145741 16 IX 25980 20 51M * 0 0 ATACCACAGAGGTGGATGAGCTGTTGCTTGAAATTGAGGTGCTATTGATAT BE9IGIIIIIJJIJJJJJIIGFIJJIHJJJJIJIIJIJGGHFHEDDB=@@@ HI:i:1 NH:i:1 NM:i:1
## ERR458493.1046703 16 IX 25987 40 51M * 0 0 AGAGGTGGATGAGCTGTTGCTTGAAATTGAGGTGCTATTGCTATAGTACTG HGHDDD9GG@DIIIHEF?FIGIGHEGGIH@GGGFC<G>DHFHFFED;DC@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.1018210 0 IX 25993 40 51M * 0 0 GGATGAGCTGTTGCTTGAAATTGAGGTGCTATTGCTATAGTACTGACCCAA ?<?BDBBDHFDFFBFGFG@?IGBCFHCFGHE@FAGGGGCG?FHIIE??BGH HI:i:1 NH:i:1 NM:i:0
## ERR458493.398571 16 IX 25995 40 51M * 0 0 ATGAGCTGTTGCTTGAAATTGAGGTGCTATTGCTATAGTACTGACCCAATG EJIJIJJJJIJJJJJJIJJJIJJJJIJJIJJJIJJJJIHHDHHDFFDB:@B HI:i:1 NH:i:1 NM:i:0
## ERR458493.679571 16 IX 26004 40 51M * 0 0 TGCTTGAAATTGAGGTGCTATTGCTATAGTACTGACCCAATGCACTTTGAG IJJJJJJJJIJJJJJJJJJJJJJJJJJJJJIJIHFJJJHHGHFBFFFF@CB HI:i:1 NH:i:1 NM:i:0
## ERR458493.144725 16 IX 26011 40 51M * 0 0 AATTGAGGTGCTATTGCTATAGTACTGACCCAATGCACTTTGAGAGTAAAG JIJJIJJJJJIJJIGIGIIIIGHBIHE?JJFJJJHGBJHHGHGFFFD=CCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.3582 16 IX 26022 40 51M * 0 0 TATTGCTATAGTACTGACCCAATGCACTTTGAGAGTAAAGGGCCAAAGCTA JJJJJJJJJJIJIJJIIJJJJJJIIIJJJJJJJIJJJJHHHHGFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.938541 16 IX 26074 40 51M * 0 0 AGCAGCTTGGAATTTGTTTAAACGATTGAACATCTCTTTCAATTAGGCGAG FGGHEEGB:000:?1<FGGIGGIGHGHGHHA>ABBFEEAFFBF?DD=?@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.226291 16 IX 26150 40 51M * 0 0 GCTTGCTGTTTCTAAATTTACTCGAAATAAAAACAACACATACTGTCCAAA BCDDDBB3EIEFEDFEEEEEEF9C<+:DDC<;C;C2:,D?42>DBDD=??: HI:i:1 NH:i:1 NM:i:0
## ERR458493.694701 0 IX 26225 40 51M * 0 0 CAAAGCCTTTTATATCATTGTAGGCGCTGTTGCAGTTATTTGTGGGCTTCC CCCFFEFFFHHHHJJJJJJJIJJJJJJJJIIJJJJHHIJJJJIHIJJJJJJ HI:i:1 NH:i:1 NM:i:0
## ERR458493.98040 16 IX 26317 40 51M * 0 0 TTTGTACGGCAATTCCACAGGACTCTTCAAAGCCTGCGAGGTCTATGCGCG HGHHEIIIGFHF?EGGHGGED?EIIIHGIGHEDG@BGHGFHHFFFFFD@?; HI:i:1 NH:i:1 NM:i:0
## ERR458493.73195 16 IX 26325 40 51M * 0 0 GCAATTCCACAGGACTCTTCAAAGCCTGCGAGGTCTATGCGCGCGCGTTTG CCDEEDFC?FFFHEDHFFGEIIGGGBEGDJJIJJIJGGGHFHHFDDDF@CC HI:i:1 NH:i:1 NM:i:0
## ERR458493.223243 16 IX 26333 40 51M * 0 0 ACAGGACTCTTCAAAGCCTGCGAGGTCTATGCGCGCGCGTTTGGTAAGTAG A>;6(A>;7;;)?)(<AB;FDACFEFFD76D:C6E:C8AFF?BBDA=D?@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.46431 0 IX 26345 40 51M * 0 0 AAAGCCTGCGAGGTCTATGCGCGCGCGTTTGGTAAGTAGGAAATTGGGTAG @@<+=BDDHHFHHHHIIIIIGIIIIGIG>GHG=EHHA7?E@BCEE@CBACC HI:i:1 NH:i:1 NM:i:0
## ERR458493.100539 0 IX 26356 40 51M * 0 0 GGTCTATGCGCGCGCGTTTGGTAAGTAGGAAATTGGGTAGTAATTTTATGC B@@+=BDDFHHHHIJIHIJJJHIIIFIIIGIIJJIJJCFGIIJIIHHHHHE HI:i:1 NH:i:1 NM:i:0
## ERR458493.936873 16 IX 26359 40 51M * 0 0 CTATGCGCGCGCGTTTGGTAAGTAGGAAATTGGGTAGTAATTTTATGCAAA FFHBGHBGHBIHIIJIJJGHGHGJIHJJIJIIJGIIGJHHFHHFFFFFC@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.566528 16 IX 26387 40 51M * 0 0 ATTGGGTAGTAATTTTATGCAAACGCAATGGTTGCCGTGAATAACAACACG IIJJJIIIHEIGJJJIJJJIIIGJJIJIJJJJJJIJJIDHHGHFFFDD@;@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.697162 16 IX 26422 40 51M * 0 0 CGTGAATAACAACACGATCCCTTTAATATTATACGAAAACGTGCAGAACGG JJIIJJIIHIGHGDIGGHJIGJIJJJJJJIFHFJJJJIHHHHHFFFFFC@B HI:i:1 NH:i:1 NM:i:0
## ERR458493.842154 16 IX 26433 40 51M * 0 0 ACACGATCCCTTTAATATTATACGAAAACGTGCAGAACGGCACAAACCAAT FC?IGGHCEEHGGIF@HFBFB=IIIHHFGGHEGEGIIIFBDA?DDBFFC@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.228435 16 IX 26466 40 51M * 0 0 AGAACGGCACAAACCAATAGGCTATGACCAAAATGGTTCATGAGAAAGAAC 1GHF?BBDB7ECBCCBBEIGIIIIHIGIIIHGHBHHEIFHHDBDDDDD@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.762001 0 IX 26480 20 51M * 0 0 CATTAGGCTATGACCAAAATGGTTCATGAGAAAGAACCAGGAGTATTATGA @@+=ADEDFBFFHGIIIIGIIGICHIHICHGHIIIEGHHFHBD:?GHIIIG HI:i:1 NH:i:1 NM:i:1
## ERR458493.111525 16 IX 26502 40 51M * 0 0 TTCATGAGAAAGAACCAGGAGTATTATGAGTAGTAACCCAATTTTTTTTTG DAGGD@IIIIGGF=CCB?>BIHEBIHBEGIGHGGFDIHFFFAHDDDBA@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.216349 16 IX 26506 40 51M * 0 0 TGAGAAAGAACCAGGAGTATTATGAGTAGTAACCCAATTTTTTTTTGAGAC BIHGJJJIGHJJJIHDIJIJIIIJHGIJJJJJIJIJGGHHGHHFEDDD@@B HI:i:1 NH:i:1 NM:i:0
## ERR458493.808342 16 IX 26509 20 51M * 0 0 GAAAGAACCAGGAGTATTATGAGTAGTAACCCAATTTTTTTTTGAGACCTA IIIIGGHGF<IIHHGCHHECF@IHGGFF?GEBDFAIIIFHHDFFFFFF@@@ HI:i:1 NH:i:1 NM:i:1
## ERR458493.628949 16 IX 26510 40 51M * 0 0 AAAGAACCAGGAGTATTATGAGTAGTAACCCAATTTTTTTTTGAGACTTAC HD===8HAHFCGB?GB@@DHGEGB??D?F@FIHBIHFDHFHHFDDBDB@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.359876 16 IX 26531 40 51M * 0 0 GTAGTAACCCAATTTTTTTTTGAGACTTACGCTACTTTCGCTCAATCCAAC CD@;AE?HGCHGIDIGGIHGGIHIIHF@IHHBF<BHEGFHHFBDDDDD@@? HI:i:1 NH:i:1 NM:i:0
## ERR458493.507364 16 IX 26536 40 51M * 0 0 AACCCAATTTTTTTTTGAGACTTACGCTACTTTCGCTCAATCCAACACGAG BEEDCFEFIIIIGGGIGIGHGGC?6GEC:GE<GF<FDGBHFFHDDDA?11@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.816228 0 IX 26556 40 51M * 0 0 CTTACGCTACTTTCGCTCAATCCAACACGAGCAATTACAGTAATCCGTACT BB:BDFFFHHHHHJJJJJJJJJJJJJJJJJJJIJJJJJJJIIJJJJJJJJI HI:i:1 NH:i:1 NM:i:0
## ERR458493.947816 0 IX 26584 40 51M * 0 0 GAGCAATTACAGTAATCCGTACTCAAAATGAGCCGTTTGGACTCACTACAG @<;DDDEEHFHHFHIIJIGHIEHHIJIIGHGIEIGJJFIJIIJJIIIGGEH HI:i:1 NH:i:1 NM:i:0
## ERR458493.480041 16 IX 26643 13 49M2S * 0 0 ACGATTAAAAGATAGAAATAGGCATCCCCCGCCTTTGGAGCAAAAAAGATT BBB@@3@>@8(;3A@>>::::(83'9'FGFCFCBEAECA<:@FB?=44111 HI:i:1 NH:i:1 NM:i:2
## ERR458493.108031 0 IX 26672 40 51M * 0 0 CGCCTTTGGAGCAAAAAAAAATGGTGGTTGCTCTTGGTCGCGAAGTGAATC @<<DDFFFDACFFIIIGIIIEHIF9BAGGGGGIGIHIC@BGEEDD?@CCA; HI:i:1 NH:i:1 NM:i:0
## ERR458493.610951 0 IX 26687 40 51M * 0 0 AAAAAATGGTGGTTGCTCTTGGTCGCGAAGTGAATCTTCAACATTGCTGAG @1?DFFFFHFHHHIJIJJJJJICGEHJIIJBDGIIJIJIJJIJJJJJJJJH HI:i:1 NH:i:1 NM:i:0
## ERR458493.62314 0 IX 26689 20 51M * 0 0 AAATTGGTGGTTGCTCTTGGTCGCGAAGTGAATCTTCAACATTGCTGAGGG ?@B+ADFAFHHHHIIIJJJJHJJJIJJJIJJJIIIJIIJJJJJJIIIGIII HI:i:1 NH:i:1 NM:i:1
## ERR458493.575185 0 IX 26710 20 51M * 0 0 CGCGAAGTGAATCTTCAACATTGCTGAGGGGACGGCTGCCTGCCTTCACTG ??@=D?D2AF>DFGEBGGEHGIHGHIIEHDGBEHEG<FCAF@(..=@CEAE HI:i:1 NH:i:1 NM:i:1
## ERR458493.41892 0 IX 26717 40 51M * 0 0 TGAATCTTCAACATTGCTGAGGGGACGGCTGCCTGACTTCACTGTGACGCG @+:ADFFFHHHHHJJIJJIHGIIGDFEHIEIEGDDCGEHIIJJIIIJGIIH HI:i:1 NH:i:1 NM:i:0
## ERR458493.87156 16 IX 26720 20 51M * 0 0 ATCTTCAACATTGCTGAGGGGACGGCTGCCTGACTTCACTGTGACGCGTAC GGGIIIGIIIIIIIIHCIIIIHGIGGDFGIEIHIIIIIHGHHFDDA=41+: HI:i:1 NH:i:1 NM:i:1
## ERR458493.416161 0 IX 26742 40 51M * 0 0 CGGCTGCCTGACTTCACTGTGACGCGAACGCGATAATTTTTGCGTTAGCGA @8@=DDFFHDDFHIIIHCHHHIIIIIGIIIIGIIIIIIIIIBHHDBFA>AC HI:i:1 NH:i:1 NM:i:0
## ERR458493.762349 0 IX 26748 20 51M * 0 0 CCAGACTTCACTGTGACGCGAACGCGATAATTTTTGCGTTAGCGAGTGGCT ?@14ADDFHHHHHEGEEHGIIJJJIJJJIJIIJIJIIJIIEEEFEF8=A@D HI:i:1 NH:i:1 NM:i:1
## ERR458493.236104 0 IX 26748 40 51M * 0 0 CCTGACTTCACTGTGACGCGAACGCGATAATTTTTGCGTTAGCGAGTGGCT @@C1BDFFHHHHHIJJJIJJIIIJIJJJJIJJJJJJIJIJHHHHDFBCDDD HI:i:1 NH:i:1 NM:i:0
## ERR458493.765110 0 IX 26748 40 51M * 0 0 CCTGACTTCACTGTGACGCGAACGCGATAATTTTTGCGTTAGCGAGTGGCT BB1ADEFFHHHHHHIIIIJJJJJJJJJJJJJJJJJFIIIIEECEFFABCDD HI:i:1 NH:i:1 NM:i:0
## ERR458493.192541 16 IX 26796 40 51M * 0 0 GCTAGTAAGTCCCTTCAAAAGCCACTCCCTGTAAAATACCGCACCTTCCGG HHIIIIIIGDIIGIIIIIIIEIGDGGC9EIIFIIIHE@HHFCAFFFDF@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.176025 16 IX 26796 40 51M * 0 0 GCTAGTAAGTCCCTTCAAAAGCCACTCCCTGTAAAATACCGCACCTTCCGG ;;;4C==8<(0B<<D?B;GCEGD?ABC?9GGIIIHEC21:C3<1ADDD@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.816618 16 IX 26836 40 51M * 0 0 GCACCTTCCGGGTTGCTTTCATAATCAAACAAAGTTTAGGACTCCTCCCAA <GIHF?@EIIGEIGIIHGGHGGEGIGFF9CHGGHDDGHDFHFFDB;D?@@? HI:i:1 NH:i:1 NM:i:0
## ERR458493.146520 16 IX 26842 20 51M * 0 0 TCCGGGTTGCTTTCATAATCAAACAAAGTTTAGGACTCCTCCCAAGAATCG IIIIHGIHDIGGF@GIGD4HF;?FEHE>IHEF?HGGHFB?DHFDAB=+@?= HI:i:1 NH:i:1 NM:i:1
## ERR458493.386 16 IX 26860 40 51M * 0 0 TCAAACAAAGTTTAGGACTCCTCCCAAGACTCGTCCCTTGGCAAAAAAAAT FEC=.C>>GEHE@7=78GH>F?HD@FDBBABEIGBBEGFFDHDFDDDB@@= HI:i:1 NH:i:1 NM:i:0
## ERR458493.549952 0 IX 26958 40 51M * 0 0 TTGAAATGCCACCCAATTTGCAAACATGATTCATGGATCTTGGAATTATCA ??@BBDD+AD?B;E7:CGCEH:A??EBBAEFFE>FIFC*9?F<39??G### HI:i:1 NH:i:1 NM:i:0
## ERR458493.733001 16 IX 27015 20 51M * 0 0 TTAACGTTATTTGCACCTATCCACGAGAGAAAGCAAGGAACATAATTCAGT IIHFIJJIGIHGGFCIGJIGIHEGGJJIIIIJJGIIIIHFHHHFFFDB+C@ HI:i:1 NH:i:1 NM:i:1
## ERR458493.120271 16 IX 27101 40 48M3S * 0 0 CTAGCGGATTCATATAACCCGCTAGATAAGCCATATATTACTAATTGCAAG AIHGGF@GIGEEGIIHFIIIGIIGIIGIIHGIIIIIGC4FHBHDDA==:1? HI:i:1 NH:i:1 NM:i:0
## ERR458493.829049 16 IX 27113 40 51M * 0 0 TATAACCCGCTAGATAAGCCATATATTACTAATTGCCTTTCGTGTCTATTT GHGG@JJJJGIIIIIJJJJIJIJIIIHEJJJIJJIJJJHHHGHFDDA4?@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.908219 0 IX 27231 40 51M * 0 0 ATCAGTGTCCGAGTTACACTCAATGAAAAACTTTTGGTGCAATTACGGGAG @@@DDDFFGHHHHJJIGGGIGIIIGJJGIEHGIGI@GDHGGIIIGIHII<F HI:i:1 NH:i:1 NM:i:0
## ERR458493.679776 16 IX 27330 20 49M2S * 0 0 CAAGTTCAAAAGCAGGCATTCGCTTCGAAAGGACTGAACCTTGGTCTACCG IIIIGIIIIIIIIIIIIIIIIIIIIIGIIIIIHIIIIHHDFHDDDD=+@@< HI:i:1 NH:i:1 NM:i:1
## ERR458493.816041 0 IX 27336 40 51M * 0 0 CAAAAGCAGGCATTCGCTTCGAAAGGACTGAACCTTGGTCTTCTGCGACAG ?:8BDDDAF?<DBFEHI@FHFIEEHI<GABGHIIIGEGGBBFGGCDAG>FH HI:i:1 NH:i:1 NM:i:0
## ERR458493.225683 0 IX 27342 40 51M * 0 0 CAGGCATTCGCTTCGAAAGGACTGAACCTTGGTCTTCTGCGACAGCAATTT ?<@BAD8:<DF?F?CGEF:3+A9:C?*?3?FFGBBDBGE############ HI:i:1 NH:i:1 NM:i:0
## ERR458493.491425 16 IX 27429 40 51M * 0 0 TTGTAGTTTTAAGGCGGGCATTATTTTGCTGAAATGGAATTTTGATAGCGG GHGHIJIGIEIJJIJIHBJJIGJJIGJHGIIHGGGJIJHHGHHFFDB+@@B HI:i:1 NH:i:1 NM:i:0
## ERR458493.511097 0 IX 27444 20 51M * 0 0 GGGCATTGTTTTGCTGAAATGGAATTTTGATAGCGGTAGTTCGATCATATG ?@<1BDD2ACDFFIFB@FIICD3BHFCEFEHIIDEGCFGFIFE89D@8BDF HI:i:1 NH:i:1 NM:i:1
## ERR458493.597055 0 IX 27447 40 51M * 0 0 CATTATTTTGCTGAAATGGAATTTTGATAGCGGTAGTTCGATCATATGGCA @@1ABDDFHGHHHIIIJJJIHIIJJJJJJJIIJIJJGHGIIIJJJJJIJJI HI:i:1 NH:i:1 NM:i:0
## ERR458493.1022937 0 IX 27451 20 51M * 0 0 ATTGTGCTGAAATGGAATTTTGATAGCGGTAGTTCGATCATATGGCAATGC :+14AABDF<=A<2:AFFIFB9+AEB@FA)?DGIEIACGFFG@DFCCEFFB HI:i:1 NH:i:1 NM:i:1
## ERR458493.479849 16 IX 27470 40 51M * 0 0 TTGATAGCGGTAGTTCGATCATATGGCAATGCAGTACGATTTGTATATCGA JIGJJIJJIIGIIJJJIIJIJJIJJJIJJIIIJHIHJJHGHGHFFFDD@?B HI:i:1 NH:i:1 NM:i:0
## ERR458493.606005 0 IX 27478 40 51M * 0 0 GGTAGTTCGATCATATGGCAATGCAGTACGATTTGTATATCGAGGCATATG ?+1BDDBDHDBFHGGFBAFHGIAFHGFHIIIGHHHDAG@HCHHABHIIEHG HI:i:1 NH:i:1 NM:i:0
## ERR458493.981978 0 IX 27482 40 51M * 0 0 GTTCGATCATATGGCAATGCAGTACGATTTGTATATCGAGGCATATGCGTG CCCFFFFFGHGHHJIJIIJJJJIJJJGIJJJHIJJJJGJJJIJJJJJJJHH HI:i:1 NH:i:1 NM:i:0
## ERR458493.705272 16 IX 27488 40 51M * 0 0 TCATATGGCAATGCAGTACGATTTGTATATCGAGGCATATGCGTGTTGTTT IJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJIJJJHHHHHFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.336141 0 IX 27491 20 51M * 0 0 TAGGGCAATGCAGTACGATTTGTATATCGAGGCATATGCGTGTTGTTTCTA ?++4ADDF?FFDHHEHI@BHGEFHIIIIIIGIIIIIIIFIAFGGGGGIIIG HI:i:1 NH:i:1 NM:i:1
## ERR458493.444372 0 IX 27511 40 51M * 0 0 TGTATATCGAGGCATATGCGTGTTGTTTCTACACACATCTAAAATATGTTT ;@1+4=A?C;F<FFIIGBGF?FDEIIIIEGIGIIIEHHIIIIFIBHIIBFG HI:i:1 NH:i:1 NM:i:0
## ERR458493.907930 0 IX 27535 40 51M * 0 0 GTTTCTACACACATCTAAAATATGTTTAGCTACTCCGTCGATAGAGTTGCG CBCFFFFFHHHHHJJJJJJJJJJJGIJJJIFIIJJIIIJJJJJJJJHIIJI HI:i:1 NH:i:1 NM:i:0
## ERR458493.861730 0 IX 27536 40 51M * 0 0 TTTCTACACACATCTAAAATATGTTTAGCTACTCCGTCGATAGAGTTGCGT CCCFFFFFHHHHHJJJGIIIIIJHJJIJIJJJJIJJHIIJFGHIIGHIIIF HI:i:1 NH:i:1 NM:i:0
## ERR458493.175066 16 IX 27579 40 51M * 0 0 AGTTGCGTATCATAAGTTATTTTGACAGCATCTCATTGTAATTTTACGTTA IIJIIGEIJJJIIJIIIGIIEGIIHIIIGGHHGGIIHGDFHFDDAFDF@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.488134 16 IX 27606 40 51M * 0 0 GCATCTCATTGTAATTTTACGTTACACCTTCTCTCTGGATGTAGTTCCAGG IHHIGIGIJJIHIGHGCDDBGGCCFAIHEGHAGGIIIJHGHHFBFDDF@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.828480 16 IX 27616 40 51M * 0 0 GTAATTTTACGTTACACCTTCTCTCTGGATGTAGTTCCAGGAGCTTGAGGC 9CFED?B90;?B1:?*C?:FC9FAEEA,@FFHFB>F>A2+8F8?<DDD@@? HI:i:1 NH:i:1 NM:i:0
## ERR458493.603548 16 IX 27616 40 51M * 0 0 GTAATTTTACGTTACACCTTCTCTCTGGATGTAGTTCCAGGAGCTTGAGGC >HGGIJIGGIIIGEFCIJJJJIGFIJIIIIGHEIGIIJFFHHHDBB4411= HI:i:1 NH:i:1 NM:i:0
## ERR458493.1067868 16 IX 27619 40 51M * 0 0 ATTTTACGTTACACCTTCTCTCTGGATGTAGTTCCAGGAGCTTGAGGCGAC IHCCGGGJIHFGDIJJJIGJIJJJJJIJJJJJIJIIJJHHHHHFFFDD:1B HI:i:1 NH:i:1 NM:i:0
## ERR458493.467943 16 IX 27673 40 51M * 0 0 ATCAACTAAAGAAGTTGGTCAGAAGAAGGATAATGCACCATCAAGTATTTG IGGGHIIJIHGIFIJGHHGGJJIJJIJIIJIJIGIHGJHHHHHFFFFFC@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.688649 16 IX 27682 20 51M * 0 0 AGAAGTTGGTCAGAAGAAGGATAATGCACCATCAAGTATTTGCTATAGACC JJJJJJJIIJIJJJJJJJJJJJJJJIHFJJIJJJJJJIHHHHHFFFDB1CB HI:i:1 NH:i:1 NM:i:1
## ERR458493.700938 16 IX 27684 20 51M * 0 0 AAGTTGGTCAGAAGAAGGATAATGCACCATCAAGTATTTGCTATAGTCTGG 99DCED?EFD<E?:@11+@EBA<:A3BFDE@CEE9AIED8DDDDDDD:+;< HI:i:1 NH:i:1 NM:i:1
## ERR458493.371281 0 IX 27731 40 51M * 0 0 CCGGTATGTCATCAAATGTCGAATGTTCCCCGTCAGGGTATTTTCTTCTAG @@<DDDFFFHHHHJIIJJHHJJJIIHIIIIJJHGIJJJ?FHEIJGJIJJFC HI:i:1 NH:i:1 NM:i:0
## ERR458493.847987 0 IX 27831 40 50M1S * 0 0 ATTACTTCTTACTTCCACTGAAGGCTAGAATAACACGCAGGACTGAAAGGG ;??DD8;DH>?BFB?<CFGEIDFIGBG>DGGHI9@GF:@FF@GCGHC3?D# HI:i:1 NH:i:1 NM:i:0
## ERR458493.684096 16 IX 27847 40 51M * 0 0 ACTGAAGGCTAGAATAACACGCAGGACTGAAAGGTGAAATGGCCTGTCTTG HJJJIJJJJJJJJJJIHIHJIIJJIIJJJJJJJJJJJJHHHHHFFFFFCCC HI:i:1 NH:i:1 NM:i:0
## ERR458493.594255 0 IX 27880 40 51M * 0 0 GTGAAATGGCCTGTCTTGTGCTCTCTTCATTTTATTTTTTATTAGCAGAAC @1:ADDFFHHHHHIIIIBCBHGHCEHHIGHHIIHGIIIIIIGIGIGIHGG# HI:i:1 NH:i:1 NM:i:0
## ERR458493.463187 0 IX 27883 20 51M * 0 0 AAACGGCCTGTCTTGTGCTCTCTTCATTTTATTTTTTATTAGCAGAACATC ?@@4=BDDFHHHHJIFHIGIIGJIHGIGIJIJJJJIJJIJIIIJJIJJGCI HI:i:1 NH:i:1 NM:i:1
## ERR458493.427935 0 IX 27894 40 4S47M * 0 0 CATCCTTGTGCTCTCTTCATTTTATTTTTTATTAGCAGAACATCCTGCTGT :1++4=DDFFDDFFGIFEHFIIIFHHIHIIFGIDFIFFFIIIIIEIIIIIF HI:i:1 NH:i:1 NM:i:0
## ERR458493.1037229 0 IX 27924 40 51M * 0 0 GCAGAACATCCTGCTGTGATGAAATTTTATGGTCAACTTCTTTTACGATGA CCCFFFFFHHHHHJJJGIIJJJJJJJJJIJJJHIJIJJJJJJJJJJJJJJJ HI:i:1 NH:i:1 NM:i:0
## ERR458493.469242 0 IX 27947 20 51M * 0 0 ATGTTATGGTCAACTTCTTTTACGATGACACCCTTTCTTGCACTTCCTTTA @B+ADEFFGHGDHJIJIJJJJIGIEHIJJJJJJJJJJJJJIJIGGIGIJJI HI:i:1 NH:i:1 NM:i:1
## ERR458493.200236 0 IX 27955 40 51M * 0 0 GTCAACTTCTTTTACGATGACACCCTTTCTTGCACTTCCTTTAGGACACTG :1:DADDDHHFHBHGE?EABF@>EGCDGBFH*9CF<CFFHCFC<F>?DAB< HI:i:1 NH:i:1 NM:i:0
## ERR458493.1035207 16 IX 27998 40 51M * 0 0 GGACACTGAATTGGGATGCCACACGTACAAATTTGTCTTCTATGTGCGGAA IIFBBIIHF?;EIHEEJIGGEJIFIGHGJIIIIIHGJIFFDGHFFFDD@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.1004377 16 IX 28003 20 51M * 0 0 CTGAATTGGGATGCCACACGTACAAATTTGTCTTCTATGTGCGGAAATAAG AB?*CFA@D?9D?1B:;?AE?1@GHFHADCCA4FA?A8DADA=,A=44:11 HI:i:1 NH:i:1 NM:i:1
## ERR458493.789981 0 IX 28128 40 51M * 0 0 AAAATTCTTCCTTATGCAGTGGAAACTTTTAATGACATGTGGTACACTTTT @?@DFFFFHFHGHJJJJJGHFGIIIJJJIIIJJIIJIJJEHGIHIJJIJJJ HI:i:1 NH:i:1 NM:i:0
## ERR458493.503366 0 IX 28195 20 26M1I24M * 0 0 GTTGACCATGCCTCTTCTTTTTTTTTTATGTCTTGCTTGTAACGCAGTTGT ?:+4ADDDAHGDFHGGHGIIHIGGII>GB89?B@@4=C38@@7@EBHHEH? HI:i:1 NH:i:1 NM:i:2
## ERR458493.4650 0 IX 28265 40 51M * 0 0 CGAGCTTCCACCTATATTGCTTTCAACACAAAGAACCTTTTGACGGTCTTG @@@FFFFFFHHHHIHHIFIGHIJJGHIIIGIJJIGIIFHIIIGIIJIEBHI HI:i:1 NH:i:1 NM:i:0
## ERR458493.347324 0 IX 28288 40 51M * 0 0 CAACACAAAGAACCTTTTGACGGTCTTGGGCTAGCTCCTATTTTTGGACCC CCCFFFFFHHHHHJJJJJJJJJJIJJJJJJJIJJJJJJJJJJJJJJIJJJI HI:i:1 NH:i:1 NM:i:0
## ERR458493.987486 16 IX 28291 40 51M * 0 0 CACAAAGAACCTTTTGACGGTCTTGGGCTAGCTCCTATTTTTGGACCCTTG GJJJJJIHDJJJJIGJIJJJJJJJJJIJJJJJIIIJJJHHHHHFDFFD?@B HI:i:1 NH:i:1 NM:i:0
## ERR458493.960602 0 IX 28305 40 51M * 0 0 TGACGGTCTTGGGCTAGCTCCTATTTTTGGACCCTTGCCCATCACAGAACT ?@@D?D82CFFDD@E;>CH>?E3AC@9??D391:DFB0B3DG;F8DF?F>B HI:i:1 NH:i:1 NM:i:0
## ERR458493.366806 0 IX 28349 40 51M * 0 0 CAGAACTATGAGAAAAGCAGTACGTCTTGACTTTAGTGCATGCTGTACATC CCCFFFFFHFGHHJJJGIJIHHIIHHIJIIIGIJGJDBFFIIJIJDGHGIJ HI:i:1 NH:i:1 NM:i:0
## ERR458493.605917 16 IX 28426 40 51M * 0 0 CCTTGCCACATAAAACTACAGATATCCCCGTTTTCTTGGTTTACGGTTGGC IIIIIIIIIIIIIIGIHFIIIHIHGIIIIIIIIIIIIIHHHHHDDDDD=<@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.1077252 0 IX 28633 20 51M * 0 0 GTTCGGGCGGAATCATATCCTCAAATTCTTGACATTGCATAGCGTAGTATA :??4:@?)@F???FGGEHDDDGHIDFFFEHFFB8.=C=FGGGGBEA8?B@7 HI:i:1 NH:i:1 NM:i:1
## ERR458493.19700 0 IX 29143 40 51M * 0 0 TTCAAATCGAGAGGAATCGGGCATTTGATAAGGAAGAGTAATGAAGAAGCG ???BBDDDFFDDDG@EGHIIE::?FHHBF<FCGA@BGGBFAFHG@HGGG>B HI:i:1 NH:i:1 NM:i:0
## ERR458493.202474 0 IX 29209 40 51M * 0 0 GGGAAGCTTGCTGTATTTTCTAGCTCTGGGGGAAATGCTGGTTTAGCAGCA CCCFFFFFHHHHHHJJJJJJJJJJJJJJJJJJHJJJJJJJJFHIIJJJJJJ HI:i:1 NH:i:1 NM:i:0
## ERR458493.408358 0 IX 29214 40 51M * 0 0 GCTTGCTGTATTTTCTAGCTCTGGGGGAAATGCTGGTTTAGCAGCAGCAAC @@@FFFFFHHHGHJJJJIFIJIIIIJIHIJJJIJJJIJJJJJIJJJIJJII HI:i:1 NH:i:1 NM:i:0
## ERR458493.94957 0 IX 29284 40 51M * 0 0 CTTAATTGCAGTGTAGTGGTTCCTAAAACTACAAAACCTAGAATGGTAAAG B114ADFFHHHHHIIJIJJJJJJJJJJJJJIJJJJJJJJJJJJJIIFIIIG HI:i:1 NH:i:1 NM:i:0
## ERR458493.605421 16 IX 29291 20 51M * 0 0 GCAGTGTAGTGGTTCCTAAAACTACAAAACCTAGAATGGTAAAGAAAATAC IJJJJJJJJJJIIGIGGIJIHJIHJIIJIJIIJJJJJIHHHGHFDD=411@ HI:i:1 NH:i:1 NM:i:1
## ERR458493.655728 0 IX 29309 40 51M * 0 0 AAACTACAAAACCTAGAATGGTAAAGAAAATTCAAAGTGCAGGAGCCAAAG @C;BDFFFGHHFDHEHGIGIIHIEGHIEFIJJEEBGEFGIEIG;DHFDDHH HI:i:1 NH:i:1 NM:i:0
## ERR458493.8733 16 IX 29342 40 51M * 0 0 AAAGTGCAGGAGCCAAAGTCATTATCCATGGTGATCATTGGGGGGAAGCAG IHGHGCEGIIGIGGIIFJJIIHGEGJJIIIIJIIIJJJHHGGHFFFDD<@B HI:i:1 NH:i:1 NM:i:0
## ERR458493.81631 0 IX 29414 20 51M * 0 0 GTATGGCGCAAGAAAGCCAACATGGTTCGAAGACACTATATGTGCACCCGT ?+:4ADDFGHHGHJJGHJJJJIJJIGHIJJJIGIIJJJIJJJGIIGJGHIJ HI:i:1 NH:i:1 NM:i:1
## ERR458493.618584 16 IX 29419 40 51M * 0 0 GCGCAAGAAAGCCAACATGGTTCGAAGACACTATATGTGCACCCGTTTGAT >BHD?FF?93EAGGDFD??)8?:>GCHFFCH>IGEEAGF@<DHA>DDB+1; HI:i:1 NH:i:1 NM:i:0
## ERR458493.230542 16 IX 29456 20 51M * 0 0 TGCACCCGTTTGATAACGAGACAATTTGGGAAGGTCATTCTACGATTGTAG GFF@IFIGHGFEEBDDIIGBAIHGGGHIFDEIHECIIIHHHDFDDDDB1+@ HI:i:1 NH:i:1 NM:i:1
## ERR458493.955355 0 IX 29592 20 51M * 0 0 ATATAGTGGCATAATTAAAGGCCTAGATAGGAATCAGCTTGCTGAAAAAAT =?14A?2=2A??FGFE@G?B<E@>3<CFH<ECBBFF>?9C;B@@?FFEEBE HI:i:1 NH:i:1 NM:i:1
## ERR458493.39522 0 IX 29619 40 51M * 0 0 TAGGAATCAGCTTGCTGAAAAAATTCCGGTCGTTGCTGTAGAAACTGCCGG :==+B?DD?FHHHIIHIBGGHIIAHHII>:D@?FCGGHDGG<?FDHIBF;8 HI:i:1 NH:i:1 NM:i:0
## ERR458493.1047607 0 IX 29621 40 51M * 0 0 GGAATCAGCTTGCTGAAAAAATTCCGGTCGTTGCTGTAGAAACTGCCGGTT @@@DD;DDFFFDAFFIGIIII;GH4AEE@D?DGIIIFCCFCCFIIIEEE5; HI:i:1 NH:i:1 NM:i:0
## ERR458493.54940 16 IX 29694 40 51M * 0 0 CAAAAAAGGTAGTCCAGTTACTCTTGAAAAATTGACAAGTGTTGCAACTTC IJJJJJJJJJJJIJIJIIHEJJJJJJJJJIJJJIHJJJHFHHHFFFFFC@B HI:i:1 NH:i:1 NM:i:0
## ERR458493.264118 16 IX 29722 40 51M * 0 0 AAATTGACAAGTGTTGCAACTTCTTTGGCCTCCCCATACATAGCATCATTC IJIGJJJJIIIDJGIHGF<EHHEIHHH?GHEGJJIHJIHHDHHFFDFF@@@ HI:i:1 NH:i:1 NM:i:0
## ERR458493.977810 0 IX 29754 40 51M * 0 0 CCCATACATAGCATCATTCGCGTTTGAGAGTTTTAACAAGTATGGATGTAA @@CFFFFFHDHHHJJIJJJIJJGHIJJJHIDGHIIGJIJIHJIJJIJJJJF HI:i:1 NH:i:1 NM:i:0
## ERR458493.640121 16 IX 29777 20 51M * 0 0 TTGAGAGTTTTAACAAGTATGGATGTAAGTCTGTAGTTTTATCAGATCAGG IIIIIIIGCIGHG@GHCIIIGGIIIIIIGIIIIIGBIGHGHFFEDB4411= HI:i:1 NH:i:1 NM:i:1
## ERR458493.116853 0 IX 29812 40 51M * 0 0 GTTTTATCAGATCAAGACGTTCTGGCAACATGCTTGAGATATGCCGATGAC CCCFFFFFHHHHHJJHIJJJJJIJJJJJJJJJJJJJJJIJJJJJJJJJJJJ HI:i:1 NH:i:1 NM:i:0
## ERR458493.778548 16 IX 29859 20 51M * 0 0 TGACTACAATTTTATAGTGGAACCAGCCTGTGGAGCATCCTTACATTTCTG DHFBEBEGGJIGFD@GF@HGGCIHGGIJIIJJJIIIJIHHHGHFFDB41BB HI:i:1 NH:i:1 NM:i:1
## ERR458493.1030927 16 IX 29887 40 51M * 0 0 TGTGGAGCATCCTTACATTTATGTTATCATCCAGAGATTCTTGAAGACATT 9B<???@GDFCE?E9EBC9@C<@DF>CFCIGI?>IIHF<DBFB;AD=A:@? HI:i:1 NH:i:1 NM:i:0
## ERR458493.521639 16 IX 29910 20 51M * 0 0 TTATCATCCAGAGATTCTTGAAGCCATTCTGGAACAAAAAATATATGAGGA GJJJJJJJIGIJJIIIIJJJIGE?IJJIIJJJJJJJJJHHHHHFFFFFCCC HI:i:1 NH:i:1 NM:i:1
## ERR458493.906587 16 IX 29942 40 51M * 0 0 AACAAAAAATATATGAGGATGATATCGTTATTATAATCGCATGCGGTGGAT IIJJJJJJJJJJJJJJJJJJJJJJGHHIJJJJJJJJJJGHHHHFFFFFCCC HI:i:1 NH:i:1 NM:i:0In addition, you can save the output as another BAM file if you want. You need -b when piping to another command.
source $HOME/miniconda3/bin/activate
conda activate variant
samtools view -b sorted.bam "IX:20000-30000" > sorted_chrIX_20000_30000.bam
samtools view sorted_chrIX_20000_30000.bam| head -n 1## ERR458493.829354 0 IX 19964 40 51M * 0 0 GTTGTAGGCATCCTAATCCAAATAACCTCCATAAACAAGTGGTACCAATAC ?@++=BDDHHHHHIIIIIIIIIIIIIIIIIIIIIIIIIIHIIHIIIIIIII HI:i:1 NH:i:1 NM:i:020.8 stats
mpileup The mpileup command produces a pileup format (or BCF) file giving, for each genomic coordinate, the overlapping read bases and indels at that position in the input BAM files(s). This can be used for SNP calling for example.
source $HOME/miniconda3/bin/activate
conda activate variant
samtools mpileup -f FASTQ/r64/Saccharomyces_cerevisiae.R64-1-1.dna.toplevel.fa -s sorted.bam | head## [E::fai_build3_core] Failed to open the file FASTQ/r64/Saccharomyces_cerevisiae.R64-1-1.dna.toplevel.faflagstat Counts the number of alignments for each FLAG type. Looking at the content of the file using samtools view:
samtools flagstat will indicate that some reads (about 10%) mapped to different chromosomes.
source $HOME/miniconda3/bin/activate
conda activate variant
samtools flagstat aln.bam## 1093957 + 0 in total (QC-passed reads + QC-failed reads)
## 1093957 + 0 primary
## 0 + 0 secondary
## 0 + 0 supplementary
## 0 + 0 duplicates
## 0 + 0 primary duplicates
## 1062112 + 0 mapped (97.09% : N/A)
## 1062112 + 0 primary mapped (97.09% : N/A)
## 0 + 0 paired in sequencing
## 0 + 0 read1
## 0 + 0 read2
## 0 + 0 properly paired (N/A : N/A)
## 0 + 0 with itself and mate mapped
## 0 + 0 singletons (N/A : N/A)
## 0 + 0 with mate mapped to a different chr
## 0 + 0 with mate mapped to a different chr (mapQ>=5)SAMtools will make the data readeble, this data is then piped through head to show the first 10 lines of the file.
20.9 BAM flags
The second column in a SAM/BAM file is the flag column; use the flags subcommand to understand specific flags.
They may seem confusing at first but the encoding allows details about a read to be stored by just using a few digits.
The trick is to convert the numerical digit into binary, and then use the table to interpret the binary numbers, where 1 = true and 0 = false.
source $HOME/miniconda3/bin/activate
conda activate variant
samtools flags## About: Convert between textual and numeric flag representation
## Usage: samtools flags FLAGS...
##
## Each FLAGS argument is either an INT (in decimal/hexadecimal/octal) representing
## a combination of the following numeric flag values, or a comma-separated string
## NAME,...,NAME representing a combination of the following flag names:
##
## 0x1 1 PAIRED paired-end / multiple-segment sequencing technology
## 0x2 2 PROPER_PAIR each segment properly aligned according to aligner
## 0x4 4 UNMAP segment unmapped
## 0x8 8 MUNMAP next segment in the template unmapped
## 0x10 16 REVERSE SEQ is reverse complemented
## 0x20 32 MREVERSE SEQ of next segment in template is rev.complemented
## 0x40 64 READ1 the first segment in the template
## 0x80 128 READ2 the last segment in the template
## 0x100 256 SECONDARY secondary alignment
## 0x200 512 QCFAIL not passing quality controls or other filters
## 0x400 1024 DUP PCR or optical duplicate
## 0x800 2048 SUPPLEMENTARY supplementary alignmentFind out about a 73 flag.
source $HOME/miniconda3/bin/activate
conda activate variant
samtools flags 73## 0x49 73 PAIRED,MUNMAP,READ120.10 Counting
SAMtools view can be used to filter the alignment based on characters like mapping quality, chromosome, orientation etc. When the -c option is added the filtered selection is counted.
Count the total number of records:
source $HOME/miniconda3/bin/activate
conda activate variant
samtools view -c sorted.bam## 1093957Count the reads that align to the reverse strand:
source $HOME/miniconda3/bin/activate
conda activate variant
samtools view -f 20 -c sorted.bam## 0Count the reads that do align:
20.11 Filtering unmapped reads
Use -F 4 to filter out unmapped reads.
source $HOME/miniconda3/bin/activate
conda activate variant
samtools view -F 4 -b sorted.bam > mapped.bamUse -f 4 to keep only unmapped reads.
source $HOME/miniconda3/bin/activate
conda activate variant
samtools view -f 4 -b sorted.bam > unmapped.bamHere -F is used to exclude reads that fit samflag 4, read unmapped. Everything else is included.
Question Sometimes you will see that this number of alignments is higher then the number of sequences in your fastq file. How can this be?
When a read multimaps (aligned to multiple positions in the genome), each of these positions is included as a separate alignment.
Count the reads that align to the forward strand:
source $HOME/miniconda3/bin/activate
conda activate variant
samtools view -F 20 -c sorted.bam## 517816Use -F 20 to exclude “read reverse strand” and “read unmapped”.
20.12 Converting a BAM file to a SAM file
A BAM file is just a SAM file but stored in binary format; you should always convert your SAM files into BAM format since they are smaller in size and are faster to manipulate.
Note: remember to use -h to ensure the SAM file contains the sequence header information.
Generally, I recommend storing only sorted BAM files as they use even less disk space and are faster to process.
source $HOME/miniconda3/bin/activate
conda activate variant
samtools view -h aln.bam > aln.samNotice that the SAM file is much larger than the BAM file.
Size of SAM file.
source $HOME/miniconda3/bin/activate
conda activate variant
ls -lh aln.sam## -rw-r--r-- 2 ggiaever staff 174M Dec 4 21:47 aln.samSize of BAM file.
source $HOME/miniconda3/bin/activate
conda activate variant
ls -lh aln.bam## -rw-r--r-- 2 ggiaever staff 52M Dec 4 21:47 aln.bamWe can use head to view a SAM file.
source $HOME/miniconda3/bin/activate
conda activate variant
head aln.sam## @HD VN:1.0 SO:coordinate
## @SQ SN:I LN:230218
## @SQ SN:II LN:813184
## @SQ SN:III LN:316620
## @SQ SN:IV LN:1531933
## @SQ SN:V LN:576874
## @SQ SN:VI LN:270161
## @SQ SN:VII LN:1090940
## @SQ SN:VIII LN:562643
## @SQ SN:IX LN:439888The lines starting with an “@” symbol contains the header information.
The @SQ tag is the reference sequence dictionary; SN refers to the reference sequence name and LN refers to the reference sequence length.
If you don’t see lines starting with the “@” symbol, the header information is probably missing.
You can generate this information again by running the command below, where ref.fa is the reference FASTA file used to map the reads.
20.13 Stats
For simple statistics use samtools flagstat.
source $HOME/miniconda3/bin/activate
conda activate variant
samtools flagstat sorted.bam## 1093957 + 0 in total (QC-passed reads + QC-failed reads)
## 1093957 + 0 primary
## 0 + 0 secondary
## 0 + 0 supplementary
## 0 + 0 duplicates
## 0 + 0 primary duplicates
## 1062112 + 0 mapped (97.09% : N/A)
## 1062112 + 0 primary mapped (97.09% : N/A)
## 0 + 0 paired in sequencing
## 0 + 0 read1
## 0 + 0 read2
## 0 + 0 properly paired (N/A : N/A)
## 0 + 0 with itself and mate mapped
## 0 + 0 singletons (N/A : N/A)
## 0 + 0 with mate mapped to a different chr
## 0 + 0 with mate mapped to a different chr (mapQ>=5)For more stats, use samtools stats.
source $HOME/miniconda3/bin/activate
conda activate variant
samtools stats sorted.bam | grep ^SN
# The columns are: GC%, unique sequence percentiles, 10th, 25th, 50th, 75th and 90th depth percentile
samtools stats sorted.bam | grep ^GCD | cut -f 2-## SN raw total sequences: 1093957 # excluding supplementary and secondary reads
## SN filtered sequences: 0
## SN sequences: 1093957
## SN is sorted: 1
## SN 1st fragments: 1093957
## SN last fragments: 0
## SN reads mapped: 1062112
## SN reads mapped and paired: 0 # paired-end technology bit set + both mates mapped
## SN reads unmapped: 31845
## SN reads properly paired: 0 # proper-pair bit set
## SN reads paired: 0 # paired-end technology bit set
## SN reads duplicated: 0 # PCR or optical duplicate bit set
## SN reads MQ0: 120777 # mapped and MQ=0
## SN reads QC failed: 0
## SN non-primary alignments: 0
## SN supplementary alignments: 0
## SN total length: 55791807 # ignores clipping
## SN total first fragment length: 55791807 # ignores clipping
## SN total last fragment length: 0 # ignores clipping
## SN bases mapped: 54167712 # ignores clipping
## SN bases mapped (cigar): 53842249 # more accurate
## SN bases trimmed: 0
## SN bases duplicated: 0
## SN mismatches: 254063 # from NM fields
## SN error rate: 4.718655e-03 # mismatches / bases mapped (cigar)
## SN average length: 51
## SN average first fragment length: 51
## SN average last fragment length: 0
## SN maximum length: 51
## SN maximum first fragment length: 51
## SN maximum last fragment length: 0
## SN average quality: 36.1
## SN insert size average: 0.0
## SN insert size standard deviation: 0.0
## SN inward oriented pairs: 0
## SN outward oriented pairs: 0
## SN pairs with other orientation: 0
## SN pairs on different chromosomes: 0
## SN percentage of properly paired reads (%): 0.0
## 0.0 0.328 0.000 0.000 0.000 0.000 0.000
## 1.2 0.493 0.008 0.008 0.008 0.008 0.008
## 32.0 0.657 0.005 0.005 0.005 0.005 0.005
## 36.0 0.821 0.502 0.502 0.502 0.502 0.502
## 37.0 0.985 2.604 2.604 2.604 2.604 2.604
## 38.0 1.806 0.201 0.240 0.523 4.900 7.558
## 39.0 6.897 0.676 0.778 1.255 1.831 2.437
## 40.0 18.227 0.643 1.170 1.614 2.480 3.907
## 41.0 38.259 1.109 1.443 2.201 3.283 5.073
## 42.0 58.621 1.476 2.035 2.732 4.126 6.226
## 43.0 73.235 1.734 2.439 3.307 5.577 8.132
## 44.0 84.729 1.273 2.573 4.251 6.293 8.256
## 45.0 91.461 2.073 3.323 4.738 7.496 16.540
## 46.0 95.567 0.965 2.593 7.964 15.410 37.320
## 47.0 97.537 0.281 1.232 6.635 30.123 69.642
## 48.0 99.015 0.490 1.466 3.815 9.868 15.542
## 49.0 99.672 2.950 3.999 7.723 19.704 23.506
## 50.0 100.000 10.962 10.962 21.675 32.388 32.388Use samtools idxstats to print stats on a BAM file; this requires an index file which is created by running samtools index.
The output of idxstats is a file with four tab-delimited columns:
- Reference name
- Sequence length of reference
- Number of mapped reads
- Number of unmapped reads
source $HOME/miniconda3/bin/activate
conda activate variant
samtools idxstats sorted.bam## I 230218 25070 0
## II 813184 66524 0
## III 316620 32445 0
## IV 1531933 129452 0
## V 576874 41845 0
## VI 270161 19243 0
## VII 1090940 114003 0
## VIII 562643 61620 0
## IX 439888 25180 0
## X 745751 68057 0
## XI 666816 68452 0
## XII 1078177 125507 0
## XIII 924431 70007 0
## XIV 784333 45479 0
## XV 1091291 87987 0
## XVI 948066 81236 0
## Mito 85779 5 0
## * 0 0 31845We can use this with awk to calculate:
The number of mapped reads by summing the third column.
source $HOME/miniconda3/bin/activate
conda activate variant
samtools idxstats sorted.bam | awk '{s+=$3} END {print s}'## 1062112The number of reads, which is the sum of mapped and unmapped reads.
source $HOME/miniconda3/bin/activate
conda activate variant
samtools idxstats sorted.bam | awk '{s+=$3+$4} END {print s}'
samtools view -F 4 -b sorted.bam > mapped.bam
samtools view -f 4 -b sorted.bam > unmapped.bam## 1093957This number can never be bigger then the number of reads in the fastq file, as all reads in the output give a single alignment.
20.14 samtools calmd/fillmd
The calmd or fillmd tool is useful for visualising mismatches and insertions in an alignment of a read to a reference genome. The -e argument changes identical bases between the read and reference into =.
Note the need for a second -b inn the second command:
source $HOME/miniconda3/bin/activate
conda activate variant
samtools view -b sorted.bam | samtools fillmd -b -e - FASTQ/r64/Saccharomyces_cerevisiae.R64-1-1.dna.toplevel.fa > fillmd.bam
samtools view fillmd.bam | head -n 1## [E::fai_build3_core] Failed to open the file FASTQ/r64/Saccharomyces_cerevisiae.R64-1-1.dna.toplevel.fa
## samtools calmd: Failed to open reference file 'FASTQ/r64/Saccharomyces_cerevisiae.R64-1-1.dna.toplevel.fa': No such file or directory20.15 Coverage
Coverage can mean the:
- average depth of each covered base
- percentage of bases covered
samtools depth and samtools mpileup can be used to indicate the depth of each covered base (and used to calculate the average depth. samtools coverage will provide both the average depth and percentage of bases covered per chromosome/reference sequence.
samtools depth will return three columns: reference, position, and coverage.
source $HOME/miniconda3/bin/activate
conda activate variant
samtools depth sorted.bam | head## I 1944 2
## I 1945 2
## I 1946 2
## I 1947 2
## I 1948 2
## I 1949 2
## I 1950 2
## I 1951 2
## I 1952 2
## I 1953 2
source $HOME/miniconda3/bin/activate
conda activate variant
samtools coverage sorted.bam## #rname startpos endpos numreads covbases coverage meandepth meanbaseq meanmapq
## I 1 230218 25070 111244 48.3212 5.52756 36 32.6
## II 1 813184 66524 447050 54.9753 4.14514 36.2 29.3
## III 1 316620 32445 157084 49.6128 5.197 36.1 35.1
## IV 1 1531933 129452 849956 55.4826 4.28245 36.2 30.5
## V 1 576874 41845 315652 54.7177 3.67724 36.2 33
## VI 1 270161 19243 135184 50.0383 3.61209 36.1 33.2
## VII 1 1090940 114003 616106 56.4748 5.29869 36.2 32.2
## VIII 1 562643 61620 296106 52.6277 5.5527 36.1 29.4
## IX 1 439888 25180 228498 51.9446 2.90019 36.1 32.7
## X 1 745751 68057 392098 52.5776 4.62022 36.1 28.3
## XI 1 666816 68452 363227 54.4718 5.20452 36.2 34.7
## XII 1 1078177 125507 568012 52.6826 5.90682 36.1 29.1
## XIII 1 924431 70007 519134 56.1571 3.83906 36.2 32.8
## XIV 1 784333 45479 437142 55.7342 2.93979 36.2 32.6
## XV 1 1091291 87987 589305 54.0007 4.08567 36.1 33.4
## XVI 1 948066 81236 502238 52.975 4.34222 36.2 25.3
## Mito 1 85779 5 209 0.243649 0.0029611 38.2 36If you want reads mapped to a single reference (e.g. chromosome), just specify the ref and leave out the start and end values.
source $HOME/miniconda3/bin/activate
conda activate variant
samtools view sorted.bam X | head -n 1## ERR458493.264522 0 X 396 0 51M * 0 0 AAAAAAATTTTCTTACAATATACTATACTACACAATACATAATCAGTGACT ???:BDDDDDDDDIEEIEIEA92??BEEEIECDDDECB?DEEDDD>?BBDE HI:i:1 NH:i:1 NM:i:0Returning to our coverage definition at the start of this section:
- average depth of each covered base =
meandepth - percentage of bases covered =
covbases
With SAMtools it is also possible to select alignments with a minimal mapping quality. Alignments with a maximal score (60 for hisat2 output files and 255 for STAR output files) are truly unique:
source $HOME/miniconda3/bin/activate
conda activate variant
samtools view -q 60 -c sorted.bam## 020.16
featureCounts -O -t gene -g ID -a ath_annotation.gff3 -o counts.txt mapped/*.bam Here is an explanation of the different arguments used:
-a
Key Points The SAM/BAM format is the end-result of a read alignment to a reference genome.
source $HOME/miniconda3/bin/activate
conda activate variant
samtools view -bT FASTQ/r64/Saccharomyces_cerevisiae.R64-1-1.dna.toplevel.fa.gz aln.sam > aln.bam## [E::fai_build3_core] Cannot index files compressed with gzip, please use bgzip
## [E::fai_path] Failed to build index file for reference file 'FASTQ/r64/Saccharomyces_cerevisiae.R64-1-1.dna.toplevel.fa.gz'Note that Rsamtools can do all these functions as well
20.17 Creating the counts file
For downstream applications e.g. differential expression analysis, the number of reads that maps within a gene has to be determined for each sample.
The featureCounts program from the Subread package can do this.
featureCounts can count the number of reads that map within a feature. The Saccharomyces genome annotation in the GTF format contain three different features to choose from.
Depending on the downstream applications the choice is gene, transcript or exon. In this study we are just looking for differentially expressed genes so our feature of interest specified by the -t will be gene.
source $HOME/miniconda3/bin/activate
conda activate variant
featureCounts -O -t gene -g gene_id -a FASTQ/r64/Saccharomyces_cerevisiae.R64-1-1.96.gtf.gz -o counts.txt BAM/*.bamHere is an explanation of the different arguments used:
-a
We can put all of the above into a shell script.