Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	887,891	G→A	R529H (CGC→CAC)	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	887,891	0	G	A	100.0%	172.0 / NA	60	R529H (CGC→CAC)	unknown	tRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/26);  total (34/26)

TCACTGCGCATTTCCCGTATCGCCGTGCAACCCAACTGGCAGCAGCGGGGACTGGGGCAACGTTTAATCGAACGAG‑TTAAACAACAACGTGCGGTGGATTTCCTTTCCGTCAGTTTCGGTTACACCCCCGAGCTGCTGGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAGGCAAAGATGCAGGCAGCGGCTGTTATTCCGTAGTGGCATTGTGCCCGCTGACGGAGCGCGGACGGGCGTTGGTGCAACAGGCGGAGCAACAATTTAAACGTAATCTGCCCCTTTCGTTGCACCCGCTGGCGACGCAATTTGTCGGAACCGAAACGGATTGGTCGCTCAATGAAACCGATTGGCAATGCTTGGCGGATTTTGCCGATTTCTTCCTGCCGTTGGCGTTTGCCTTGCCCGCTATTCGGCGCTTGATGAAATTAAGCGGGGAGCGGGATTTTCCTTTATTGACAGCCTATTGCAA  >  1/887650‑888134                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    tctttGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCAcgggcggg                                                                                                                                                                                                                                             >  1:369461/5‑250 (MQ=255)  ctttGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGggcgggcg                                                                                                                                                                                                                                           >  1:289964/4‑251 (MQ=255)                      gCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGaaca                                                                                                                                                                                                                       >  1:602337/1‑248 (MQ=255)                      gCCGTGCAGCCCAACTAGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGaaca                                                                                                                                                                                                                       >  1:562036/1‑248 (MQ=255)                              gcccAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTa                                                                                                                                                                                                               >  1:635730/2‑251 (MQ=255)                                      ggCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAAt                                                                                                                                                                                                       >  1:168332/1‑251 (MQ=255)                                           cagcGAGGGCTGGGGCAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGc                                                                                                                                                                                                   >  1:819576/1‑250 (MQ=255)                                              cgagggcTGGGGCAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGcccc                                                                                                                                                                                                <  1:465959/244‑1 (MQ=255)                                                  ggcTGGGGCAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTc                                                                                                                                                                                            <  1:447361/248‑1 (MQ=255)                                                       gggCAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGc                                                                                                                                                                                       <  1:704049/250‑1 (MQ=255)                                                       gggCAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCa                                                                                                                                                                                      >  1:476885/1‑251 (MQ=255)                                                          cAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGTAcc                                                                                                                                                                                    <  1:709743/250‑1 (MQ=255)                                                            aCGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACggg                                                                                                                                                                                                                                                 >  1:528975/1‑186 (MQ=255)                                                            aCGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTg                                                                                                                                                                                        <  1:881949/244‑1 (MQ=255)                                                            aCGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCg                                                                                                                                                                                  >  1:810302/1‑250 (MQ=255)                                                                            gattAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTg                                                                                                                                                                  >  1:559844/3‑251 (MQ=255)                                                                            gattAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTg                                                                                                                                                                  >  1:140263/3‑251 (MQ=255)                                                                            gattAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTg                                                                                                                                                                  >  1:807808/3‑251 (MQ=255)                                                                              ttAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGc                                                                                                                                                                                       >  1:77605/1‑228 (MQ=255)                                                                              ttAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTgcc                                                                                                                                                                <  1:422809/251‑3 (MQ=255)                                                                                  acaacaacGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGaa                                                                                                                                                            >  1:394151/1‑251 (MQ=255)                                                                                               gTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCaa                                                                                                                                                                                                                                  >  1:315291/1‑168 (MQ=255)                                                                                                  gATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGt                                                                                                                                             <  1:28006/250‑1 (MQ=255)                                                                                                           tccGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGa                                                                                                                                   <  1:25636/251‑1 (MQ=255)                                                                                                             cGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGaag                                                                                                                                 <  1:8584/251‑2 (MQ=255)                                                                                                                  gTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCAcc                                                                                                                                                                                    >  1:497447/1‑195 (MQ=255)                                                                                                                  gTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGa                                                                                                                             >  1:622352/1‑250 (MQ=255)                                                                                                                              cccccGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCgg                                                                                                                                                              <  1:780018/205‑1 (MQ=255)                                                                                                                              cccccGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGaagc                                                                                                                                <  1:571006/235‑3 (MQ=255)                                                                                                                              cccccGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGaccggccaggtttg                                                                                                                >  1:855603/1‑238 (MQ=255)                                                                                                                              cccccGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGaccggcaaggtttg                                                                                                                >  1:766329/1‑238 (MQ=255)                                                                                                                              cccccGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGaccggcaaggtttg                                                                                                                >  1:497213/1‑238 (MQ=255)                                                                                                                              cccccGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGaccggcaaggtttg                                                                                                                >  1:38965/1‑238 (MQ=255)                                                                                                                              cccccGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGaccggcaaggtttg                                                                                                                >  1:271732/1‑238 (MQ=255)                                                                                                                                cccGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGaccggcaaggtttgg                                                                                                               >  1:748746/1‑236 (MQ=255)                                                                                                                                 ccGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCg                                                                                                             <  1:41675/251‑1 (MQ=255)                                                                                                                                                    gCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATttcttc                                                                                           <  1:383121/250‑1 (MQ=255)                                                                                                                                                       aaaGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCtac                                                                                       <  1:35151/251‑3 (MQ=255)                                                                                                                                                           tGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCg                                                                                     <  1:672984/249‑1 (MQ=255)                                                                                                                                                           tGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCg                                                                                     <  1:639961/249‑1 (MQ=255)                                                                                                                                                           tGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCCAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCg                                                                                     <  1:429372/249‑1 (MQ=255)                                                                                                                                                               gTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCg                                                                               >  1:139986/1‑251 (MQ=255)                                                                                                                                                               gTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCg                                                                               >  1:550499/1‑251 (MQ=255)                                                                                                                                                                                                 tgcgggcagCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCGTTTGTCTTGCCCGCTATTCGGCGCTTGATGaaag                                             <  1:88031/247‑2 (MQ=255)                                                                                                                                                                                                   cgggcagCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCGTTTGTCTTGCCCGCTATTCGGCGCTTGATGaaagta                                           <  1:458958/249‑4 (MQ=255)                                                                                                                                                                                                         gCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCGTTTGTCTTGCCCGCTATTCGGCGCTTGATGAAAGTAAGCggg                                     >  1:323976/1‑251 (MQ=255)                                                                                                                                                                                                             cTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCGTTTGTCTTGCCCGCTATTCGGCGCTTGATGAAAGTAAGCGGGGaac                                 <  1:205980/251‑3 (MQ=255)                                                                                                                                                                                                             cTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCGTTTGTCTTGCCCGCTATTCGGCGCTTGATGAAAGTAAGCGGGGaa                                  >  1:23021/1‑249 (MQ=255)                                                                                                                                                                                                             cTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCGTTTGTCTTGCCCGCTATTCGGCGCTTGATGAAAGTAAGCGGGGaa                                  >  1:548670/1‑249 (MQ=255)                                                                                                                                                                                                             cTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCGTTTGTCTTGCCCGCTATTCGGCGCTTGATGAAAGTAAGCGGGGaa                                  >  1:858142/1‑249 (MQ=255)                                                                                                                                                                                                                   ttCCGTGGTGGCGCTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCGTTTGTCTTGCCCGCTATTCGGCGCTTGATGAAAGTAAGCGGGGAACAGGAtt                           >  1:728512/1‑251 (MQ=255)                                                                                                                                                                                                                       gtggTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCGTTTGTCTTGCCCGCTATTCGGCGCTTGATGAAAGTAAGCGGGGAACAGGATTTTc                        >  1:214156/4‑250 (MQ=255)                                                                                                                                                                                                                       gtggTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCGTTTGTCTTGCCCGCTATTCGGCGCTTGATGAAAGTAAGCGGGGAACAGGATTTTc                        >  1:457608/4‑250 (MQ=255)                                                                                                                                                                                                                        tggTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCGTTTGTCTTGCCCGCTATTCGGCGCTTGATGAAAGTAAGCGGGGAACAGGATTTTCCt                      <  1:616559/249‑1 (MQ=255)                                                                                                                                                                                                                          gTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCGTTTGTCTTGCCCGCTATTCGGCGCTTGATGAAAGTAAGc                                        <  1:752862/231‑1 (MQ=255)                                                                                                                                                                                                                          gTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCGTTTGTCTTGCCCGCTATTCGGCGCTTGATGAAAGTAAGCGGGGAACAGGATTTTCCttt                    >  1:332777/1‑251 (MQ=255)                                                                                                                                                                                                                                   tgtccGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCGTTTGTCTTGCCCGCTATTCGGCGCTTGATGAAAGTAAGCGGGGAACAGGATTTTCCTTTATTGGCAGc           <  1:410682/248‑1 (MQ=255)                                                                                                                                                                                                                                          taaCGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCGTTTGTCTTGCCCGCTATTCGGCGCTTGATGAAAGTAAGCGGGGAACAGGATTTTCCTTTATTGGCAGCCTATTGc    <  1:664681/249‑1 (MQ=255)                                                                                                                                                                                                                                           aaCGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCGTTTGTCTTGCCCGCTATTCGGCGCTTGATGAAAGTAAGCGGGGAACAGGATTTTCCTTTATTGGCAGCCTATTGCa   <  1:551660/250‑1 (MQ=255)                                                                                                                                                                                                                                            aCGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCGTTTGTCTTGCCCGCTATTCGGCGCTTGATGAAAGTAAGCGGGGAACAGGATTTTCCTTTATTGGCAGCCTATTGCaa  <  1:475604/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    TCACTGCGCATTTCCCGTATCGCCGTGCAACCCAACTGGCAGCAGCGGGGACTGGGGCAACGTTTAATCGAACGAG‑TTAAACAACAACGTGCGGTGGATTTCCTTTCCGTCAGTTTCGGTTACACCCCCGAGCTGCTGGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAGGCAAAGATGCAGGCAGCGGCTGTTATTCCGTAGTGGCATTGTGCCCGCTGACGGAGCGCGGACGGGCGTTGGTGCAACAGGCGGAGCAACAATTTAAACGTAATCTGCCCCTTTCGTTGCACCCGCTGGCGACGCAATTTGTCGGAACCGAAACGGATTGGTCGCTCAATGAAACCGATTGGCAATGCTTGGCGGATTTTGCCGATTTCTTCCTGCCGTTGGCGTTTGCCTTGCCCGCTATTCGGCGCTTGATGAAATTAAGCGGGGAGCGGGATTTTCCTTTATTGACAGCCTATTGCAA  >  1/887650‑888134

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