Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,475,880	T→C	Q66Q (CAA→CAG)	unknown ←	tRNA cytosine(34) acetyltransferase (EC 2.3.1.193)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,475,880	0	T	C	100.0%	172.2 / NA	63	Q66Q (CAA→CAG)	unknown	tRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/42);  total (21/42)

TTAGGGAAAACGAGCGCAGACGGCAAATGTTGATAAATGGGGAAGCTGTATTTCCGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACTGCTCCAACGCAGGCTGTCGCCGTCGGGTTGTGCCGCCAAATCCTGCCACTCATTGAATAACAGCAACAATGTGCCGCCTTCCGACAACGCACCGGCGGCAATAGCGAGGGCATCTAAATGAATGCCTTGGCGCCCATCATAAATGAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAATGAAGTTGAGCGGAAAAGTGCGGCAAAACACCGCCAATCCACACCGCACTTTTCAGGGCAGGAAAGTATTCCCATTGGCGCGCTAACGTGGCATCGTCGCCAACAAGGAGAATTAGCTGCCGCGCCATGCTTTACCAACCCCGCGAATACGGATCCATCC  >  1/1475641‑1476109                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                      ctAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCcat                                                                                                                                                                                                                             >  1:703132/2‑250 (MQ=255)          accagcgcgggcGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAAta                                                                                                                                                                                                                     <  1:466534/238‑2 (MQ=255)          accagcgcgggcGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAAta                                                                                                                                                                                                                     <  1:467746/238‑2 (MQ=255)          accagcgcgggcGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAAta                                                                                                                                                                                                                     <  1:51859/238‑2 (MQ=255)          accagcgcgggcGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAAta                                                                                                                                                                                                                     <  1:734236/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)                                     aTGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGAtt                                                                                                                                                                                        >  1:36715/1‑251 (MQ=255)                                                                           ggCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGaaa                                                                                                                                                  >  1:163458/1‑251 (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)                                                                               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)                                                                                                                                tgtcaccatcaggttgcaccgccAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGaaaatgt                                                                                             <  1:309866/233‑5 (MQ=255)                                                                                                                                    cccatcaggttgcaccgccAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTccc                                                                                         <  1:253168/237‑1 (MQ=255)                                                                                                                                     ccatcaggttgcaccgccAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGAtttt                                                                                                                                                                                      >  1:485204/14‑157 (MQ=255)                                                                                                                                     ccatcaggttgcaccgccAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCa                                                                                        >  1:63206/14‑251 (MQ=255)                                                                                                                                           ggttgcaccgccAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAACAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCAc                                                                                                                 >  1:893088/8‑220 (MQ=255)                                                                                                                                            gttgcaccgccAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTgg                                                                                    <  1:756796/242‑1 (MQ=255)                                                                                                                                            gttgcaccgccAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTgg                                                                                    <  1:2743/242‑1 (MQ=255)                                                                                                                                                    gccAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGt                                                                         >  1:29227/1‑251 (MQ=255)                                                                                                                                                     ccAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGt                                                                         <  1:630816/250‑1 (MQ=255)                                                                                                                                                     ccAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGt                                                                         <  1:59991/250‑1 (MQ=255)                                                                                                                                                                  atCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACcaa                                                           <  1:307769/250‑1 (MQ=255)                                                                                                                                                                   tCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACcaat                                                          >  1:867890/1‑250 (MQ=255)                                                                                                                                                                   tCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACcaa                                                           >  1:14211/1‑250 (MQ=255)                                                                                                                                                                      ttGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCggtga                                                                                                                                        <  1:382846/170‑4 (MQ=255)                                                                                                                                                                        gAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACcaatgagaa                                                     >  1:104009/1‑245 (MQ=255)                                                                                                                                                                          aTAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACcaatgagaaga                                                   <  1:726572/251‑9 (MQ=255)                                                                                                                                                                                 caacaaTATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTaact                                            <  1:870287/251‑4 (MQ=255)                                                                                                                                                                                 caacaaTATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTaact                                            <  1:470535/251‑4 (MQ=255)                                                                                                                                                                                 caacaaTATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTaa                                              >  1:56150/1‑248 (MQ=255)                                                                                                                                                                                   acaaTATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCAtcgtcacc                                                             <  1:411686/232‑4 (MQ=255)                                                                                                                                                                                   acaaTATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTaac                                             >  1:271834/1‑246 (MQ=255)                                                                                                                                                                                   acaaTATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTg                                           >  1:300262/1‑250 (MQ=255)                                                                                                                                                                                   acaaTATGCCGCCTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTg                                           >  1:193059/1‑250 (MQ=255)                                                                                                                                                                                              ccTTCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCTGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCAtag                                <  1:212219/250‑3 (MQ=255)                                                                                                                                                                                                ttCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCAtagttt                             >  1:248402/1‑246 (MQ=255)                                                                                                                                                                                                 tCCGACAGCGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCcat                                                                                                                        <  1:287215/159‑2 (MQ=255)                                                                                                                                                                                                     acagcGCACCGGCGGCAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAggg                                                                                                        <  1:71496/167‑1 (MQ=255)                                                                                                                                                                                                                    cAATGGCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGa         <  1:179552/251‑1 (MQ=255)                                                                                                                                                                                                                      atggCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGatc       <  1:238826/248‑1 (MQ=255)                                                                                                                                                                                                                        ggCGAGAGCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCAtagt                               <  1:592869/224‑4 (MQ=255)                                                                                                                                                                                                                             gagCATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGatccgtcc  <  1:28446/247‑1 (MQ=255)                                                                                                                                                                                                                                cATCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGatc       <  1:301001/241‑1 (MQ=255)                                                                                                                                                                                                                                 aTCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGc                                                                                                       <  1:27031/144‑1 (MQ=255)                                                                                                                                                                                                                                 aTCTAAATGAATGCCCTGACGCCCATCATAAATAAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCAtagttt                             <  1:492585/218‑6 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                      TTAGGGAAAACGAGCGCAGACGGCAAATGTTGATAAATGGGGAAGCTGTATTTCCGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACTGCTCCAACGCAGGCTGTCGCCGTCGGGTTGTGCCGCCAAATCCTGCCACTCATTGAATAACAGCAACAATGTGCCGCCTTCCGACAACGCACCGGCGGCAATAGCGAGGGCATCTAAATGAATGCCTTGGCGCCCATCATAAATGAGCGCGGGGAATTCGCTGCCGAGCAGATTTTTACTTTTGCTAAAAGGGAAATGAAGTTGAGCGGAAAAGTGCGGCAAAACACCGCCAATCCACACCGCACTTTTCAGGGCAGGAAAGTATTCCCATTGGCGCGCTAACGTGGCATCGTCGCCAACAAGGAGAATTAGCTGCCGCGCCATGCTTTACCAACCCCGCGAATACGGATCCATCC  >  1/1475641‑1476109

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