Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	3	82,004	G→T	A477A (GCG→GCT)	unknown →	Recombinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	3	82,004	0	G	T	100.0%	177.8 / NA	66	A477A (GCG→GCT)	unknown	Recombinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (26/40);  total (26/40)

TAGGTGGACGTGTATCTATCGCTCCTAATGATCTTCTAAAAATCAATGAAGAATTTACCGCCTCTCTGGTTGTAAGCCCTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCGGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGATCTTGTAGAAAATATTGCGTCTGTTGGTTTAAAACGACCAATTACAGTCTCAAAGACGTCTGATGGTTATAAC  >  3/81760‑82245                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  tAGGTGGACGTGTATCTATCGCTCCTAATGATCTTCTAAAAATCAATGAAGAATTTACCGCCTCTCTGGTTGTAAGCCCTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGc                                                                                                                                                                                                                                             <  1:234149/251‑1 (MQ=255)  aGGTGGACGTGTATCTATCGCTCCTAATGATCTTCTAAAAATCAATGAAGAATTTACCGCCTCTCTGGTTGTAAGCCCTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCa                                                                                                                                                                                                                                            >  1:399102/1‑251 (MQ=255)       gACGTGTATCTATCGCTCCTAATGATCTTCTAAAAATCAATGAAGAATTTACCGCCTCTCTGGTTGTAAGCCCTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCg                                                                                                                                                                                                                                       >  1:527444/1‑251 (MQ=255)            gtATCTATCGCTCCTAATGATCTTCTAAAAATCAATGAAGAATTTACCGCCTCTCTGGTTGTAAGCCCTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGa                                                                                                                                                                                                                                  >  1:787091/1‑251 (MQ=255)                      cTCCTAATGATCTTCTAAAAATCAATGAAGAATTTACCGCCTCTCTGGTTGTAAGCCCTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATAtt                                                                                                                                                                                                                        >  1:781800/1‑251 (MQ=255)                        ccTAATGATCTTCTAAAAATCAATGAAGAATTTACCGCCTCTCTGGTTGTAAGCCCTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTa                                                                                                                                                                                                                      >  1:523055/1‑251 (MQ=255)                                                    aaTTTACCGCCTCTCTGGTTGTAAGCCCTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTttcttc                                                                                                                                                                                          <  1:348929/251‑1 (MQ=255)                                                    aaTTTACCGCCTCTCTGGTTGTAAGCCCTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTttcttc                                                                                                                                                                                          <  1:683709/251‑1 (MQ=255)                                                         aCCGCCTCTCTGGTTGTAAGCCCTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTAt                                                                                                                                                                                     <  1:103186/251‑1 (MQ=255)                                                              ctctctGGTTGTAAGCCCTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAg                                                                                                                                                                                <  1:210039/251‑1 (MQ=255)                                                                  ctGGTTGTAAGCCCTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACg                                                                                                                                                                            <  1:471317/251‑1 (MQ=255)                                                                  ctGGTTGTAAGCCCTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACg                                                                                                                                                                            <  1:361988/251‑1 (MQ=255)                                                                      ttGTAAGCCCTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATc                                                                                                                                                                        <  1:370329/251‑1 (MQ=255)                                                                           aGCCCTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCt                                                                                                                                                                       <  1:649336/247‑1 (MQ=255)                                                                            gCCCTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCttatta                                                                                                                                                                  >  1:362826/1‑251 (MQ=255)                                                                              ccTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAgg                                                                                                                                                                <  1:211974/251‑1 (MQ=255)                                                                              ccTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAgg                                                                                                                                                                <  1:488045/251‑1 (MQ=255)                                                                              ccTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAgg                                                                                                                                                                <  1:520481/251‑1 (MQ=255)                                                                                           cTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGt                                                                                                                                                                           >  1:648857/1‑227 (MQ=255)                                                                                           cTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAgg                                                                                                                                                                >  1:622798/1‑238 (MQ=255)                                                                                           cTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTAt                                                                                                                                                      <  1:338378/248‑1 (MQ=255)                                                                                           cTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGa                                                                                                                                                   >  1:463523/1‑251 (MQ=255)                                                                                                  gggTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGTAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGatat                                                                                                                                                                                                                         <  1:307966/174‑1 (MQ=255)                                                                                                  gggTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACg                                                                                                                                                                            >  1:586564/1‑219 (MQ=255)                                                                                                  gggTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAg                                                                                                                                                              >  1:804075/1‑233 (MQ=255)                                                                                                  gggTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTAt                                                                                                                                                      >  1:561026/1‑241 (MQ=255)                                                                                                            cGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGa                                                                                                                                  <  1:711733/251‑1 (MQ=255)                                                                                                              cTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATa                                                                                                                                >  1:305576/1‑251 (MQ=255)                                                                                                                          ttttGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCgg                                                                                                                    >  1:106234/1‑251 (MQ=255)                                                                                                                                tACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAg                                                                                                                                                                                >  1:614363/1‑185 (MQ=255)                                                                                                                                        ttGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAg                                                                                                      <  1:782251/251‑1 (MQ=255)                                                                                                                                           aaCCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGaa                                                                                                    <  1:116322/250‑1 (MQ=255)                                                                                                                                           aaCCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGaa                                                                                                    <  1:229972/250‑1 (MQ=255)                                                                                                                                                gATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGaa                                                                                                                                                                              >  1:17189/1‑171 (MQ=255)                                                                                                                                                     aaCTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTAt                                                                                                                                                                                     <  1:625784/159‑1 (MQ=255)                                                                                                                                                         aTTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGaaaaa                                                                                     <  1:747728/251‑1 (MQ=255)                                                                                                                                                                   gTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGa                                                                           >  1:777025/1‑251 (MQ=255)                                                                                                                                                                    tCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGAt                                                                          >  1:361201/1‑251 (MQ=255)                                                                                                                                                                            aaGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGaa                                                                                                                                                                              <  1:75989/143‑1 (MQ=255)                                                                                                                                                                                  cAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGAt                                                                                                                                                                                                   <  1:364970/116‑1 (MQ=255)                                                                                                                                                                                  cAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGaa                                                                                                    <  1:4815/211‑1 (MQ=255)                                                                                                                                                                                        aCATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGATCTTGTAg                                                                   <  1:392912/238‑1 (MQ=255)                                                                                                                                                                                          attattGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGaa                                                                                                    <  1:490128/203‑1 (MQ=255)                                                                                                                                                                                             attGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGATCTTGTAGAAAATATTGCGTCTGTTg                                                 >  1:339315/1‑251 (MQ=255)                                                                                                                                                                                                    attattTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCCGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAgg                                                                                                                                                                <  1:337958/133‑1 (MQ=255)                                                                                                                                                                                                     ttatttatTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGATCTTGTAGAAAATATTGCGTCTGTTGGTTTaaaa                                         <  1:363661/251‑1 (MQ=255)                                                                                                                                                                                                       atttatttCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGaa                                                                                                                                                                              <  1:160164/116‑1 (MQ=255)                                                                                                                                                                                                       atttatttCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCgg                                                                                                                    <  1:497947/174‑1 (MQ=255)                                                                                                                                                                                                       atttatttCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGaa                                                                                                    >  1:104689/1‑190 (MQ=255)                                                                                                                                                                                                            tttCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATcc                                                                                                         >  1:279995/1‑180 (MQ=255)                                                                                                                                                                                                            tttCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGATCTTGTAGAAAATATTGCGt                                                       <  1:676005/230‑1 (MQ=255)                                                                                                                                                                                                               ccAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAAt                                                                                    <  1:201752/198‑1 (MQ=255)                                                                                                                                                                                                                cAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGATCTTGTAGAAAATATTGCGTCTGTTGGTTTAAAACGACCAATTAc                              <  1:7673/251‑1 (MQ=255)                                                                                                                                                                                                                    gAAGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGa                                                                                                                                  >  1:450039/1‑147 (MQ=255)                                                                                                                                                                                                                      aGGGGAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGATCTTGTAGAAAATATTGCGTCTGTTGGTTTAAAACGACCAATTAc                              <  1:704795/245‑1 (MQ=255)                                                                                                                                                                                                                       ggggAGAGAAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGATCTTGTAGAAAATATTGCGTCTGTTGGTTTAAAACGACCAATTACAGTCTCa                       >  1:574296/1‑251 (MQ=255)                                                                                                                                                                                                                           agagaAAGTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGATCTTGTAGAAAATATTGCGTCTGTTGGTTTAAAACGACCAATTACAGTCTCAAAg                    <  1:39861/250‑1 (MQ=255)                                                                                                                                                                                                                                  gTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAgg                                                                                                                                                                <  1:569766/103‑1 (MQ=255)                                                                                                                                                                                                                                  gTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGATCTTGTAGAAAATATTGCGTCTGTTGGTTTAAAACGACCAATTACAGTCTCAAAGACGTCTGa            >  1:408949/1‑251 (MQ=255)                                                                                                                                                                                                                                  gTAATAAATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGATCTTGTAGAAAATATTGCGTCTGTTGGTTTAAAACGACCAATTACAGTCTCAAAGACGTCTGa            >  1:68202/1‑251 (MQ=255)                                                                                                                                                                                                                                    aataaATTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGATCTTGTAGAAAATATTGCGTCTGTTGGTTTAAAACGACCAATTAc                              <  1:228178/231‑1 (MQ=255)                                                                                                                                                                                                                                       aaaTTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGATCTTGTAGAAAATATTGCGTCTGTTGGTTTAAAACGACCAATTACAGtc                          <  1:158484/232‑1 (MQ=255)                                                                                                                                                                                                                                       aaaTTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGATCTTGTAGAAAATATTGCGTCTGTTGGTTTAAAACGACCAATTACAGTCTCAAAGAc                  <  1:628949/240‑1 (MQ=255)                                                                                                                                                                                                                                        aaTTGAGATTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGCTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGATCTTGCAGAAAATATTGCGTCTGTTGGTTTAAAACGACCAATTACAGTCTCAAAGACGACTGATGGTta      >  1:78919/1‑251 (MQ=255)                                                                                                                                                                                                                                             agaTTGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGATCTTGTAGAAAATATTGCGTCTGTTGGTTTAAAACGACCAATTACAGTCTCAAAGACGTCTGATGGTTATAAc  <  1:79100/250‑1 (MQ=255)                                                                                                                                                                                                                                                 tGGCTGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGATCTTGTAGAAAATATTGCGTCTGTTGGTTTAAAACGACCAATTACAGTCTCAAAGACGTc               <  1:869675/233‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  TAGGTGGACGTGTATCTATCGCTCCTAATGATCTTCTAAAAATCAATGAAGAATTTACCGCCTCTCTGGTTGTAAGCCCTTGCCGACCAACTACAAGGGGTAATAAACGCTGGATGATTCGTTTTGATACCGGGCTTGAACCTGATTTAACTATTGCAGTTCGTCTTGACGAAGATGCAGAAAACATTATTGATTATTATTTATTTCCAATGAAGGGGAGAGAAAGTAATAAATTGAGATTGGCGGAAAGCAACCCGGCTGAACTTGATATTTATCGATTTTCTAACTTGGATCGTTTCTTCATTATGGTAGAACGTATCTTATTAGGAGTTCATTATGGATATTCAAACTATTCCGATAGAAGATATTCGGATTATCAATCCAAGAACGCGTAATCGAAAAATCTTCGCAGATCTTGTAGAAAATATTGCGTCTGTTGGTTTAAAACGACCAATTACAGTCTCAAAGACGTCTGATGGTTATAAC  >  3/81760‑82245

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: ‑