Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	887,844	A→G	A513A (GCA→GCG)	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,844	0	A	G	100.0%	148.5 / NA	52	A513A (GCA→GCG)	unknown	tRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/19);  total (33/19)

GGATTGGAACAGGCGTGTAAATTGCGTTCACTGCGCATTTCCCGTATCGCCGTGCAACCCAACTGGCAGCAGCGGGGACTGGGGCAACGTTTAATCGAACGAG‑TTAAACAACAACGTGCGGTGGATTTCCTTTCCGTCAGTTTCGGTTACACCCCCGAGCTGCTGGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAGGCAAAGATGCAGGCAGCGGCTGTTATTCCGTAGTGGCATTGTGCCCGCTGACGGAGCGCGGACGGGCGTTGGTGCAACAGGCGGAGCAACAATTTAAACGTAATCTGCCCCTTTCGTTGCACCCGCTGGCGACGCAATTTGTCGGAACCGAAACGGATTGGTCGCTCAATGAAACCGATTGGCAATGCTTGGCGGATTTTGCCGATTTCTTCCTGCCGTTGGCG  >  1/887623‑888057                                                                                                                                                                                                                               |                                                                                                                                                                                                                      ggATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCg                                                                                                                                                                                                    >  1:808127/1‑241 (MQ=255)  gATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTtg                                                                                                                                                                                         >  1:85014/1‑251 (MQ=255)  gATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTtg                                                                                                                                                                                         >  1:827358/1‑251 (MQ=255)   aTTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTtg                                                                                                                                                                                         <  1:687759/250‑1 (MQ=255)            ggCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCtaa                                                                                                                                                                               >  1:505795/1‑249 (MQ=255)                            tctttGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCAcgggcggg                                                                                                                                                                >  1:369461/5‑250 (MQ=255)                             ctttGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCt                                                                                                                                                                                 >  1:256412/4‑232 (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)                                                                            gggcTGGGGCAACGTTTAATCGAAC‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCtaac                                                                                                                                                                              >  1:888813/4‑185 (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)                                                                                                       gattAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCtaac                                                                                                                                                                              <  1:287239/160‑4 (MQ=255)                                                                                                       gattAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTg                                                                                     >  1:807808/3‑251 (MQ=255)                                                                                                       gattAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTg                                                                                     >  1:559844/3‑251 (MQ=255)                                                                                                       gattAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTg                                                                                     >  1:140263/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)                                                                                                                          gTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTtgtgtcc                                                                                                                                                                                    >  1:592848/1‑134 (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:497213/1‑238 (MQ=255)                                                                                                                                                         cccccGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGaccggcaaggtttg                                   >  1:271732/1‑238 (MQ=255)                                                                                                                                                         cccccGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGaccggcaaggtttg                                   >  1:766329/1‑238 (MQ=255)                                                                                                                                                         cccccGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGaccggcaaggtttg                                   >  1:38965/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:550499/1‑251 (MQ=255)                                                                                                                                                                                          gTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGAACAACAATTTAAACGCAATTTGCCCCTTTCATTGCACCCGCTGGCGACACAATTTGCCGGAAGCGAAACGGATTGGTCGCTCAATGAAGCCGACCGGCAAGGTTTGGCGGATTTTGCCGATTTCTTCCTACCGTTGGCg  >  1:139986/1‑251 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                      GGATTGGAACAGGCGTGTAAATTGCGTTCACTGCGCATTTCCCGTATCGCCGTGCAACCCAACTGGCAGCAGCGGGGACTGGGGCAACGTTTAATCGAACGAG‑TTAAACAACAACGTGCGGTGGATTTCCTTTCCGTCAGTTTCGGTTACACCCCCGAGCTGCTGGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAGGCAAAGATGCAGGCAGCGGCTGTTATTCCGTAGTGGCATTGTGCCCGCTGACGGAGCGCGGACGGGCGTTGGTGCAACAGGCGGAGCAACAATTTAAACGTAATCTGCCCCTTTCGTTGCACCCGCTGGCGACGCAATTTGTCGGAACCGAAACGGATTGGTCGCTCAATGAAACCGATTGGCAATGCTTGGCGGATTTTGCCGATTTCTTCCTGCCGTTGGCG  >  1/887623‑888057

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