Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	248,853	T→C	R159R (CGT→CGC)	unknown →	Pyruvate formate‑lyase (EC 2.3.1.54)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	248,853	0	T	C	100.0%	191.0 / NA	70	R159R (CGT→CGC)	unknown	Pyruvate formate‑lyase (EC 2.3.1.54)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/45);  total (25/45)

TCTACAATCACTTCTCACAAACCGGGTTATATCAATAAAGAATTGGAAAAAATCGTCGGTTTACAAACCGATGCACCGTTAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGTTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGCCACCATTCAATTACGTGAAGAAATTGCCGAACAACATCGCGCATTGGGTA  >  1/248614‑249073                                                                                                                                                                                                                                                |                                                                                                                                                                                                                             tCTACAATCACTTCTCACAAACCGGGTTATATCAATAAAGAATTGGAAAAAATCGTCGGTTTACAAACCGATGCACCGTTAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAAt                                                                                                                                                                                                                    >  1:659450/1‑250 (MQ=255)  cTACAATCACTTCTCACAAACCGGGTTATATCAATAAAGAATTGGAAAAAATCGTCGGTTTACAAACCGATGCACCGTTAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATcc                                                                                                                                                                                                                  <  1:396779/251‑1 (MQ=255)     cAATCACTTCTCACAAACCGGGTTATATCAATAAAGAATTGGAAAAAATCGTCGGTTTACAAACCGATGCACCGTTAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCgg                                                                                                                                                                                                                >  1:705795/1‑250 (MQ=255)     cAATCACTTCTCACAAACCGGGTTATATCAATAAAGAATTGGAAAAAATCGTCGGTTTACAAACCGATGCACCGTTAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCgg                                                                                                                                                                                                                >  1:304037/1‑250 (MQ=255)           cTTCTCACAAACCGGGTTATATCAATAAAGAATTGGAAAAAATCGTCGGTTTACAAACCGATGCACCGTTAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTa                                                                                                                                                                                                         >  1:201176/1‑251 (MQ=255)              ctcACAAACCGGGTTATATCAATAAAGAATTGGAAAAAATCGTCGGTTTACAAACCGATGCACCGTTAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAAcc                                                                                                                                                                                                      <  1:672838/251‑1 (MQ=255)                        gggTTATATCAATAAAGAATTGGAAAAAATCGTCGGTTTACAAACCGATGCACCGTTAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATc                                                                                                                                                                                                                   >  1:599627/1‑228 (MQ=255)                          gTTATATCAATAAAGAATTGGAAAAAATCGTCGGTTTACAAACCGATGCACCGTTAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGAt                                                                                                                                                                                          >  1:425455/1‑251 (MQ=255)                                 cAATAAAGAATTGGAAAAAATCGTCGGTTTACAAACCGATGCACCGTTAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAAt                                                                                                                                                                                                                    <  1:119160/218‑1 (MQ=255)                                   aTAAAGAATTGGAAAAAATCGTCGGTTTACAAACCGATGCACCGTTAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACggt                                                                                                                                                                                 >  1:634855/1‑251 (MQ=255)                                          aTTGGAAAAAATCGTCGGTTTACAAACCGATGCACCGTTAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTTTTAACCGGTTTGCCGGATGCTTACggtcgtggtc                                                                                                                                                                          <  1:569048/251‑1 (MQ=255)                                               aaaaaaTCGTCGGTTTACAAACCGATGCACCGTTAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCAGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGtat                                                                                                                                                                      <  1:367268/250‑1 (MQ=255)                                                                  aaCCGATGCACCGTTAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTAtcgtcg                                                                                                                                                    <  1:1117/249‑1 (MQ=255)                                                                    ccGATGCACCGTTAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATc                                                                                                                                                                  <  1:606386/233‑1 (MQ=255)                                                                       aTGCACCGTTAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACggtcgtggt                                                                                                                                                                           <  1:631161/221‑1 (MQ=255)                                                                            ccGTTAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTa                                                                                                                                        >  1:344367/1‑251 (MQ=255)                                                                               ttAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATTCACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATgg                                                                                                                                     <  1:806980/251‑1 (MQ=255)                                                                               ttAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGGCATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATgg                                                                                                                                     <  1:883246/251‑1 (MQ=255)                                                                               ttAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATgg                                                                                                                                     <  1:549267/251‑1 (MQ=255)                                                                                                   gTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATc                                                                                                                                                                                                                   >  1:298940/1‑153 (MQ=255)                                                                                                        gtggtATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAg                                                                                                            >  1:113834/1‑251 (MQ=255)                                                                                                          ggtATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGaa                                                                                                          >  1:898918/1‑251 (MQ=255)                                                                                                            tATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGc                                                                                                                   <  1:707301/240‑1 (MQ=255)                                                                                                               cAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAAc                                                                                                     <  1:395763/251‑1 (MQ=255)                                                                                                                aaaaTGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTCCGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACa                                                                                                    >  1:445915/1‑251 (MQ=255)                                                                                                                    tGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTc                                                                                                <  1:575122/251‑1 (MQ=255)                                                                                                                        gAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGc                                                                                                                                              <  1:237031/201‑1 (MQ=255)                                                                                                                             gATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCttt                                                                                       >  1:425264/1‑250 (MQ=255)                                                                                                                              aTCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTAc                                                                                                                                                                                    <  1:715266/157‑1 (MQ=255)                                                                                                                                  tGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAAc                                                                                   >  1:791661/1‑250 (MQ=255)                                                                                                                                        gTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAAtt                                                                            >  1:594091/1‑251 (MQ=255)                                                                                                                                            aTGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGaa                                                                        <  1:622950/251‑1 (MQ=255)                                                                                                                                                  gtcAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAAc                                                                                                     <  1:562073/216‑1 (MQ=255)                                                                                                                                                   tcAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACgtggtg                                                                 <  1:116688/251‑1 (MQ=255)                                                                                                                                                      aTTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAAc                                                                                                     <  1:369605/212‑1 (MQ=255)                                                                                                                                                      aTTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAAc                                                                       <  1:727310/242‑1 (MQ=255)                                                                                                                                                           aTCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTAtcgtcg                                                                                                                                                    <  1:23448/160‑1 (MQ=255)                                                                                                                                                           aTCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACa                                                           <  1:876736/249‑1 (MQ=255)                                                                                                                                                           aTCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACAtt                                                         >  1:577665/1‑251 (MQ=255)                                                                                                                                                              cTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACaa                                                                                    >  1:830554/1‑221 (MQ=255)                                                                                                                                                                gAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAAcc                                                                                  <  1:755843/221‑1 (MQ=255)                                                                                                                                                                gAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAg                                                     >  1:90304/1‑250 (MQ=255)                                                                                                                                                                 aaGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATgg                                                                                                                                     <  1:777786/169‑1 (MQ=255)                                                                                                                                                                 aaGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGcc                                                   <  1:318825/251‑1 (MQ=255)                                                                                                                                                                 aaGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGc                                                    >  1:395849/1‑250 (MQ=255)                                                                                                                                                                   gTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAg                                                              <  1:99414/238‑1 (MQ=255)                                                                                                                                                                   gTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGcca                                                  >  1:324453/1‑250 (MQ=255)                                                                                                                                                                      gAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTAtcgtcg                                                                                                                                                    >  1:577163/1‑149 (MQ=255)                                                                                                                                                                       aaCATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGAc                                                            <  1:79805/236‑1 (MQ=255)                                                                                                                                                                           tatCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGa                                                             <  1:569283/231‑1 (MQ=255)                                                                                                                                                                           tatCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGCCACCATTCAAt                                         <  1:606025/251‑1 (MQ=255)                                                                                                                                                                            atCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCT‑AATTAGAACGTGGTGAAGACATTCAAGCCACCATTCAATTa                                       <  1:880327/251‑1 (MQ=255)                                                                                                                                                                             tCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGTCATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGCCACCATTCAATTa                                       <  1:229926/251‑1 (MQ=255)                                                                                                                                                                             tCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGCCACCATTCAATTa                                       <  1:666871/251‑1 (MQ=255)                                                                                                                                                                             tCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGCCACCATTCAATTa                                       <  1:48401/251‑1 (MQ=255)                                                                                                                                                                                          tCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGaa                                                                        <  1:179845/205‑1 (MQ=255)                                                                                                                                                                                          tCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGCCACCATTCAATTa                                       <  1:724379/238‑1 (MQ=255)                                                                                                                                                                                          tCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGCCACCATTCAATTa                                       <  1:790699/238‑1 (MQ=255)                                                                                                                                                                                            gTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGCCACCATTCAAt                                         <  1:386112/234‑1 (MQ=255)                                                                                                                                                                                                 aCCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGCCACCATTCAATTACGTGAAGAAATTGCCGaaca                   >  1:751988/1‑251 (MQ=255)                                                                                                                                                                                                  cccACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGCCACCATTCAATTACGTGAAGAAATTGCCGaacaa                  <  1:61593/251‑1 (MQ=255)                                                                                                                                                                                                     ccaACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGCCACCATTCAATTACGTGAAGAAATTGCCGAACAACAt               <  1:49715/250‑1 (MQ=255)                                                                                                                                                                                                     acaACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGCCACCATTCAATTACGTGAAGAAATTGCCGAACAACAt               <  1:652628/251‑1 (MQ=255)                                                                                                                                                                                                      caACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGAc                                                            <  1:380102/205‑1 (MQ=255)                                                                                                                                                                                                            aGGTGTTTTCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGCCACCATTCAAt                                         >  1:40787/1‑218 (MQ=255)                                                                                                                                                                                                                 ttttCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGCCACCATTCAATTACGTGAAGAAATTGCCGAACAACATCGCGCATTGGGt   <  1:167275/251‑1 (MQ=255)                                                                                                                                                                                                                  tttCGACGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGCCACCATTCAATTACGTGAAGAAATTGCCGAACAACATCGCGCATTGGGTa  >  1:589893/1‑251 (MQ=255)                                                                                                                                                                                                                        cGTTTATACACCGGACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGcc                                                   <  1:366426/196‑1 (MQ=255)                                                                                                                                                                                                                                     gACATTCTTCGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGCCACCATTCAATTa                                       >  1:537282/1‑195 (MQ=255)                                                                                                                                                                                                                                         ttcttcGCTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAg                                                              <  1:799519/168‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                             TCTACAATCACTTCTCACAAACCGGGTTATATCAATAAAGAATTGGAAAAAATCGTCGGTTTACAAACCGATGCACCGTTAAAACGTGCAATTATGCCGTTCGGTGGTATCAAAATGGTGAAAGGATCTTGCGAAGTTTATGGTCGTCAATTAGATCCTGAAGTAGAACATATCTTTACCGAATATCGTAAAACCCACAACCAAGGTGTTTTCGACGTTTATACACCGGACATTCTTCGTTGCCGTAAATCCGGCGTGTTAACCGGTTTGCCGGATGCTTACGGTCGTGGTCGTATTATCGGTGACTATCGTCGTTTAGCGGTTTATGGTGTCGATTTCTTGATGGCGGATAAGAAAAAACAATTCGCTTCCTTACAACCTAAATTAGAACGTGGTGAAGACATTCAAGCCACCATTCAATTACGTGAAGAAATTGCCGAACAACATCGCGCATTGGGTA  >  1/248614‑249073

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