Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,524,596	C→T	D139N (GAT→AAT)	yjhH ←	putative 2‑dehydro‑3‑deoxy‑D‑pentonate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,524,596	0	C	T	100.0%	171.6 / NA	57	D139N (GAT→AAT)	yjhH	putative 2‑dehydro‑3‑deoxy‑D‑pentonate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/33);  total (24/33)

AGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCCCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTG  >  NC_000913/4524355‑4524842                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       agTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCCAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGGTGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCCGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGtt                                                                                                                                                                                                                                               >  1:253267/1‑251 (MQ=255)                  aGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCAGCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCTCGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGt                                                                                                                                                                                                                               <  2:136422/250‑1 (MQ=255)                       gggCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTGAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTg                                                                                                                                                                                                                         <  2:41392/250‑1 (MQ=255)                              tATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCTCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTTTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCATCGATATTCTCGTTTTGCATAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTa                                                                                                                                                                                                                   <  1:402082/250‑1 (MQ=255)                                 cGCACCGTCGCCGCTCAGCAGCATCGTAGTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTGTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACggg                                                                                                                                                                                                               <  2:376529/251‑1 (MQ=255)                                  gCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACTCTACGGGc                                                                                                                                                                                                              <  2:380059/251‑1 (MQ=255)                                       gTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATct                                                                                                                                                                                                         >  1:348776/1‑251 (MQ=255)                                         cgccgcTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATctg                                                                                                                                                                                                        <  1:214175/250‑1 (MQ=255)                                          gccgctCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGCATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTATTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCCCGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATctgct                                                                                                                                                                                                      >  1:55404/1‑251 (MQ=255)                                               tcagcagCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCTGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTa                                                                                                                                                                                                  >  1:212494/2‑250 (MQ=255)                                                cagcagCATCTTATCCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACTGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTaa                                                                                                                                                                                                 <  1:78325/250‑1 (MQ=255)                                                 agcagcATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCTGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAAt                                                                                                                                                                                                >  2:286510/1‑250 (MQ=255)                                                   cagcaTCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTATGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAgt                                                                                                                                                                                             <  1:90738/251‑1 (MQ=255)                                                    agcaTCGTATTCAGCAAATGATACTCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTCCGCAAGTGACCAACGATGTCGATGGTGTCTTTGATGCCAACGATATTCCCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCCGATTCGGAAAGTTGTACAGGATCACAGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAgtc                                                                                                                                                                                            >  2:255592/1‑251 (MQ=255)                                                      caTCGTATTCAGCAAATGATCATCGTAACCGCAGAATCACGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCACGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCCGACGCGTCACGGTTTCCGGGGTTCAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTCATAgtcgg                                                                                                                                                                                          >  2:384694/1‑250 (MQ=255)                                                      caTCGTATTAATCTAATGATCATCGTCACCGCAGAATAACGAAAACGACGGGCGTGCCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCGTTGATGCCAACGATAGTCTCGTTTTGCAGAGCCGGACGCTTAACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGGACAGGATCACCGGTAGGTTGTCGCTACGGGCTATCTGCTGGTAATAgtcgt                                                                                                                                                                                          <  1:134861/251‑1 (MQ=255)                                                       aTCGTATTCAGCAAATGATCATCGTAACCGCAGAATCCCGAAAACGCCGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGACAGTTGAACAGGGTCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAgtcgtc                                                                                                                                                                                         >  2:95391/1‑251 (MQ=255)                                                        tCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGCGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTCAGTCCTGACCCGTCAGATTCGGAAAGTTGCACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAgtcgtc                                                                                                                                                                                         >  1:64009/1‑250 (MQ=255)                                                         cGTAGTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCa                                                                                                                                                                                        <  1:89537/250‑1 (MQ=255)                                                               cAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTc                                                                                                                                                                                 <  1:258356/251‑1 (MQ=255)                                                                    aaTGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCgtggt                                                                                                                                                                            <  1:97981/251‑1 (MQ=255)                                                                    aaTGATCAACGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCgtggt                                                                                                                                                                            <  2:211441/251‑1 (MQ=255)                                                                       gATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGc                                                                                                                                                                          <  1:127501/250‑1 (MQ=255)                                                                          catcGTAACCGCAGAATACCGAAAACGACGGGCTTACCGACTTATCTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGTGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGAc                                                                                                                                                                       <  2:225504/250‑1 (MQ=255)                                                                           atcGTAACCGCAGAATACCGAAAACGACGGTCGTACCGTCTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCAAGACGCGTCACGGTTTCCGGGGTTACGTCCTGACCCGTCCGATTCGGAAAGTTGTACCGAATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTCATAGTCGTCAAGATTTCGTGGTGCGACtt                                                                                                                                                                     >  2:358899/1‑251 (MQ=255)                                                                           atcGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGGGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGGCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCGGGTAATAGTCGTCAAGATTTCGTGGTGCGACtt                                                                                                                                                                     <  2:53475/251‑1 (MQ=255)                                                                             cGTATCCGCAGAATCCCGAAAACGACGGTCGTCTCGACTTAACTGTGTTGATCATCGTTCTCAAGTGTCCAACGCTGTAGATTGTGTCTTTGATGCCATCGATATTCTCGTTTTGCAGATCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGACCACCGGTAGGGTGACTCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGTCTTTc                                                                                                                                                                   <  2:17632/251‑1 (MQ=255)                                                                                  ccGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGa                                                                                                                                                                                                            >  1:105074/1‑205 (MQ=255)                                                                                  ccGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGa                                                                                                                                                                                                            <  2:105074/205‑1 (MQ=255)                                                                                   cGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCagtag                                                                                                                                                             <  1:267628/251‑1 (MQ=255)                                                                                      aGAATACCGAAAACGACGGGCGTCACGACTTAACTGTGTTTATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGtagtag                                                                                                                                                          <  2:336609/251‑1 (MQ=255)                                                                                         aTACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTg                                                                                                                                                                                                                         <  1:266095/185‑1 (MQ=255)                                                                                         aTACCGAAAACGACGGGCGTACAGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAACGTTGTACAGGATCACCGGTAGGGTg                                                                                                                                                                                                                         >  2:266095/1‑185 (MQ=255)                                                                                            ccGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGGGCGACTTTCCAGTAGTAgggggtg                                                                                                                                                    >  2:30857/1‑248 (MQ=255)                                                                                              gAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGCTCACCGGTAGGGAGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGAt                                                                                                                                                  >  2:357568/1‑251 (MQ=255)                                                                                                aaaCGACGGGCGTACCGCCTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGCTTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATgg                                                                                                                                                >  2:260795/1‑251 (MQ=255)                                                                                                   cgacgGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATACTCGTTTTGCAGAGCCAGACGCGGCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACATGATCACCGGTAGGGTGACGCTACGGGCTATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGGAGGGGTTGATGGcta                                                                                                                                             <  2:18879/251‑3 (MQ=255)                                                                                                     acgGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGCAGGGGTTGATGGcgacg                                                                                                                                           <  2:183729/251‑1 (MQ=255)                                                                                                               cGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATc                                                                                                                                   >  1:152175/1‑249 (MQ=255)                                                                                                                    tAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGAGTCGGAAAGTTGTACAGGATCACCTGTAGGGTGACGCTACGGGCTATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGTCGACGATACCATCAgcgc                                                                                                                              <  2:64009/249‑1 (MQ=255)                                                                                                                             tGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTg                                                                                                                    <  1:358899/250‑1 (MQ=255)                                                                                                                                 catcGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCAt                                                                                                               >  1:362569/1‑251 (MQ=255)                                                                                                                                    cGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCAtgct                                                                                                            <  2:39299/251‑1 (MQ=255)                                                                                                                                    cGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTGGATGCCAACGATAGTCTCGTTGTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCGACGGGCGATCTGCTGGGAATAGTCGTCAATATTTCGTGGTGCGACCTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCTCAtgct                                                                                                            <  2:149186/251‑1 (MQ=255)                                                                                                                                           aaGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGcc                                                                                                       >  2:313756/1‑249 (MQ=255)                                                                                                                                             gTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGACAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGcc                                                                                                       >  1:177002/1‑247 (MQ=255)                                                                                                                                                    aCGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGc                                                                                            >  1:244004/1‑251 (MQ=255)                                                                                                                                                                                          gTTTTGCAGAGCAAGACGCGGCACGGTTTCCTGGGTTAAGTCCGGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCg                                                      <  2:348776/251‑1 (MQ=255)                                                                                                                                                                                                       aGACGCGTTACGTTTTCCTGGGATAAGTGCTGACCCGGCAGATTCGAAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGGAATAGTCGTCCAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGAAGGATTCCGCATGCTGCGCCAGATTAACCGCTTCGTCAGGGGAAGGGGAACAGACGCCAATCAATACCGGCACTCGCCCggc                                         >  2:154043/1‑249 (MQ=255)                                                                                                                                                                                                          cgcgTAACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGAc                                      >  1:84381/1‑251 (MQ=255)                                                                                                                                                                                                               cACGGTTTCCGGGGTTAAGTCCTGACACGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACattg                                  <  1:384694/250‑4 (MQ=255)                                                                                                                                                                                                                              tAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGCTTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCgg                  >  1:253870/1‑251 (MQ=255)                                                                                                                                                                                                                                 gtcctgacaaggcAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTCGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGGCGACAATGGTTACAGCTTCTTCGGCgg               >  2:156614/13‑250 (MQ=255)                                                                                                                                                                                                                                      gACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTg            <  1:271943/249‑1 (MQ=255)                                                                                                                                                                                                                                        cccGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCAt        <  1:25805/251‑1 (MQ=255)                                                                                                                                                                                                                                          cGTCAGATTCGGAAATTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGGAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATgc      <  2:212494/251‑1 (MQ=255)                                                                                                                                                                                                                                              aGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTGCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAATGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCGTCGGCGAGTTCCATGCGCTg  <  2:253267/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       AGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCCCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTG  >  NC_000913/4524355‑4524842

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