Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,475,772	G→A	G102G (GGC→GGT)	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	1,475,772	0	G	A	100.0%	171.0 / NA	64	G102G (GGC→GGT)	unknown	tRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/42);  total (22/42)

ACGATTAACATCTCCGCGTCAGATGTCGCCATTTGCCGCAATAAATGGGCTTGTTCCAAGGTCGGTTGGTGTTGTTCTGAAAGCGACTTTTCAATCTTAGGGAAAACGAGCGCAGACGGCAAATGTTGATAAATGGGGAAGCTGTATTTCCGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACTGCTCCAACGCAGGCTGTCGCCGTCGGGTTGTGCCGCCAAATCCTGCCACTCATTGAATAACAGCAACAATGTGCCGCCTTCCGACAACGCACCGGCGGCAATAGCGAGGGCATCTAAATGAATGCCTTGGCGCCCATCATAAATGAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAATGAAGTTGAGCGGAAAAGTGCGGCAAAACACCGCCAATCCACAC  >  1/1475545‑1475993                                                                                                                                                                                                                                    |                                                                                                                                                                                                                              aCGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATc                                                                                                                                                                                                        <  1:780503/251‑1 (MQ=255)   gATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATc                                                                                                                                                                                                        >  1:707442/1‑249 (MQ=255)                    cggATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAAc                                                                                                                                                                                     <  1:775271/249‑1 (MQ=255)                       aTGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAg                                                                                                                                                                                   <  1:661606/250‑1 (MQ=255)                                  tGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATgcc                                                                                                                                                                       <  1:195416/251‑1 (MQ=255)                                       cAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCtt                                                                                                                                                                  >  1:581539/1‑251 (MQ=255)                                         aTAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTc                                                                                                                                                                 >  1:168614/1‑250 (MQ=255)                                           aaaTGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGcaacaa                                                                                                                                                                             >  1:181389/1‑236 (MQ=255)                                           aaaTGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATgc                                                                                                                                                                        >  1:755147/1‑241 (MQ=255)                                           aaaTGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGa                                                                                                                                                              >  1:58554/1‑251 (MQ=255)                                            aaTGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGAc                                                                                                                                                             <  1:408960/251‑1 (MQ=255)                                            aaTGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGAc                                                                                                                                                             <  1:380027/251‑1 (MQ=255)                                                ggtttGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACagcg                                                                                                                                                         >  1:200952/4‑248 (MQ=255)                                                     gTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAAc                                                                                                                                                                                     <  1:49114/218‑1 (MQ=255)                                                     gTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACcg                                                                                                                                                    >  1:840272/1‑251 (MQ=255)                                                              tCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAAtg                                                                                                                                           <  1:49567/251‑2 (MQ=255)                                                              tCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAAt                                                                                                                                            <  1:578240/250‑1 (MQ=255)                                                               cGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTTTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAAtgg                                                                                                                                          <  1:739493/251‑3 (MQ=255)                                                                 gTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGAc                                                                                                                                                             <  1:120324/230‑1 (MQ=255)                                                                 gTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACcggcg                                                                                                                                                 <  1:182502/242‑1 (MQ=255)                                                                 gTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCatattg                                                                                                                                         <  1:650412/250‑6 (MQ=255)                                                                      tgtagttCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAgagc                                                                                                                                   <  1:320599/247‑4 (MQ=255)                                                                              tGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAaat                                                                                                                           <  1:7795/251‑1 (MQ=255)                                                                                aaaGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAaatg                                                                                                                          <  1:373382/250‑1 (MQ=255)                                                                                aaaGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAaatg                                                                                                                          <  1:726986/250‑1 (MQ=255)                                                                                aaaGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAaatg                                                                                                                          <  1:132686/250‑1 (MQ=255)                                                                                      acgtttCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGccc                                                                                                                   <  1:722041/248‑2 (MQ=255)                                                                                         tttCAATCCTAGGGAAGGCCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAaatg                                                                                                                          <  1:774487/241‑1 (MQ=255)                                                                                         tttCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAAc                                                                                                                                                                                     <  1:263503/182‑1 (MQ=255)                                                                                            cAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGccctgacgc                                                                                                             <  1:7940/251‑8 (MQ=255)                                                                                                 ctAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCcat                                                                                                         >  1:703132/2‑250 (MQ=255)                                                                                                         gaccagcgcgggcGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCGACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAaatg                                                                                                                          >  1:424225/13‑225 (MQ=255)                                                                                                          accagcgcgggcGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAAta                                                                                                 <  1:466534/238‑2 (MQ=255)                                                                                                          accagcgcgggcGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAAta                                                                                                 <  1:467746/238‑2 (MQ=255)                                                                                                          accagcgcgggcGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAAta                                                                                                 <  1:734236/238‑2 (MQ=255)                                                                                                          accagcgcgggcGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAAta                                                                                                 <  1:51859/238‑2 (MQ=255)                                                                                                          accagcgcgggcGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAAta                                                                                                 <  1:545766/238‑2 (MQ=255)                                                                                                          acaagcgcgggcGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGACAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAAta                                                                                                 <  1:582163/238‑2 (MQ=255)                                                                                                               cgcgggcGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAgcgc                                                                                            <  1:591353/243‑1 (MQ=255)                                                                                                                 cgggcGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCgggg                                                                                        <  1:802091/247‑1 (MQ=255)                                                                                                                           atattGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGc                                                                                 <  1:473493/245‑1 (MQ=255)                                                                                                                                 aTAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATgcc                                                                                                                                                                       >  1:245815/1‑156 (MQ=255)                                                                                                                                     aTGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGAtt                                                                    >  1:36715/1‑251 (MQ=255)                                                                                                                                                                           ggCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAaatg                                                                                                                          >  1:116375/1‑159 (MQ=255)                                                                                                                                                                           ggCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGaaa                              >  1:163458/1‑251 (MQ=255)                                                                                                                                                                            gCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGAc                                                                                                                                                             >  1:83046/1‑123 (MQ=255)                                                                                                                                                                            gCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAgcgc                                                                                            <  1:335624/188‑1 (MQ=255)                                                                                                                                                                            gCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGaaaa                             >  1:25093/1‑251 (MQ=255)                                                                                                                                                                            gCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGaaaa                             >  1:322993/1‑251 (MQ=255)                                                                                                                                                                            gCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGaaa                              <  1:151475/250‑1 (MQ=255)                                                                                                                                                                            gCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGaaa                              <  1:239994/250‑1 (MQ=255)                                                                                                                                                                            gCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGaaa                              <  1:552108/250‑1 (MQ=255)                                                                                                                                                                            gCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGATCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGaaa                              <  1:546230/250‑1 (MQ=255)                                                                                                                                                                             cAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAg                            <  1:436916/251‑1 (MQ=255)                                                                                                                                                                               aTAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAaatg                                                                                                                          >  1:782028/1‑155 (MQ=255)                                                                                                                                                                               aTAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTg                          >  1:37136/1‑251 (MQ=255)                                                                                                                                                                                         gtgtgtTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCggtgaaacg                >  1:409167/1‑244 (MQ=255)                                                                                                                                                                                             gtTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCggtgaaacgccg             >  1:719897/1‑240 (MQ=255)                                                                                                                                                                                             gtTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCggtgaaacgc               <  1:71257/248‑9 (MQ=255)                                                                                                                                                                                                 aTCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAAt        >  1:642971/1‑251 (MQ=255)                                                                                                                                                                                                  tCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATc       <  1:210987/251‑1 (MQ=255)                                                                                                                                                                                                   cGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATcc      <  1:212676/251‑1 (MQ=255)                                                                                                                                                                                                       tgtTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCcatac  <  1:235089/251‑4 (MQ=255)                                                                                                                                                                                                                  gCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCggtatggattggcggcgtttcac  <  1:335772/240‑22 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                              ACGATTAACATCTCCGCGTCAGATGTCGCCATTTGCCGCAATAAATGGGCTTGTTCCAAGGTCGGTTGGTGTTGTTCTGAAAGCGACTTTTCAATCTTAGGGAAAACGAGCGCAGACGGCAAATGTTGATAAATGGGGAAGCTGTATTTCCGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACTGCTCCAACGCAGGCTGTCGCCGTCGGGTTGTGCCGCCAAATCCTGCCACTCATTGAATAACAGCAACAATGTGCCGCCTTCCGACAACGCACCGGCGGCAATAGCGAGGGCATCTAAATGAATGCCTTGGCGCCCATCATAAATGAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAATGAAGTTGAGCGGAAAAGTGCGGCAAAACACCGCCAATCCACAC  >  1/1475545‑1475993

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