Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,320,785	C→T	E155K (GAG→AAG)	phnH ←	carbon‑phosphorus lyase core complex subunit PhnH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,320,785	0	C	T	100.0%	238.9 / NA	69	E155K (GAG→AAG)	phnH	carbon‑phosphorus lyase core complex subunit PhnH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/39);  total (30/39)

CTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGG  >  NC_000913/4320544‑4321026                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  ctgctgTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTg                                                                                                                                                                                                                                          >  2:161062/1‑251 (MQ=255)   gctgTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGc                                                                                                                                                                                                                                         <  2:174267/250‑1 (MQ=255)   gctgTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGATCGTGCAGAATGCACTTCGGCAGCTGCg                                                                                                                                                                                                                                        >  1:261873/1‑251 (MQ=255)               gggCGACGCTCAGTTCGGGCAAATCTGCTTCGCCGCTGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCTCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTAAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCTGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTc                                                                                                                                                                                                                            <  2:11930/251‑1 (MQ=255)               gggCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCAtt                                                                                                                                                                                                                             <  2:80015/250‑1 (MQ=255)               gggCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTc                                                                                                                                                                                                                            >  1:130141/1‑251 (MQ=255)                 gCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTCCATCAGCACACCTCCACCTGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCCCAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGt                                                                                                                                                                                                                          >  2:250319/1‑251 (MQ=255)                   gACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTGAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGGGATCATTCGttc                                                                                                                                                                                                                        <  1:192986/251‑1 (MQ=255)                   gACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGttc                                                                                                                                                                                                                        <  1:193004/251‑1 (MQ=255)                     cGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCGTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGttctt                                                                                                                                                                                                                      <  1:15464/251‑1 (MQ=255)                                  cAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGc                                                                                                                                                                                                         >  1:185495/1‑251 (MQ=255)                                   aaaTCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGc                                                                                                                                                                                                         >  1:54743/1‑250 (MQ=255)                                    aaTCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCTTTCTTCGGCAAt                                                                                                                                                                                                               <  2:34868/243‑1 (MQ=255)                                    aaTCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAAt                                                                                                                                                                                                               >  1:34868/1‑243 (MQ=255)                                        ggTATCGCATCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCATTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGGTGGT‑CGCGGAATAGGCAGCCGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTGGTGCAGAATGCACTTAGGCAGATGCGGAGCGATACTTCGTTCTTCGGCAATACCCGCAcccgt                                                                                                                                                                                                   >  2:78239/1‑248 (MQ=255)                                         gTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTa                                                                                                                                                                                                  <  2:249393/251‑1 (MQ=255)                                             cGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAgg                                                                                                                                                                                               <  2:44060/250‑1 (MQ=255)                                             cGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAgg                                                                                                                                                                                               <  1:161062/250‑1 (MQ=255)                                                     cGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATg                                                                                                                                                                                      <  2:78148/251‑1 (MQ=255)                                                          gCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGTCGCTCGGTGAGTTCGTTCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTcc                                                                                                                                                                                 <  1:138431/251‑1 (MQ=255)                                                          gCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTcc                                                                                                                                                                                 <  1:71632/251‑1 (MQ=255)                                                                   cAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAg                                                                                                                                                                        <  1:238193/251‑1 (MQ=255)                                                                    aGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAgg                                                                                                                                                                       <  1:232392/251‑1 (MQ=255)                                                                       gCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTg                                                                                                                                                                    <  1:307823/251‑1 (MQ=255)                                                                         gTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGAAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGc                                                                                                                                                                  >  2:149562/1‑251 (MQ=255)                                                                               ggcgTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTg                                                                                                                                                            <  1:292667/251‑1 (MQ=255)                                                                                 cgTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTa                                                                                                                                                          <  1:178851/251‑1 (MQ=255)                                                                                  gTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTaa                                                                                                                                                         <  2:122333/251‑1 (MQ=255)                                                                                          ccTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCg                                                                                                                                                 >  2:189130/1‑251 (MQ=255)                                                                                                         tCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACATGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATCCCTGCACCAGTAAGGCGCAACATGCTTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTcc                                                                                                                                   <  1:15741/250‑1 (MQ=255)                                                                                                          cACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCgg                                                                                                                                 >  1:17351/1‑251 (MQ=255)                                                                                                           aCGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGgc                                                                                                                                >  1:159716/1‑251 (MQ=255)                                                                                                            cGGCAACGTACATCAGCACACCTCCACATGAG‑TCGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGACGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGgcg                                                                                                                               >  1:256744/1‑251 (MQ=255)                                                                                                                aaCGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAAc                                                                                                                           >  1:98548/1‑251 (MQ=255)                                                                                                                aaCGTACATCAGCACACCTCCAAATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACAAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAAc                                                                                                                           >  2:268166/1‑251 (MQ=255)                                                                                                                      cATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGt                                                                                                                     >  1:52/1‑251 (MQ=255)                                                                                                                       aTCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTg                                                                                                                    <  2:219108/251‑1 (MQ=255)                                                                                                                       aTCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCGCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTg                                                                                                                    <  1:201488/251‑1 (MQ=255)                                                                                                                          aGCACACCTCCACATGAG‑TGGCTCGCGGAATAGCCAGCAGGCGCTCGACACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGAGCTCGGTGAGTTCGTGCAGAATGAACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCACTACCCCCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGCCCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGAGCCg                                                                                                                 >  1:102036/1‑251 (MQ=255)                                                                                                                                 cTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATACACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGaa                                                                                                           >  1:3730/1‑250 (MQ=255)                                                                                                                                                 gcgcgGAATAGCCAGCAGGCGCGCGCCAGAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGtt                                                                                           <  2:33041/250‑1 (MQ=255)                                                                                                                                                  cgcgGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTc                                                                                          <  1:120445/251‑1 (MQ=255)                                                                                                                                                  cgcgGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGGGGAAAGGGCGTTgggcgctgc                                                                                          >  2:223258/1‑243 (MQ=255)                                                                                                                                                   gcgGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCg                                                                                         <  1:99827/251‑1 (MQ=255)                                                                                                                                                           gCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAAt                                                                                 <  1:14696/251‑1 (MQ=255)                                                                                                                                                            ccggcagGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGATTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTTAGCTGTTCGCTGGAAATc                                                                                <  2:214413/248‑1 (MQ=255)                                                                                                                                                            ccAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGGCGGGGGAAAGGGCGGTTAGCTGTTCGCTGGAAATc                                                                                >  1:85039/1‑251 (MQ=255)                                                                                                                                                             cagcggGCGCTCACCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCTGTGAGTTCTTGCAGAATGCCCTTTGGCAGCTGAGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCATTAAGGCTCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAATATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAGTCg                                                                               <  1:206515/251‑1 (MQ=255)                                                                                                                                                             cagcagGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCCGTGGAAAggggggtgggcgtgtcgcttgacctcg                                                                               >  2:224571/1‑227 (MQ=255)                                                                                                                                                                    cgcTCGCCACAGGTCAGGATCAGGTCGCTGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGCAAGGGCGTTGGGCTGTTCGCTGGAAATCGCCTCATc                                                                        >  1:236678/1‑251 (MQ=255)                                                                                                                                                                      cTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGGGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCCCATCCg                                                                      >  1:235151/1‑251 (MQ=255)                                                                                                                                                                         gCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGGGCCGGGGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTcc                                                                   >  2:134747/1‑250 (MQ=255)                                                                                                                                                                           cacaGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCAcc                                                                 <  1:19443/251‑1 (MQ=255)                                                                                                                                                                             caGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACcgc                                                               <  1:166736/251‑1 (MQ=255)                                                                                                                                                                              aGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACcgcg                                                              <  2:263849/251‑1 (MQ=255)                                                                                                                                                                                       tCAGGTGGATCCCCAGCGGGTACGGATGCGGGCGCTCGGTGAGTTCGGCCAGCATGCCCTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGCAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCTACGCCAGCGCGCGCGGAAAGGGGGTTGTGCTTTTCTCTGGAAATCCCCCCATCCGCCCCCGCGAAggggcgc                                                     >  2:284565/1‑246 (MQ=255)                                                                                                                                                                                              gATGCCGGGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCg                                              <  2:17378/251‑1 (MQ=255)                                                                                                                                                                                                   cGAGCGGAAACGGAGTCGGGGGCGCGGTGAGTGTGTGCAGAATGCACTTCGGCAGCTGCGGAGCTATCATGCGTGTTTCGGCAATACCCGCTCCAGTAAGGCGCAACAGGCTGCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCATCCGCTGCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGCGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCggctgt                                         <  2:160309/251‑2 (MQ=255)                                                                                                                                                                                                              ggATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAgcgg                              <  1:208004/251‑1 (MQ=255)                                                                                                                                                                                                                   cGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATCCCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGt                         <  2:58591/251‑1 (MQ=255)                                                                                                                                                                                                                    gggCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGGAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCCTCCGTCCCCCCGACGGGCGCCTGTTCCGGCCGGCTCAcccgctgcccgtt                        >  2:145719/1‑240 (MQ=255)                                                                                                                                                                                                                           cGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCCTCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATg                  >  1:147564/1‑250 (MQ=255)                                                                                                                                                                                                                             gTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGaaa               <  1:95907/251‑1 (MQ=255)                                                                                                                                                                                                                             gTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGaa                >  1:71026/1‑250 (MQ=255)                                                                                                                                                                                                                                   tCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCa          <  2:304982/250‑1 (MQ=255)                                                                                                                                                                                                                                     gTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGGTGGTATGAAAACGCCGGc       >  1:182071/1‑251 (MQ=255)                                                                                                                                                                                                                                        cAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGctct    <  2:235151/251‑1 (MQ=255)                                                                                                                                                                                                                                        cAGAATGCACTTCGGCAGCTGCGGAGCGATCAGTCGTTCTTCGGCAATCCCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGctct    <  2:145014/251‑1 (MQ=255)                                                                                                                                                                                                                                          gAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGGGGAAAGGGCGTTGAGCTGTTCGCTGGGAATCGCCTCATCCGTCACCGCGAAGGGCGCCTGTTCCGGCTGGCTGACCAGCGGCACGGTGGGATGAAAACGCAGGCTCTgg  >  1:57132/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  CTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAG‑TGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGG  >  NC_000913/4320544‑4321026

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