Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,660,102	C→T	N678N (AAC→AAT)	ynfE →	putative selenate reductase YnfE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,660,102	0	C	T	100.0%	146.7 / NA	45	N678N (AAC→AAT)	ynfE	putative selenate reductase YnfE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/21);  total (24/21)

ACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCT  >  NC_000913/1659876‑1660345                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    acgaACTTAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTCCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAAGTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCa                                                                                                                                                                                                                             <  2:406320/251‑1 (MQ=255)        gAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAAc                                                                                                                                                                                                                      >  1:256361/1‑251 (MQ=255)                 gggTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGtt                                                                                                                                                                                                             >  1:6305/1‑251 (MQ=255)                       cTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCAtt                                                                                                                                                                                                       >  2:186280/1‑251 (MQ=255)                           taGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACa                                                                                                                                                                                                    <  2:394371/249‑1 (MQ=255)                                        ccAATGGTCATTTTGGCGCATACAAATCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATAGAAATCTATTACAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCAGTACAAACCCCTTtctcac                                                                                                                                                                                      <  2:194463/251‑1 (MQ=255)                                                           cTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTg                                                                                                                                                                   >  1:361966/1‑251 (MQ=255)                                                             aCAAAGCATTTCGTGACGACCCCGATGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGAt                                                                                                                                                                 <  2:55099/251‑1 (MQ=255)                                                                aaGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACCAATCCCGTACTCACTCGACCTACGGCAATATTGAtctc                                                                                                                                                              >  2:466887/1‑251 (MQ=255)                                                                 aGCATTTCGTGACTGCCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCAGTACCTGGGAACTGGAAAAAGATTAAGGGATAAGCCCATTGCCAGTTTATGCCTCAAACATTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCCACTGTTTGGTTTCCATTACAAATCCCGTACTCAATCGCCCTACGGCAATATTGATCTcc                                                                                                                                                             >  2:290423/1‑251 (MQ=255)                                                                  gcaattCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAAGTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCt                                                                                                                                                            <  2:299272/247‑1 (MQ=255)                                                                         gTGACGACCCCGAGGCAAATCCACTTAAACCGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATAAAGTGATAAGCCCATTGCCAGTTTCTGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGGTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGct                                                                                                                                                     >  2:360072/1‑251 (MQ=255)                                                                                      ggCAAATCCACTTAAAACGCCTTCCGGTAAGATTGCAATTTATTCCAGAAGGCTCAAGGAAATTGCCCGTACCTGGGAACTGGAAACAGATGAAGTGCTAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGACTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGATTGGTTTCCATTACAAATCCCGTACTCACTCGACTTACGCCTATCTTGATCTACTGTAGGCTGCCAGCCGTcggg                                                                                                                                        >  2:132500/1‑248 (MQ=255)                                                                                          aaTCCACTTAAAACGCCTTCCGGAAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGTTAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCACCTGCACGTAGAACCTTCCCACTGAAACTGTTTGGTTTCCATTACAAATACCGTACTCACTCGACCTACGGCAAACTTGATCACCTGAAGGCTGCCTGCCGTCaggagg                                                                                                                                     >  2:320398/1‑250 (MQ=255)                                                                                           aaccACTTAAAACGCCTGCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCGCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAat                                                                                                                                                                       <  2:136503/213‑1 (MQ=255)                                                                                           aTCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAat                                                                                                                                                                       >  1:136503/1‑215 (MQ=255)                                                                                                      aCGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAAccc                                                                                                                        <  2:344995/251‑1 (MQ=255)                                                                                                                  aaGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTCCCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGGTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCa                                                                                                            <  2:389600/251‑1 (MQ=255)                                                                                                                              tATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATAAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGa                                                                                                 <  1:360072/250‑1 (MQ=255)                                                                                                                               aTTCCAGCAGGCTTGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTGCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCGATAGATGCGCAGAAACGttggctt                                                                                               <  2:331293/251‑7 (MQ=255)                                                                                                                               aTTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGAtt                                                                                               >  1:151734/1‑251 (MQ=255)                                                                                                                               aTTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGAt                                                                                                <  1:119767/250‑1 (MQ=255)                                                                                                                                   cagcagGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAt                                                                                                                              >  2:239087/1‑216 (MQ=255)                                                                                                                                   cagcagGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACGCACTCGTCCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGAt                                                                                                                              <  1:239087/216‑1 (MQ=255)                                                                                                                                     gcagGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAAc                                                                                         <  2:6305/251‑1 (MQ=255)                                                                                                                                      cagGCTGGCGGAAATTGCCCGTTCCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAAGCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACg                                                                                        <  2:256361/251‑1 (MQ=255)                                                                                                                                       agGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCGGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACgg                                                                                       <  1:293049/251‑1 (MQ=255)                                                                                                                                       agGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCGCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACgg                                                                                       <  2:216544/251‑1 (MQ=255)                                                                                                                                              cGGAAATTGCCCGTTCCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTAGGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATg                                                                                <  2:151734/251‑1 (MQ=255)                                                                                                                                               ggAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCACCCCTATAGATGCGCAGAAACGTGGGATTGCCAACGACGATATgg                                                                               >  1:351291/1‑251 (MQ=255)                                                                                                                                               agAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATACGCCAATTGCCAGTTTATGACTCAACCTTTGAAGGCTGGAATTCCCCTGCACGTAGAACCTTCCCACTGCCACTGTTTGGTTTCCATTACAACTCCCGTACTCACTCGACTTACGGCCATATTGATCTCATGACGGCTGCCTGCCGTCCGGCGGTGTGGCTCAACCCTCTAGATGCGCAGAAACGTGGGATTCCCAACGGCGATATgg                                                                               >  2:177679/2‑251 (MQ=255)                                                                                                                                                   aTTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCg                                                                           <  1:262711/251‑1 (MQ=255)                                                                                                                                                      gCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGgt                                                                        <  1:285055/251‑1 (MQ=255)                                                                                                                                                                            aaGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGtt                                                  >  1:403944/1‑251 (MQ=255)                                                                                                                                                                             aGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTc                                                 <  1:104198/251‑1 (MQ=255)                                                                                                                                                                                 gAAGTGATAAGCCCATTGCCAGTTTATGCATCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGCGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCggat                                             >  1:79662/1‑249 (MQ=255)                                                                                                                                                                                               aTTGCCAGTTTATGCCTCATCCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGCTCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATCGATGCGCAGAAACGTGGGGTTGCCAACGGCGATATGGTGCGGGTGTTTACACATCGTGGCGAAGTTCGGCTCCCCGCCAAAGTa                                >  2:30973/1‑250 (MQ=255)                                                                                                                                                                                                ttGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAg                                  >  1:94579/1‑247 (MQ=255)                                                                                                                                                                                                ttGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAg                                  <  2:94579/247‑1 (MQ=255)                                                                                                                                                                                                        ttATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGACCCTTCCCACTTCAACTGTTTGGCTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGCTTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGCAGTTCGGCTACCACCCAAAGTAACACCACGtc                      >  1:169730/1‑250 (MQ=255)                                                                                                                                                                                                          aTGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTAtt                    >  2:203063/1‑251 (MQ=255)                                                                                                                                                                                                             ccTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATtctc                 >  1:168625/1‑251 (MQ=255)                                                                                                                                                                                                                 aaCCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATGCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTg             >  1:339514/1‑251 (MQ=255)                                                                                                                                                                                                                            gCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCt  >  2:168829/1‑251 (MQ=255)                                                                                                                                                                                                                             cTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACAACGTATTCTCCCTGGAGTTAGCGCt  >  1:58232/1‑250 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    ACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCT  >  NC_000913/1659876‑1660345

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