Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,216,473	A→G	Y7H (TAC→CAC)	nfeF ←	NADPH‑dependent ferric‑chelate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,216,473	0	A	G	100.0%	238.6 / NA	84	Y7H (TAC→CAC)	nfeF	NADPH‑dependent ferric‑chelate reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/54);  total (30/54)

ATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTG‑TGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCA  >  NC_000913/3216228‑3216721                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                           ataGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACACAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCggg                                                                                                                                                                                                                                                       >  2:423583/1‑251 (MQ=255)      cACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGgt                                                                                                                                                                                                                                                   <  1:94817/250‑1 (MQ=255)       aCGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGgtgt                                                                                                                                                                                                                                                 <  1:487030/251‑1 (MQ=255)        cgcgACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGtt                                                                                                                                                                                                                                                <  1:235742/251‑1 (MQ=255)            ccGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTAATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTc                                                                                                                                                                                                                                            <  1:506001/250‑1 (MQ=255)              ggtgggcgtgtgcagtCCGGCCAGACGATGCCCTCTTCCGTTACCGGTGGCGGCACAAAGTGAGCGTCAGGTTGATGAAAGAAGAGTTTGCTGTGATCTTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTTAAAACCGGCGCTGATGCGCTCAACGCGTGACACAGTCAGTTCACGGAAGCGCATATCATTGCGAACGCGCTGAGGGTGGCGGGTGGTGGTATTCAt                                                                                                                                                                                                                                          <  2:250051/236‑1 (MQ=255)              ggTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCAt                                                                                                                                                                                                                                          >  1:131804/1‑251 (MQ=255)               gTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACAGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGCGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCAtt                                                                                                                                                                                                                                         >  1:97184/1‑251 (MQ=255)                tGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTg                                                                                                                                                                                                                                        >  1:577721/1‑251 (MQ=255)                 gggCGTTGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTTCGAACGCGCTGCGGGGGGCGGGGGGTGTTATTCATTGt                                                                                                                                                                                                                                       <  1:537334/251‑1 (MQ=255)                 gggCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCTTTGt                                                                                                                                                                                                                                       >  2:288058/1‑251 (MQ=255)                  ggCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGtt                                                                                                                                                                                                                                      <  1:577844/251‑1 (MQ=255)                     gTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATc                                                                                                                                                                                                                                   >  1:662577/1‑251 (MQ=255)                       ggTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCCACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGACCGCGCTGCGGGTGGCGGGGGGGGGTTTTCATTGTTTTCGc                                                                                                                                                                                                                                 >  2:4006/1‑251 (MQ=255)                       ggTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCg                                                                                                                                                                                                                                  <  1:459103/250‑1 (MQ=255)                        gTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGcc                                                                                                                                                                                                                                >  2:312198/1‑251 (MQ=255)                        gTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTAGTCATTGTTATCGcc                                                                                                                                                                                                                                <  1:176834/251‑1 (MQ=255)                        gTCCTTCCGGCCAGACGATGCCCTCTGCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGcc                                                                                                                                                                                                                                <  2:446406/251‑1 (MQ=255)                         cccTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCATATCATTGCGATCGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCt                                                                                                                                                                                                                               <  1:588854/250‑1 (MQ=255)                          ccTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGTCGGGGGGGGGTATTCATTGTTATCGCCtt                                                                                                                                                                                                                              <  2:608601/251‑1 (MQ=255)                          ccTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCt                                                                                                                                                                                                                               <  1:540032/250‑1 (MQ=255)                          ccTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCGCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGGTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCt                                                                                                                                                                                                                               <  2:387738/250‑1 (MQ=255)                          ccTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGGGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCtt                                                                                                                                                                                                                              <  2:221145/251‑1 (MQ=255)                             tCCGGCCAGACGATGCCCTCTTTCGTTACCGGTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATTTAAAACCTTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACTGAAGCGCAGATCAGTTCGAACGCGCTCCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTc                                                                                                                                                                                                                             <  1:646141/249‑1 (MQ=255)                             tCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTc                                                                                                                                                                                                                             <  1:481097/249‑1 (MQ=255)                                ggCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACCAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTCTCGCCTTCGTGAt                                                                                                                                                                                                                        >  2:544458/1‑251 (MQ=255)                                ggCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACACAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCCGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAt                                                                                                                                                                                                                        >  1:388021/1‑251 (MQ=255)                                   cAGACGAAGACCTCGTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGGCAAAGCCACGCGATGTAAAACCGTCCCGCGCCTCGCCGCCGAGGGCAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGGAACACAGTCAGTTCACGGAGGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGt                                                                                                                                                                                                                    <  2:410549/251‑1 (MQ=255)                                    aGACGATGCCCTCTTCCGTGACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTa                                                                                                                                                                                                                   <  2:27302/251‑1 (MQ=255)                                      aCGAAGCCCTCTTCCGTTACCGTTGGCGGCACAAATTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGGAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAAt                                                                                                                                                                                                                 <  1:20063/251‑1 (MQ=255)                                       cgaggCCCTCGTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTAGTCATTGTTAGCGCCTTCTTGAT‑GGTAATc                                                                                                                                                                                                                <  2:653629/247‑1 (MQ=255)                                                ctTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATc                                                                                                                                                                                                       >  1:330748/1‑251 (MQ=255)                                                        aCCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCtaaataaa                                                                                                                                                                                               <  1:645954/251‑1 (MQ=255)                                                        aCCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCtaaataa                                                                                                                                                                                                <  1:304180/250‑1 (MQ=255)                                                               gcggcACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGc                                                                                                                                                                                         >  1:100037/1‑250 (MQ=255)                                                                      aaaGTGAGCGTCAGGTTGCGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCa                                                                                                                                                                                                               >  1:474622/1‑221 (MQ=255)                                                                      aaaGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAaatgata                                                                                                                                                                                 <  1:81332/251‑1 (MQ=255)                                                                      aaaGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTAGCGCCTTCGTGAT‑GGTAATCa                                                                                                                                                                                                               <  2:474622/221‑1 (MQ=255)                                                                           gAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAaatgataatgat                                                                                                                                                                            >  2:436816/1‑251 (MQ=255)                                                                             gCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCTCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCTAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GTTAAGCAGATATATCTAAATAAAACTCGCAAATGATAATGATTg                                                                                                                                                                          <  2:63868/251‑1 (MQ=255)                                                                              cGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGt                                                                                                                                                                         <  2:458957/251‑1 (MQ=255)                                                                                    gTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCATTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGTCTTCGTGAT‑TGTAATCATATAGAGCTAAATAAAACTCGCAAATGATAATGATTGTTAATCa                                                                                                                                                                   <  1:608466/251‑1 (MQ=255)                                                                                    gTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATc                                                                                                                                                                    >  1:510045/1‑250 (MQ=255)                                                                                            aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAAt                                                                                                                                                           >  1:525119/1‑251 (MQ=255)                                                                                            aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAAt                                                                                                                                                           <  1:423583/251‑1 (MQ=255)                                                                                                  ggTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGc                                                                                                                                                     <  2:169900/250‑1 (MQ=255)                                                                                                  agTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGAGGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCATATATATCGAAATAAAACTCGCAAATGAGAATGAGTGTTAATCATTATAAATGCAAGc                                                                                                                                                     <  2:401325/251‑1 (MQ=255)                                                                                                                gTCAAAGCCACGCGATGTAAACCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACGCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTg                                                                                                                                       <  2:52188/251‑1 (MQ=255)                                                                                                                gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGCTCATTGCGAACGCGCTGCGGGTGGCGGGGGGCGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTg                                                                                                                                       >  1:465247/1‑251 (MQ=255)                                                                                                                  cAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAAGGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGat                                                                                                                                     <  2:299495/251‑1 (MQ=255)                                                                                                                   aaaGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGata                                                                                                                                    <  2:149508/251‑1 (MQ=255)                                                                                                                   aaaGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGat                                                                                                                                     >  2:126038/1‑250 (MQ=255)                                                                                                                    aaGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑TGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATg                                                                                                                                  <  1:562966/251‑1 (MQ=255)                                                                                                                    aaGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCGCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAGCGCGCTGCGGGTGGCGGGTGGTGTTATTCATTGTTATCGCGTGCGTTAG‑GGTAATCAGATATAGCTAATTAGAACTCGCTGATGTTAGTGAGGGTTGATCATGATAAATGCAAGCGATTTGTAGAACTGatat                                                                                                                                   <  1:637626/251‑1 (MQ=255)                                                                                                                        cACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACGGATATGTCt                                                                                                                               <  2:408893/251‑1 (MQ=255)                                                                                                                         aCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCACCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCta                                                                                                                              <  1:629609/251‑1 (MQ=255)                                                                                                                            cgATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAg                                                                                                                           <  1:232743/251‑1 (MQ=255)                                                                                                                                gTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTg                                                                                                                       >  2:149855/1‑251 (MQ=255)                                                                                                                                     aCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGGAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATa                                                                                                                    >  1:420758/1‑249 (MQ=255)                                                                                                                                          ccAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTAGCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTTTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGa                                                                                                              <  1:348481/250‑1 (MQ=255)                                                                                                                                            aGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAAcc                                                                                                           >  1:230965/1‑251 (MQ=255)                                                                                                                                             gcgcCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCg                                                                                                          <  1:436816/251‑1 (MQ=255)                                                                                                                                                               aagggcgcTAAAAACCGGCGTTGAGGCGCCCAACGCGTAACACAGTAAGTTCCGGGAAGCGCAGATCATTGGGACCGCGCTGGGGGGGGCGGGGGGTTTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGAAATATCAAAATAAACCCCGCAAATGATAAGGTTTGTTAATCATGAAAAAGGCAAGCGATTTGTAGAACTGAAATGTTAATAGTCTAAAAAGACGAACCGCCTCTTCTCAGGcatcat                                                                                        <  2:662577/247‑1 (MQ=255)                                                                                                                                                                aaggcgcTAAAAACCGGCGCTGAGGCGCACAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGACCGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTGCGTGAT‑GGTAATCAAATATATCAAAAAAAACCTCGCAAATGATAATGATTGTTAATCAGGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGCCTGATAAGACGAACCGCCTCTTCTCAGGCATCAtt                                                                                       <  2:100037/248‑1 (MQ=255)                                                                                                                                                                aaTGCGCTGAAAACCGGAGCTGATGCGCTCAACGCGTAACAAAGTCCGTTCACGGAAGCGCAGATCATTCCGAACGCGCTGCGCGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTAGTGAT‑GGTAATCAGATATATCTCAATAAAACTAGAACATGCTAATGCTTGTTAATAATGATAAATGCAAGCGCTTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCCTCAtt                                                                                       >  1:592061/1‑251 (MQ=255)                                                                                                                                                                           aacgggCGCGGATGCGCTCAACGGGTAACCCAGTCAGTTCGCGGAGGCGAAGATCATTGGGAACGCGCTGGGGGTGGGGGGGGGTGTTATTAATTGTTACCGCCTGCGTGAT‑GGTAATGAGATATATCTAAAAAAACCTCGAAAAGGATAATGATTGTTAATCATGATAAAGGCAGCCGATTTGTAGAACGGATATGTCTAGAGGCTGATAAGACGAGCCGCCTCTTCTCAGGCATCATTACTCAACGCCg                                                                            <  2:493142/247‑1 (MQ=255)                                                                                                                                                                                    ggATGCGCTCAACGCGTAACACGGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCAAAATAAAACTCGCAAATAATAATGATTGTTAATCATAATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATgcggcg                                                                   <  2:230965/250‑1 (MQ=255)                                                                                                                                                                                          gcTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAAcc                                                                                                           <  2:296276/205‑1 (MQ=255)                                                                                                                                                                                          gcTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAAcc                                                                                                           >  1:296276/1‑205 (MQ=255)                                                                                                                                                                                                cgcgTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCGGAAAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATgcggc                                                                    <  1:242085/238‑1 (MQ=255)                                                                                                                                                                                                cgcgTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCCTTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATgcggc                                                                    >  2:242085/1‑238 (MQ=255)                                                                                                                                                                                                 gcgTAACACCGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATc                                                      >  1:491094/1‑251 (MQ=255)                                                                                                                                                                                                 gagTAACACGGCCAGTCCCCGGAAGCGCAGTTCATGGCGAACGCGCGGGGGGTGGGGGGGGTTTTTTTTCATTTTTATCGCCTTCGTGAG‑GGTAATCAGATATATCAAAAAAAAACTCGAAAATGATAATGATGGTTAAACATGATAAAGGCAAGCGTTTTGTAAACCGGATTGGTCAATAGTCTGAAAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATc                                                      <  2:606521/249‑1 (MQ=255)                                                                                                                                                                                                          agtcagtTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACg                                             >  1:105969/1‑251 (MQ=255)                                                                                                                                                                                                           gtgatttCACGGAAGCCCAGATCTTGGCTACCGCGCGGCGGGTGGGGGGGGGTGTTATTCATTGTTATCGCCTTC‑TGTTTGTTAATCAAATATATCAAAATAAAACGCGCAAATGATAATGATTGTTATGCATGATAAATTCAAGCTATGGGTGGAATTGATATGTTTATGTTCTGAATAGACGGACCGCCTCTTCTCAGGCATCATTACGCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACgt                                            <  1:620501/246‑1 (MQ=255)                                                                                                                                                                                                                  cACGGAAGCGCAGATCTTTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAACCTCGCAAATGATAATGTTTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTG‑TGAGAt                                    <  1:361716/251‑1 (MQ=255)                                                                                                                                                                                                                     ggAAGCGCAGATCATGGCGAGCGCGCTGCGGGTGGCGGGGGGTTTTATTCATTGTTATCGCCTTCGTGAT‑GGTAAACAGATATATCTAAATAAAACTCGCAAATGAAAATGATTGTTAATCTTGATAAATGCACGCGATTTGTAGAACTAATATGTCTATAGTCTGTTAAGACGGACCGCCTCTTCTCAGGCATCATTGCTCAACGCCGGATGCGGCGTGACCGCCTTATCCGTCCTACGTG‑TGAGATGAg                                 <  1:497492/251‑1 (MQ=255)                                                                                                                                                                                                                          cgcAGATCATGGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTTTTATCGCTTTCGTGAT‑GGTAATCAGTTATATCTAAATAAACCTCGCAAATGATAATGATGGTTAATCATGATAAATTCAAGCGATTTGTAGATCTGATATTTCAATAGTGTGATAAGCCGATCCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACGTG‑TGAGATGAGTCCAc                            <  1:294564/251‑1 (MQ=255)                                                                                                                                                                                                                             aGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATACTGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCGTTATCCGGCCTAC‑TGCTGAGATGAGTACACGtt                         >  2:37444/1‑251 (MQ=255)                                                                                                                                                                                                                               aTCATTGCGAACGCGCTGCGGGGGGCGGGGGGTGTTATTCATTTTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGAAAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTG‑TGAGATGAGTCCACGTTTa                       <  2:525119/251‑1 (MQ=255)                                                                                                                                                                                                                                  aTTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCACTGTTATCACCTTCGAGAT‑GGTACTCAGATATATCTAAATAAAACTCGCAAATGATCATGATTGTTACTCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGCACGCCTTATCCGGCCTACGTG‑AGAGATGAGTCCACGTTTaaa                     >  1:614797/1‑250 (MQ=255)                                                                                                                                                                                                                                    ggCGAACGCGCTGCGGGGGGCGGGGGGTTTTATTAATTGTTATCGCCTTCGTGAG‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTG‑TGAGATGAGTCCACGTTTaaaaa                   <  1:517516/249‑1 (MQ=255)                                                                                                                                                                                                                                           gcgcTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTG‑TGAGATGAGTCCACGTTTAAAAATCAttttt           >  2:252466/1‑251 (MQ=255)                                                                                                                                                                                                                                             gcTGCGGGGGGCGGGGGGGGTTATTCATTATTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATACTGATTGTTATTCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTCTAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCTTTACTCAACTCCGGTGGCGGCGTGACCGCCTTTTCCGGCCTACGCG‑TGAGTTGGGGCCACGTTTAAAACTCGttttttc         >  2:4510/1‑249 (MQ=255)                                                                                                                                                                                                                                                      aGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTG‑TGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCa  >  1:500744/1‑249 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                           ATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAT‑GGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTG‑TGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCA  >  NC_000913/3216228‑3216721

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑