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%	176.8 / NA	54	E155K (GAG→AAG)	phnH	carbon‑phosphorus lyase core complex subunit PhnH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/28);  total (26/28)

GCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGCCACACCGGCGTGTCGT  >  NC_000913/4320500‑4321075                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                   gCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAg                                                                                                                                                                                                                                                                                      <  1:281476/300‑1 (MQ=255) gCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAg                                                                                                                                                                                                                                                                                      <  2:399472/300‑1 (MQ=255) gCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCACTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAg                                                                                                                                                                                                                                                                                      <  2:399410/300‑1 (MQ=255)                             cTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGNATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGc                                                                                                                                                                                                                                                          >  1:189529/1‑300 (MQ=255)                                     aGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGACGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTaa                                                                                                                                                                                                                                                  <  1:478447/300‑1 (MQ=255)                                         tAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGgcg                                                                                                                                                                                                                                              <  2:720989/300‑1 (MQ=255)                                          aaGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGTGCAAATCGGTATTGCCTCGTCGTCGGCTCTCTTGCAGGGCGTGGGCTGCGTCGATCGCCTCCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGgcgc                                                                                                                                                                                                                                             <  2:144246/300‑1 (MQ=255)                                          aaGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTCCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATTAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGgcgc                                                                                                                                                                                                                                             <  2:142480/300‑1 (MQ=255)                                                  gTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGc                                                                                                                                                                                                                                      >  2:489456/1‑299 (MQ=255)                                                            ggCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTc                                                                                                                                                                                                                           <  1:827071/300‑1 (MQ=255)                                                                 cGCTCAGTTCGGGCAAATCGGTATCGCCTAGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCt                                                                                                                                                                                                                      >  1:407109/1‑300 (MQ=255)                                                                       gTTCGGGCAAATCGGTATCGCCTAGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGNCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCAACAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGa                                                                                                                                                                                                                 >  2:603814/1‑299 (MQ=255)                                                                               aaaTCGGTATCGCCTCGTCTGCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCNCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTaaaa                                                                                                                                                                                                        <  2:251448/300‑1 (MQ=255)                                                                                                             tGCAGGGCGTGGGCGGCGTCGATCGTCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAAcgg                                                                                                                                                                          <  2:515159/300‑1 (MQ=255)                                                                                                                  ggCGTGGGCGGCGTCGATTGCCTTCTCGCCCCCTTTCACTGCAACGTACATCAGCATACCTCCACATGAGGGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGGCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCAGGCGCTCGGTGAGTTCGTGCAGAATGCACGTCGGCAGCTGCGGAGCGATCAGTCGTTCTGCGGCAATAccc                                                                                                                                                                                                                                                            <  2:149487/213‑1 (MQ=255)                                                                                                                  ggCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACa                                                                                                                                                                                                                                         <  1:154903/232‑1 (MQ=255)                                                                                                                  ggCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACa                                                                                                                                                                                                                                         >  2:154903/1‑232 (MQ=255)                                                                                                                  ggCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACACCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAGATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGt                                                                                                                                                                      >  2:240193/1‑299 (MQ=255)                                                                                                                            gcgTCGATCGCCTTCTCGCCCCCTTTAACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGACTGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAACGTCAGCGTCGCACCCGCTTCCGgcgc                                                                                                                                                                               >  2:456098/1‑280 (MQ=255)                                                                                                                                  aTCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGt                                                                                                                                                     <  1:765355/300‑1 (MQ=255)                                                                                                                                      ccTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCNGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGa                                                                                                                                                  >  1:598005/1‑299 (MQ=255)                                                                                                                                        ttCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGc                                                                                                                                                <  1:318928/299‑1 (MQ=255)                                                                                                                                        ttCTCGCCCCCCTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGTCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGc                                                                                                                                                <  1:315252/299‑1 (MQ=255)                                                                                                                                              cccccTTTCACGGCAACGTACATCAGCACACCTCCCCATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGNTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGc                                                                                                                                         >  1:164002/1‑300 (MQ=255)                                                                                                                                              cccccTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGNTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGc                                                                                                                                         >  1:600102/1‑300 (MQ=255)                                                                                                                                              cccccTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGc                                                                                                                                         >  1:397068/1‑300 (MQ=255)                                                                                                                                                   tttCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTNTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGaa                                                                                                                                    >  2:277999/1‑300 (MQ=255)                                                                                                                                                             aagTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCAt                                                                                                                          >  1:33566/3‑300 (MQ=255)                                                                                                                                                                  cATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGNGTCGCACCTGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTc                                                                                                                     >  1:173955/1‑300 (MQ=255)                                                                                                                                                                    tCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCAc                                                                                                                   <  2:794716/300‑1 (MQ=255)                                                                                                                                                                                                ggAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCggctggct                                                                                        >  2:487305/1‑299 (MQ=255)                                                                                                                                                                                                   aTAGCCGGCATGCGCTCGCCACAGGTCAGGATCAGGTCGATTCCTAGCGGGAACGAATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCTGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTTTAAAATCAGCGTCGCATCCGCTTCCGGCGCAACGGCGGTGCCGTTGGAAAGGGCGTTGAGCGGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGAcc                                                                                    <  2:180522/300‑1 (MQ=255)                                                                                                                                                                                                   aTAGCCGGCAGGCACTCGCCACAGTTCAAGATCAGGTCGATGCCGAGCGGGAACGGATGCTGGCGCTCGGTGAGTTCGTGCAGGATGCACTTCGGCAGCTGCGGAGCAATCATTCGTTCTTCTGCAATACCCGTACCAGTAAGGCGCAACATGCGTCCGCCGCTCATGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAGAGGGCTTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGAcc                                                                                    <  2:180407/300‑1 (MQ=255)                                                                                                                                                                                                   aTAGCCAGCAGGCACTCGCCACAAGTCAGGATCAGGTCGATGCCGAGCGGGAACAGATGCGTGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGTCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGAcc                                                                                    <  2:181225/300‑1 (MQ=255)                                                                                                                                                                                                                  tCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGTGAACGGATGCGGGCGCTCGTTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTa                                                                     <  2:33566/300‑1 (MQ=255)                                                                                                                                                                                                                   cGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCggct                                                                                            >  2:52652/1‑276 (MQ=255)                                                                                                                                                                                                                   cGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCggct                                                                                            <  1:52652/276‑1 (MQ=255)                                                                                                                                                                                                                     ccACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTNCGGCAATACCCGCACCAGTAAGGCGCAACATGNGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCAGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGa                                                                  >  2:632658/1‑300 (MQ=255)                                                                                                                                                                                                                     ccACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTNCGGCAATACCCGCACCAGTAAGGCGCAACATGNGTCCGCAGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGa                                                                  >  2:631756/1‑300 (MQ=255)                                                                                                                                                                                                                          ggTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAAc                                                              >  2:763148/1‑299 (MQ=255)                                                                                                                                                                                                                                 aTCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAATGCGCAACATGCGTCCGCCGCGCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGctc                                                      <  1:433995/300‑1 (MQ=255)                                                                                                                                                                                                                                 aTCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGctc                                                      <  1:433444/300‑1 (MQ=255)                                                                                                                                                                                                                                     ggTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAGGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGAACCCGCTTCCGGCGCAACGGCGGTGCAGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCATCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGt                                                  >  2:438781/1‑300 (MQ=255)                                                                                                                                                                                                                                     ggTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTAGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTAATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGGATGACAACGCAGGCTCTGGt                                                  >  2:439128/1‑300 (MQ=255)                                                                                                                                                                                                                                      gTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGtt                                                 >  1:453931/1‑300 (MQ=255)                                                                                                                                                                                                                                               gAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATc                                        >  2:837553/1‑300 (MQ=255)                                                                                                                                                                                                                                                         gtatgcgtgtgcTCGGTTAGTTCGTGCAGAGTGCACTTCGGCGGCTGCGGAGCGGTCGTTCGTTCTTCGGCAATACTCGTACCAGTAATGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTTTAAAATTAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCTGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAAt                              <  2:125506/290‑1 (MQ=255)                                                                                                                                                                                                                                                                    cTCGGTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCCATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGTGGTGCCGGTGGAAAGGGCGTTGAGCTGCTCGCTGGAAATCGCCTAATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGaa                     >  2:43798/1‑298 (MQ=255)                                                                                                                                                                                                                                                                      cGGTGAGTTCGTGCAGAATGCACTTAGGCAGCTGTGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCAGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGATGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGcc                 >  2:364043/1‑300 (MQ=255)                                                                                                                                                                                                                                                                        gTGAGTTCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGCcac               <  1:53222/300‑1 (MQ=255)                                                                                                                                                                                                                                                                             ttCGTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCTTCCGCCGCTCAGGCTGGCGACCGGTAAAATCAGCGTCGCACCCGCTTCCTGCGCAACTGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGCCACACCgg          <  2:355585/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                gTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCNGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGCCACACCGGCg        <  2:600102/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                gTGCAGAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCACAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGCCACACCGGCg        <  2:164002/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                     gAATGCACTTCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCATGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTTTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGCCACACCGGCGTGTCGt  >  2:16053/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                   GCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGGCAACGTACATCAGCACACCTCCACATGAGTGGTTCGCGGAATAGCCAGCAGGCGCTCGCCACAGGTCAGGATCAGGTCGATGCCGAGCGGGAACGGATGCGGGCGCTCGGTGAGTTCGTGCAGAATGCACTCCGGCAGCTGCGGAGCGATCATTCGTTCTTCGGCAATACCCGCACCAGTAAGGCGCAACATGCGTCCGCCGCTCAGGCTGGCGACCTGTAAAATCAGCGTCGCACCCGCTTCCGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGCCACACCGGCGTGTCGT  >  NC_000913/4320500‑4321075

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