Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,132,227	C→A	A798E (GCG→GAG)	metL →	Bifunctional aspartokinase/homoserine dehydrogenase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,132,227	0	C	A	100.0%	241.2 / NA	62	A798E (GCG→GAG)	metL	Bifunctional aspartokinase/homoserine dehydrogenase 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (28/34);  total (28/34)

GATGGTGCAACGGCTGGAAGCGGCCCGCGAAATGGGGCTGGTGCTGCGCTACGTGGCGCGTTTCGATGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGA  >  NC_000913/4131997‑4132440                                                                                                                                                                                                                                       |                                                                                                                                                                                                                      gATGGTGCAACGGCTGGAAGCGGCCCGCGAAATGGGGCTGGTGCTGCGCTACGTGGCGCGTTTCGATGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGt                                                                                                                                                                                                               >  1:642806/1‑239 (MQ=255)  atgcggCAAGGGCTGGAACCGGCCCGCGAAAGGGGGCTGGTGCTGCGCTACGGGGCGCGTTTCGATGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCCCTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGt                                                                                                                                                                                                               <  2:15445/233‑1 (MQ=255)    ggTGCAACGGCTGGAAGCGGCCCGCGAAATGGGGCTGGTGCTGCGCTACGTGGCGCGTTTCGATGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCgg                                                                                                                                                                                                            >  2:264895/1‑239 (MQ=255)       gCAACGGCTGGAAGCGGCCCGCGAAATGGGGCTGGTGCTGCGCTACGTGGCGCGTTTCGATGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGat                                                                                                                                                                                                          <  2:690453/238‑1 (MQ=255)               tGGAAGCGGCCCGCGAAATGGGGCTGGTGCTGCGCTACGTGGCGCGTTTCGATGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCaa                                                                                                                                                                                                     <  1:264895/235‑1 (MQ=255)                      gtcccGCGAAATGGGGCTGGTGCTGCGCTACGTGGCGCGTTTCGATGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGc                                                                                                                                                                                            <  1:281931/235‑1 (MQ=255)                        cccGCGAAATGGGGCTGGTGCTGCGCTACGTGGCGCGTTTCGATGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGc                                                                                                                                                                                            >  1:328676/1‑235 (MQ=255)                               aaTGGGGCTGGTGCTGCGCTACGTGGCGCGTTTCGATGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCgg                                                                                                                                                                                                 <  1:199167/223‑1 (MQ=255)                                               cgcTACGTGGCGCGTTTCGATGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCaa                                                                                                                                                                                                     >  2:145565/1‑203 (MQ=255)                                               cgcTACGTGGCGCGTTTCGATGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCaa                                                                                                                                                                                                     <  1:145565/203‑1 (MQ=255)                                                ccTACGGTGCGCTTTTCGATGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGTTCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTtaa                                                                                                                                                                 <  2:642806/237‑1 (MQ=255)                                                 tttctgggCGCGTTTCGATGCAAACGTTAAGGCGCGTTTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTtaata                                                                                                                                                               <  2:647098/233‑1 (MQ=255)                                                       ggCGCGTTTCGATGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTGATAATGcc                                                                                                                                                          <  2:77028/238‑1 (MQ=255)                                                       ggCGCGTTTCGATGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGcc                                                                                                                                                          <  2:328676/238‑1 (MQ=255)                                                                      cAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGtttt                                                                                                                                           <  1:161012/238‑1 (MQ=255)                                                                       aaCGGTAAAGCGCGGTTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCg                                                                                                                                         <  1:649331/239‑1 (MQ=255)                                                                          ggTAAAGCGAGGGAAGGCGTGGAAGCGGTGGTGGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAGGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGGGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTc                                                                                                                                          <  2:365009/235‑1 (MQ=255)                                                                                    gtgtAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTtaata                                                                                                                                                               >  1:493443/1‑204 (MQ=255)                                                                                    gtgtAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTtaata                                                                                                                                                               <  2:493443/204‑1 (MQ=255)                                                                                    gtgtAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGtt                                                                                                                              >  1:484962/1‑237 (MQ=255)                                                                                     tgtAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTACTGTTTTCGGGGTTTGGGGTTCt                                                                                                                           >  2:659968/1‑239 (MQ=255)                                                                                            gTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGaa                                                                                                                     <  1:659968/238‑1 (MQ=255)                                                                                               gAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAATCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGCCGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGGAATTTAATAAAGCCCGGTACTCATGCTTTCGGGATTATGGTTTCTAATGAAatat                                                                                                                 >  2:216771/1‑239 (MQ=255)                                                                                               gAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATGACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGGTTCTAATGAAatat                                                                                                                 >  2:564367/1‑239 (MQ=255)                                                                                                aaGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTc                                                                                                                                                  >  2:413076/1‑205 (MQ=255)                                                                                                aaGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTc                                                                                                                                                  <  1:413076/205‑1 (MQ=255)                                                                                                      tGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGCTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGaa                                                                                                           <  1:128188/238‑1 (MQ=255)                                                                                                                  atcCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGaa                                                                                                           <  2:138651/226‑1 (MQ=255)                                                                                                                  atcCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGaa                                                                                                           >  1:138651/1‑226 (MQ=255)                                                                                                                         ggCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGcc                                                                                                                                                          <  1:292236/172‑1 (MQ=255)                                                                                                                         ggCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGcc                                                                                                                                                          >  2:292236/1‑172 (MQ=255)                                                                                                                         ggCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGa                                                                                                >  1:83256/1‑230 (MQ=255)                                                                                                                         ggCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCgg                                                                                       >  1:598094/1‑239 (MQ=255)                                                                                                                         ggCATCACTGCTGCCGTGCGATAACGTCTTGGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGa                                                                                                <  2:83256/230‑1 (MQ=255)                                                                                                                                   ctgcCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGtt                                                                             >  2:613896/1‑239 (MQ=255)                                                                                                                                    tgcCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAAt                                                                                                          >  1:132597/1‑209 (MQ=255)                                                                                                                                    tgcCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAAt                                                                                                          <  2:132597/209‑1 (MQ=255)                                                                                                                                     gcCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAAtt                                                                                                                                                                    >  2:514419/1‑150 (MQ=255)                                                                                                                                     gcCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAAtt                                                                                                                                                                    <  1:514419/150‑1 (MQ=255)                                                                                                                                        gTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCaa                                                                                                                                                                                                     <  2:153572/114‑1 (MQ=255)                                                                                                                                        gTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCaa                                                                                                                                                                                                     >  1:153572/1‑114 (MQ=255)                                                                                                                                         tGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGaa                                                                                                                     >  2:279598/1‑193 (MQ=255)                                                                                                                                         tGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGaa                                                                                                                     <  1:279598/193‑1 (MQ=255)                                                                                                                                                       tGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaagaa                                                         >  1:87743/1‑239 (MQ=255)                                                                                                                                                         ccATCGGAAGCCTGTGGTAGCGCGATACCCCTATGGTGATCCGCGGACCGGGCGCTGGGCGGCGCGTCACCGCAGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGc                                                       <  2:415420/239‑1 (MQ=255)                                                                                                                                                            tCGAAAGCCGTGGGTATCGCGATAACCTTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAAt                                                    <  2:598094/239‑1 (MQ=255)                                                                                                                                                             cGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGc                                                                                                                                                                                                >  2:47593/1‑98 (MQ=255)                                                                                                                                                             cGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGc                                                                                                                                                                                                <  1:47593/98‑1 (MQ=255)                                                                                                                                                                   ccGCTGGTATCGCGATACCCCTCTGTGGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTa                                               <  2:484962/237‑1 (MQ=255)                                                                                                                                                                      cTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGaa                                                                                                                     >  2:380865/1‑164 (MQ=255)                                                                                                                                                                      cTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGaa                                                                                                                     <  1:380865/164‑1 (MQ=255)                                                                                                                                                                                         ctGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAg                                                                                                                                                                         >  2:563518/1‑93 (MQ=255)                                                                                                                                                                                         ctGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAg                                                                                                                                                                         <  1:563518/93‑1 (MQ=255)                                                                                                                                                                                          tGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTcc                      >  1:236286/1‑239 (MQ=255)                                                                                                                                                                                                tCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTtaata                                                                                                                                                               <  2:214187/96‑1 (MQ=255)                                                                                                                                                                                                tCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTtaata                                                                                                                                                               >  1:214187/1‑96 (MQ=255)                                                                                                                                                                                                    cgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgcttt            >  2:59828/1‑239 (MQ=255)                                                                                                                                                                                                       aCCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCa         >  2:726738/1‑239 (MQ=255)                                                                                                                                                                                                          tGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATAc      <  2:87743/239‑1 (MQ=255)                                                                                                                                                                                                               cTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGa  <  1:613896/238‑1 (MQ=255)                                                                                                                                                                                                                       gACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTAc                                              >  1:340682/1‑186 (MQ=255)                                                                                                                                                                                                                       gACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTAc                                              <  2:340682/186‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                      GATGGTGCAACGGCTGGAAGCGGCCCGCGAAATGGGGCTGGTGCTGCGCTACGTGGCGCGTTTCGATGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGA  >  NC_000913/4131997‑4132440

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