Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,510,711	C→T	R249R (AGG→AGA)	fecE ←	ferric citrate ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,510,711	0	C	T	100.0%	176.6 / NA	52	R249R (AGG→AGA)	fecE	ferric citrate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/27);  total (25/27)

GATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCA‑CAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAG  >  NC_000913/4510475‑4510953                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    gatgtccggtctcccGACTGGTTTGAAATCAGCGCCCAAGCCTGTACAGCACTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCATTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCATGTTCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCg                                                                                                                                                                                                                                        <  2:447905/236‑1 (MQ=255)  attgCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACNGGc                                                                                                                                                                                                                                      >  2:593249/1‑251 (MQ=255)      ccGTGATCTCGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCGGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCGAGATACCGGCTCgg                                                                                                                                                                                                                                  <  2:68311/251‑1 (MQ=255)       cGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCggg                                                                                                                                                                                                                                 >  1:267323/1‑251 (MQ=255)                       ttGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGGGTATTTCCGCTTcc                                                                                                                                                                                                                 >  1:327481/1‑251 (MQ=255)                                 cgcCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACAGCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATa                                                                                                                                                                                                       <  1:495825/251‑1 (MQ=255)                                      aaGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGtt                                                                                                                                                                                                  <  1:495730/251‑1 (MQ=255)                                          cTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTTACTACTGTtctc                                                                                                                                                                                               >  1:592359/1‑250 (MQ=255)                                          cTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCAGCTTCCACGCTGAATACTGTtctc                                                                                                                                                                                               >  1:574102/1‑250 (MQ=255)                                                     tGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTggg                                                                                                                                                                                   >  1:200245/1‑251 (MQ=255)                                                           gCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTCTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGCTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGtcatc                                                                                                                                                                             <  2:65738/251‑1 (MQ=255)                                                                gATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTTAATACTTTTCTCAGCAATCCTGGGGTCATCACct                                                                                                                                                                         <  1:499250/250‑1 (MQ=255)                                                                   tATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAAGCCTGGGGTCATCACCtctt                                                                                                                                                                      <  1:578151/250‑1 (MQ=255)                                                                    aTGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCtcttct                                                                                                                                                                    <  1:305875/251‑1 (MQ=255)                                                                     tGAACCTGTCAGCCTACAGCAAGCGGATGGACGATGTGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATCCCGGCTCGGGGGGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCCATCCTTGGGTCATCACCTCTTCTg                                                                                                                                                                   >  1:571615/1‑251 (MQ=255)                                                                         ccTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATACTGGGGGCATCCCCTCTTCTTgtgt                                                                                                                                                               >  2:582275/1‑251 (MQ=255)                                                                               aGCNTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACAGCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTg                                                                                                                                                         <  2:608646/251‑1 (MQ=255)                                                                                     aagcaagcGGATGGACGATGAGTATTGGTAATCTTTCTGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCAt                                                                                                                                                   <  2:221701/251‑1 (MQ=255)                                                                                     aagcaagcGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCAGGAAAAGTTCAGCACCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTAAGCAGTGGATTGGTTCTTGAGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCCGGGCGTATTTCCGCTTCCCCGCAGAATACTGTTCTACGCAATCCTGGGGTCACTACCtcttctagactggcttacgccat                                                                                                                                                   >  1:91673/1‑234 (MQ=255)                                                                                                tGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGGCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGt                                                                                                                                        >  1:221383/1‑251 (MQ=255)                                                                                                     gATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGGGTGCCTTGCGCCATAACATGTCCGTTTGcc                                                                                                                                   >  1:406280/1‑251 (MQ=255)                                                                                                     gATGAGTATTGGTAATCTGTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATCCCGGCTCGGGGTGTAGTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGcc                                                                                                                                   <  2:561288/251‑1 (MQ=255)                                                                                                        gAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCAtt                                                                                                                                <  2:192471/251‑1 (MQ=255)                                                                                                          gTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTac                                                                                                                              <  2:327481/251‑1 (MQ=255)                                                                                                          gTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGGCGTGGTTGGGTTCTTGCGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCGCGGGGTGTATTTCCGCTTCCTCGCTGAATACTGTTCTCAGCAATCCTGGTGTCATCCCCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTac                                                                                                                              <  2:285332/251‑1 (MQ=255)                                                                                                              tGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGGGTGCCTTGCGCCATAACATGTCCGTTTGCCATTaccacc                                                                                                                          >  1:178458/1‑251 (MQ=255)                                                                                                                 tAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGgtgt                                                                                                                                                                                                                             >  2:447979/1‑149 (MQ=255)                                                                                                                 tAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGgtgt                                                                                                                                                                                                                             <  1:447979/149‑1 (MQ=255)                                                                                                                      tttCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCAAATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACCGTTCTCAGCAATCCCGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCTTAACATGTCCGTTTGCCATTAACACCAGTTGATc                                                                                                                  >  1:443162/1‑251 (MQ=255)                                                                                                                           gagTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACCCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTTCTCCATAACATGGCCGTTTGCCATTACCACCCGTTGATCGCAGt                                                                                                             >  2:467470/1‑251 (MQ=255)                                                                                                                              tCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTAACTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCCGGGGTCATCACCTCTTCTTGGGTGCCTTGCGCCATAACATGTCCGGTTGCCATTACAACCAGTTGATCGCAGTAcc                                                                                                          >  2:479120/1‑251 (MQ=255)                                                                                                                                     aaGTGCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGc                                                                                                   <  2:428846/251‑1 (MQ=255)                                                                                                                                         tCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTg                                                                                                <  2:406280/250‑1 (MQ=255)                                                                                                                                             ccccAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGTTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTaa                                                                                           <  2:101348/251‑1 (MQ=255)                                                                                                                                                 aGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTGCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATGTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGTTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTc                                                                                       <  2:658504/251‑1 (MQ=255)                                                                                                                                                       aaCAGGCTTGCTGGCGCCATTCCAGTGTCTTTCAGTCGTTGTTTGGTTCTTACGGCCTGTGCACTCTACCTCATTAGGCACATCGGTCTGCCAGATGCCTGCTCGGGGTGTATTTCCGCTTCCACTCTTAATACTGTTCTCAGCAATCCTGGTGTCATCCCCTCTTCTTGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGCTAAAGGTCGTgca                                                                                  <  2:227222/250‑1 (MQ=255)                                                                                                                                                       aaCAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTgcag                                                                                 <  1:433369/251‑1 (MQ=255)                                                                                                                                                       aaCAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTgcag                                                                                 <  1:658896/251‑1 (MQ=255)                                                                                                                                                       aaCAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGGCGTgc                                                                                   >  1:129057/1‑249 (MQ=255)                                                                                                                                                        aCAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGGATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTgcagc                                                                                >  1:361390/1‑251 (MQ=255)                                                                                                                                                                  gggNGCCAGTCCAGTGTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCA‑CAGCGacca                                                                     <  2:592359/250‑1 (MQ=255)                                                                                                                                                                           tccagtTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACCGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACCTGTCCGTTTGCCATTACCCCCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCA‑CAGCGGCCACCGTTTTcc                                                            >  1:63908/1‑251 (MQ=255)                                                                                                                                                                            ccagtTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTTCGACTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCCCCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTTCCATTTCCCCCAGTTGATCTCAGTACCGGCTAGCCTGATTAATGTCGGGCCGCA‑ACGCGACCACCGTTTTccc                                                           >  2:192160/1‑251 (MQ=255)                                                                                                                                                                                   cattcaGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGCTTTAGGTCGTGCCGCC‑CAGCGACCACCGTTTTCCCCTGAGTc                                                    >  2:554697/1‑251 (MQ=255)                                                                                                                                                                                   cattcaGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCAACA‑CGaccacc                                                                   >  2:286033/1‑236 (MQ=255)                                                                                                                                                                                   cattcaGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCA‑CAGCGaccacc                                                                   <  1:286033/236‑1 (MQ=255)                                                                                                                                                                                           cGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCA‑CAGCGACCACCGTTTTCCCCTGAGTCCGGAGtt                                             <  1:108315/250‑1 (MQ=255)                                                                                                                                                                                           cGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCA‑CAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTc                                            <  2:200245/251‑1 (MQ=255)                                                                                                                                                                                                              ccTGTGCAATCTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCA‑CAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAgg                         >  1:47698/1‑251 (MQ=255)                                                                                                                                                                                                                        cTACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCCCGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCC‑CCGCGACCACCGTTTTCCCCCGAGGCCGGAGTTCCCCCATCAACCGCACCAGGTCCACCTggg               >  2:542876/1‑250 (MQ=255)                                                                                                                                                                                                                         tACCTCATTAGGCACATCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCCGTACCGGCTAGCCTGATTAAGGTCGGGCAGCA‑CAGCGACCACCGCTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTgggg              >  2:239786/1‑249 (MQ=255)                                                                                                                                                                                                                                     cacaTCGGTCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCA‑CAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAg  <  2:221383/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    GATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCA‑CAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAG  >  NC_000913/4510475‑4510953

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: ‑