Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	887,580	T→C	R425R (CGT→CGC)	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,580	0	T	C	100.0%	190.7 / NA	68	R425R (CGT→CGC)	unknown	tRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/35);  total (33/35)

AAGAACAGATTCCGCCTTATCAAATCGAATCGGTTTATGGCTTGCTGACCTTGGCGCATTACCGAACCTCGCCGTTGGATTTACGTCGATTATTGGATGCGCCGCGACAACAATTTTATCTTGCGCAGACGCAGGAGGCCTTGTTGGGTTGTGTGTGGGCTGTGCCTGAGGGCGAATTTGGTGACGACGCATTAATTCGTAAAATTCGGCGCGGACAACGTCGTCCACGTGGTAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCACTGCGCATTTCCCGTATCGCCGTGCAACCCAACTGGCAGCAGCGGGGACTGGGGCAACGTTTAATCGAA‑CGAG‑TTAAACAACAACGTGCGGTGGATTTCCTTTCCGTCAGTTTCGGTTACACCCCCGAGCTGCTGGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGT  >  1/887357‑887819                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  aagaaCGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAAtg                                                                                                                                                                                                                        <  1:798858/251‑1 (MQ=255)    aacggATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAAtgttg                                                                                                                                                                                                                     <  1:809807/247‑1 (MQ=255)      cggATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTtgtg                                                                                                                                                                                                                   <  1:390128/249‑1 (MQ=255)                               cgGTTTATGGCTTGCTGACTTTGGCGCATTACCTCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGGCGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCg                                                                                                                                                                                          <  1:506287/251‑1 (MQ=255)                                 gTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAAc                                                                                                                                                                                               <  1:882176/244‑1 (MQ=255)                                    tATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTaa                                                                                                                                                                                     <  1:733341/251‑1 (MQ=255)                                            gCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAAtg                                                                                                                                                                                                                        >  1:390954/1‑208 (MQ=255)                                                           ttacagcaCCTCGCCGTTGGATTTACGACGATTATTGGATGCACCGAGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGAGTGTGGGCGGAGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGGAAATTGGGTCCTTTGCGCATTTcc                                                                                                                                                               >  1:210964/8‑250 (MQ=255)                                                           ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTtgtg                                                                                                                                                                                                                   <  1:701921/198‑1 (MQ=255)                                                           ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAAc                                                                                                                                                                                               >  1:728742/1‑218 (MQ=255)                                                           ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTcc                                                                                                                                                               <  1:582402/250‑1 (MQ=255)                                                           ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTcc                                                                                                                                                               >  1:51847/1‑250 (MQ=255)                                                           ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTcc                                                                                                                                                               >  1:698893/1‑250 (MQ=255)                                                           ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTcc                                                                                                                                                               >  1:209876/1‑250 (MQ=255)                                                             aCCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTaa                                                                                                                                                                                     >  1:632915/1‑226 (MQ=255)                                                                            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)                                                                                                        gcgACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTaa                                                                                                                                                                                     <  1:283615/183‑1 (MQ=255)                                                                                                               caaTTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAAt                                                                                                          <  1:143521/251‑1 (MQ=255)                                                                                                                 attttattttGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGTCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGGATCGCCGTGc                                                                                                                                                  >  1:615681/1‑209 (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)                                                                                                                                             ctgttgGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAAc                                                                                                                                                                                               >  1:701442/2‑136 (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)                                                                                                                                                       gtgtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTa                                                                                                                                                           >  1:833883/1‑162 (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)                                                                                                                                                           gtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGc                                                                                                                                                       >  1:830742/1‑162 (MQ=255)                                                                                                                                                                        gaaggaggatttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTaca                                                >  1:338311/22‑250 (MQ=255)                                                                                                                                                                               gatttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGGTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaa                                        >  1:770770/15‑250 (MQ=255)                                                                                                                                                                                atttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGc                                                                                                                                                  <  1:147719/133‑1 (MQ=255)                                                                                                                                                                                atttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac                                       <  1:543915/238‑3 (MQ=255)                                                                                                                                                                                atttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac                                       <  1:836462/238‑3 (MQ=255)                                                                                                                                                                                atttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac                                       <  1:149617/238‑3 (MQ=255)                                                                                                                                                                                atttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac                                       <  1:217141/238‑3 (MQ=255)                                                                                                                                                                                atttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac                                       <  1:578174/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)                                                                                                                                                                                           atgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAAt                                                                                                                                                                                   <  1:473155/100‑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)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  AAGAACAGATTCCGCCTTATCAAATCGAATCGGTTTATGGCTTGCTGACCTTGGCGCATTACCGAACCTCGCCGTTGGATTTACGTCGATTATTGGATGCGCCGCGACAACAATTTTATCTTGCGCAGACGCAGGAGGCCTTGTTGGGTTGTGTGTGGGCTGTGCCTGAGGGCGAATTTGGTGACGACGCATTAATTCGTAAAATTCGGCGCGGACAACGTCGTCCACGTGGTAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCACTGCGCATTTCCCGTATCGCCGTGCAACCCAACTGGCAGCAGCGGGGACTGGGGCAACGTTTAATCGAA‑CGAG‑TTAAACAACAACGTGCGGTGGATTTCCTTTCCGTCAGTTTCGGTTACACCCCCGAGCTGCTGGCGTTTTGGCAAAAGTGCGGTTTTATTTTGGT  >  1/887357‑887819

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