Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,989,036	G→A	D884D (GAC→GAT)	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,036	0	G	A	100.0%	181.8 / NA	62	D884D (GAC→GAT)	unknown	Exodeoxyribonuclease V gamma chain (EC 3.1.11.5)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/32);  total (30/32)

GAATTGCACAAAACGGCTTTGCGGCGTTTGGCTTACATAATCGGCGATCTTGCCTTTGAATTCCGCCACATCCGCCAGCAACTGATTGGCATAAATGGCGGAGAACTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCCGCTTCCACCAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTGTCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCAATTTCATTTTCCGGCTCCGCATCCGTCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTACCTGCGAATTTGCCAAAGGTAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCAGT  >  1/1988799‑1989259                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                gAATTGCAAAAAACGGCTTTGCGGCGTTTGGCTTACATAATCGGCGATATTGCCTTTAAATTCCGCCACATCCGCCAGCAACTGATTGGCATAAATGGCGGAGAACTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCttt                                                                                                                                                                                                                    >  1:407508/1‑251 (MQ=255)     tGCAAAAAACGGCTTTGCGGCGTTTGGCTTACATAATCGGCGATATTGCCTTTAAATTCCGCCACATCCGCCAGCAACTGATTGGCATAAATGGCGGAGAACTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTttcttc                                                                                                                                                                                                                <  1:643523/251‑1 (MQ=255)                           tttGGCTTACATAATCGGCGATATTGCCTTTAAATTCCGCCACATCCGCCAGCAACTGATTGGCATAAATGGCGGAGAACTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGcc                                                                                                                                                                                           <  1:824162/250‑1 (MQ=255)                              ggCTTACATAATCGGCGATATTGCCTTTAAATTCCGCCACATCCGCCAGCAACTGATTGGCATAAATGGCGGAGAACTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAgg                                                                                                                                                                                       >  1:77054/1‑251 (MQ=255)                                     aTAATCGGCGATATTGCCTTTAAATTCCGCCACATCCGCCAGCAACTGATTGGCATAAATGGCGGAGAACTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGttttt                                                                                                                                                                                <  1:493412/251‑1 (MQ=255)                                         tCGGCGATATTGCCTTTAAATTCCGCCACATCCGCCAGCAACTGATTGGCATAAATGGCGGAGAACTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCaaa                                                                                                                                                                            >  1:854633/1‑251 (MQ=255)                                            gCGATATTGCCTTTAAATTCCGCCACATCCGCCAGCAACTGATTGGCATAAATGGCGGAGAACTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTcc                                                                                                                                                                                                               <  1:520978/213‑1 (MQ=255)                                            gCGATATTGCCTTTAAATTCCGCCACATCCGCCAGCAACTGATTGGCATAAATGGCGGAGAACTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGaa                                                                                                                                                                         >  1:793173/1‑251 (MQ=255)                                                     ccTTTAAATTCCGCCACATCCGCCAGCAACTGATTGGCATAAATGGCGGAGAACTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAAt                                                                                                                                                                                                                         <  1:819985/194‑1 (MQ=255)                                                          aaaTTCCGCCACATCCGCCAGCAACTGATTGGCATAAATGGCGGAGAACTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTc                                                                                                                                                                                                          >  1:421111/2‑204 (MQ=255)                                                                       tCCGCCAGCAACTGATTGGCATAAATGGCGGAGAACTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGc                                                                                                                                              <  1:474783/251‑1 (MQ=255)                                                                               cAACTGATTGGCATAAATGGCGGAGAACTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTCTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGc                                                                                                                                      <  1:447882/251‑1 (MQ=255)                                                                                    gATTGGCATAAATGGCGGAGAACTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGaa                                                                                                                                 >  1:825363/1‑251 (MQ=255)                                                                                        ggCATAAATGGCGGAGAACTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGaaaa                                                                                                                               >  1:867202/1‑249 (MQ=255)                                                                                               aTGGCGGAGAACTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGAttt                                                                                                                       >  1:741047/1‑250 (MQ=255)                                                                                                    ggAGAACTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTcttctt                                                                                                                 >  1:143447/1‑246 (MQ=255)                                                                                                       gaACTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCg                                                                                                              >  1:269773/1‑251 (MQ=255)                                                                                                        aaCTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCg                                                                                                              <  1:197982/250‑1 (MQ=255)                                                                                                          cTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTc                                                                                                                                                                                                          >  1:476977/1‑156 (MQ=255)                                                                                                                 cgcgGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCAt                                                                                                     >  1:9434/1‑250 (MQ=255)                                                                                                                           cGCCTTTGACTGTCAGACGCGCAAAAAACGCGTCTGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTGCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACACGATTTCTTCTTCCGGCTCCGCATCCGCCatcatc                                                                                          <  1:876211/251‑1 (MQ=255)                                                                                                                            gCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCg                                                                                         >  1:721882/1‑251 (MQ=255)                                                                                                                                  gACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGAtt                                                                                                                        <  1:328339/214‑1 (MQ=255)                                                                                                                                   aCTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAAt                                                                                  >  1:800121/1‑251 (MQ=255)                                                                                                                                           aCGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAAt                                                                             <  1:889286/248‑1 (MQ=255)                                                                                                                                           aCGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATGTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTAACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAAt                                                                             <  1:86227/248‑1 (MQ=255)                                                                                                                                                 aaaaaaCGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAgg                                                                    <  1:420152/251‑1 (MQ=255)                                                                                                                                                   aaaaCGCGTCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGtt                                                                  <  1:123813/251‑1 (MQ=255)                                                                                                                                                           tCGGTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCaa                                                                                                                                        >  1:601721/1‑173 (MQ=255)                                                                                                                                                             ggTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTcttcttc                                                                                                                <  1:601988/195‑7 (MQ=255)                                                                                                                                                             ggTTTGCTGCTCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATaa                                                                                >  1:638378/1‑227 (MQ=255)                                                                                                                                                                   ctgctCGTCTGCTTCTAACAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAAt                                                  <  1:717573/251‑1 (MQ=255)                                                                                                                                                                             gcttctaacaAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATcc                                                                           <  1:117100/208‑1 (MQ=255)                                                                                                                                                                             gcttctaacaAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAgg                                        >  1:61372/9‑251 (MQ=255)                                                                                                                                                                              cttctaacaAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAgg                                        <  1:210145/243‑1 (MQ=255)                                                                                                                                                                              cttctaacaAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAgg                                        >  1:296274/8‑250 (MQ=255)                                                                                                                                                                              cttctaacaAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAgg                                        <  1:852350/243‑1 (MQ=255)                                                                                                                                                                                  taacaAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGc                                                                                                                                      >  1:720480/4‑152 (MQ=255)                                                                                                                                                                                  taacaAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAAt                                                  <  1:588551/233‑1 (MQ=255)                                                                                                                                                                                      aaaTCCGCATGAATGGTGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTg                               <  1:525848/251‑1 (MQ=255)                                                                                                                                                                                       aaTCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCaa                                                                      >  1:52211/1‑211 (MQ=255)                                                                                                                                                                                       aaTCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCTACTGt                              <  1:258587/251‑1 (MQ=255)                                                                                                                                                                                       aaTCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGt                              <  1:342022/251‑1 (MQ=255)                                                                                                                                                                                           cGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAggc                                       >  1:798453/1‑237 (MQ=255)                                                                                                                                                                                             cATGAATGGCGTATTTTTCCAAGGCACTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGc                                                                                                                                      >  1:591618/1‑141 (MQ=255)                                                                                                                                                                                                 tgacaagtgtttttCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAg                    <  1:829681/242‑1 (MQ=255)                                                                                                                                                                                                 aaTGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAg                    <  1:431620/251‑1 (MQ=255)                                                                                                                                                                                                 aaTGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAg                    <  1:191840/251‑1 (MQ=255)                                                                                                                                                                                                   tGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCTAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGAt                  <  1:440580/251‑1 (MQ=255)                                                                                                                                                                                                   tGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGAt                  <  1:521263/251‑1 (MQ=255)                                                                                                                                                                                                   tGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGAt                  <  1:623965/251‑1 (MQ=255)                                                                                                                                                                                                       gTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCCTCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCa              >  1:455530/1‑251 (MQ=255)                                                                                                                                                                                                       gTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAAt                                                                             >  1:229831/1‑188 (MQ=255)                                                                                                                                                                                                       gTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCTTGCAAGTTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTg                >  1:867636/1‑249 (MQ=255)                                                                                                                                                                                                        tATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCAt             <  1:687923/251‑1 (MQ=255)                                                                                                                                                                                                             ttCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGGTCAGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAg        <  1:383112/251‑1 (MQ=255)                                                                                                                                                                                                             ttCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAg        >  1:780869/1‑251 (MQ=255)                                                                                                                                                                                                             ttCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAg        <  1:853074/251‑1 (MQ=255)                                                                                                                                                                                                             ttCCAAGGCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTC‑AAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAGt       >  1:545110/1‑251 (MQ=255)                                                                                                                                                                                                                    gCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATaa                                                                                >  1:213294/1‑172 (MQ=255)                                                                                                                                                                                                                    gCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGcggt  >  1:396076/1‑247 (MQ=255)                                                                                                                                                                                                                    gCATTATCCAAGGCGAAATTTTCGCTATCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCGATTTCTTCTTCCGGCTCCGCATCCGCCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTGCCTGCGAATTTGCCAAAGGCAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGcggt  >  1:868613/1‑247 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                GAATTGCACAAAACGGCTTTGCGGCGTTTGGCTTACATAATCGGCGATCTTGCCTTTGAATTCCGCCACATCCGCCAGCAACTGATTGGCATAAATGGCGGAGAACTCCCCGCGCGGCAACACGCCTTTGACTTTCAGACGCGCAAAAAACGCGTCGGTTTGCTGCTCGTCCGCTTCCACCAAATCCGCATGAATGGCGTATTTTTCCAAGGCATTATCCAAGGCGAAATTTTCGCTGTCGGCAATGCTTTCTTCCTGTTCGGAGAAATTCACGCCAAGGCGTTTTTCAAAGAAAAATTTCACCGGATTTTTCACAAACGCCACCAAGCGGGAAAACTCAATTTCATTTTCCGGCTCCGCATCCGTCATCATCGGCTGAATAAAATCCTGCAAGGTTTGATTTACCTGCGAATTTGCCAAAGGTAGCCACTGTTGCGCGAAAGATTGCATAAAAGTGCAGT  >  1/1988799‑1989259

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