Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,006,782	T→C	S67G (AGC→GGC)	unknown ←	Ribulosamine/erythrulosamine 3‑kinase potentially involved in protein deglycation

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,006,782	0	T	C	100.0%	195.9 / NA	68	S67G (AGC→GGC)	unknown	Ribulosamine/erythrulosamine 3‑kinase potentially involved in protein deglycation
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/38);  total (30/38)

GAAAAACTTCGCCCAGTGGTTGCTCCATTCGTTCGGTTGATATTGCGGACCAAGCCAAGTATCAAAAGACAAACCGTAGTTTTCTGAGCCGTGATATTGGTGCAAACGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCTTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATACCATGTTCCTTTTTTTCAGTTAATCATGTTTATAGTGTAAAAGAT  >  1/1006537‑1007026                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     gAAAAACTTCGCCCAGTGGTTGCTCCATTCGTTCGGTTGATATTGCGGACCAAGCCAAGTATCAAAAGACAAACCGTAGTTTTCTGAGCCGTGATATTGGTGCAAACGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTaaaa                                                                                                                                                                                                                                                 >  1:292303/1‑251 (MQ=255)      aCTTCGCCCAGTGGTTGCTCCATTCGTTCGGTTGATATTGCGGACCAAGCCAAGTATCAAAAGACAAACCGTAGTTTTCTGAGCCGTGATATTGGTGCAAACGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGAt                                                                                                                                                                                                                                             >  1:321032/1‑250 (MQ=255)                gTGGTTGCTCCATTCGTTCGGTTGATATTGCGGACCAAGCCAAGTATCAAAAGACAAACCGTAGTTTTCTGAGCCGTGATATTGGTGCAAACGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCg                                                                                                                                                                                                                                  >  1:477713/1‑251 (MQ=255)                                    gTTGATATTGCGGACCAAGCCAAGTATCAAAAGACAAACCGTAGTTTTCTGAGCCGTGATATTGGTGCAAACGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGa                                                                                                                                                                                                               <  1:121513/250‑1 (MQ=255)                                     ttGATATTGCGGACCAAGCCAAGTATCAAAAGACAAACCGTAGTTTTCTGAGCCGTGATATTGGTGCAAACGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAAc                                                                                                                                                                                                             <  1:635506/251‑1 (MQ=255)                                        atatTGCGGACCAAGCCAAGTATCAAAAGACAAACCGTAGTTTTCTGAGCCGTGATATTGGTGCAAACGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATgg                                                                                                                                                                                                                <  1:146989/245‑1 (MQ=255)                                              cGGACCAAGCCAAGTATCAAAAGACAAACCGTAGTTTTCTGAGCCGTGATATTGGTGCAAACGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGAc                                                                                                                                                                                                     <  1:758029/250‑1 (MQ=255)                                                gACCAAGCCAAGTATCAAAAGACAAACCGTAGTTTTCTGAGCCGTGATATTGGTGCAAACGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAAc                                                                                                                                                                                                             >  1:297923/1‑240 (MQ=255)                                                         aagTATCAAAAGACAAACCGTAGTTTTCTGAGCCGTGATATTGGTGCAAACGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGt                                                                                                                                                                                         <  1:847486/251‑1 (MQ=255)                                                                    gACAAACCGTAGTTTTCTGAGCCGTGATATTGGTGCAAACGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTttc                                                                                                                                                                                                >  1:216242/1‑233 (MQ=255)                                                                    gACAAACCGTAGTTTTCTGAGCCGTGATATTGGTGCAAACGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGa                                                                                                                                                                               >  1:545788/1‑250 (MQ=255)                                                                                ttttCTGAGCCGTGATATTGGTGCAAACGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGcc                                                                                                                                                                  <  1:886182/251‑1 (MQ=255)                                                                                 tttCTGAGCCGTGATATTGGTGCAAACGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGATGTTGGTTTTGCCTAAAAg                                                                                                                                                                                                                                                <  1:584185/172‑1 (MQ=255)                                                                                                    gTGCAAACGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGAtt                                                                                                                                              >  1:694682/1‑251 (MQ=255)                                                                                                    gTGCAAACGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGAt                                                                                                                                               <  1:824843/250‑1 (MQ=255)                                                                                                     tGCAAACGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGAtt                                                                                                                                              <  1:150817/250‑1 (MQ=255)                                                                                                         aaCGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGAc                                                                                                                                                                                                     >  1:128767/1‑191 (MQ=255)                                                                                                         aaCGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATg                                                                                                                                          >  1:28793/1‑250 (MQ=255)                                                                                                                           cGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAAtt                                                                                                                       <  1:275590/251‑1 (MQ=255)                                                                                                                           cGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCAAACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAAtt                                                                                                                       <  1:157740/251‑1 (MQ=255)                                                                                                                                      atatTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATg                                                                                                                                                        >  1:332841/1‑207 (MQ=255)                                                                                                                                                 cgTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGa                                                                                                                                                                                                               <  1:738517/141‑1 (MQ=255)                                                                                                                                                  gTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGcc                                                                                                                                 >  1:184050/1‑218 (MQ=255)                                                                                                                                                  gTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGaata                                                                                                >  1:548465/1‑251 (MQ=255)                                                                                                                                                      tCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATAc                                                                                            <  1:138734/251‑1 (MQ=255)                                                                                                                                                      tCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATAc                                                                                            <  1:153910/251‑1 (MQ=255)                                                                                                                                                      tCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATAc                                                                                            <  1:626521/251‑1 (MQ=255)                                                                                                                                                               tatTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGa                                                                                                                                                                                                               <  1:389155/127‑1 (MQ=255)                                                                                                                                                               tatTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGAt                                                                                                      <  1:548428/232‑1 (MQ=255)                                                                                                                                                               tatTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATAc                                                                                            >  1:408976/1‑242 (MQ=255)                                                                                                                                                                atTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATAc                                                                                            >  1:568088/1‑241 (MQ=255)                                                                                                                                                                atTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAAtt                                                                                  >  1:220000/1‑251 (MQ=255)                                                                                                                                                                             aaTAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAAc                                                                     <  1:162283/251‑1 (MQ=255)                                                                                                                                                                                 aCAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGcc                                                                                                                                 <  1:361788/187‑1 (MQ=255)                                                                                                                                                                                       aaCTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAAt                                                            >  1:35147/1‑250 (MQ=255)                                                                                                                                                                                       aaCTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAAt                                                            >  1:571976/1‑250 (MQ=255)                                                                                                                                                                                                    cGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATAc                                              <  1:645705/251‑1 (MQ=255)                                                                                                                                                                                                        cATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTg                                                                                                                                                                                                                                       <  1:111157/62‑1 (MQ=255)                                                                                                                                                                                                        cATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGc                                                                                                                                                                   >  1:893964/1‑130 (MQ=255)                                                                                                                                                                                                        cATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATAcc                                             >  1:568932/1‑248 (MQ=255)                                                                                                                                                                                                        cATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATAAATTTTGGTTTTATTTTTGATGGAATAATACGc                                                                                          >  1:425318/1‑203 (MQ=255)                                                                                                                                                                                                         aTCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACtt                                                                                                                                                                                                   >  1:611116/1‑97 (MQ=255)                                                                                                                                                                                                         aTCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGc                                                                                                                                         >  1:892745/1‑155 (MQ=255)                                                                                                                                                                                                         aTCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATACCATGt                                         >  1:659148/1‑251 (MQ=255)                                                                                                                                                                                                          tCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTttc                                                                                                                                                                                                <  1:189544/99‑1 (MQ=255)                                                                                                                                                                                                          tCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACt                                                                    <  1:662800/223‑1 (MQ=255)                                                                                                                                                                                                          tCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATAcc                                             <  1:209406/246‑1 (MQ=255)                                                                                                                                                                                                          tCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATACCATGtt                                        >  1:217489/1‑251 (MQ=255)                                                                                                                                                                                                          tCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATACCATGtt                                        >  1:480660/1‑251 (MQ=255)                                                                                                                                                                                                            ccACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATACCATGTTcc                                      <  1:21969/251‑1 (MQ=255)                                                                                                                                                                                                              acacCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTg                                                                                                                                                                                                                                       <  1:172001/56‑1 (MQ=255)                                                                                                                                                                                                                acCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATACCATGTTCCtttt                                  >  1:480517/1‑251 (MQ=255)                                                                                                                                                                                                                 ccATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAAtt                                                                                  >  1:610679/1‑202 (MQ=255)                                                                                                                                                                                                                        cTTCAGGCACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATg                                                                                                                                          <  1:554004/139‑1 (MQ=255)                                                                                                                                                                                                                               cACGCGTACGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATACCATGTTCCTTTTTTTCAGTTAATCATg                   >  1:694091/1‑251 (MQ=255)                                                                                                                                                                                                                                      aCGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATACCATGTTCCTTTTTTTCAGTTAATCATGTTTATAg            <  1:899899/251‑1 (MQ=255)                                                                                                                                                                                                                                      aCGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATACCATGTTCCTTTTTTTCAGTTAATCATGTTTATAg            <  1:9935/251‑1 (MQ=255)                                                                                                                                                                                                                                      aCGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATACCATGTTCCTTTTTTTCAGTTAATCATGTTTATAg            <  1:416037/251‑1 (MQ=255)                                                                                                                                                                                                                                      aCGGTGTTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATACCATGTTCCTTTTTTTCAGTTAATCATGTTTATAg            <  1:298669/251‑1 (MQ=255)                                                                                                                                                                                                                                         gtgtTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATAc                                                                                            <  1:337565/168‑1 (MQ=255)                                                                                                                                                                                                                                         gtgtTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATAcc                                             <  1:241702/215‑1 (MQ=255)                                                                                                                                                                                                                                         gtgtTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATACCa                                            <  1:616089/216‑1 (MQ=255)                                                                                                                                                                                                                                         gtgtTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATACCATGTTCCTTTTTTTCAGTTAATCATGTTTATAgtgt         >  1:758506/1‑251 (MQ=255)                                                                                                                                                                                                                                         gtgtTGGTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATACCATGTTCCTTTTTTTCAGTTAATCATGTTTATAgtgt         >  1:875136/1‑251 (MQ=255)                                                                                                                                                                                                                                               gTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATACCATGTTCCTTTTTTTCAGTTAATCATGTTTATAGTGTAAAAg    <  1:5314/250‑1 (MQ=255)                                                                                                                                                                                                                                               gTTTTGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATACCATGTTCCTTTTTTTCAGTTAATCATGTTTATAGTGTAAAAg    <  1:460000/250‑1 (MQ=255)                                                                                                                                                                                                                                                ttttGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATACCATGTTCCTTTTTTTCAGTTAATCATGTTTATAGTGTAAAAGAt  <  1:39358/251‑1 (MQ=255)                                                                                                                                                                                                                                                ttttGCCTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATACCATGTTCCTTTTTTTCAGTTAATCATGTTTATAGTGTAAAAGAt  <  1:172244/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     GAAAAACTTCGCCCAGTGGTTGCTCCATTCGTTCGGTTGATATTGCGGACCAAGCCAAGTATCAAAAGACAAACCGTAGTTTTCTGAGCCGTGATATTGGTGCAAACGGGCGAGTTGCTCACCGAATTCCGCCATATTTTGCGGCGTATTCGGTTGCATATTTAAGCCTTCCAATAACAAAAAACTATGATTTTGCGAGCATCCCACACCATACACTTCAGGCACGCGTACGGTGTTGGTTTTGCTTAAAAGGATCAATTGATCCGCTTCGGCGCGGAACATGGAACGGTAGGACTTTTCATTCAGTTTTACGAAGACGGGTTGAATGCCGTCATTAATGAGCCAAGATTCATGCATTTCGCCCGTATGAATTTTGGTTTTATTTTTGATGGAATAATACGCCCCGAATTGATCAGCTAAAACTTGAGAAATCGATTTCCACATACCATGTTCCTTTTTTTCAGTTAATCATGTTTATAGTGTAAAAGAT  >  1/1006537‑1007026

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