Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	886,527	T→C	A74A (GCT→GCC)	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,527	0	T	C	100.0%	200.4 / NA	68	A74A (GCT→GCC)	unknown	tRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (44/24);  total (44/24)

TCCGTATTCGCGGGGTTGGTAAAGCATGGCGCGGCAGCTAATTCTCCTTGTTGGCGACGATGCCACGTTAGCGCGCCAATGGGAATACTTTCCTGCCCTGAAAAGTGCGGTGTGGATTGGCGGTGTTTTGCCGCACTTTTCCGCTCAACTTCATTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTCATTTATGATGGGCGCCAAGGCATTCATTTAGATGCCCTCGCTATTGCCGCCGGTGCGTTGTCGGAAGGCGGCACATTGTTGCTGTTATTCAATGAGTGGCAGGATTTGGCGGCACAACCCGACGGCGACAGCCTGCGTTGGAGCAGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCGGAAATACAGCTTCCCCATTTATCAACATTTGCCGTCTGCGCTCGTTTTCCCTAAGATTGAAAAGTCGCTTTCAGAACAACA  >  1/886281‑886769                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   tCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTg                                                                                                                                                                                                                                                >  1:179552/1‑251 (MQ=255)                     aaactaTGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGaa                                                                                                                                                                                                                            <  1:248402/246‑1 (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:59991/1‑250 (MQ=255)                                                                 aCGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTgg                                                                                                                                                                                 >  1:630816/1‑250 (MQ=255)                                                                            ccAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTggcggtgcaac                                                                                                                                                                        >  1:756796/1‑242 (MQ=255)                                                                            ccAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTggcggtgcaac                                                                                                                                                                        >  1:2743/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)                                                                                                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:151475/1‑250 (MQ=255)                                                                                                                                          tttCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGc                                                                                                        >  1:239994/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:466534/2‑246 (MQ=255)                                                                                                                                                                                                             taTTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCtggt                                      >  1:545766/2‑246 (MQ=255)                                                                                                                                                                                                             taTTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCtggt                                      >  1:467746/2‑246 (MQ=255)                                                                                                                                                                                                             taTTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCtggt                                      >  1:734236/2‑246 (MQ=255)                                                                                                                                                                                                             taTTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCtggt                                      >  1:51859/2‑246 (MQ=255)                                                                                                                                                                                                                     atgGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTag                             <  1:703132/250‑2 (MQ=255)                                                                                                                                                                                                                         gcgtcagggCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTg                        >  1:7940/8‑251 (MQ=255)                                                                                                                                                                                                                               gggCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGaaacgt                  >  1:722041/2‑248 (MQ=255)                                                                                                                                                                                                                                      cattTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTAt                                                                                                           <  1:782028/155‑1 (MQ=255)                                                                                                                                                                                                                                      cattTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGcc                                                                                                       <  1:116375/159‑1 (MQ=255)                                                                                                                                                                                                                                      cattTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCttt            >  1:726986/1‑250 (MQ=255)                                                                                                                                                                                                                                      cattTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCttt            >  1:132686/1‑250 (MQ=255)                                                                                                                                                                                                                                      cattTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCttt            >  1:373382/1‑250 (MQ=255)                                                                                                                                                                                                                                      cattTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCTGGCCTTCCCTAGGATTGaaa                     >  1:774487/1‑241 (MQ=255)                                                                                                                                                                                                                                      cattTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTCGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCtggtc                                     <  1:424225/225‑5 (MQ=255)                                                                                                                                                                                                                                       attTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCTTTCa          >  1:7795/1‑251 (MQ=255)                                                                                                                                                                                                                                               gctcTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCTTTCAGAactaca  >  1:320599/4‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   TCCGTATTCGCGGGGTTGGTAAAGCATGGCGCGGCAGCTAATTCTCCTTGTTGGCGACGATGCCACGTTAGCGCGCCAATGGGAATACTTTCCTGCCCTGAAAAGTGCGGTGTGGATTGGCGGTGTTTTGCCGCACTTTTCCGCTCAACTTCATTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTCATTTATGATGGGCGCCAAGGCATTCATTTAGATGCCCTCGCTATTGCCGCCGGTGCGTTGTCGGAAGGCGGCACATTGTTGCTGTTATTCAATGAGTGGCAGGATTTGGCGGCACAACCCGACGGCGACAGCCTGCGTTGGAGCAGTGAAAAACACCCGATAAACACACCGCACTTTATTGCCTTTCTGCAAGAAAAAATCCGGAAATACAGCTTCCCCATTTATCAACATTTGCCGTCTGCGCTCGTTTTCCCTAAGATTGAAAAGTCGCTTTCAGAACAACA  >  1/886281‑886769

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