Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	887,374	T→A	Y357N (TAT→AAT)	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,374	0	T	A	100.0%	211.1 / NA	71	Y357N (TAT→AAT)	unknown	tRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (30/41);  total (30/41)

ACGACCACCATTCACAGCTACGAAGGCACGGGCAGGGGCTTCTTGTTAAAATTCATGACGAAAACCGACCGCACTTTGCGGCATTTTGAATTGTTCAGCCCGTTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGATGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGACGCAACCGCGTTACGATGCCAACCTTGCCAAACAGATTCAGATTTCTTACGAAGAACAGATTCCGCCTTATCAAATCGAATCGGTTTATGGCTTGCTGACCTTGGCGCATTACCGAACCTCGCCGTTGGATTTACGTCGATTATTGGATGCGCCGCGACAACAATTTTATCTTGCGCAGACGCAGGAGGCCTTGTTGGGTTGTGTGTGGGCTGTGCCTGAGGGCGAATTTGGTGACGACGCATTAATTCGTAAAATTCGGCGCGGACAACGTCGTCCACGTGGTAATTTGGTGGCGCAAATGTTGTG  >  1/887128‑887612                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               accacgaccaTTCACAGTTACGAGGGTACGGGCAGGGGCTTCTTGTTAAAATTCATAGCGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATc                                                                                                                                                                                                                                             <  1:499629/244‑1 (MQ=255)       accaTTCACAGTTACGAGGGTACGGGCAGGGGCTTCTTGTTAAAATTCATAGCGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAaatcga                                                                                                                                                                                                                                      >  1:103846/1‑251 (MQ=255)               cagttACGAGGGTACGGGCAGGGGCTTCTTGTTAAAATTCATAGCGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGttt                                                                                                                                                                                                                              >  1:807094/5‑251 (MQ=255)                 gttACGAGGGTACGGGCAGGGGCTTCTTGTTAAAATTCATAGCGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTAt                                                                                                                                                                                                                            >  1:212776/3‑251 (MQ=255)                    acgagggtaCGGGCAGGGGCTTCTTGTTAAAATTCATAGCGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGc                                                                                                                                                                                                                         <  1:833517/243‑1 (MQ=255)                         ggtaCGGGCAGGGGCTTCTTGTTAAAATTCATAGCGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGGCGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTg                                                                                                                                                                                                                      >  1:443813/4‑249 (MQ=255)                         ggtaCGGGCAGGGGCTTCTTGTTAAAATTCATAGCGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTg                                                                                                                                                                                                                      >  1:573906/4‑249 (MQ=255)                          gtaCGGGCAGGGGCTTCTTGTTAAAATTCATAGCGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATc                                                                                                                                                                                                                                             <  1:358157/223‑1 (MQ=255)                           taCGGGCAGGGGCTTCTTGTTAAAATTCATAGCGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTg                                                                                                                                                                                                                   <  1:508527/249‑1 (MQ=255)                            aCGGGCAGGGGCTTCTTGTTAAAATTCATAGCGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGAc                                                                                                                                                                                                                 <  1:291446/251‑1 (MQ=255)                                    gggCTTCTTGTTAAAATTCATAGCGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGgc                                                                                                                                                                                                           <  1:810904/249‑1 (MQ=255)                                            tGTTAAAATTCATAGCGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACTTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGATTTATGGCTTGCTGACTTTGGCGCATTACCgc                                                                                                                                                                                                 <  1:386945/251‑2 (MQ=255)                                            tGTTAAAATTCATAGCGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCg                                                                                                                                                                                                  >  1:480086/1‑250 (MQ=255)                                             gTTAAAATTCATAGCGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCgca                                                                                                                                                                                                >  1:326311/1‑249 (MQ=255)                                              ttAAAATTCATAGCGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCgcac                                                                                                                                                                                               <  1:75896/251‑4 (MQ=255)                                              ttAAAATTCATAGCGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCgcac                                                                                                                                                                                               <  1:140583/251‑4 (MQ=255)                                                         agcGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGttt                                                                                                                                                                                                                              <  1:90505/207‑1 (MQ=255)                                                         agcGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTgg                                                                                                                                                                                    >  1:624709/3‑251 (MQ=255)                                                         agcGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTgg                                                                                                                                                                                    >  1:151873/3‑251 (MQ=255)                                                           cGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGAt                                                                                                                                                                                  <  1:831458/251‑1 (MQ=255)                                                           cGAACACCGACCGCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGAt                                                                                                                                                                                  <  1:12318/251‑1 (MQ=255)                                                                 ccgaccgCACTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGt                                                                                                                                                                            <  1:888508/251‑1 (MQ=255)                                                                          cTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGc                                                                                                                                                                                                                         <  1:203440/197‑1 (MQ=255)                                                                          cTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGc                                                                                                                                                                                                                         >  1:876514/1‑197 (MQ=255)                                                                          cTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCgcac                                                                                                                                                                                               >  1:514904/1‑220 (MQ=255)                                                                          cTTTTCGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTg                                                                                                                                                                   >  1:823289/1‑251 (MQ=255)                                                                            tttcGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGGGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGa                                                                                                                                                                 <  1:445058/248‑1 (MQ=255)                                                                             ttcGGCATTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGAt                                                                                                                                                                <  1:477740/249‑1 (MQ=255)                                                                                   aTTTTGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATgg                                                                                                                                                                                                                          <  1:881742/187‑1 (MQ=255)                                                                                     tttGAATTGTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACcgcg                                                                                                                                                        >  1:437937/1‑251 (MQ=255)                                                                                             gTTCAGCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGAt                                                                                                                                                                >  1:535960/1‑235 (MQ=255)                                                                                                  gCCCGCTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTAt                                                                                                                                           >  1:229495/1‑251 (MQ=255)                                                                                                       ctACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTg                                                                                                                                       >  1:11101/2‑250 (MQ=255)                                                                                                                   gCGGATGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTAt                                                                                                                                           >  1:284801/1‑234 (MQ=255)                                                                                                                        tGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCgggaaacc                                                                                                                      >  1:428935/1‑246 (MQ=255)                                                                                                                          aCAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATgcacc                                                                                                                                                           <  1:327108/211‑4 (MQ=255)                                                                                                                             aaCTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCgggaaaccctgttg                                                                                                                <  1:604455/251‑11 (MQ=255)                                                                                                                                tGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTgggc                                                                                                             <  1:367894/251‑2 (MQ=255)                                                                                                                                tGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTgggc                                                                                                             <  1:527982/251‑2 (MQ=255)                                                                                                                                       cTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCgcac                                                                                                                                                                                               <  1:189234/162‑4 (MQ=255)                                                                                                                                       cTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTgggc                                                                                                             <  1:644225/244‑2 (MQ=255)                                                                                                                                       cTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCtgtgtgt                                                                                                      <  1:39237/251‑1 (MQ=255)                                                                                                                                       cTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCtgtgtg                                                                                                       <  1:420895/250‑1 (MQ=255)                                                                                                                                          tATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTggg                                                                                                   <  1:278758/251‑1 (MQ=255)                                                                                                                                             tgacgaATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGcgg                                                                                                <  1:607949/247‑3 (MQ=255)                                                                                                                                              gacgaATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTgggctg                                                                                                           >  1:52133/4‑236 (MQ=255)                                                                                                                                                     tgttgttgtGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGaa                                                                                        <  1:175357/251‑2 (MQ=255)                                                                                                                                                       ttgttgtGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGaag                                                                                       <  1:619874/250‑3 (MQ=255)                                                                                                                                                            gtgtgATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGaaggaggat                                                                                 >  1:406502/1‑243 (MQ=255)                                                                                                                                                                  ttGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGaag                                                                                       >  1:342096/1‑237 (MQ=255)                                                                                                                                                                     cGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCTTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGaaggaggatttgaggat                                                                         <  1:443662/250‑17 (MQ=255)                                                                                                                                                                     cGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGaaggaggatttgaggat                                                                         <  1:674652/250‑17 (MQ=255)                                                                                                                                                                          aTCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGaag                                                                                       >  1:823891/1‑229 (MQ=255)                                                                                                                                                                                tgccGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAAtt                                                             <  1:819051/248‑1 (MQ=255)                                                                                                                                                                                 gccGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAAtt                                                             >  1:408787/3‑250 (MQ=255)                                                                                                                                                                                          gtgttatgacgctaaCCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGCTTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAAtt                                                    >  1:890573/14‑250 (MQ=255)                                                                                                                                                                                          gtgttatgacgctaaCCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGaa                                                                                        <  1:252535/201‑2 (MQ=255)                                                                                                                                                                                          gtgttatgacgctaaCCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAAtt                                                    >  1:283929/14‑250 (MQ=255)                                                                                                                                                                                            gttatgacgctaaCCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCga                                                 >  1:177839/12‑250 (MQ=255)                                                                                                                                                                                                  acgctaaCCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCtgtgtg                                                                                                       >  1:290264/6‑191 (MQ=255)                                                                                                                                                                                                  acgctaaCCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGaaggaggatttg                                                                              <  1:695858/211‑12 (MQ=255)                                                                                                                                                                                                                     gattcagattTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTgggctg                                                                                                           <  1:735718/168‑4 (MQ=255)                                                                                                                                                                                                                           gattTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAAtt                                                             >  1:77889/1‑208 (MQ=255)                                                                                                                                                                                                                            attTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAAcgtcgcctgcgcggcaatttgg                  <  1:28453/250‑18 (MQ=255)                                                                                                                                                                                                                            attTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAAcgtcgcctgcgcggcaatttgg                  <  1:130677/250‑18 (MQ=255)                                                                                                                                                                                                                                   acaaagaaCGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTgggc                                                                                                             >  1:81350/4‑151 (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)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               ACGACCACCATTCACAGCTACGAAGGCACGGGCAGGGGCTTCTTGTTAAAATTCATGACGAAAACCGACCGCACTTTGCGGCATTTTGAATTGTTCAGCCCGTTACGTTGGCAGGCGGATGACAAACTGGAGCGCTTTATTGATGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGACGCAACCGCGTTACGATGCCAACCTTGCCAAACAGATTCAGATTTCTTACGAAGAACAGATTCCGCCTTATCAAATCGAATCGGTTTATGGCTTGCTGACCTTGGCGCATTACCGAACCTCGCCGTTGGATTTACGTCGATTATTGGATGCGCCGCGACAACAATTTTATCTTGCGCAGACGCAGGAGGCCTTGTTGGGTTGTGTGTGGGCTGTGCCTGAGGGCGAATTTGGTGACGACGCATTAATTCGTAAAATTCGGCGCGGACAACGTCGTCCACGTGGTAATTTGGTGGCGCAAATGTTGTG  >  1/887128‑887612

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