Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	201,005	C→T	V556I (GTC→ATC)	unknown ←	Putative sulfate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	201,005	0	C	T	100.0%	263.3 / NA	75	V556I (GTC→ATC)	unknown	Putative sulfate permease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/40);  total (35/40)

CGATTATATTGAATGAGTCTATGAGCAGGCAAAAAATGTTAGGTTGTATTCTGATTTTAAGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGACTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAACACGGTCACGGCATCGCATTGTAAAACGATATGTTTGATTTCGTGATCCGTGTGCACGGTTTTATCGTGCAAATCGGCGAAGAGATTATCTGCCGCCGCGAAAAAG  >  1/200763‑201240                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                            cGATTATATTGAATGAGTCTATGAGCAGGCAAAAAATGTTAGGTTGTATTCTGATTTTAAGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCa                                                                                                                                                                                                                                     >  1:654720/1‑251 (MQ=255) cGATTATATTGAATGAGTCTATGAGCAGGCAAAAAATGTTAGGTTGTATTCTGATTTTAAGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCa                                                                                                                                                                                                                                     >  1:41957/1‑251 (MQ=255)   aTTATATTGAATGAGTCTATGAGCAGGCAAAAAATGTTAGGTTGTATTCTGATTTTAAGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAgg                                                                                                                                                                                                                                   >  1:585636/1‑251 (MQ=255)         ttGAATGAGTCTATGAGCAGGCAAAAAATGTTAGGTTGTATTCTGATTTTAAGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCg                                                                                                                                                                                                                             <  1:178722/251‑1 (MQ=255)                aGTCTATGAGCAGGCAAAAAATGTTAGGTTGTATTCTGATTTTAAGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGa                                                                                                                                                                                                                      >  1:313223/1‑251 (MQ=255)                 gTCTATGAGCAGGCAAAAAATGTTAGGTTGTATTCTGATTTTAAGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGa                                                                                                                                                                                                                      >  1:2189/1‑250 (MQ=255)                 gTCTATGAGCAGGCAAAAAATGTTAGGTTGTATTCTGATTTTAAGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGa                                                                                                                                                                                                                      >  1:746023/1‑250 (MQ=255)                                     tGTTAGGTTGTATTCTGATTTTAAGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGAc                                                                                                                                                                                                                     >  1:743400/1‑231 (MQ=255)                                      gTTAGGTTGTATTCTGATTTTAAGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGaa                                                                                                                                                                                                >  1:658983/1‑251 (MQ=255)                                       ttAGGTTGTATTCTGATTTTAAGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGaa                                                                                                                                                                                                <  1:883/250‑1 (MQ=255)                                            ttGTATTCTGATTTTAAGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCat                                                                                                                                                                                          <  1:211473/251‑1 (MQ=255)                                            ttGTATTCTGATTTTAAGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCat                                                                                                                                                                                          <  1:50126/251‑1 (MQ=255)                                            ttGTATTCTGATTTTAAGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCat                                                                                                                                                                                          <  1:116595/251‑1 (MQ=255)                                              gTATTCTGATTTTAAGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCGCCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGAc                                                                                                                                                                                                                     >  1:776607/1‑222 (MQ=255)                                                   cTGATTTTAAGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGcac                                                                                                                                                                                   <  1:61512/251‑1 (MQ=255)                                                     gATTTTAAGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACat                                                                                                                                                                                 >  1:588013/1‑251 (MQ=255)                                                           aaGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAg                                                                                                                                                                           <  1:475613/251‑1 (MQ=255)                                                           aaGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAg                                                                                                                                                                           <  1:130729/251‑1 (MQ=255)                                                                cACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAgttgt                                                                                                                                                                       <  1:418801/250‑1 (MQ=255)                                                                            attgtttTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCa                                                                                                                                                          <  1:313220/251‑1 (MQ=255)                                                                                 tttCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATAtgc                                                                                                                                                     <  1:27691/251‑1 (MQ=255)                                                                                 tttCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATAtgc                                                                                                                                                     <  1:409968/251‑1 (MQ=255)                                                                                 tttCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATAtgc                                                                                                                                                     <  1:796878/251‑1 (MQ=255)                                                                                 tttCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATAtgc                                                                                                                                                     <  1:160401/251‑1 (MQ=255)                                                                                 tttCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATAtgc                                                                                                                                                     <  1:684821/251‑1 (MQ=255)                                                                                      tGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGcac                                                                                                                                                <  1:36538/251‑1 (MQ=255)                                                                                          aTGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGc                                                                                                                                                                                     <  1:92340/210‑1 (MQ=255)                                                                                           tGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACTTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGc                                                                                                                                           <  1:451225/251‑1 (MQ=255)                                                                                           tGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGc                                                                                                                                           <  1:237136/251‑1 (MQ=255)                                                                                           tGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGc                                                                                                                                           <  1:323826/251‑1 (MQ=255)                                                                                              tAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCggg                                                                                                                                        <  1:589146/251‑1 (MQ=255)                                                                                                       gaaCAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATAtgc                                                                                                                                                     <  1:530535/229‑1 (MQ=255)                                                                                                       gaaCAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCATCCGCAGCGGTTGGAATTGCATATTGCACATATAAAg                                                                                                                                                                           >  1:450145/1‑207 (MQ=255)                                                                                                          cAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGa                                                                                                                            <  1:640150/251‑1 (MQ=255)                                                                                                                aaGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGcccc                                                                                                                      <  1:665997/251‑1 (MQ=255)                                                                                                                  ggTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAg                                                                                                                                                            >  1:521568/1‑211 (MQ=255)                                                                                                                  ggTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCg                                                                                                                    >  1:381507/1‑251 (MQ=255)                                                                                                                   gTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGa                                                                                                                            >  1:453408/1‑242 (MQ=255)                                                                                                                     aTGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCat                                                                                                                                                                                          <  1:833125/178‑1 (MQ=255)                                                                                                                      tGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTAt                                                                                                                <  1:254344/251‑1 (MQ=255)                                                                                                                       gTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGAt                                                                                                                                                                                                                          <  1:825105/144‑1 (MQ=255)                                                                                                                       gTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGt                                                                                                                  >  1:406147/1‑248 (MQ=255)                                                                                                                            attaAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCGAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATAtgc                                                                                                                                                     >  1:344226/1‑208 (MQ=255)                                                                                                                            attaAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACa                                                                                                                                                                                                                    <  1:542675/145‑1 (MQ=255)                                                                                                                            attaAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATAtgc                                                                                                                                                     <  1:479183/208‑1 (MQ=255)                                                                                                                            attaAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGc                                                                                                                                           <  1:669704/218‑1 (MQ=255)                                                                                                                            attaAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAaca                                                                                                          >  1:165155/1‑251 (MQ=255)                                                                                                                                      ctctCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAg                                                                                                                                                                           >  1:284371/1‑176 (MQ=255)                                                                                                                                      ctctCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTGGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGcc                                                                                                                        <  1:384036/227‑1 (MQ=255)                                                                                                                                          cGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCat                                                                                                                                                                                          >  1:29786/1‑157 (MQ=255)                                                                                                                                           gTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGcac                                                                                                                                                <  1:69289/198‑1 (MQ=255)                                                                                                                                           gTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGc                                                                                                                                           <  1:67167/203‑1 (MQ=255)                                                                                                                                                aaCCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGc                                                                                                                                                                                     <  1:200218/156‑1 (MQ=255)                                                                                                                                                      tCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGc                                                                                                                                                                                     <  1:834899/150‑1 (MQ=255)                                                                                                                                                       cATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAACACGGTCACGGCATCGCATTGTAAAACGa                                                                               >  1:271181/1‑251 (MQ=255)                                                                                                                                                        aTTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGAc                                                                                                                                                                                                                     >  1:142836/1‑116 (MQ=255)                                                                                                                                                        aTTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCat                                                                                                                                                         >  1:513828/1‑176 (MQ=255)                                                                                                                                                        aTTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAACACGGTCACGGCATCGCATTGTAAAACGat                                                                              >  1:638854/1‑251 (MQ=255)                                                                                                                                                        aTTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAACACGGTCACGGCATCGCATTGTAAAACGat                                                                              >  1:494812/1‑251 (MQ=255)                                                                                                                                                              ctctCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCatat                                                                                                                                                                                        <  1:362629/139‑1 (MQ=255)                                                                                                                                                              ctctCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAACACGGTCACGGCATCGCATTGTAAAACGATATGttt                                                                        >  1:446193/1‑251 (MQ=255)                                                                                                                                                                   aCTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAACACGGTCACGGCATCGCATTGTAAAACGATATGTTTGAttt                                                                   <  1:402268/251‑1 (MQ=255)                                                                                                                                                                      ctTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAACACGGTCACGGCATCGCATTGTAAAACGATATGTTTGATTTCGt                                                                >  1:782619/1‑251 (MQ=255)                                                                                                                                                                                   cccTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAACACGGTCACGGCATCGCATTGTAAAACGATATGTTTGATTTCGTGATCCGTGTGCAc                                                   <  1:39354/251‑1 (MQ=255)                                                                                                                                                                                    ccTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAACACGGTCACGGCATCGc                                                                                           >  1:708232/1‑210 (MQ=255)                                                                                                                                                                                    ccTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAACACGGTCACGGCATCGCATTGTAAAACGATATGTTTGATTTCGTGATCCGTGTGCACg                                                  >  1:372152/1‑251 (MQ=255)                                                                                                                                                                                           cTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAACACGGTCACGGCATCGCATTGTAAAACGATATGTTTGATTTCGTGATCCGTGTGCACGGTTTTAt                                           >  1:485457/1‑251 (MQ=255)                                                                                                                                                                                                  aTCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAACACGGTCACGGCATCGCATTGTaa                                                                                    <  1:630590/203‑1 (MQ=255)                                                                                                                                                                                                  aTCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAACACGGTCACGGCATCGCATTGTAAAACGATATGTTTGATTTCGTGATCCGTGTGCACGGTTTTAt                                           >  1:569904/1‑244 (MQ=255)                                                                                                                                                                                                  aTCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAACACGGTCACGGCATCGCATTGTAAAACGATATGTTTGATTTCGTGATCCGTGTGCACGGTTTTATCGTGCaa                                    >  1:157377/1‑251 (MQ=255)                                                                                                                                                                                                      tATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAACACGGTCACGGCATCGCATTGTAAAACGATATg                                                                           >  1:144675/1‑208 (MQ=255)                                                                                                                                                                                                              gCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAacac                                                                                                         >  1:110770/1‑170 (MQ=255)                                                                                                                                                                                                                   cGGCAAATCCGTGGCGAAGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGa                                                                                                                            >  1:188783/1‑146 (MQ=255)                                                                                                                                                                                                                                    aGTTGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAACACGGTCACGGCATCGCATTGTAAAACGATATGTTTGATTTCGTGATCCGTGTGCACGGTTTTATCGTGCAAATCGGCGAAGAGATTATCTGCCGCCGCGAAAAag  >  1:537844/1‑251 (MQ=255)                                                                                                                                                                                                                                      ttGATTTTTTTGATTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAACACGGTCACGGCATCGCATTGTAAAACGATATGTTTGATTTCGTGATCCGTGTGCACGGTTTTATCGTGCAAATCGGCGAAGAGATTATCTGCCGCCGCGAAAAag  <  1:94671/249‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                            CGATTATATTGAATGAGTCTATGAGCAGGCAAAAAATGTTAGGTTGTATTCTGATTTTAAGTGCACTTATTGTTTATTGTTTTCCTGGGATGTTAAAACGAAGAACAAATAAAGGTATGTTTTATTAAAAGCACTCTCGTTTTAACCAATCATTTAACTCTCACTCTTCACCTCGATTCCCTCAAACTCCCGAATCTTATTCAATGCTTCCGGCAAATCCGTGGCGAAGTTGATTTTTTTGACTTCCGGCAGGCTACCGGATTTGACAATCGTCCGCAGCGGTTGGAATTGCATATTGCACATATAAAGTTGTTGATGCGGAAGCATATGCTGCACAAAGCGGGTGAGGGCATGGATGCCCCCGGTATCTAACACGGTCACGGCATCGCATTGTAAAACGATATGTTTGATTTCGTGATCCGTGTGCACGGTTTTATCGTGCAAATCGGCGAAGAGATTATCTGCCGCCGCGAAAAAG  >  1/200763‑201240

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