Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,033,820	T→C	C298R (TGC→CGC)	unknown →	Esterase/lipase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,033,820	0	T	C	100.0%	135.5 / NA	50	C298R (TGC→CGC)	unknown	Esterase/lipase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/28);  total (22/28)

GATCAGCCTTGGGTGAACGTCAGCAAAGAAGGTGTTACCGACATTGATTGGGACGGTTATATTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGTGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCATTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATCGGTGCTATGTTAGCTATAAAATTGCGTATGGCGGGTAAGCAGGTGGATTATGAAACGCCATGGGGCGTGCCGCATGCGGGGGATTATGATTTGGATGAGTTGTTTCAGTGGGTGGATAGTATTAGCA  >  1/1033592‑1034042                                                                                                                                                                                                                                     |                                                                                                                                                                                                                               gATCAGCCTTGGGTGAACGTCAGCAAAGAAGGTGTTACCGACATTGATTGGGACGGTTATATTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGaa                                                                                                                                                                                                             >  1:325119/1‑248 (MQ=255) gATCAGCCTTGGGTGAACGTCAGCAAAGAAGGTGTTACCGACATTGATTGGGACGGTTATATTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTa                                                                                                                                                                                                          >  1:588023/1‑251 (MQ=255)      gCCTTGGGTGAACGTCAGCAAAGAAGGTGTTACCGACATTGATTGGGACGGTTATATTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTgg                                                                                                                                                                                                     >  1:506456/1‑251 (MQ=255)       ccTTGGGTGAACGTCAGCAAAGAAGGTGTTACCGACATTGATTGGGACGGTTATATTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTAtgt                                                                                                                                                                                                       <  1:504345/248‑1 (MQ=255)       ccTTGGGTGAACGTCAGCAAAGAAGGTGTTACCGACATTGATTGGGACGGTTATATTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTAtgt                                                                                                                                                                                                       <  1:43536/248‑1 (MQ=255)         ttGGGTGAACGTCAGCAAAGAAGGTGTTACCGACATTGATTGGGACGGTTATATTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATc                                                                                                                                                                                                  <  1:114959/251‑1 (MQ=255)             gTGAACGTCAGCAAAGAAGGTGTTACCGACATTGATTGGGACGGTTATATTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAAc                                                                                                                                                                                                                     <  1:494981/228‑1 (MQ=255)                   gTCAGCAAAGAAGGTGTTACCGACATTGATTGGGACGGTTATATTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTaaaa                                                                                                                                                                                         >  1:759255/1‑250 (MQ=255)                       gCAAAGAAGGTGTTACCGACATTGATTGGGACGGTTATATTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGc                                                                                                                                                                                                                   >  1:502284/1‑220 (MQ=255)                                  gtTACCGACATTGATTGGGACGGTTATATTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCa                                                                                                                                                                          >  1:751843/1‑250 (MQ=255)                                         aCATTGATTGGGACGGTTATATTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAAt                                                                                                                                                                                                            >  1:814785/1‑209 (MQ=255)                                              gattgGGACGGTTATATTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAg                                                                                                                                                                                              >  1:60073/1‑218 (MQ=255)                                                   ggACGGTTATATTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCg                                                                                                                                                        >  1:127297/1‑251 (MQ=255)                                                      cGGTTATATTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATATGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGgcg                                                                                                                                                                            >  1:454661/1‑228 (MQ=255)                                                        gTTATATTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACgc                                                                                                                                                                                      <  1:370468/216‑1 (MQ=255)                                                        gTTATATTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGt                                                                                                                                                   >  1:719169/1‑251 (MQ=255)                                                        gTTATATTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGt                                                                                                                                                   >  1:483800/1‑251 (MQ=255)                                                             atTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATAc                                                                                                                                              <  1:21785/251‑1 (MQ=255)                                                             atTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGCTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATAc                                                                                                                                              <  1:720349/251‑1 (MQ=255)                                                               tCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACtt                                                                                                                                            <  1:680738/251‑1 (MQ=255)                                                                                    aTCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGt                                                                                                                       <  1:45715/251‑1 (MQ=255)                                                                                         cACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCt                                                                                                                  <  1:387728/251‑1 (MQ=255)                                                                                          aCCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCta                                                                                                                 <  1:95481/251‑1 (MQ=255)                                                                                              gCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCTATaaa                                                                                                             >  1:145930/1‑251 (MQ=255)                                                                                               cATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATc                                                                                                                                                                                                  <  1:665882/165‑1 (MQ=255)                                                                                               cATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCTATaaaa                                                                                                            >  1:744244/1‑251 (MQ=255)                                                                                                                         ggCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCTATAAAATTGCGTATGGCGGGGAAGCaccgt                                                                                    >  1:860908/1‑245 (MQ=255)                                                                                                                             aaaaTAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAAtt                                                                                                                                                                                                           >  1:800316/1‑126 (MQ=255)                                                                                                                               aataatTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCTATAAAATTGCGTATGGCGGGTAAGCAGGTGGATTATGa                                                                            <  1:432274/251‑1 (MQ=255)                                                                                                                                 taatTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCTATAAAATTGCGTATGGCGGGTAAGCAGGTGGATTATGaaa                                                                          <  1:512607/251‑1 (MQ=255)                                                                                                                                            gTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCTATAAAATTGCGTATGGCGGGTAAGCAGGTGGATTATGAAACGCCATgggg                                                                <  1:861637/250‑1 (MQ=255)                                                                                                                                                 gAAAGCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGGACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATGCTTCGCACGCTATTGGTGCTATGTTAGCTATAAAATTGCGTATGGCGGGTAAGCAGGTGGATTATGAAACGCCATGGGGCGTGc                                                           <  1:395143/250‑1 (MQ=255)                                                                                                                                                     ccGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTg                                                                                                                             >  1:779200/1‑180 (MQ=255)                                                                                                                                                        tGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCTATAAAATTGCGTATGGCGGGTAAGCAGGTGGATTATGAAACGCCATGGGGCGTGCCGCATGc                                                    <  1:856688/250‑1 (MQ=255)                                                                                                                                                             aatCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCTATAAAATTGCGTATGGCGGGTAAGCAGGTGGATTATGAAACGCCATgggg                                                                <  1:252296/233‑1 (MQ=255)                                                                                                                                                             aatCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCTATAAAATTGCGTATGGCGGGTAAGCAGGTGGATTATGAAACGCCATGGGGCGTGCCGCATGCGGGGGa                                              <  1:35027/251‑1 (MQ=255)                                                                                                                                                                  gCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCt                                                                                                                  >  1:147507/1‑178 (MQ=255)                                                                                                                                                                  gCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCTATAAAATTGCGTATGGCGGGTAAGCAGGTGGATTATGa                                                                            <  1:347831/216‑1 (MQ=255)                                                                                                                                                                          ccGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAAt                                                                                                                                                                                                            <  1:261861/80‑1 (MQ=255)                                                                                                                                                                          ccGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCTATAAAATTGCGTATGGCGGGTAAGCAGGTGGATTATGAAACGCCATGGGGCGTGCCGCATGCGGGGGATTATGATTTGGAt                                 <  1:886790/251‑1 (MQ=255)                                                                                                                                                                                        cAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTg                                                                                                                             >  1:679117/1‑145 (MQ=255)                                                                                                                                                                                          acacAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCTATAAAATTGCGTATGGCGGGTAAGCAGGTGGATTATGAAACGc                                                                       <  1:84188/197‑1 (MQ=255)                                                                                                                                                                                               gCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCTATAAAATTGCGTATGGc                                                                                                 <  1:417289/166‑1 (MQ=255)                                                                                                                                                                                               gCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCTATAAAATTGCGTATGGCGGGTAAGAAGGTGGATTATGAAACGCCATGGGGCGTGCCGCATGCGGGGGATTATGATTTGGATGAGTTGGTTCAGTGGGTGGAt            >  1:261573/1‑251 (MQ=255)                                                                                                                                                                                                 aCGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCTATAAAATTGCGTATGGCGGGTAAGCAGGTGGATTATGAAACGCCATGGGGCGTGCCGCATGCGGGGGATTATGATTTGGATGAGTTGTTTCAGTGGGTGGATAg          <  1:111146/251‑1 (MQ=255)                                                                                                                                                                                                    gAAAAGGGAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCt                                                                                                                  >  1:787782/1‑144 (MQ=255)                                                                                                                                                                                                         gggAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTg                                                                                                                             <  1:854424/128‑1 (MQ=255)                                                                                                                                                                                                          ggAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCTATAAAATTGCGTATGGCGGGTAAGCAGGTGGATTATGAAACGCCATGGGGCGTGCCGCATGCGGGGGATTATGATTTGGATGAGTTGTTTCAGTGGGTGGATAGTATTAGCa  <  1:732043/250‑1 (MQ=255)                                                                                                                                                                                                          ggAAAGATGGCGGATAAAAATGTGGTGCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCTATAAAATTGCGTATGGCGGGTAAGCAGGTGGATTATGAAACGCCATGGGGCGTGCCGCATGCGGGGGATTATGATTTGGATGAGTTGTTTCAGTGGGTGGATAGTATTAGCa  <  1:87170/250‑1 (MQ=255)                                                                                                                                                                                                                                 ggtgCGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCACTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATTGGTGCTATGTTAGCTATAAAATTGCGTATGGCGGGTAAGCAGGGGGATTATGAAACGCCATGGGGCGTGCCGCATGCGGGGGATTATGATTTGGATGAGGTGTTTCAgtg                   >  1:479642/1‑210 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                               GATCAGCCTTGGGTGAACGTCAGCAAAGAAGGTGTTACCGACATTGATTGGGACGGTTATATTCACGCGGAAAAACGAATGAAATCGCCACCGGCATTTGATGCGCTGAATTTGAGTTCGGGCGAAAATAATTTATTTGGTACGGAAACCGTGAATAATCAGCATTTCACCGTTTTTTCTTGGCAACACAGCACGGAAAAGGGAAAGATGGCGGATAAAAATGTGGTGTGCTTAATGAACGCCATGAATTATGTGGATCAAAGTAAAACGCAGCATTGGCGCATTCGGGTCGGCACCTCCGACCGTGATACTTCGCACGCTATCGGTGCTATGTTAGCTATAAAATTGCGTATGGCGGGTAAGCAGGTGGATTATGAAACGCCATGGGGCGTGCCGCATGCGGGGGATTATGATTTGGATGAGTTGTTTCAGTGGGTGGATAGTATTAGCA  >  1/1033592‑1034042

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