Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	887,541	C→T	D412D (GAC→GAT)	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,541	0	C	T	100.0%	171.1 / NA	60	D412D (GAC→GAT)	unknown	tRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/26);  total (34/26)

TGACGCAACCGCGTTACGATGCCAACCTTGCCAAACAGATTCAGATTTCTTACGAAGAACAGATTCCGCCTTATCAAATCGAATCGGTTTATGGCTTGCTGACCTTGGCGCATTACCGAACCTCGCCGTTGGATTTACGTCGATTATTGGATGCGCCGCGACAACAATTTTATCTTGCGCAGACGCAGGAGGCCTTGTTGGGTTGTGTGTGGGCTGTGCCTGAGGGCGAATTTGGTGACGACGCATTAATTCGTAAAATTCGGCGCGGACAACGTCGTCCACGTGGTAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCACTGCGCATTTCCCGTATCGCCGTGCAACCCAACTGGCAGCAGCGGGGACTGGGGCAACGTTTAATCGAA‑CGAG‑TTAAACAACAACGTGCGGTGGATTTCCTTTCCGT  >  1/887303‑887759                                                                                                                                                                                                                                               |                                                                                                                                                                                                                             tgccGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAAtt                                                                                                                                                                                                                  <  1:819051/248‑1 (MQ=255)  gccGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAAtt                                                                                                                                                                                                                  >  1:408787/3‑250 (MQ=255)           gtgttatgacgctaaCCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGCTTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAAtt                                                                                                                                                                                                         >  1:890573/14‑250 (MQ=255)           gtgttatgacgctaaCCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAAtt                                                                                                                                                                                                         >  1:283929/14‑250 (MQ=255)             gttatgacgctaaCCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCga                                                                                                                                                                                                      >  1:177839/12‑250 (MQ=255)                                            gattTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAAtt                                                                                                                                                                                                                  >  1:77889/1‑208 (MQ=255)                                             attTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAAcgtcgcctgcgcggcaatttgg                                                                                                                                                                       <  1:130677/250‑18 (MQ=255)                                             attTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAAcgtcgcctgcgcggcaatttgg                                                                                                                                                                       <  1:28453/250‑18 (MQ=255)                                                     caaagaaCGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTc                                                                                                                                                                                                        >  1:803115/3‑209 (MQ=255)                                                       aagaaCGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAAtg                                                                                                                                                            <  1:798858/251‑1 (MQ=255)                                                         gaaCGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAAtt                                                                                                                                                                                                                  <  1:775515/195‑1 (MQ=255)                                                          aacggATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAAtgttg                                                                                                                                                         <  1:809807/247‑1 (MQ=255)                                                            cggATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTtgtg                                                                                                                                                       <  1:390128/249‑1 (MQ=255)                                                                                  aatcgGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAAcgtcgcctg                                                                                                                                                                                    >  1:114989/1‑196 (MQ=255)                                                                                     cgGTTTATGGCTTGCTGACTTTGGCGCATTACCTCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGGCGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCg                                                                                                                              <  1:506287/251‑1 (MQ=255)                                                                                      ggTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAAcgtc                                                                                                                                                                                         >  1:739850/1‑191 (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)                                                                                                                 ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGAc                                                                                                                                                                                               <  1:97440/158‑1 (MQ=255)                                                                                                                 ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTtgtg                                                                                                                                                       <  1:701921/198‑1 (MQ=255)                                                                                                                 ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAAc                                                                                                                                   >  1:728742/1‑218 (MQ=255)                                                                                                                 ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTcc                                                                                                   >  1:209876/1‑250 (MQ=255)                                                                                                                 ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTcc                                                                                                   >  1:51847/1‑250 (MQ=255)                                                                                                                 ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTcc                                                                                                   <  1:582402/250‑1 (MQ=255)                                                                                                                 ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTcc                                                                                                   >  1:698893/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:308059/2‑250 (MQ=255)                                                                                                                                                                                                   ctgttgGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTggg              >  1:596187/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)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                             TGACGCAACCGCGTTACGATGCCAACCTTGCCAAACAGATTCAGATTTCTTACGAAGAACAGATTCCGCCTTATCAAATCGAATCGGTTTATGGCTTGCTGACCTTGGCGCATTACCGAACCTCGCCGTTGGATTTACGTCGATTATTGGATGCGCCGCGACAACAATTTTATCTTGCGCAGACGCAGGAGGCCTTGTTGGGTTGTGTGTGGGCTGTGCCTGAGGGCGAATTTGGTGACGACGCATTAATTCGTAAAATTCGGCGCGGACAACGTCGTCCACGTGGTAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCACTGCGCATTTCCCGTATCGCCGTGCAACCCAACTGGCAGCAGCGGGGACTGGGGCAACGTTTAATCGAA‑CGAG‑TTAAACAACAACGTGCGGTGGATTTCCTTTCCGT  >  1/887303‑887759

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