Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,786,535	G→A	E426K (GAA→AAA)	gpmM →	2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,786,535	0	G	A	100.0%	214.1 / NA	62	E426K (GAA→AAA)	gpmM	2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/40);  total (22/40)

GAAAAATATGCCCACGTTACTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGG  >  NC_000913/3786262‑3786821                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                               gaaaAATATGCCCACGTTACTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCTAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTc                                                                                                                                                                                                                                                                       <  1:459063/299‑1 (MQ=255) gaaaAATATGCCCACGTTACTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTcc                                                                                                                                                                                                                                                                      <  1:573480/300‑1 (MQ=255)      atatGCCCACGTTACTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCgtggg                                                                                                                                                                                                                                                                 <  1:11686/300‑1 (MQ=255)       tatGCCCACGTTACTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCgtggg                                                                                                                                                                                                                                                                 <  2:422481/299‑1 (MQ=255)             cACGTTACTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATTAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTACGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGTGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGGCaa                                                                                                                                                                                                                                                          >  2:147070/1‑300 (MQ=255)              aCGTTACTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGAcgac                                                                                                                                                                                                                                                         <  1:649271/300‑4 (MQ=255)              aCGTTACTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAAc                                                                                                                                                                                                                                                         <  1:649241/300‑1 (MQ=255)                                acgttGGCGTATAATAGTCGTTCAAATGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAACTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGAcc                                                                                                                                                                                                                                       <  2:134051/296‑1 (MQ=255)                                aCGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGAGCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGAcc                                                                                                                                                                                                                                       <  2:370727/300‑1 (MQ=255)                                  ggTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCAc                                                                                                                                                                                                                                     <  2:283379/300‑1 (MQ=255)                                     ggCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGt                                                                                                                                                                                                                                  >  1:430055/1‑300 (MQ=255)                                         tAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTACTGACCACCGCTGACCACGGTaa                                                                                                                                                                                                                                >  2:274491/1‑298 (MQ=255)                                         tAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGAAATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGATAAc                                                                                                                                                                                                                               >  2:274391/1‑299 (MQ=255)                                          agaagaGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGc                                                                                                                                                                                                                             <  1:242013/300‑1 (MQ=255)                                                gTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAGGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACCGCTGATCACCGCTGACCACGGTAACGCTGAgtc                                                                                                                                                                                                                       >  2:445525/1‑298 (MQ=255)                                                    ttCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGAta                                                                                                                                                                                                                   <  1:371189/300‑2 (MQ=255)                                                    ttCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCGACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGACATCAAGAGCGGCAGATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATAGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAATCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTAATCACCGCCGACCACGGTAACGCTGAGCAGAta                                                                                                                                                                                                                   <  2:435610/300‑2 (MQ=255)                                                       aaaGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGagacg                                                                                                                                                                                                                 <  1:724922/299‑5 (MQ=255)                                                       aaaGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGAtacg                                                                                                                                                                                                                 <  1:251096/299‑4 (MQ=255)                                                         aGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGa                                                                                                                                                                                                              <  1:458238/300‑1 (MQ=255)                                                         aGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGa                                                                                                                                                                                                              <  2:704306/300‑1 (MQ=255)                                                                          tGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAgg                                                                                                                                                                                             <  1:369434/300‑1 (MQ=255)                                                                             tCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCAGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGTTAGGTCACACCTGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCTTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAATGCTGAGCAGATACGCGATCCGGCAACGGGTCAggc                                                                                                                                                                                            <  2:69309/298‑1 (MQ=255)                                                                               aaCTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAggcacac                                                                                                                                                                                        <  1:573935/300‑1 (MQ=255)                                                                                  tCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATNAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAggcacacac                                                                                                                                                                                      <  1:525227/299‑1 (MQ=255)                                                                                  tCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAggcacacac                                                                                                                                                                                      <  1:413503/299‑1 (MQ=255)                                                                                  tCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTAGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGCCACACCTGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAggcacacac                                                                                                                                                                                      <  2:343862/299‑1 (MQ=255)                                                                                       gAAAGTGGCTACCTACGATCTGCAACCGTAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGGAACTATCCGAACGGCGACATGGTAGGTCANACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAggcacacacggcac                                                                                                                                                                                 <  2:235110/299‑1 (MQ=255)                                                                                       gAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAggcacacacggca                                                                                                                                                                                  <  1:510863/298‑1 (MQ=255)                                                                                           gTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGANGCGCTGGATCACTGCGTGAAAGAAGTCGCGAANGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACAcc                                                                                                                                                                            <  2:544649/300‑1 (MQ=255)                                                                                           gTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACAcc                                                                                                                                                                            <  1:405452/300‑1 (MQ=255)                                                                                             ggCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCAAACACCaa                                                                                                                                                                          >  2:445590/1‑300 (MQ=255)                                                                                                           gaacacaggAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAggcacacac                                                                                                                                                                                      <  2:246660/267‑1 (MQ=255)                                                                                                                   gAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTa                                                                                                                                                     >  2:608145/1‑299 (MQ=255)                                                                                                                          gCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACAACAACCTGACAGTTCCGCTGATTTACGTTGGt                                                                                                                                              >  2:605924/1‑299 (MQ=255)                                                                                                                               gtagagCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACTACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTGAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAg                                                                                                                                        <  2:137340/298‑1 (MQ=255)                                                                                                                                  gagCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAAc                                                                                                                                     >  1:106109/1‑300 (MQ=255)                                                                                                                                                    gTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAggc                                                                                                                   <  1:632025/300‑1 (MQ=255)                                                                                                                                                                           aTACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCgcgc                                                                                             >  1:272937/1‑299 (MQ=255)                                                                                                                                                                                acacCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCa                                                                                       >  1:444764/1‑300 (MQ=255)                                                                                                                                                                                   ccATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCAtgt                                                                                    >  2:422078/1‑300 (MQ=255)                                                                                                                                                                                   ccATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAAAGGGTCAGGCACACACGGCACACACCAACCTGCAAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCAtgt                                                                                    >  2:376226/1‑300 (MQ=255)                                                                                                                                                                                          cTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGAACATGTTGTCGCt                                                                             >  2:679261/1‑300 (MQ=255)                                                                                                                                                                                                cTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGAt                                                                          >  1:717207/1‑297 (MQ=255)                                                                                                                                                                                                      gAACTGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCTGTTAAAGCGGTTGAAGCGCTGGATCACTGCTTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGa                                                                 <  2:157454/300‑1 (MQ=255)                                                                                                                                                                                                       aaCGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCCACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGaa                                                                <  2:505280/300‑1 (MQ=255)                                                                                                                                                                                                       aaCGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGaa                                                                <  1:657177/300‑1 (MQ=255)                                                                                                                                                                                                                     aGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGAt                                                  >  2:728333/1‑300 (MQ=255)                                                                                                                                                                                                                     aGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGCTGTCGCTGATGGGTATGGAAATCCCGCAAGAGAt                                                  >  1:571049/1‑300 (MQ=255)                                                                                                                                                                                                                           cacCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACt                                              >  1:547594/1‑298 (MQ=255)                                                                                                                                                                                                                                gggTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCTTGAAAGAAGTCGCTAAAGCTGTTTAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCTGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGc                                       <  2:161400/300‑1 (MQ=255)                                                                                                                                                                                                                                        gAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTc                               <  2:180520/300‑1 (MQ=255)                                                                                                                                                                                                                                        gAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTc                               <  1:588748/300‑1 (MQ=255)                                                                                                                                                                                                                                          aGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCAt                             <  1:728333/300‑1 (MQ=255)                                                                                                                                                                                                                                                ggTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCATACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGa                       >  1:55624/1‑300 (MQ=255)                                                                                                                                                                                                                                                ggTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACTTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCAGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGa                       >  1:55471/1‑300 (MQ=255)                                                                                                                                                                                                                                                ggTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGa                       >  2:424261/1‑300 (MQ=255)                                                                                                                                                                                                                                                       gCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCNCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAAtcc                <  1:147070/300‑1 (MQ=255)                                                                                                                                                                                                                                                       gCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAAtcc                <  1:404128/300‑1 (MQ=255)                                                                                                                                                                                                                                                       gCGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGGTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCATGCACACACTGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAAtcc                <  2:700444/300‑1 (MQ=255)                                                                                                                                                                                                                                                        cGGTTGAAGCGCTGGATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACCACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAACCTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATAGTGGAATAAtccc               >  2:505889/1‑300 (MQ=255)                                                                                                                                                                                                                                                                     ggATCACTGCGTGAAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATACGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATTAgggg  <  2:106109/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                               GAAAAATATGCCCACGTTACTTTCTTCTTCAACGGTGGCGTAGAAGAGTCGTTCAAAGGCGAAGATCGCATTCTGATCAACTCACCGAAAGTGGCTACCTACGATCTGCAACCGGAAATGAGCTCCGCAGAGCTGACCGAAAAACTGGTTGCGGCCATCAAGAGCGGCAAATACGACACCATCATCTGTAACTATCCGAACGGCGACATGGTAGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGG  >  NC_000913/3786262‑3786821

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