Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	705,013	T→C	L357S (TTG→TCG)	nagE →	N‑acetylglucosamine‑specific PTS enzyme IIABC component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	705,013	0	T	C	100.0%	138.6 / NA	42	L357S (TTG→TCG)	nagE	N‑acetylglucosamine‑specific PTS enzyme IIABC component
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/21);  total (21/21)

ACTGAGCCGCTGGAATTCCTGTTCATGTTCCTTGCTCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTTGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTATATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGTGCTTCTGGG  >  NC_000913/704772‑705257                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     aCTGAGCCGCTGGAATTCCTGTTCATGTTCCTTGCTCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCg                                                                                                                                                                                                                                             <  1:293805/251‑1 (MQ=255)    gAGCCGCTGGAATTCCTGTTCATGTGCCTTGCTCCGTGGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCGCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGc                                                                                                                                                                                                                                            <  1:30885/249‑1 (MQ=255)     aGCCGCTGGAATTCCTGTTCATGTTCCTTGCTCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATg                                                                                                                                                                                                                                         >  2:353656/1‑251 (MQ=255)           gggAATTCCTGTTCATGTTCCTTGCTCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAAc                                                                                                                                                                                                                                   >  2:152636/2‑251 (MQ=255)              aaTTCCTGTTCATGTTCCTTGCTCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTg                                                                                                                                                                                                                                <  1:62746/251‑1 (MQ=255)                  ccTGTTCATGTTCCTTGCTCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGaaaa                                                                                                                                                                                                                            <  1:125246/251‑1 (MQ=255)                                       gctgctGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaag                                                                                                                                                                                                       <  1:124683/251‑1 (MQ=255)                                           ctgTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaaga                                                                                                                                                                                                      >  1:405314/1‑248 (MQ=255)                                                                     cGGTATCAGCCTGTTTGTTGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCTTGGTGTTCATTTCGGTTATCCGCATGTTCAACCTGAAAACTCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACa                                                                                                                                                                         <  1:397997/251‑1 (MQ=255)                                                                               cTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTg                                                                                                                                                                >  2:3850/1‑250 (MQ=255)                                                                                     gtgGCAACTCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGGGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAATACGAGATCTTTACTGAAGAATCCAACATCAACACTGaagaag                                                                                                                                                          <  2:380591/250‑1 (MQ=255)                                                                                     gtgGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGaagaag                                                                                                                                                          <  1:3850/250‑1 (MQ=255)                                                                                             gctgctGGGTCTCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTgaag                                                                                                                                                                                  >  2:122266/1‑218 (MQ=255)                                                                                             gctgctGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTgaag                                                                                                                                                                                  <  1:122266/218‑1 (MQ=255)                                                                                                 ctgGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACt                                                                                                                                             >  2:319105/1‑251 (MQ=255)                                                                                                           aCGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGCCTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGcaaccaac                                                                                                                                   >  1:207407/1‑251 (MQ=255)                                                                                                                               gCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTATATTGCTGCGGTTggcgg                                                                                                               <  1:319105/251‑1 (MQ=255)                                                                                                                                               aCGCGTTGATGGATAACCTGCCGGGCGCCAGCCAGAACGTCTGGATGCTGCTGGTTATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTTGCAACCAACTATATTGCTGCGGTTGGCGGCACTGACAACCTGaaa                                                                                               <  1:40757/251‑1 (MQ=255)                                                                                                                                                        tGTATAACCTGCCGGCCGCCAGCCCGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTATATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGAc                                                                                      >  2:204770/1‑251 (MQ=255)                                                                                                                                                                     ggccgCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTATATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCAccc                                                                         <  2:302812/251‑1 (MQ=255)                                                                                                                                                                       ccgccAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTATATTGCTGCGGTTGGCGGCACTGACAACCTGCAAGCGATTGACGCCTGTATCACCCGt                                                                       >  2:17533/1‑251 (MQ=255)                                                                                                                                                                                 gAACGTCTGGATGCTGCTGGTGATGGGCGTTATTTTCTTCGCTATATCCCTCGTGGTGTTAAGTTCGGTTATCCGAATGTTCACCCTGAAACCGCCGGGTCGTGCAGATAAAGCAGCCGAGATCGTTACTGAAGAAGCCACCAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACAAACTATATTGCTGCGGTTGGCGGCACTGACAAcc                                                                                                    >  2:202876/1‑212 (MQ=255)                                                                                                                                                                                 gAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGGCTGACTCAACTGGCAACCAGCTATATTGCTGCGGTTGGCGGCACTGACAAcc                                                                                                    <  1:202876/212‑1 (MQ=255)                                                                                                                                                                                         ggATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTATATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGc                                                      <  1:315224/250‑1 (MQ=255)                                                                                                                                                                                           aTGCTGCTGGTGATGGGAGTTATCTTCTTCGCTATCTACTTCGTGGAGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTATATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTg                                                    >  1:264674/1‑250 (MQ=255)                                                                                                                                                                                           aTGCTGCTGGGGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTGCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAATAAGCCAACAGCAACACTGAAGAATGTCTGACTCAACTGGCAGCCAACGATAGTGCTGCGTTTGGCGGCACTGATAACCTGAAAGCGATTGACGCCTGTAGCACCCGTCTTCTCCTTACAGTGGCTGa                                                   <  2:27269/251‑1 (MQ=255)                                                                                                                                                                                             gctgctcattaggagcgacaTCTTCTTCGCTTTCTACTACGTGATGTTCCGTTCGGTTATCCGCATGTTCAACCTGACAACGACGGGTCGTCAAGCTACAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGACGAAGGTCTGACTCAACTGGAAACCAACTATATTGCTGCGGTTGGCGGCACTGACCACCTGAAAGCGATTGAAGCCTGTATCCCCCGTCTGCGCCTTACAGTGGCTGAct                                                 >  2:60432/20‑251 (MQ=255)                                                                                                                                                                                                gctgGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTATATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTg                                              >  1:330374/1‑251 (MQ=255)                                                                                                                                                                                                  tgGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAACCGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTATATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGcc                                            >  2:250320/1‑251 (MQ=255)                                                                                                                                                                                                   ggTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTATATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGcc                                            <  1:204770/250‑1 (MQ=255)                                                                                                                                                                                                    gTGATGGGCGTTATCTTCTTCGATATCTACTTCGTGGTGTTCCGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGCAGATAAAGAAGACGAGATCGTTACTGAAGCAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTCTATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCcg                                          >  2:168772/1‑251 (MQ=255)                                                                                                                                                                                                        tGGGCGTTATCTTCTTCGCTATCTACTTCGTAGTGATAAGTTCGGTTATCAGCATGATAAACTTGAACACGCCGGGTCGTGAAGATAAAGAAGACAAGATCGTTACCGAAGAAGACAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACAATATTGATGCGGTTGCCGGCACTGACAACCTGAATGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTgctcttgcccgagtt                                      >  1:276405/1‑237 (MQ=255)                                                                                                                                                                                                           gCGTTATCTTCTTCGCTATCGACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTATATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTAac                                   <  1:367893/251‑1 (MQ=255)                                                                                                                                                                                                             gTTATCTTCTTCGATATCTGCTTAGAGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAACACGCCGGGTCGTGAAGATAAAGCAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGACGAAGGTCTGACTCAACTGGAAACCAACTATATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTAacg                                  >  2:382235/1‑250 (MQ=255)                                                                                                                                                                                                                            aTCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTATATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTTGCTGACTCTGCCCGCGTTAACGATACGATGTGTaaa                    <  1:24076/249‑1 (MQ=255)                                                                                                                                                                                                                               tACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTATATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCt               >  1:360577/1‑251 (MQ=255)                                                                                                                                                                                                                                aCTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGATATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTATATTGCTGCGGTTTGCGTCACTGACAACCTGAAAGCGATTGACTCCTGTATCTCCCGGCTGCGTCTTACAGTGGGTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGt                 <  1:183965/248‑1 (MQ=255)                                                                                                                                                                                                                                 cTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTATATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTgg             >  1:244155/1‑251 (MQ=255)                                                                                                                                                                                                                                    cGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCGGGGTCGTGAAGATAAATAAGACGAGATCGGTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTTACTCAACTGGGAACCAACTATAGTGCTGCGGTTGGCTGCACTTACAACCTGAAAGCGATTGACGCCGGTATCACCCGTCTGCGCTTTACATTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGt           <  1:166747/250‑1 (MQ=255)                                                                                                                                                                                                                                    cGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTATATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGTg          >  1:211976/1‑251 (MQ=255)                                                                                                                                                                                                                                          gtTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTATATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGGCGCCTGTATCACCCGTCTGCGCCTTACAGTTGCTGACTGTGACCGCGGTAACGATACGATGTGTAAACGTCTGGGTGCTTCt     <  1:60432/250‑1 (MQ=255)                                                                                                                                                                                                                                            tCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTATATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGTGCTTCTggg  >  1:325875/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     ACTGAGCCGCTGGAATTCCTGTTCATGTTCCTTGCTCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTTGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTGACTCAACTGGCAACCAACTATATTGCTGCGGTTGGCGGCACTGACAACCTGAAAGCGATTGACGCCTGTATCACCCGTCTGCGCCTTACAGTGGCTGACTCTGCCCGCGTTAACGATACGATGTGTAAACGTCTGGGTGCTTCTGGG  >  NC_000913/704772‑705257

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