Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,989,414	C→T	E758E (GAG→GAA)	unknown ←	Exodeoxyribonuclease V gamma chain (EC 3.1.11.5)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,989,414	0	C	T	100.0%	223.8 / NA	78	E758E (GAG→GAA)	unknown	Exodeoxyribonuclease V gamma chain (EC 3.1.11.5)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (48/30);  total (48/30)

GCTGAATAAAATCCTGCAAGGTTTGATTTACCTGCGAATTTGCCAAAGGTAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCAGTCGTTTTTTTCAAAGTTTTTACGACTGAACACGGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGATGGGGTATTTTCCGCCGGCAGATTTTCCACGATGTAATCCAGCAACTGACTCACCAGCACCGACGGCTCCCGCTCCTGATTATCGGCAATAGAACGACCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCGCGACGGGCGCGATCGCCTTTTAAGCGGTGATATTGCATTAAATCGAAACTGTTCGGCGCTTGTTGGCGCGGGTATTCGCCGTCGTTCATACCGAGCAGACAAACCACTTTGAACGGAATGGAGCGCATCGGCAACAATGTACAAAAATT  >  1/1989173‑1989659                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      gCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTa                                                                                                                                                                                                                                              >  1:5541/1‑251 (MQ=255)    gAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCg                                                                                                                                                                                                                                           >  1:534587/1‑251 (MQ=255)     aaTAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCgg                                                                                                                                                                                                                                          >  1:570521/1‑251 (MQ=255)      aTAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCgg                                                                                                                                                                                                                                          >  1:579779/1‑250 (MQ=255)      aTAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCgg                                                                                                                                                                                                                                          >  1:27562/1‑250 (MQ=255)           aTCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGCATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGcgatgg                                                                                                                                                                                                                                    >  1:241492/1‑246 (MQ=255)                   aGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCg                                                                                                                                                                                                                            >  1:527006/1‑251 (MQ=255)                    ggTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGcgatggaa                                                                                                                                                                                                                                  >  1:666502/1‑237 (MQ=255)                    ggTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGa                                                                                                                                                                                                                           >  1:704603/1‑251 (MQ=255)                    ggTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGa                                                                                                                                                                                                                           >  1:847548/1‑251 (MQ=255)                          aTTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGc                                                                                                                                                                                                                     <  1:879774/251‑1 (MQ=255)                            ttgccTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTg                                                                                                                                                                                                                   <  1:40145/248‑1 (MQ=255)                               ccTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATg                                                                                                                                                                                                                <  1:41794/251‑1 (MQ=255)                                       tttGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAaata                                                                                                                                                                                                        <  1:155836/251‑1 (MQ=255)                                          gCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGAATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAaataat                                                                                                                                                                                                      >  1:898826/1‑250 (MQ=255)                                             aaaGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATc                                                                                                                                                                                                     >  1:631029/1‑248 (MQ=255)                                                ggcaGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTgcg                                                                                                                                                                                                >  1:478886/4‑250 (MQ=255)                                                  caGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCg                                                                                                                                                                                             >  1:98473/2‑251 (MQ=255)                                                   aGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGt                                                                                                                                                                                                                        <  1:532898/223‑1 (MQ=255)                                                     ccACTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACa                                                                                                                                                                                          >  1:231105/1‑251 (MQ=255)                                                        cTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACaaa                                                                                                                                                                                                          <  1:497541/232‑1 (MQ=255)                                                        cTGTTGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCa                                                                                                                                                                                       >  1:257666/1‑251 (MQ=255)                                                            tGCGCGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGc                                                                                                                                                                                   <  1:862253/251‑1 (MQ=255)                                                             gcgcGAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCTATGTTGTTGAATTAGTCGATCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGGCTCACCAGCACCGACGGCCCCCGTTCCTGATTATCGGCGATCGCATGTCCGATGTAGCTGATGTACAAATAATcc                                                                                                                                                                                                    >  1:485305/1‑233 (MQ=255)                                                              cgcgAAAGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATcc                                                                                                                                                                                                    >  1:848777/1‑232 (MQ=255)                                                                    aGATTGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAgg                                                                                                                                                                            >  1:249722/1‑250 (MQ=255)                                                                       ttGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGAtt                                                                                                                                                                                                                                               >  1:669061/1‑180 (MQ=255)                                                                       ttGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAg                                                                                                                                                                                         <  1:562584/234‑1 (MQ=255)                                                                       ttGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACa                                                                                                                                                                        >  1:109521/1‑251 (MQ=255)                                                                        tGCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAg                                                                                                                                                                       <  1:309185/251‑1 (MQ=255)                                                                         gCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGCACAGa                                                                                                                                                                      >  1:604735/1‑251 (MQ=255)                                                                         gCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGa                                                                                                                                                                      >  1:309648/1‑251 (MQ=255)                                                                         gCATAAAAGTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGa                                                                                                                                                                      >  1:876884/1‑251 (MQ=255)                                                                                 gTGCGGTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGAt                                                                                                                                                               >  1:300565/1‑250 (MQ=255)                                                                                  tgcggTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGAt                                                                                                                                                               <  1:791166/245‑1 (MQ=255)                                                                                  tgcggTCGTTTTTTTCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGAt                                                                                                                                                               <  1:183703/245‑1 (MQ=255)                                                                                            tttttCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGAc                                                                                                                                                   <  1:48330/251‑1 (MQ=255)                                                                                               ttCAAAGTTTTTGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACggg                                                                                                                                                <  1:448391/251‑1 (MQ=255)                                                                                                      tttttGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATc                                                                                                                                         <  1:470865/251‑1 (MQ=255)                                                                                                       ttttGCGACTGAACGCCGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCg                                                                                                                                        >  1:613964/1‑251 (MQ=255)                                                                                                                tgaacgccgTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTaaa                                                                                                                               <  1:328640/243‑2 (MQ=255)                                                                                                                 gaacgccgTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAAc                                                                                                                                                                         >  1:562185/8‑208 (MQ=255)                                                                                                                        gTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCGCTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAgg                                                                                                                                                                            <  1:662272/198‑1 (MQ=255)                                                                                                                        gTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGata                                                                                                                       >  1:727960/1‑251 (MQ=255)                                                                                                                         tCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGatat                                                                                                                      <  1:810107/251‑1 (MQ=255)                                                                                                                           aTGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGAt                                                                                                                                                               >  1:526049/1‑208 (MQ=255)                                                                                                                           aTGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTg                                                                                                                    >  1:512700/1‑251 (MQ=255)                                                                                                                              ccATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCAt                                                                                                                 <  1:532992/251‑1 (MQ=255)                                                                                                                                aTGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGgcg                                                                                                                                              >  1:10409/1‑220 (MQ=255)                                                                                                                                  gttgttGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTaaa                                                                                                             >  1:666731/1‑251 (MQ=255)                                                                                                                                   ttgttgAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAAt                                                                                                            <  1:685163/251‑1 (MQ=255)                                                                                                                                     gttgAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAgg                                                                                                                                                                            >  1:538290/1‑185 (MQ=255)                                                                                                                                        gAATCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGaaa                                                                                                       >  1:696121/1‑251 (MQ=255)                                                                                                                                         aaTCAGTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTg                                                                                                                                                                                                                   <  1:835099/142‑1 (MQ=255)                                                                                                                                              gTCGTTCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGatat                                                                                                                      >  1:73391/1‑230 (MQ=255)                                                                                                                                                   tCGCGCCACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGgcgc                                                                                            <  1:272560/251‑1 (MQ=255)                                                                                                                                                        ccACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGGCGCttgtt                                                                                       >  1:100251/1‑251 (MQ=255)                                                                                                                                                         cACTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATcc                                                                                                                                                                                                    <  1:743503/141‑1 (MQ=255)                                                                                                                                                          aCTTGAGATTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGGCGCttgttga                                                                                     >  1:242367/1‑250 (MQ=255)                                                                                                                                                               agaTTTCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCg                                                                                                                                                             <  1:799944/174‑1 (MQ=255)                                                                                                                                                                   ttCCCGCTGGGGTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGGCGCTTGTTGACGCGGGTAt                                                                            <  1:310672/251‑1 (MQ=255)                                                                                                                                                                           ggggTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGGCGCTTGTTGACGCGGGTATTCGCcgt                                                                     >  1:534461/1‑250 (MQ=255)                                                                                                                                                                           ggggTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGGCGCTTGTTGACGCGGGTATTCGCcgt                                                                     >  1:256208/1‑250 (MQ=255)                                                                                                                                                                            gggTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGGCGCTTGTTGACGCGGGTATTCGCcgtcg                                                                   >  1:158779/1‑251 (MQ=255)                                                                                                                                                                              gTATTTTCTGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCgg                                                                                                                                                                                                                                          >  1:767695/1‑82 (MQ=255)                                                                                                                                                                                   ttCCGCCGGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGGCGCTTGTTGACGCGGGTATTCGCCGTCGTTCATAc                                                            <  1:761603/251‑1 (MQ=255)                                                                                                                                                                                        ccgGCAGATTTTCCACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGGCGCTTGTTGACGCGGGTATTCGCCGTCGTTCATACCGAGc                                                       <  1:33677/251‑1 (MQ=255)                                                                                                                                                                                                    ccACGATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGGCGCTTGTTGACGCGGGTATTCGCCGTCGTTCATACCGAGCAGACAAActa                                             >  1:724864/1‑247 (MQ=255)                                                                                                                                                                                                        gATGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGGCGCTTGTTGACGCGGGTATTCGCCGTCGTTCATACCGAGCAGACAAACTACTCTGa                                       >  1:27081/1‑251 (MQ=255)                                                                                                                                                                                                         aTGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGGCGCTTGTTGACGCGGGTATTCGCCGTCGTTCATACCGAGCAGACAAACTACTCTGaa                                      >  1:372735/1‑251 (MQ=255)                                                                                                                                                                                                         aTGTAATCCAACAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGGCGCTTGTTGACGCGGGTATTCGCCGTCGTTCATACCGAGCAGACAAACTACTCTGaa                                      >  1:523683/1‑251 (MQ=255)                                                                                                                                                                                                                    cAACTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGGCGCTTGTTGACGCGGGTATTCGCCGTCGTTCATACCGAGCAGACAAACTACTCTGAACGGAATAGAgc                           <  1:891037/251‑1 (MQ=255)                                                                                                                                                                                                                       cTGACTCACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGGCGCttgttgacg                                                                                   <  1:315327/192‑4 (MQ=255)                                                                                                                                                                                                                             cACCAGCACCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGGCGCTTGTTGACGCGGGTATTCGCCGTCGTTCATACCGAGCAGACAAAc                                               >  1:529927/1‑222 (MQ=255)                                                                                                                                                                                                                                    aCCGACGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGgcgc                                                                                            <  1:779839/170‑1 (MQ=255)                                                                                                                                                                                                                                         cGGCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGGCGCTTGTTGACGCGGGTATTCGCCGTCGTTCATACCGAGCAGACAAACTACTCTGAACGGAATAGAGCGCATCGGCAGCAATGTACaaa      >  1:816794/1‑251 (MQ=255)                                                                                                                                                                                                                                           gCTCCCGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGGCGCTTGTTGACGCGGGTATTCGCCGTCGTTCATACCGAGCAGACAAACTACTCTGAACGGAATAGAGCGCATCGGCAGCAATGTACa        <  1:528612/247‑1 (MQ=255)                                                                                                                                                                                                                                              cccGTTCCTGATTATCGGCGATGGAACGTCCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCACGACGGGCGCGATCGCCTTTTAAACGATGATATTGCATTAAATCGAAACTGTTCGGCGCTTGTTGACGCGGGTATTCGCCGTCGTTCATACCGAGCAGACAAACTACTCTGAACGGAATAGAGCGCATCGGCAGCAATGTACAAAAAtt  >  1:671082/1‑250 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      GCTGAATAAAATCCTGCAAGGTTTGATTTACCTGCGAATTTGCCAAAGGTAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCAGTCGTTTTTTTCAAAGTTTTTACGACTGAACACGGTCATGCCATGTTGTTGAATCAGTCGTTCGCGCCACTTGAGATTTCCCGATGGGGTATTTTCCGCCGGCAGATTTTCCACGATGTAATCCAGCAACTGACTCACCAGCACCGACGGCTCCCGCTCCTGATTATCGGCAATAGAACGACCGATGTAGCTGATGTACAAATAATCCTGCGCCGACAGCAAGGCTTCAAGGAACAGATAACGATCGTCATCGCGACGGGCGCGATCGCCTTTTAAGCGGTGATATTGCATTAAATCGAAACTGTTCGGCGCTTGTTGGCGCGGGTATTCGCCGTCGTTCATACCGAGCAGACAAACCACTTTGAACGGAATGGAGCGCATCGGCAACAATGTACAAAAATT  >  1/1989173‑1989659

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