Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,785,550	G→A	P97P (CCG→CCA)	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,785,550	0	G	A	100.0%	178.4 / NA	53	P97P (CCG→CCA)	gpmM	2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/23);  total (30/23)

ATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGC  >  NC_000913/3785306‑3785790                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 atggctatCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGc                                                                                                                                                                                                                                              <  1:394352/249‑1 (MQ=255)       atCGTGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGgc                                                                                                                                                                                                                                      >  1:208729/1‑251 (MQ=255)              agaaCAGCCGGATAACGCCATTTTTAGTACTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGACATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGCCTGGCGCGGt                                                                                                                                                                                                                                 >  2:359389/1‑249 (MQ=255)                aaCAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGt                                                                                                                                                                                                                                 >  1:426328/1‑247 (MQ=255)                            aCGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAAc                                                                                                                                                                                                                 <  1:337192/251‑1 (MQ=255)                               ccATTTTTAGTGCTCAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGc                                                                                                                                                                                                               >  1:117355/1‑250 (MQ=255)                                     ttAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGGCGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCGGACTGGCTCGGTAGATAAAGCGAAAAACGCAggtaa                                                                                                                                                                                                         <  2:255202/250‑4 (MQ=255)                                       aGTGCTAAAAACCCAGTAAAGTATGCACTATGGGACAATCGTCCGCATACCCTGCTCACCGCTTCAAGTATGGAAGTAAGTCTGCCTGAAAGTCAGATGGGAAACTAAGAAATCGACCAAGATAACATGGGTGCCGGAAGCATCGTGTATCAGGAACTATCTCGTCTGGACGTTGCAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAACGAAAAACGCAggtaa                                                                                                                                                                                                         >  2:182789/1‑245 (MQ=255)                                        gTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGc                                                                                                                                                                                                      >  2:316068/1‑250 (MQ=255)                                                aCCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGTAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCCGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATCAAACGAAAAACGCAGGTAAAGCATTACACAt                                                                                                                                                                                             >  1:128666/1‑251 (MQ=255)                                                 ccgcGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGCCGCTTCCGGTATGGAAGTCGGTCTGCCTGACCGTCAGATTGTTAAATCCGCAGTAGGAAACGTTAACCTGGGCGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGCCGTTGAAGTCAATGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAGGCGAAATACGccagtaaatcagtacacgtt                                                                                                                                                                                            >  2:233891/4‑232 (MQ=255)                                                      gTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGTAATTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATggg                                                                                                                                                                                       <  1:289656/251‑1 (MQ=255)                                                       tAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGt                                                                                                                                                                                      >  2:205320/1‑251 (MQ=255)                                                         aTGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTCACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCt                                                                                                                                                                                    >  2:139444/1‑251 (MQ=255)                                                          tGGATGCACTGTGGGCCACTCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCCTCGTGTATCAGGACCTGACTCGTCTGGACGTTGAACTCAAAGATCGCGCTTTCTTTGCTAATCTAGTGCTGACTGGCGCGGTAGATAAAGCGAACAACGCAGGTAAAGCGGTACACATTATGGGTCTg                                                                                                                                                                                   >  2:58041/1‑251 (MQ=255)                                                           ggATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTg                                                                                                                                                                                   <  1:132069/250‑1 (MQ=255)                                                           ggATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGc                                                                                                                                                                                  >  2:406285/1‑251 (MQ=255)                                                                       ggCCAATCGTCCGCATGCCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAggcg                                                                                                                                                                      <  1:184790/251‑1 (MQ=255)                                                                                                 aCGCTTCCGGTCTGGAAGTCGGTCTGACTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGACACGAAGATCAcatc                                                                                                                                            >  2:59556/1‑251 (MQ=255)                                                                                                         ggTCTGGAAGTCGGTCTGCCTGCCCGTCCGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACGCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTGTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGAACACAGCCACGAAGa                                                                                                                                                   <  2:361375/236‑1 (MQ=255)                                                                                                         ggTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAGGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGa                                                                                                                                                   >  1:361375/1‑236 (MQ=255)                                                                                                             tGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTACAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGt                                                                                                                                 >  1:21572/1‑250 (MQ=255)                                                                                                                 aGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAAc                                                                                                                            <  1:306634/251‑1 (MQ=255)                                                                                                                     ggTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGc                                                                                                                        >  1:278536/1‑251 (MQ=255)                                                                                                                      gTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGc                                                                                                                        >  1:216167/1‑250 (MQ=255)                                                                                                                                 cGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAggc                                                                                                                                                                       >  1:138723/1‑192 (MQ=255)                                                                                                                                 cGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAggc                                                                                                                                                                       <  2:138723/192‑1 (MQ=255)                                                                                                                                 cGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACgcgg                                                                                                            >  1:122191/1‑251 (MQ=255)                                                                                                                                  gTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACgcggc                                                                                                           >  1:329258/1‑251 (MQ=255)                                                                                                                                    cAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGgcgc                                                                                                         >  2:406125/1‑251 (MQ=255)                                                                                                                                                      gTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTAGACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCa                                                                                       <  2:121398/251‑1 (MQ=255)                                                                                                                                                              cGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGAGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTCGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGCACGCGGCGCAGAAAAAATCTACCTGCCCGCATTcc                                                                               >  1:262234/1‑251 (MQ=255)                                                                                                                                                                          ggccggccgCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTcgcg                                                                   <  2:21572/250‑1 (MQ=255)                                                                                                                                                                            ccggccgCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGac                                                                 >  1:123910/1‑251 (MQ=255)                                                                                                                                                                                         gtATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCa                                                    <  1:395276/251‑1 (MQ=255)                                                                                                                                                                                           aTCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGt                                                  >  2:321877/1‑251 (MQ=255)                                                                                                                                                                                               ggACCTGCCTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCTCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATTGGGCTGCGCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCTATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTTCCCGCACGCATTCCGTGACGGGCTCGACACCCCGCCGCGGAGTGCTg                                              <  2:432369/251‑1 (MQ=255)                                                                                                                                                                                               ggACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAAGCCAGTGCTGACTGGCGCTGTAGATAAAGCGAAAAACGCAGGTAAATCAGTACACATTATGGTTCTGCTCGCCGCAGGCGGCGTACACAGCCACGAAGAGCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCATAAAAAATCTACCTGCACGCATTCCTTGACGGGCGCGACACCCCGCTGCGCAGTGCTg                                              <  2:432372/251‑1 (MQ=255)                                                                                                                                                                                                gACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAACAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTg                                              >  1:101811/1‑250 (MQ=255)                                                                                                                                                                                                    tGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATcc                                         >  1:275578/1‑251 (MQ=255)                                                                                                                                                                                                      acgcGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTc                                       <  2:156849/248‑1 (MQ=255)                                                                                                                                                                                                         cGTCTGGACGTTGAAATCAAAGATCGCTCTATCTTTGCTAATCCAGTTCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATTGGTCTGCTCTCCGCAGGCGTCGTACACAGCCACGAAGATCACATCATGTCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGc                                     <  1:47706/250‑1 (MQ=255)                                                                                                                                                                                                         cGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGc                                     <  1:316068/250‑1 (MQ=255)                                                                                                                                                                                                          gtgtGGACGTTGAAATAAAAGATCGCTCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATGAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTCCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCAGGGCGATGGTAGAACTTGCAGCTGAACGCGGCTCAGAAAAAATCTACCTGCACTCATTCCTGGACGGTCGCGACACCTCGCCGCGCAGTGCTGAATCCTCGCTg                                   <  2:376139/248‑1 (MQ=255)                                                                                                                                                                                                            cTGGACGTTGAAATCAAAGCTCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTCCACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTc                                                                                >  2:32580/1‑204 (MQ=255)                                                                                                                                                                                                            cTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACAGTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACAGCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCa                                                                                                        <  2:390532/180‑1 (MQ=255)                                                                                                                                                                                                            cTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACGGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTc                                                                                <  1:32580/204‑1 (MQ=255)                                                                                                                                                                                                            cTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGAATGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCa                                                                                                        >  1:390532/1‑180 (MQ=255)                                                                                                                                                                                                              ggACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGaaa                                >  1:373766/1‑250 (MQ=255)                                                                                                                                                                                                                        aTCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTgaaga                     <  1:406125/251‑1 (MQ=255)                                                                                                                                                                                                                           aaaGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCCGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGAGACACCCCGCCGCGCAGTGCTGAATACTCGCTGAAAAAATTTgaagaga                   >  1:391970/1‑248 (MQ=255)                                                                                                                                                                                                                            aaGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGagaaa                 <  2:34558/251‑1 (MQ=255)                                                                                                                                                                                                                                            tGCTAAACCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAgggt  <  2:77194/250‑2 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 ATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGC  >  NC_000913/3785306‑3785790

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