Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	886,428	A→G	Q41Q (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,428	0	A	G	100.0%	197.0 / NA	66	Q41Q (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 (46/20);  total (46/20)

CAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGCGAACCGCTCTATAAAAGCAATAAAAATATCAATGTAGACTGGATGGATCCGTATTCGCGGGGTTGGTAAAGCATGGCGCGGCAGCTAATTCTCCTTGTTGGCGACGATGCCACGTTAGCGCGCCAATGGGAATACTTTCCTGCCCTGAAAAGTGCGGTGTGGATTGGCGGTGTTTTGCCGCACTTTTCCGCTCAACTTCATTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTCATTTATGATGGGCGCCAAGGCATTCATTTAGATGCCCTCGCTATTGCCGCCGGTGCGTTGTCGGAAGGCGGCACATTGTTGCTGTTATTCAATGAGTGGCAGGATTTGGCGGCACAACCCGACGGCGACAGCCTGCGTTGGAGCAGTGAAAAACACCCGATAAACACACCGCACTTTATTGCC  >  1/886195‑886668                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 cAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAg                                                                                                                                                                                                                                 >  1:572754/1‑251 (MQ=255)   aCGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCa                                                                                                                                                                                                                               >  1:198051/1‑251 (MQ=255)              ttGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAAt                                                                                                                                                                                                                    <  1:108532/251‑1 (MQ=255)                                              ctctATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCTCGCGCTTATTTa                                                                                                                                                                                    >  1:467338/1‑251 (MQ=255)                                                ctATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCttatttat                                                                                                                                                                                   >  1:154465/1‑250 (MQ=255)                                                ctATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTatg                                                                                                                                                                                  >  1:11930/1‑251 (MQ=255)                                                ctATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTatg                                                                                                                                                                                  >  1:590757/1‑251 (MQ=255)                                                ctATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTatg                                                                                                                                                                                  >  1:823761/1‑251 (MQ=255)                                                           aTAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGgcgtcag                                                                                                                                                                       >  1:160657/1‑247 (MQ=255)                                                                  tattaaTGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTccc                                                                                                                                                                                                >  1:577465/5‑219 (MQ=255)                                                                      aaTGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGgcgtcag                                                                                                                                                                       >  1:558484/1‑236 (MQ=255)                                                                                ggacggatCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGctc                                                                                                                                                    >  1:28446/5‑247 (MQ=255)                                                                                     gatCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTatga                                                                                                                                                                                 <  1:373470/215‑1 (MQ=255)                                                                                     gatCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATg                                                                                                                                                       >  1:301001/1‑241 (MQ=255)                                                                                     gatCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGccat                                                                                                                                             >  1:238826/1‑248 (MQ=255)                                                                                       tCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAGGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTatga                                                                                                                                                                                 >  1:437524/1‑213 (MQ=255)                                                                                       tCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTc                                                                                                                                                                                                  >  1:647836/1‑196 (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)                                                                                                                                       cattgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTatga                                                                                                                                                                                 >  1:23565/3‑165 (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:59991/1‑250 (MQ=255)                                                                                                                                                       aCGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTgg                                                                            >  1:630816/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)                                                                                                                                                                                               gTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGa                                                                                                                                                                                                      >  1:575030/1‑88 (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:322993/251‑1 (MQ=255)                                                                                                                                                                                                                               ttttCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGc   <  1:25093/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:552108/1‑250 (MQ=255)                                                                                                                                                                                                                                tttCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGc   >  1:151475/1‑250 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 CAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGCGAACCGCTCTATAAAAGCAATAAAAATATCAATGTAGACTGGATGGATCCGTATTCGCGGGGTTGGTAAAGCATGGCGCGGCAGCTAATTCTCCTTGTTGGCGACGATGCCACGTTAGCGCGCCAATGGGAATACTTTCCTGCCCTGAAAAGTGCGGTGTGGATTGGCGGTGTTTTGCCGCACTTTTCCGCTCAACTTCATTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTCATTTATGATGGGCGCCAAGGCATTCATTTAGATGCCCTCGCTATTGCCGCCGGTGCGTTGTCGGAAGGCGGCACATTGTTGCTGTTATTCAATGAGTGGCAGGATTTGGCGGCACAACCCGACGGCGACAGCCTGCGTTGGAGCAGTGAAAAACACCCGATAAACACACCGCACTTTATTGCC  >  1/886195‑886668

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: ‑