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%	128.1 / NA	40	N678N (AAC→AAT)	ynfE	putative selenate reductase YnfE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/18);  total (22/18)

CGTTACCGTCTTACGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTAT  >  NC_000913/1659861‑1660347                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      cGTTACCGTCTTACGACGAACTGAAAAAAATGGGTAACTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAATCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAAGTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGaa                                                                                                                                                                                                                                              <  1:248445/251‑1 (MQ=255)  gTTACCGTCTTACGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAAc                                                                                                                                                                                                                                             >  1:208450/1‑251 (MQ=255)    tACCGTCTTACGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACg                                                                                                                                                                                                                                            <  2:392732/250‑1 (MQ=255)                   aaCTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTg                                                                                                                                                                                                                            <  2:167803/251‑1 (MQ=255)                                                   gATCCCAATTGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTac                                                                                                                                                                                            <  2:123242/251‑1 (MQ=255)                                                              ccATTTAGTCGCATACCACGCATTTGGTGCCGACCCCGAGCCAAATCCACTTAAAACGCCTGCCGTCAAGATTGAAATTTATTCCAGCAGGCTCGCGGAAATAGCCAGTACCGGGGAACTGGAAAAAGATGAAGAGATAAGCCCATTGCAAGTTCATGCCTCAACCTTTGAAGGCTGGAATTCCCATGAACGTAGAACCTTCCCACTGCAACTGTTTAGTTTCCATTACAAATCCCGTACTCACTCGACCt                                                                                                                                                                                 >  2:243268/1‑251 (MQ=255)                                                                aTTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCTATTACAAATCCCGTACTCACTCGACCTAc                                                                                                                                                                               >  2:500501/1‑251 (MQ=255)                                                                     gTCGCCTACAAAGCCTTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGCAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCAATCGACCTACGGCa                                                                                                                                                                           >  2:434334/1‑250 (MQ=255)                                                                       cGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAat                                                                                                                                                                         <  1:385875/250‑1 (MQ=255)                                                                              aaaGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGAtct                                                                                                                                                                 >  2:190719/1‑251 (MQ=255)                                                                                             gacCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGACCGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGCAGGctgcctg                                                                                                                                                  >  1:333969/1‑251 (MQ=255)                                                                                                    aGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCag                                                                                                                                           <  1:274925/251‑1 (MQ=255)                                                                                                     ggCAAATCCACTTAAAACGCCTTCCGGTACGATTGAAATTTCTTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGCAAACGATGAAGAGATAAGCCCATTGCCAGTTTATGACTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCagg                                                                                                                                          >  2:86105/1‑251 (MQ=255)                                                                                                       cAAATCCACTTAAAACGCCTTCCGGAAAGATTGAAATTTATTCCAGAAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCACTATTGATCTCCTGACGGCTGCCTGCCGTCagg                                                                                                                                          >  1:115309/1‑249 (MQ=255)                                                                                                            ccACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGtgtg                                                                                                                                   <  1:157813/251‑1 (MQ=255)                                                                                                                     aCGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAAccc                                                                                                                          >  1:537113/1‑251 (MQ=255)                                                                                                                       gCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACt                                                                                                                                                                                                                             >  1:442185/1‑150 (MQ=255)                                                                                                                       gCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACt                                                                                                                                                                                                                             <  2:442185/150‑1 (MQ=255)                                                                                                                       gCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCGTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAGTCCCGTACTCACTCGACCTACGGCAATATGGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCta                                                                                                                        <  1:260583/251‑1 (MQ=255)                                                                                                                             cGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGCAGGctgc                                                                                                                                                     >  1:277357/1‑216 (MQ=255)                                                                                                                             cGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCGCCTGAAGGctgc                                                                                                                                                     <  2:277357/216‑1 (MQ=255)                                                                                                                             cGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTGTCCATTACAATTCCCTTACGCACTCGACCTGCGGCAATAGTGATCTTCTGAAGGCTGCCTGCCTTGAGGAGGTGGGGATCAGCCCTATAGATg                                                                                                                  <  2:364410/251‑1 (MQ=255)                                                                                                                              ggTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACATGGGAACTGGAAAAAGATGACGTGATAAGCCAATTGCCAGTTTATGCATCAACCTTTGACGGCTGGACTTCCCCTGAACGTAGCACCTTCCAACTGCAAATGTTTGGTTTCCATTACAAATCCAGTACTAAATAGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCAGTCAGGAGGCGTGGATCAACCCTATAGATgc                                                                                                                 >  1:30789/1‑251 (MQ=255)                                                                                                                                        aaaTTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAATGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGAGGCGCATAAACGt                                                                                                       <  2:51782/251‑1 (MQ=255)                                                                                                                                                    gcagGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCCTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGTATTGCCAAc                                                                                           <  2:515947/251‑1 (MQ=255)                                                                                                                                                                       ccGTTCCTGGGAACTTGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCGTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTGCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTTCTCACTCGACCGACGGCAATATTGATCTCCGGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGCTATGGTGCGGgtgt                                                                        <  2:452370/251‑1 (MQ=255)                                                                                                                                                                        cGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGtt                                                                       <  2:208450/251‑1 (MQ=255)                                                                                                                                                                        cGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGtt                                                                       <  1:190719/251‑1 (MQ=255)                                                                                                                                                                              tGGGAACTGGAAAAAGATGAAGTGATAAGCCCCTTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGAAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAACTCCCGTAATCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAAcc                                                                  >  1:463524/1‑250 (MQ=255)                                                                                                                                                                                 gAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCGCCTGAAGGCTGCCTGCCGTCAGGAGGTGTTGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCg                                                              <  1:86105/251‑1 (MQ=255)                                                                                                                                                                                           aaGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGtt                                                    <  1:159612/251‑1 (MQ=255)                                                                                                                                                                                                           cccATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCaa                                      >  1:78509/1‑249 (MQ=255)                                                                                                                                                                                                                   cAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAAACTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGCCCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAAcacc                             >  2:388362/1‑250 (MQ=255)                                                                                                                                                                                                                     gTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGAGAAGTTCGGCTACAAGCCAAAGTAACACCACg                          >  1:338396/1‑251 (MQ=255)                                                                                                                                                                                                                                        aaGGCTGGATTTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAgc       >  2:270133/1‑251 (MQ=255)                                                                                                                                                                                                                                        aaGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAgc       >  2:338000/1‑251 (MQ=255)                                                                                                                                                                                                                                        aaGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCAGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAgc       >  1:351493/1‑251 (MQ=255)                                                                                                                                                                                                                                        aaGGCTGGAATTCCCCTGAAAGTAGAACCTTCCCACGGCAAATGTTTGGTTTCAATTACAAATCCAGTACTCACTCGACCTACGGCAATCTTGTTCTCCTGAAGGCTGCATGCCGTCAGGTGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGAGTGCCAACGGCGATATGGTGCGGGTGCTTACCCATCGTGGCGAAGTTCTGCTACCAGACAAAGTAACACCACGTATTATCACTGGAGTTAgc       >  2:402070/1‑251 (MQ=255)                                                                                                                                                                                                                                            cTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCt    >  1:496524/1‑250 (MQ=255)                                                                                                                                                                                                                                             tGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTAt  >  2:489027/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      CGTTACCGTCTTACGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTAT  >  NC_000913/1659861‑1660347

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: ‑