Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	2,958,523	C→T	A428V (GCA→GTA)	PP_2589 →	Aldehyde dehydrogenase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	2,958,523	0	C	T	100.0%	232.1 / NA	65	A428V (GCA→GTA)	PP_2589	Aldehyde dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/33);  total (32/33)

GTAACCCGCTGGATCCGGCGAGCCGCGCCGGGGCCATTGTCGATGCCGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGCAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGT  >  AE015451/2958229‑2958811                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                  gTAACCCGCTGGATCCGGCGAGCCGCGCCGGGGCCATTGTCGATGCCGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTc                                                                                                                                                                                                                                                                                             >  4:437509/1‑301 (MQ=255)           gttcagacgtgtgctctgccgatctaTTGTCGATGCCGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGc                                                                                                                                                                                                                                                                                   <  4:401562/276‑1 (MQ=255)           ggATCCGGCGAGCCGCGCCGGGGCCATTGTCGATGCCGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGCACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAgg                                                                                                                                                                                                                                                                                   >  4:413976/1‑300 (MQ=255)                     ggCCGCGCCGGGGCCATTGTCGATGCCGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCa                                                                                                                                                                                                                                                                        <  3:249583/300‑1 (MQ=255)                                  ccATTGTCGATGCCGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCg                                                                                                                                                                                                                                                           >  4:538293/1‑301 (MQ=255)                                   cATTGTCGATGCCGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTTCACCgg                                                                                                                                                                                                                                                          >  3:490355/1‑301 (MQ=255)                                     tcttccgatctgggCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCAAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGt                                                                                                                                                                                                                                                        <  4:406079/290‑1 (MQ=255)                                                gggCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCAAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCCCCGGGTAgatcggaagg                                                                                                                                                                                                                                             >  2:405996/1‑292 (MQ=255)                                                 ggcagaacaccggcacggAGGCCGCCATGGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGtt                                                                                                                                                                                                                                            <  4:725712/284‑1 (MQ=255)                                                 ggCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGGAGCCGGTCGCCTGGCCACCGACAGCCTCTCTTGCCGGGTAGCCTCGGTGTGGTACGAACGACTACACCCAGGTGCCCCGCGTAGCCCggcgtg                                                                                                                                                                                                                                             >  2:725519/1‑297 (MQ=255)                                                    aGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACAGGGTCGCGCGGGCGTTGaa                                                                                                                                                                                                                                         >  1:265694/1‑301 (MQ=255)                                                             gCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCACGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTAc                                                                                                                                                                                                                                >  4:386520/1‑301 (MQ=255)                                                             caaTCGGGGCGGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTAc                                                                                                                                                                                                                                <  1:640108/299‑1 (MQ=255)                                                                  gAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTAt                                                                                                                                                                                                                           <  4:718722/301‑1 (MQ=255)                                                                              gtgcgtgcCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAg                                                                                                                                                                                                               <  4:47829/301‑1 (MQ=255)                                                                                    gcggggCAAGGGGGCAGGCGATGGGTGGGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATg                                                                                                                                                                                                         <  2:388689/298‑1 (MQ=255)                                                                                      cGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCAACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATgcg                                                                                                                                                                                                       >  3:613124/1‑301 (MQ=255)                                                                                          cAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGTATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑GCGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGGTTAAAGCCGGTTCCGTAACGGTGTacaccttccgcccgccgg                                                                                                                                                                                                   >  2:658087/1‑287 (MQ=255)                                                                                          cAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGAG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTTAAACCAGTCGATGCGCTgg                                                                                                                                                                                                   >  3:399531/1‑301 (MQ=255)                                                                                                    cgcgATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGt                                                                                                                                                                                         <  1:490341/301‑1 (MQ=255)                                                                                                       gATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGAGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATAGGTGAACACAGTCGATGCGCTGGATGTGAAGGTGcc                                                                                                                                                                                      >  1:420820/1‑301 (MQ=255)                                                                                                           ggTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGCCAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTc                                                                                                                                                                                  <  3:553260/301‑1 (MQ=255)                                                                                                           ggTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCggtc                                                                                                                                                                                  >  2:526867/1‑298 (MQ=255)                                                                                                                       ccGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAg                                                                                                                                                                      <  1:24605/301‑1 (MQ=255)                                                                                                                       ccGACGCCTGACCATCGAAGGTTCGGACAACGACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAg                                                                                                                                                                      <  1:31233/301‑1 (MQ=255)                                                                                                                        cGACGCCTGACCATCGAAGTTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGc                                                                                                                                                                     <  3:265699/301‑1 (MQ=255)                                                                                                                           cgccggaccagagAAGGTTCTGACAACGACATCGCGCCGACGATTTTTGCCGGTGTCGGGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGt                                                                                                                                                                  <  2:413893/289‑1 (MQ=255)                                                                                                                            gcttGACCGTCGAAGGTTCGGACAACTACATCGAGCCGACCAGTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGAACGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTc                                                                                                                                                                 <  2:386444/298‑1 (MQ=255)                                                                                                                                     tCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTgg                                                                                                                                                        <  4:527000/301‑1 (MQ=255)                                                                                                                                     gcGAAGGTTCGGAAAACTACATCGAGCCGACCATTTTTGCCGGTGTGGAGGGCCGAATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCAACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTgg                                                                                                                                                        <  1:613110/300‑1 (MQ=255)                                                                                                                                             tCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATc                                                                                                                                                <  3:552576/301‑1 (MQ=255)                                                                                                                                              cGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCt                                                                                                                                               >  1:581674/1‑301 (MQ=255)                                                                                                                                                aaCAACTGCATCTAGCCGACCATTTTGGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGTTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCTTGGCCAACGACAGCATCTATTGCCTGTTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCTTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGt                                                                                                                                             <  1:399521/300‑1 (MQ=255)                                                                                                                                                    aCTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACt                                                                                                                                         <  3:420831/301‑1 (MQ=255)                                                                                                                                                       aCATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCa                                                                                                                                      >  4:213317/1‑301 (MQ=255)                                                                                                                                                         aTCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCAtt                                                                                                                                    >  3:538421/1‑301 (MQ=255)                                                                                                                                                            gAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCAttcgt                                                                                                                                 >  4:363989/1‑301 (MQ=255)                                                                                                                                                                 gACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGACCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGCACACAGTCGATGCGCTGGCTGTGACGGTGCCGTTCGGTGGGGCCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGGCTTCGTTCGAc                                                                                                                            >  2:532756/1‑301 (MQ=255)                                                                                                                                                                    cATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAg                                                                                                                         <  1:538407/301‑1 (MQ=255)                                                                                                                                                                     aTTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTTCCGTTCGGGTGGGGCAAGCAGTCGGGGTTTGGCCGCGGTCTGTCACTGCAttcgttcgcccatg                                                                                                                        >  2:70290/1‑295 (MQ=255)                                                                                                                                                                       tttggCCGGGGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTAt                                                                                                                      <  2:213279/297‑1 (MQ=255)                                                                                                                                                                           gCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGc                                                                                                                  >  4:317829/1‑301 (MQ=255)                                                                                                                                                                              ggTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAAc                                                                                                               >  1:53318/1‑301 (MQ=255)                                                                                                                                                                                 gtCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCa                                                                                                            <  3:581689/301‑1 (MQ=255)                                                                                                                                                                                        ggcgcaggggCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGaccac                                                                                                     <  2:363916/292‑1 (MQ=255)                                                                                                                                                                                          ggaaTGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGAAACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCt                                                                                                   <  1:207750/298‑1 (MQ=255)                                                                                                                                                                                                    tGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGTAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGtt                                                                                         >  3:325129/1‑301 (MQ=255)                                                                                                                                                                                                           cgAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCggg                                                                                  >  1:276645/1‑300 (MQ=255)                                                                                                                                                                                                           cgAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCggg                                                                                  >  1:414451/1‑300 (MQ=255)                                                                                                                                                                                                           cgAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGc                                                                                  >  1:7425/1‑301 (MQ=255)                                                                                                                                                                                                                  gtgtTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTa                                                                           <  3:276650/301‑1 (MQ=255)                                                                                                                                                                                                                     ttCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTg                                                                        <  3:7425/301‑1 (MQ=255)                                                                                                                                                                                                                          cccGGTACTGGCGGTAAGCGCCGTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTAGGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCAt                                                                   <  2:538160/301‑1 (MQ=255)                                                                                                                                                                                                                                       gTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCgtcgtc                                                      >  4:489649/1‑301 (MQ=255)                                                                                                                                                                                                                                                 gcGACACGGAAGAAGAAGCCGTGCGCCTGGCGAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCgg                                            <  2:317765/300‑1 (MQ=255)                                                                                                                                                                                                                                                     cacGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTAACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACAACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCAt                                        >  1:139313/1‑301 (MQ=255)                                                                                                                                                                                                                                                        ggAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGcc                                     <  3:139319/301‑1 (MQ=255)                                                                                                                                                                                                                                                           agaagaagCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGc                                  >  4:563042/1‑301 (MQ=255)                                                                                                                                                                                                                                                                          ccTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGCCAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCCCGGGCCTGCCCGCTCCCCCAGGTTCacc                   >  2:266778/1‑299 (MQ=255)                                                                                                                                                                                                                                                                            tGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCCCGGGCATGCCCGCGCCCACATGTTCAgacc                 >  2:574113/1‑298 (MQ=255)                                                                                                                                                                                                                                                                                cAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGt             <  3:414461/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                   cGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCCCAAGTCtt          >  1:289633/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                    gACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGa         <  4:658257/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                     aCAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAg        <  3:289639/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                          aTCTATGGCCTGGTAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGT‑GGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtgg  >  4:540875/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                  GTAACCCGCTGGATCCGGCGAGCCGCGCCGGGGCCATTGTCGATGCCGGGCAGACCGGCCGCATCGAGGCCGCCATTGTGCGTGCCGGGCAAGAGGGCGCGCGATTGGTGTGTGGCGGCCGACGCCTGACCATCGAAGGTTCGGACAACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGCAGCCTCGGTGTGGAGCG‑ACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGT  >  AE015451/2958229‑2958811

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑