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

CCGGCGTTACCGTCTTACGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGC  >  NC_000913/1659857‑1660344                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   ccGGCGTTACCGTCTGACGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACCAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAAATGGAAAAAGATGAAGTGATAAGCCCATTGCAAGTTTATGCCTCAACCTTTGAAGGCTGGAATTcccc                                                                                                                                                                                                                                               >  2:514822/1‑251 (MQ=255)   ggCGTTACCGTCTTACGACGAACTGAAAAAACTGGGTATCTATAAGCGTAACGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAACTTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAACAAGCTGAAGTGCTAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCtc                                                                                                                                                                                                                                             >  2:78458/1‑250 (MQ=255)                   gacgaACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTccc                                                                                                                                                                                                                             >  2:131540/1‑251 (MQ=255)                            aaaaaaaTGGGTATCTATAAGCGTAAAGATCCCACTGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACt                                                                                                                                                                                                                    >  2:71889/1‑251 (MQ=255)                                    gggTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGtt                                                                                                                                                                                                            <  1:178305/251‑1 (MQ=255)                                     ggTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGttt                                                                                                                                                                                                           >  1:472017/1‑251 (MQ=255)                                            ataAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGCACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTAc                                                                                                                                                                                                    >  2:550233/1‑251 (MQ=255)                                                 cGTAAAGATCCCAATGGCCATTTTGTCGCCTAAAAAGCATTTCGTGACGACCCCGAGGCAAATACACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAAACTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTCGGTTTCCATTACaactc                                                                                                                                                                                               >  2:451296/1‑248 (MQ=255)                                                   tAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATccc                                                                                                                                                                                             >  1:7515/1‑251 (MQ=255)                                                                 gCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAATGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGAcc                                                                                                                                                                               >  1:376963/1‑251 (MQ=255)                                                                    aTTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTAc                                                                                                                                                                            >  1:367057/1‑251 (MQ=255)                                                                      tttGTCGCCTACAAAGCATGTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGGTTGGTTGCCATTACAAATCCCGTACTCACTCGACCTACg                                                                                                                                                                           <  2:160873/250‑1 (MQ=255)                                                                       ttGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAAGCCCGTACTCACTCGACCTACgg                                                                                                                                                                          <  1:408195/250‑1 (MQ=255)                                                                                         ttCGTGACGCCCCCGAGGCAAATCCACTTCAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCACCCTTTGAAGGCTGGAATTCCCCTGACCGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCCATATTGATCTCCTGa                                                                                                                                                         >  1:503934/1‑249 (MQ=255)                                                                                            gTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTAGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGct                                                                                                                                                    >  1:39919/1‑251 (MQ=255)                                                                                                 gacCCAGAGGCAAATACACTTAACACGCCTTACGGTAAGCTTGCAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCCGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCACTATTGATCTCCTGAAGGctgcctg                                                                                                                                               >  2:407045/1‑251 (MQ=255)                                                                                                        aGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTc                                                                                                                                          >  1:192187/1‑249 (MQ=255)                                                                                                                                      aGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAg                                                                                                          <  1:134193/251‑1 (MQ=255)                                                                                                                                       gATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAATGCTGCCTGCCGTCATGAGGTGTGGATCAACCCTATAGATGCGCAGa                                                                                                         <  1:550233/251‑1 (MQ=255)                                                                                                                                                ttATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTggg                                                                                                 <  2:176799/250‑1 (MQ=255)                                                                                                                                                  aTTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGAt                                                                                               <  2:472017/250‑1 (MQ=255)                                                                                                                                                          agGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCCTTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGCGGTGTGGATCAACCCTATAGCTGCGCCGAAACGTGGGATTGCCAACgg                                                                                      >  2:554402/1‑251 (MQ=255)                                                                                                                                                                     aaTTGCCCGTACCTGGGAACTGGAAAAAGATGCAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCACCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGc                                                                           >  1:510200/1‑251 (MQ=255)                                                                                                                                                                      aTTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCg                                                                          <  1:252455/251‑1 (MQ=255)                                                                                                                                                                               aCCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTaa                                                                 <  1:78458/251‑1 (MQ=255)                                                                                                                                                                                ccTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAAc                                                                <  1:71889/251‑1 (MQ=255)                                                                                                                                                                                   gggAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCGTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAAcc                                                               <  2:488972/249‑1 (MQ=255)                                                                                                                                                                                     gAACTGGTAAAAGATGAAGAGCTCAGACCATTGCCCGTTTATGCCACCACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAAAATTCCCACTGCCACAGTTTGGTTTCCATTAAAAATCCCGTACTCACTCGACATACGGCAATATAGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGCGTGGATAAACCCTATAGATGCGCAGAAACGCGGGATTGCCCACGGCGATATGGTGCGGGTGTTTAACCCTCg                                                           >  2:565484/1‑251 (MQ=255)                                                                                                                                                                                             aaaaGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACATTTGAAGGCTGGAATTCCACTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGCGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGagg                                                   >  2:560150/1‑249 (MQ=255)                                                                                                                                                                                              aaaGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCACCTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGt                                                  >  2:292018/1‑251 (MQ=255)                                                                                                                                                                                                     aaGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTa                                           <  1:131540/251‑1 (MQ=255)                                                                                                                                                                                                         gATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCa                                        <  1:560150/250‑1 (MQ=255)                                                                                                                                                                                                                  aTTGCCAGTTTATGCCTCGACCTTTGAAGGCTGGCATTCCCCTGAACGTAGAACCTTCCCACTGCCACTGTTTGGCTTCCATTACAAATCCCGTACTCACTCGACCTACGGCACTATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGCTGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGGGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTa                               >  2:408349/1‑250 (MQ=255)                                                                                                                                                                                                                            tATGCCTCAACCTTTGAAGGCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCTTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTTCCAACGGCGATATGGTGCGGGTGTTTCACCATCGTGGCGAAGTTCGGCTACCAGCCAAACTAACACCACGTa                     <  1:514822/250‑1 (MQ=255)                                                                                                                                                                                                                             aTGCCTCAACCTATGAAGGCTGGAATTCCCCTGAACGTAGAGCCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAACGTCGGCTGCCATCCACAGTGGCATCACGTAt                    <  1:179608/250‑1 (MQ=255)                                                                                                                                                                                                                                    aaCCGTTGAAGGCTGGTATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGAGCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTg            <  1:6893/251‑1 (MQ=255)                                                                                                                                                                                                                                              ggCTGGAATTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAgcgc  >  2:509188/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   CCGGCGTTACCGTCTTACGACGAACTGAAAAAAATGGGTATCTATAAGCGTAAAGATCCCAATGGCCATTTTGTCGCCTACAAAGCATTTCGTGACGACCCCGAGGCAAATCCACTTAAAACGCCTTCCGGTAAGATTGAAATTTATTCCAGCAGGCTGGCGGAAATTGCCCGTACCTGGGAACTGGAAAAAGATGAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGC  >  NC_000913/1659857‑1660344

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