Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,648,617	T→C	D1224G (GAT→GGT)	yfhM ←	alpha2‑macroglobulin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,648,617	0	T	C	100.0%	113.1 / NA	38	D1224G (GAT→GGT)	yfhM	alpha2‑macroglobulin
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/21);  total (17/21)

ACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCC‑GTC‑AGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAACTCTTTGATATAACGTGCGATGTTCAGCGGTGGTTTGCCGCTCAACAACAATTGCCCTTCCAGCGTAACAGGCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCT  >  NC_000913/2648376‑2648863                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       ccTGGCTTTGATATTATCCGTGTAGGGGACCGACATCAGGCCCGGATGTAGAAAAGAGCCGAGCAGACTCTCTTTACCACGGTAAAGGGCGTCTGGCGGCACGCGGAAACCCTTTCCGCAGGCGGGGACGAGGAAACCAATCACG‑AAGCCCTTC‑AGCTATTTTTCTTCATCACCGTTTTTACCCCAAAGCGAAAGGCCGCCGTTATCGCGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGAt                                                                                                                                                                                                                                                <  2:16900/249‑1 (MQ=255)    ggCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCC‑GTC‑AGCCAGTACTCTTCGTCAACGCTTTTATACCACAGCGCAACGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCAc                                                                                                                                                                                                                                            >  1:211118/1‑251 (MQ=255)               tATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCC‑GTC‑AGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTtctc                                                                                                                                                                                                                                  >  1:55301/1‑250 (MQ=255)                 tCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCC‑GTC‑AGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCtcatc                                                                                                                                                                                                                               >  2:341765/1‑251 (MQ=255)                            gttcgtcatcatcccagatgcTTGTAAATAGCGCAGCAGACGCTCATTCCCCCGGTCAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCACGCGCGGACCAGGAAATCCATCCCGTAAGC‑‑GTCCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGAAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATGCCGACCGATGCACGGCGTTTCTCATCACTGTCGCCtt                                                                                                                                                                                                                    <  2:330672/231‑1 (MQ=255)                             atagactccatgcctggttctgttAAATAGCCCAGCAGTCGTTCATTCCCCCGGTTAATGGCTTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGTAAATCCATCACGTAATCC‑GCC‑AGCCAGTACTCTTCGTCCCCGTTTTTACCCCACAGCGCAAAGCCGCCGTTATCACGGTGCATGTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCttt                                                                                                                                                                                                                   <  2:395680/228‑1 (MQ=255)                               cGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTACCCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCC‑GTC‑AGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGa                                                                                                                                                                                                                 <  1:121574/251‑1 (MQ=255)                                     catGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGCTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCC‑GTC‑AGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAACGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCTCTGTCGCCTTTGCTGccc                                                                                                                                                                                                            >  2:298338/1‑250 (MQ=255)                                          ccGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGCCCAGGAAATCCATCACGTAAGCC‑GTC‑AGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGcc                                                                                                                                                                                                      <  1:324383/251‑1 (MQ=255)                                            ggATCTTAAAAATAGCGCGGCAGACGCTCATTCCCCCGGTTAATGGCGTGTGTCGGCCCGCTGTAACCCTGTTCGCCTGCGGGGCCAAGGAAACCCAGCCCGTAAGCC‑GCC‑AGCCAGTACTCTTCGCCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTg                                                                                                                                                                                                    <  2:329908/251‑1 (MQ=255)                                              aTCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGCAATCCATCACGTAAGCC‑GTC‑AGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAACTGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGcc                                                                                                                                                                                                  >  2:155030/1‑250 (MQ=255)                                                   gTAAATAGCGCAGCAGCCGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCC‑GTC‑AGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTg                                                                                                                                                                                              <  1:5857/250‑1 (MQ=255)                                                            gcagcaGCCGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCC‑GTC‑AGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTg                                                                                                                                                                                    <  1:99751/251‑1 (MQ=255)                                                                    cGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAAACATGTTCGCATGCGCGGACCAGGAAATACATCACGTAAGCC‑GTC‑AGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCAGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGACGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGACCAACGCCTGCAGTTGGGCTGCGTTGGaatac                                                                                                                                                                              >  1:264385/1‑248 (MQ=255)                                                                                      aggggtgcgcgcagCACGTTGTACCCGTGTTCGCCTCGGCGGACCAGAAAATCCATCCCGTAAGCC‑GTC‑ACCCAGTCATCTGCGTCCCCGTTTTTCGCCCACCGCGCAAGGCCGCCGTTATCACGGTGCATTTGCAGCAAACGGGAAATGCCGATTCCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCGGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCt                                                                                                                                                          <  2:102771/238‑1 (MQ=255)                                                                                      aTGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCC‑GTC‑AGCCAGTACTCTTCGTCACCGTTTTTATCCAACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTAGAAAGTGACGGAAACAGGCCGCt                                                                                                                                                          <  1:106697/251‑1 (MQ=255)                                                                                             cTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCC‑GTC‑AGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGtt                                                                                                                                                   >  1:329783/1‑251 (MQ=255)                                                                                                              cccTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCC‑GTC‑AGCCAGTCCTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTa                                                                                                                                  <  1:178083/251‑1 (MQ=255)                                                                                                                  gTTCGCCGGCGCGGACCAGGAAACCCATCACGTAAGCC‑GCC‑AGCCAGTCCTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACgg                                                                                                                               <  1:135007/250‑1 (MQ=255)                                                                                                                       ccTGCGCGGACCAGGAAATCCATCACGTAAGCC‑GTC‑AGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGGTTGCTCAAGACAGCCGTACGGATACGc                                                                                                                         >  1:371272/1‑251 (MQ=255)                                                                                                                             cgGACCAGGAAATCCATCACGAAAGCC‑GTC‑AGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTaa                                                                                                                   <  1:254424/251‑1 (MQ=255)                                                                                                                              ggACCAGGAAATCCATCACGTAAGCC‑GTC‑AGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAAc                                                                                                                  >  1:373476/1‑251 (MQ=255)                                                                                                                                      aaaTCCATCACGTAGGCC‑GGC‑AGCCAGTCCTCTTCGTCGCCGTTGTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAACTCTTTGat                                                                                                          <  1:405558/251‑1 (MQ=255)                                                                                                                                                   aaGCC‑GTC‑AGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGAAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAACTCTTTGATATAACGTGCGAt                                                                                              <  1:277091/250‑1 (MQ=255)                                                                                                                                                             aGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAACTCTTTGATATAACGTGCGATGTTCAGCgg                                                                                     >  1:177155/1‑251 (MQ=255)                                                                                                                                                               ccAGTCCTCTTCGTCACCGTTTTTATCCCACAGCGCACAGCCGCCGTTATCACGTTGCATTTGCATCAAACGGGAAATGCCGCTACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAACTCTTTGATTTAACGTGCGCTGTTCAGCgggg                                                                                   >  2:182981/1‑249 (MQ=255)                                                                                                                                                                            tCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAACCGGGAAATGCCGATACCGACCGATGAACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAACTCTTTGATATAACGGGCGATGTACAGCGGTGGTTTGCCGCTca                                                                       >  2:240297/1‑250 (MQ=255)                                                                                                                                                                             ccccttttttGTCCCACAGAGAAAGGCCGCCGTTATCACGTTGAATTTGCAGCAATCGTGAAATGCCGATACCGACCGAGGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTGGGGCTGCGTTGGTATAAAGTGGCGGAAACAGGCCGCTGGCGGTTTGGGTAAGACAGCCTTACGGATACGCTTTTAACGCTTTGATATTACGTGCGATTTTCAGCGGTGGTTTGCCGCTcaac                                                                     <  2:272199/246‑1 (MQ=255)                                                                                                                                                                                        cccccacgcgaacaaccgccgTTATCACGTTGCATTTGCAGCAAACTGGAAATGCAGAGTCCGACCGATGCACGGCGTTTCTCATCACTGTTGCCTTTGATGCCCAACGCCTGCTGTGGGGCTGCGTTGTTATAAAGTGACGGAAACAGGCCGCGGGCGGTTTGCTCAAGACAGTCGGACGGTTACGCGGGTAACTCGTGGATATAGCGTGGGATGTTCAGCGGTGGTTTTCCGCTCAACAACAATTGccc                                                          <  2:155214/236‑1 (MQ=255)                                                                                                                                                                                         cccACAGAGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACCTCCGTACGGATACGCTTTTAACTCTTTGATAAAACGTGCGAAGGTAAGCGGAGGTTTGCCGCTCAACAACAATTGACCt                                                         >  1:439310/1‑251 (MQ=255)                                                                                                                                                                                          ccACAGCGCAAAGCCGCCGTTATCACTTTGCATTTGCAGCAACCGGTAAATGCCTATCCCGACCGATGCACGGCGTTTCGCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGATGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAACTCTTTGATATAACGTGCGATGTTCAGCGGTGGTTTGCCGCTCAACAACAATTGCCCtt                                                        <  2:384860/251‑1 (MQ=255)                                                                                                                                                                                               ccgcAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACATCCGTACGGATACGCTTTTAACTCTTTGATATAACGGGCGATGTTCAGCGGTGGTTTGCCGCTCAACAACAATTGCCCTTCCAGc                                                   <  1:211997/250‑1 (MQ=255)                                                                                                                                                                                                                     ttGCATTTGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACATCCGAACGGATACGCTTTTAACTCTTTGATATAACGTGCGATGGTCAGCGGTGGTTTGCCGCTCAACAACAATTGCCCTTCCAGCGTAACAGGCGAGCAGTTTTGcc                             >  2:241624/1‑250 (MQ=255)                                                                                                                                                                                                                        cATTTGCAGCAAACGGGAAATGCCGATACCGACCGATTCACGGCGTTTCTCATCACTTTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAACTCTTTGATATAACGTGCGATGTTCAGCGGTGGTTTGCCGCTCAACAACAATTGCCCTTCCAGCGTAACAGGCGAGAAGTTTTGCAatc                          <  2:62713/251‑1 (MQ=255)                                                                                                                                                                                                                            tGCAGCAAACGGGAAATGCCGATACCGACCGATGCACGGCGTCTCTCATCCCTGTCGCATTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGACACAGGCCGATGGCGGTTTGCTAACGCCAGCCGTACGGATACGCTTTTAACTCTTTAATATAACGTGCGATGTTCAGCGGTGGTTTGACGCTCAACAACAATTGCCCTTCCAGCGTAACAGGCGCGAAGTTTTGCAatccatc                      >  1:369606/1‑251 (MQ=255)                                                                                                                                                                                                                                  aaaCGGGAAATGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAACTCTTTGATATAACGTGCGATGTTCAGCGGTGGTTTGCCGCTCAACAACAATTGCCCTTCCAGCGTAACAGGCGAGAAGTTTTGCAATCCATCCGCCgg                <  1:317829/251‑1 (MQ=255)                                                                                                                                                                                                                                         aaaTGCCGATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAACTCTTTGATATAACGTGCGATGTTCAGCGGTGGTTTGCCGCTCAACAACAATTGCCCTTCCAGCGTAACAGGCGAGAAGTTTTGCAATCCATCCGCCGGAATCGcc         >  1:431763/1‑251 (MQ=255)                                                                                                                                                                                                                                                gATACCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAACTCTTTGATATAACGTGCGATGTTCAGCGGTGGTTTGCCGCTCAACAACAATTGCCCTTCCAGCGTAACAGGCGAGACGTTTTGCAATCCATCCGCCGGAATCGCCCATGtct  >  1:234801/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       ACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCC‑GTC‑AGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATACGCTTTTAACTCTTTGATATAACGTGCGATGTTCAGCGGTGGTTTGCCGCTCAACAACAATTGCCCTTCCAGCGTAACAGGCGAGAAGTTTTGCAATCCATCCGCCGGAATCGCCCATGTCT  >  NC_000913/2648376‑2648863

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