Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	570,274	G→A	P411S (CCG→TCG)	unknown ←	5‑methyltetrahydropteroyltriglutamate‑‑ homocysteine methyltransferase (EC 2.1.1.14)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	570,274	0	G	A	100.0%	213.2 / NA	76	P411S (CCG→TCG)	unknown	5‑methyltetrahydropteroyltriglutamate‑‑ homocysteine methyltransferase (EC 2.1.1.14)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/43);  total (33/43)

CCATGTCGTTACGTTCCGCTTCGCCGTGAACCAGCACGTCTAAATCCAATTTTTCTTGCTCACGCACCACCAATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGGGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGACGCGTAATTCTTGAATTTTTTGTAAACAGAATGCAAGCCAGCTGTAAAGTTTAGGTTTGTTTACTTGTAACTGTGTTTCCAC  >  1/570028‑570517                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    ccaTGTCGTTACGTTCCGCTTCGCCGTGAACCAGCACGTCTAAATCCAATTTTTCTTGCTCACGCACCACCAATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGAt                                                                                                                                                                                                                                                 <  1:421193/251‑1 (MQ=255)   aTGTCGTTACGTTCCGCTTCGCCGTGAACCAGCACGTCTAAATCCAATTTTTCTTGCTCACGCACCACCAATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGAttt                                                                                                                                                                                                                                               >  1:718250/1‑251 (MQ=255)     gTCGTTACGTTCCGCTTCGCCGTGAACCAGCACGTCTAAATCCAATTTTTCTTGCTCACGCACCACCAATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTa                                                                                                                                                                                                                                              >  1:636218/1‑250 (MQ=255)     gTCGTTACGTTCCGCTTCGCCGTGAACCAGCACGTCTAAATCCAATTTTTCTTGCTCACGCACCACCAATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTa                                                                                                                                                                                                                                              >  1:163362/1‑250 (MQ=255)            cGTTCCGCTTCGCCGTGAACCAGCACGTCTAAATCCAATTTTTCTTGCTCACGCACCACCAATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATc                                                                                                                                                                                                                                      <  1:539216/251‑1 (MQ=255)            cGTTCCGCTTCGCCGTGAACCAGCACGTCTAAATCCAATTTTTCTTGCTCACGCACCACCAATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATc                                                                                                                                                                                                                                      <  1:78408/251‑1 (MQ=255)                   cTTCGCCGTGAACCAGCACGTCTAAATCCAATTTTTCTTGCTCACGCACCACCAATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCt                                                                                                                                                                                                                               <  1:90675/251‑1 (MQ=255)                             aaCCAGCACGTCTAAATCCAATTTTTCTTGCTCACGCACCACCAATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAg                                                                                                                                                                                                                      <  1:178292/250‑1 (MQ=255)                             aaCCAGCACGTCTAAATCCAATTTTTCTTGCTCACGCACCACCAATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGt                                                                                                                                                                                                                     >  1:274660/1‑251 (MQ=255)                                     cGTCTAAATCCAATTTTTCTTGCTCACGCACCACCAATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCaaa                                                                                                                                                                                                             >  1:117463/1‑251 (MQ=255)                                     cGTCTAAATCCAATTTTTCTTGCTCACGCACCACCAATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCaaa                                                                                                                                                                                                             <  1:422788/251‑1 (MQ=255)                                       tCTAAATCCAATTTTTCTTGCTCACGCACCACCAATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACg                                                                                                                                                                                                           >  1:659676/1‑251 (MQ=255)                                            aTCCAATTTTTCTTGCTCACGCACCACCAATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGAt                                                                                                                                                                                                                                                 <  1:95654/208‑1 (MQ=255)                                                       cTTGCTCCCGCACCACCAATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAgag                                                                                                                                                                                           <  1:707907/251‑3 (MQ=255)                                                        ttGCTCACGCACCACCAATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGc                                                                                                                                                                                          <  1:413889/251‑1 (MQ=255)                                                         tGCTCACGCACCACCAATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGCTTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAggg                                                                                                                                                                                           >  1:700650/1‑249 (MQ=255)                                                             cACGCACCACCAATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGa                                                                                                                                                                                     <  1:195502/251‑1 (MQ=255)                                                                   ccaccaATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAggg                                                                                                                                                                                           >  1:547952/1‑239 (MQ=255)                                                                   ccaccaATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAggg                                                                                                                                                                                           >  1:503765/1‑239 (MQ=255)                                                                      ccaATTCGATTTCTTTTTGCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCtttt                                                                                                                                                                            <  1:823800/251‑1 (MQ=255)                                                                       caATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTg                                                                                                                                                                           <  1:821137/251‑1 (MQ=255)                                                                         aTTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAg                                                                                                                                                                         <  1:709518/251‑1 (MQ=255)                                                                                   ttttttCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTGTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACggg                                                                                                                                                               <  1:422993/251‑1 (MQ=255)                                                                                      tttCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCAt                                                                                                                                                            <  1:277819/251‑1 (MQ=255)                                                                                          aTTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATcggc                                                                                                                                                        <  1:136972/251‑1 (MQ=255)                                                                                                   tCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAatat                                                                                                                                                                                 <  1:288667/217‑1 (MQ=255)                                                                                                   tCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTtg                                                                                                                                               >  1:603333/1‑251 (MQ=255)                                                                                                        aTCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTtggctgg                                                                                                                                          >  1:791939/1‑251 (MQ=255)                                                                                                        aTCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTtggctgg                                                                                                                                          >  1:649663/1‑251 (MQ=255)                                                                                                                 ggACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGtt                                                                                                                                 >  1:322816/1‑251 (MQ=255)                                                                                                                 ggACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTCTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGtt                                                                                                                                 >  1:405156/1‑251 (MQ=255)                                                                                                                   aCAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTc                                                                                                                                <  1:426366/250‑1 (MQ=255)                                                                                                                    cAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGt                                                                                                                                                                       >  1:334821/1‑210 (MQ=255)                                                                                                                           ccTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTa                                                                                                                        >  1:153591/1‑250 (MQ=255)                                                                                                                           ccTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTAc                                                                                                                       <  1:569716/251‑1 (MQ=255)                                                                                                                              tttttGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCGTTTTGGCAAGTTCGCCAAAGGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTAcggc                                                                                                                    <  1:851475/251‑1 (MQ=255)                                                                                                                                      ggcggcGCGCGCGTGTCGAAGTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAg                                                                                                                                                                         <  1:455853/190‑1 (MQ=255)                                                                                                                                        cggcgCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTtggctgg                                                                                                                                          >  1:44537/1‑219 (MQ=255)                                                                                                                                                 cgTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGt                                                                                                 >  1:579113/1‑251 (MQ=255)                                                                                                                                                   tgtCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTg                                                                                                >  1:539238/1‑250 (MQ=255)                                                                                                                                                   tgtCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTg                                                                                                <  1:559170/250‑1 (MQ=255)                                                                                                                                                    gtCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGc                                                                                               <  1:47313/250‑1 (MQ=255)                                                                                                                                                        aaTTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGt                                                                                            >  1:167559/1‑249 (MQ=255)                                                                                                                                                        aaTTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCt                                                                                          <  1:431704/251‑1 (MQ=255)                                                                                                                                                          tttCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTa                                                                                        >  1:887176/1‑251 (MQ=255)                                                                                                                                                          tttCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAAGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTa                                                                                        >  1:453822/1‑251 (MQ=255)                                                                                                                                                                tagtTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGa                                                                                    <  1:164288/249‑1 (MQ=255)                                                                                                                                                                tagtTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGa                                                                                    <  1:563657/249‑1 (MQ=255)                                                                                                                                                                tagtTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGa                                                                                    <  1:741013/249‑1 (MQ=255)                                                                                                                                                                      gCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGGCTTAATGACGCGTAAt                                                                            >  1:131905/1‑251 (MQ=255)                                                                                                                                                                                          gTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGa                                                                                    <  1:640965/223‑1 (MQ=255)                                                                                                                                                                                          gTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACCGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGACGCGTAATTCTTGAATTTTTTGTAAACa                                                        >  1:202162/1‑251 (MQ=255)                                                                                                                                                                                                ggCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGa                                                                                                                                                                                     >  1:367377/1‑120 (MQ=255)                                                                                                                                                                                                ggCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATcggc                                                                                                                                                        >  1:410727/1‑149 (MQ=255)                                                                                                                                                                                                   aaTAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTAc                                                                                                                       <  1:685879/179‑1 (MQ=255)                                                                                                                                                                                                   aaTAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGACGCGTAATTCTTGAATTTTTTGTAAACAGAATGCAAg                                               <  1:118751/251‑1 (MQ=255)                                                                                                                                                                                                     tAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATg                                                                                     <  1:709866/211‑1 (MQ=255)                                                                                                                                                                                                     tAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGa                                                                                    <  1:740109/212‑1 (MQ=255)                                                                                                                                                                                                     tAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGACGCGTAATTCTTGAATTTTTTGTAAACAGAATGCAAGcc                                             <  1:858097/251‑1 (MQ=255)                                                                                                                                                                                                        cGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGa                                                                                    >  1:282394/1‑209 (MQ=255)                                                                                                                                                                                                        cGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGACGCGTAATTCTTGAATTTTTTGTAAACAGAATGCAAGCCAg                                           >  1:880608/1‑250 (MQ=255)                                                                                                                                                                                                              aTTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGACGCGTAATTCTTGAATTTTTTGTAAACAGAATGCAAGCCAGCTGTaaa                                    >  1:187002/1‑251 (MQ=255)                                                                                                                                                                                                               tttAACCAGGCATTTTGTAATTTGATACGTTCCGCAAGCGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCTGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGACGCGTAATTCTTGAATTTTTTGTAAACAGAATGCAAGCCAGCTGTAAAg                                   <  1:481818/251‑1 (MQ=255)                                                                                                                                                                                                               tttAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGACGCGTAATTCTTGAATTTTTTGTAAACAGAATGCAAGCCAGCTGTAAAg                                   <  1:461423/251‑1 (MQ=255)                                                                                                                                                                                                                       ggCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGACGCGTAATTCTTGAATTTTTTGTAAACAGAATGCAAGCCAGCTGTAAAGTTTAGgtt                           >  1:230127/1‑251 (MQ=255)                                                                                                                                                                                                                        gCATTTTGTAATTTGATGCGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTCAATGACGCGTAATTCTTGAATTTTTTGTAAACAGAATGCAAGCCAGCTGTAAAGTTTAGgtt                           <  1:517938/250‑1 (MQ=255)                                                                                                                                                                                                                        gCATTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGACGCGTAATTCTTGAATTTTTTGTAAACAGAATGCAAGCCAGCTGTAAAGTTTAGgttt                          >  1:662959/1‑251 (MQ=255)                                                                                                                                                                                                                          aTTTTGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGACGCGTAATTCTTGAATTTTTTGTAAACAGAATGCAAGCCAGCTGTAAAGTTTAGgtttgt                        <  1:711447/251‑1 (MQ=255)                                                                                                                                                                                                                           ttttGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTtggctgg                                                                                                                                          >  1:571966/1‑136 (MQ=255)                                                                                                                                                                                                                           ttttGTAATTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGACGCGTAATTCTTGAATTTTTTg                                                              <  1:871208/212‑1 (MQ=255)                                                                                                                                                                                                                                 aaTTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGc                                                                                               <  1:228895/173‑1 (MQ=255)                                                                                                                                                                                                                                  aTTTGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGACGCGTAATTCTTGAATTTTTTGTAAACAGAATGCAAGCCAGCTGTAAAGTTTAGGTTTGTTTACTTGt                >  1:356011/1‑251 (MQ=255)                                                                                                                                                                                                                                     tGATACGTTCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGACGCGTAATTCTTGAATTTTTTGTAAACAGAATGCAAGCCAGCTGTAAAGTTTAGGTTTGTTTACTTGTAAc             <  1:177559/251‑1 (MQ=255)                                                                                                                                                                                                                                            ttCCGCAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGACGCGTAATTCTTGAATTTTTTGTAAACAGAATGCAAGCCAGCTGTAAAGTTTAGGTTTGTTTACTTGTAACTGTGtt       <  1:380226/250‑1 (MQ=255)                                                                                                                                                                                                                                                 cAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGACGCGTAATTCTTGAATTTTTTGTAAACAGAATGCAAGCCAGCTGTAAAGTTTAGGTTTGTTTACTTGTAACTGTGTTTCCAc  <  1:231612/250‑1 (MQ=255)                                                                                                                                                                                                                                                 cAAACGAGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGACGCGTAATTCTTGAATTTTTTGTAAACAGAATGCAAGCCAGCTGTAAAGTTTAGGTTTGTTTACTTGTAACTGTGTTTCCAc  <  1:624371/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    CCATGTCGTTACGTTCCGCTTCGCCGTGAACCAGCACGTCTAAATCCAATTTTTCTTGCTCACGCACCACCAATTCGATTTCTTTTTTCATTGCCGCTTCATAATCCGCAAGGGACAATGTGCCTTTTTTGAAGGCGGCGCGCGCGTGTCGAATTTCAGTAGTTTGCGGGAAAGAACCGATATTCGTGGTTGGCAATAACGGTAGATTTAACCAGGCATTTTGTAATTTGATACGTTCCGCAAACGGGGATTTACGTTGATCCGCGCCTTTTGGCAAGTTCGCCAAACGTTCCGCCACTTCAGGGCGATGAATATCTTTTGAGGTCGCACGGGCATCGGCGGCAACTTGGCTGGCATCCAGTTCCGCTTGTACGGCGGCTCTGCCTTGCTCAAGTGCGGTCTTAATGACGCGTAATTCTTGAATTTTTTGTAAACAGAATGCAAGCCAGCTGTAAAGTTTAGGTTTGTTTACTTGTAACTGTGTTTCCAC  >  1/570028‑570517

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