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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	887,679	0	A	G	100.0%	189.1 / NA	64	Q458Q (CAA→CAG)	unknown	tRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/29);  total (35/29)

TGGATTTACGTCGATTATTGGATGCGCCGCGACAACAATTTTATCTTGCGCAGACGCAGGAGGCCTTGTTGGGTTGTGTGTGGGCTGTGCCTGAGGGCGAATTTGGTGACGACGCATTAATTCGTAAAATTCGGCGCGGACAACGTCGTCCACGTGGTAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCACTGCGCATTTCCCGTATCGCCGTGCAACCCAACTGGCAGCAGCGGGGACTGGGGCAACGTTTAATCGAA‑CGAG‑TTAAACAACAACGTGCGGTGGATTTCCTTTCCGTCAGTTTCGGTTACACCCCCGAGCTGCTGGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAGGCAAAGATGCAGGCAGCGGCTGTTATTCCGTAGTGGCATTGTGCCCGCTGACGGAGCGCGGACGGGCGTTGGTGCAACAGGCGGAGCA  >  1/887432‑887921                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     tGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGccc                                                                                                                                                                                                                                                   <  1:158901/251‑1 (MQ=255)              attattGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCagcagc                                                                                                                                                                                                                                      <  1:69320/251‑1 (MQ=255)                 attGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCgag                                                                                                                                                                                                                                   >  1:483325/1‑249 (MQ=255)                 attGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCgag                                                                                                                                                                                                                                   >  1:232337/1‑249 (MQ=255)                  ttGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGcagcag                                                                                                                                                                                                                                       >  1:789812/1‑246 (MQ=255)                  ttGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCgagg                                                                                                                                                                                                                                  >  1:781539/1‑248 (MQ=255)                         caccgcGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTgggg                                                                                                                                                                                                                           <  1:759887/249‑1 (MQ=255)                                    caaTTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAAt                                                                                                                                                                                                                <  1:143521/251‑1 (MQ=255)                                      attttattttGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCg                                                                                                                                                                                                              <  1:445071/251‑1 (MQ=255)                                           attttgcacaaacgcgggaaaccctgttgGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGAAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAAt                                                                                                                                                                                                                <  1:9417/220‑1 (MQ=255)                                                     aacgcgggaaaccctgttgGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAaacaac                                                                                                                                                                                             <  1:144890/237‑1 (MQ=255)                                                       cgcgggaaaccctgttgGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAAt                                                                                                                                                                                                                >  1:790218/13‑232 (MQ=255)                                                       cgcgggaaaccctgttgGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAaacaacaa                                                                                                                                                                                           <  1:585054/239‑1 (MQ=255)                                                           ggaaaccctgttgGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGt                                                                                                                                                                                        <  1:91252/242‑1 (MQ=255)                                                                  ctgttgGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTggg                                                                                                                                                                                >  1:596187/2‑250 (MQ=255)                                                                  ctgttgGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTggg                                                                                                                                                                                >  1:308059/2‑250 (MQ=255)                                                                    gttgGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGAtt                                                                                                                                                                              >  1:386119/1‑251 (MQ=255)                                                                            gtgtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑cagattaaac                                                                                                                                                                                                 <  1:539036/225‑10 (MQ=255)                                                                            gtgtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAaacaacaa                                                                                                                                                                                           <  1:471552/230‑1 (MQ=255)                                                                            gtgtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGt                                                                                                                                                                                        >  1:297852/1‑233 (MQ=255)                                                                            gtgtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGAttttctttc                                                                                                                                                                       >  1:892622/1‑250 (MQ=255)                                                                            gtgtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGGGGAttttctttc                                                                                                                                                                       >  1:30271/1‑250 (MQ=255)                                                                              gtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTTCCGTATCGCCGTGCAGCCCAACTGGCCGCAGCGAGGGCTGGGGCCACGTTTAATCGAACCGA‑‑TTAAACAACAACCTGCGGTGGATTTTCTTTCCGt                                                                                                                                                                    >  1:16077/1‑251 (MQ=255)                                                                              gtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAaacaacaa                                                                                                                                                                                           <  1:437368/228‑1 (MQ=255)                                                                              gtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACg                                                                                                                                                                                         >  1:323252/1‑230 (MQ=255)                                                                              gtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGt                                                                                                                                                                    >  1:598662/1‑251 (MQ=255)                                                                               tgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGATCTTTGCGCATTTCCCGCATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCCGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGACTTTCTTTCCGt                                                                                                                                                                    >  1:85052/1‑250 (MQ=255)                                                                                             gaaggaggatttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTaca                                                                                                                                                      >  1:338311/22‑250 (MQ=255)                                                                                                    gatttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGGTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaa                                                                                                                                              >  1:770770/15‑250 (MQ=255)                                                                                                     atttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac                                                                                                                                             <  1:836462/238‑3 (MQ=255)                                                                                                     atttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac                                                                                                                                             <  1:543915/238‑3 (MQ=255)                                                                                                     atttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac                                                                                                                                             <  1:578174/238‑3 (MQ=255)                                                                                                     atttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac                                                                                                                                             <  1:149617/238‑3 (MQ=255)                                                                                                     atttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac                                                                                                                                             <  1:217141/238‑3 (MQ=255)                                                                                                     atttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGGTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac                                                                                                                                             <  1:446276/238‑3 (MQ=255)                                                                                                       ttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAactg                                                                                                                                           <  1:305880/240‑1 (MQ=255)                                                                                                           ggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑cagattaaac                                                                                                                                                                                                 >  1:487161/8‑185 (MQ=255)                                                                                                              tgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCgagggc                                                                                                                                                                                                                                <  1:227295/157‑6 (MQ=255)                                                                                                              tgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCg                                                                                                                                                                     <  1:390200/214‑1 (MQ=255)                                                                                                                atgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑cagattaaaca                                                                                                                                                                                                >  1:325522/3‑180 (MQ=255)                                                                                                                atgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACAcccc                                                                                                                                                  <  1:848980/233‑1 (MQ=255)                                                                                                                atgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac                                                                                                                                             >  1:366912/3‑238 (MQ=255)                                                                                                                            gtcaaaTTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac                                                                                                                                             <  1:143746/225‑3 (MQ=255)                                                                                                                                      acgcgGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTAttt                                                                                                            <  1:135587/250‑1 (MQ=255)                                                                                                                                          ggACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAaacaacaa                                                                                                                                                                                           >  1:489061/1‑168 (MQ=255)                                                                                                                                          ggACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTggg                                                                                                        >  1:101494/1‑250 (MQ=255)                                                                                                                                                   cgcctgcgcggcaaTTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAg                                                                                               >  1:27586/13‑251 (MQ=255)                                                                                                                                                             gcaaTTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGc                                                                                                                                                                                      >  1:569874/3‑154 (MQ=255)                                                                                                                                                                 tttGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATg                                                                                 <  1:159577/251‑1 (MQ=255)                                                                                                                                                                              aTGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac                                                                                                                                             >  1:340791/1‑176 (MQ=255)                                                                                                                                                                                     gtTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTAttt                                                                                                            <  1:127976/204‑1 (MQ=255)                                                                                                                                                                                                ggATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCg                                                            >  1:808127/1‑241 (MQ=255)                                                                                                                                                                                                 gATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAactg                                                                                                                                           >  1:287907/1‑161 (MQ=255)                                                                                                                                                                                                 gATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTtg                                                 >  1:827358/1‑251 (MQ=255)                                                                                                                                                                                                 gATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTtg                                                 >  1:85014/1‑251 (MQ=255)                                                                                                                                                                                                  aTTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTtg                                                 <  1:687759/250‑1 (MQ=255)                                                                                                                                                                                                           ggCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCtaa                                       >  1:505795/1‑249 (MQ=255)                                                                                                                                                                                                                           tctttGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCAcgggcggg                        >  1:369461/5‑250 (MQ=255)                                                                                                                                                                                                                            ctttGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCt                                         >  1:256412/4‑232 (MQ=255)                                                                                                                                                                                                                            ctttGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGggcgggcg                      >  1:289964/4‑251 (MQ=255)                                                                                                                                                                                                                                           gTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGttt                                                                                                                 <  1:423762/145‑1 (MQ=255)                                                                                                                                                                                                                                                gCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGaaca  >  1:602337/1‑248 (MQ=255)                                                                                                                                                                                                                                                gCCGTGCAGCCCAACTAGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAACTGCTCGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAAGCAAAGATGCGGGCAGCGGCTGTTATTCCGTGGTGGCGTTGTGTCCGCTAACGGAGCACGGGCGGGCGTTGGTGCAACGGGCGGaaca  >  1:562036/1‑248 (MQ=255)                                                                                                                                                                                                                                                   gTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac                                                                                                                                             <  1:704660/109‑3 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     TGGATTTACGTCGATTATTGGATGCGCCGCGACAACAATTTTATCTTGCGCAGACGCAGGAGGCCTTGTTGGGTTGTGTGTGGGCTGTGCCTGAGGGCGAATTTGGTGACGACGCATTAATTCGTAAAATTCGGCGCGGACAACGTCGTCCACGTGGTAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCACTGCGCATTTCCCGTATCGCCGTGCAACCCAACTGGCAGCAGCGGGGACTGGGGCAACGTTTAATCGAA‑CGAG‑TTAAACAACAACGTGCGGTGGATTTCCTTTCCGTCAGTTTCGGTTACACCCCCGAGCTGCTGGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGTGCATTTCAGCGAAGGCAAAGATGCAGGCAGCGGCTGTTATTCCGTAGTGGCATTGTGCCCGCTGACGGAGCGCGGACGGGCGTTGGTGCAACAGGCGGAGCA  >  1/887432‑887921

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