Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	886,503	A→G	Q66Q (CAA→CAG)	unknown →	tRNA cytosine(34) acetyltransferase (EC 2.3.1.193)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	886,503	0	A	G	100.0%	180.2 / NA	63	Q66Q (CAA→CAG)	unknown	tRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (42/21);  total (42/21)

GGATGGATCCGTATTCGCGGGGTTGGTAAAGCATGGCGCGGCAGCTAATTCTCCTTGTTGGCGACGATGCCACGTTAGCGCGCCAATGGGAATACTTTCCTGCCCTGAAAAGTGCGGTGTGGATTGGCGGTGTTTTGCCGCACTTTTCCGCTCAACTTCATTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTCATTTATGATGGGCGCCAAGGCATTCATTTAGATGCCCTCGCTATTGCCGCCGGTGCGTTGTCGGAAGGCGGCACATTGTTGCTGTTATTCAATGAGTGGCAGGATTTGGCGGCACAACCCGACGGCGACAGCCTGCGTTGGAGCAGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCGGAAATACAGCTTCCCCATTTATCAACATTTGCCGTCTGCGCTCGTTTTCCCTAA  >  1/886274‑886742                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                ggacggatCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGctc                                                                                                                                                                                                                              >  1:28446/5‑247 (MQ=255)      gatCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATg                                                                                                                                                                                                                                 >  1:301001/1‑241 (MQ=255)      gatCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGccat                                                                                                                                                                                                                       >  1:238826/1‑248 (MQ=255)        tCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTg                                                                                                                                                                                                                     >  1:179552/1‑251 (MQ=255)                            aaactaTGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGAt                                                                                                                                                                                                                                  >  1:492585/6‑218 (MQ=255)                            aaactaTGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGaa                                                                                                                                                                                                 <  1:248402/246‑1 (MQ=255)                              actaTGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGcc                                                                                                                                                                                                                         >  1:592869/4‑224 (MQ=255)                               ctaTGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCAGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAgg                                                                                                                                                                                               >  1:212219/3‑250 (MQ=255)                                          cagttaattcttctcattgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATAttgt                                                                                                                                                                                    <  1:300262/234‑1 (MQ=255)                                          cagttaattcttctcattgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATAttgt                                                                                                                                                                                    <  1:193059/234‑1 (MQ=255)                                           agttaattcttctcattgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATAttgttg                                                                                                                                                                                  >  1:470535/16‑251 (MQ=255)                                           agttaattcttctcattgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATAttgttg                                                                                                                                                                                  >  1:870287/16‑251 (MQ=255)                                            gttaattcttctcattgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATAttgt                                                                                                                                                                                    <  1:271834/234‑1 (MQ=255)                                             ttaattcttctcattgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATAttgttg                                                                                                                                                                                  <  1:56150/236‑1 (MQ=255)                                                  tcttctcattgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTAt                                                                                                                                                                           >  1:726572/9‑251 (MQ=255)                                                    ttctcattgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTc                                                                                                                                                                         <  1:104009/245‑1 (MQ=255)                                                         attgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGa                                                                                                                                                                    <  1:867890/250‑1 (MQ=255)                                                          ttgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGat                                                                                                                                                                   >  1:307769/1‑250 (MQ=255)                                                          ttgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGa                                                                                                                                                                    <  1:14211/250‑1 (MQ=255)                                                            ggtgacgaTGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATAttgt                                                                                                                                                                                    >  1:411686/4‑232 (MQ=255)                                                                        aCGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTggc                                                                                                                                                     <  1:29227/251‑1 (MQ=255)                                                                        aCGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTgg                                                                                                                                                      >  1:630816/1‑250 (MQ=255)                                                                        aCGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTgg                                                                                                                                                      >  1:59991/1‑250 (MQ=255)                                                                                   ccAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTggcggtgcaac                                                                                                                                             >  1:2743/1‑242 (MQ=255)                                                                                   ccAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTggcggtgcaac                                                                                                                                             >  1:756796/1‑242 (MQ=255)                                                                                       tGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTggcggtgcaacctgatgg                                                                                                                                      <  1:63206/251‑14 (MQ=255)                                                                                        gggaacattttCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTggcggtgcaacctgatggg                                                                                                                                     >  1:253168/9‑237 (MQ=255)                                                                                            acattttCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTggcggtgcaacctgatggtgaca                                                                                                                                 >  1:309866/5‑233 (MQ=255)                                                                                                      gCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGAt                                                                                                                                                                                                                                  >  1:27031/1‑144 (MQ=255)                                                                                                       cccTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGt                                                                                                                                                                                                      >  1:71496/1‑171 (MQ=255)                                                                                                                gTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTGTTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTggcggtgcaacc                                                                                                                                            <  1:893088/220‑8 (MQ=255)                                                                                                                     gtgaaacgccgccaatccataccGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGc                                                                                                                   >  1:335772/22‑240 (MQ=255)                                                                                                                     gtatgGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAaca                                                                                                        >  1:235089/4‑251 (MQ=255)                                                                                                                       atgGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGa                                                                                                                                                                                                  >  1:287215/2‑159 (MQ=255)                                                                                                                         ggATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCg                                                                                                    >  1:212676/1‑251 (MQ=255)                                                                                                                          gATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGa                                                                                                   >  1:210987/1‑251 (MQ=255)                                                                                                                           aTTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGAt                                                                                                  <  1:642971/251‑1 (MQ=255)                                                                                                                                cggcgtttcaccGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAac                                                                                              <  1:719897/240‑1 (MQ=255)                                                                                                                                  gcgtttcaccGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAac                                                                                              >  1:71257/9‑248 (MQ=255)                                                                                                                                   cgtttcaccGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAacacac                                                                                          <  1:409167/244‑1 (MQ=255)                                                                                                                                       tcaccGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCaa                                                                                                                                                                       >  1:382846/4‑170 (MQ=255)                                                                                                                                             cACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTAt                                                                                <  1:37136/251‑1 (MQ=255)                                                                                                                                               cTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTg                                                                              >  1:436916/1‑251 (MQ=255)                                                                                                                                                ttttCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGc                                                                             <  1:25093/251‑1 (MQ=255)                                                                                                                                                ttttCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGc                                                                             <  1:322993/251‑1 (MQ=255)                                                                                                                                                 tttCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGATCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGc                                                                             >  1:546230/1‑250 (MQ=255)                                                                                                                                                 tttCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGcc                                                                            <  1:163458/251‑1 (MQ=255)                                                                                                                                                 tttCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGc                                                                             >  1:239994/1‑250 (MQ=255)                                                                                                                                                 tttCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGc                                                                             >  1:151475/1‑250 (MQ=255)                                                                                                                                                 tttCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGc                                                                             >  1:552108/1‑250 (MQ=255)                                                                                                                                                                                     aaaaTCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTggcggtgcaacctgatgg                                                                                                                                      <  1:485204/157‑14 (MQ=255)                                                                                                                                                                                       aaTCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCAt                                      <  1:36715/251‑1 (MQ=255)                                                                                                                                                                                                    gCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCaatat                            >  1:473493/1‑245 (MQ=255)                                                                                                                                                                                                           ccccGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCgcccg                  >  1:802091/1‑247 (MQ=255)                                                                                                                                                                                                               gcgcTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGc                                                                             >  1:335624/1‑188 (MQ=255)                                                                                                                                                                                                               gcgcTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCgcccgcg                >  1:591353/1‑243 (MQ=255)                                                                                                                                                                                                                    taTTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGTCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCttgt           >  1:582163/2‑246 (MQ=255)                                                                                                                                                                                                                    taTTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCtggt           >  1:545766/2‑246 (MQ=255)                                                                                                                                                                                                                    taTTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCtggt           >  1:51859/2‑246 (MQ=255)                                                                                                                                                                                                                    taTTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCtggt           >  1:734236/2‑246 (MQ=255)                                                                                                                                                                                                                    taTTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCtggt           >  1:467746/2‑246 (MQ=255)                                                                                                                                                                                                                    taTTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCtggt           >  1:466534/2‑246 (MQ=255)                                                                                                                                                                                                                            atgGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTag  <  1:703132/250‑2 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                GGATGGATCCGTATTCGCGGGGTTGGTAAAGCATGGCGCGGCAGCTAATTCTCCTTGTTGGCGACGATGCCACGTTAGCGCGCCAATGGGAATACTTTCCTGCCCTGAAAAGTGCGGTGTGGATTGGCGGTGTTTTGCCGCACTTTTCCGCTCAACTTCATTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTCATTTATGATGGGCGCCAAGGCATTCATTTAGATGCCCTCGCTATTGCCGCCGGTGCGTTGTCGGAAGGCGGCACATTGTTGCTGTTATTCAATGAGTGGCAGGATTTGGCGGCACAACCCGACGGCGACAGCCTGCGTTGGAGCAGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCGGAAATACAGCTTCCCCATTTATCAACATTTGCCGTCTGCGCTCGTTTTCCCTAA  >  1/886274‑886742

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