Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	23,038	A→G	A216A (GCA→GCG)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	23,038	0	A	G	93.1%	187.7 / 2.8	73	A216A (GCA→GCG)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base G (41/26);  total (44/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

AACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTT‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCAC  >  NC_000913/22800‑23274                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             aaCGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGAt                                                                                                                                                                                                                                  >  1:340027/1‑251 (MQ=255)       aaGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTg                                                                                                                                                                                                                            <  1:208830/251‑1 (MQ=255)                     cTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGc                                                                                                                                                                                                              >  1:380799/1‑251 (MQ=255)                     cTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCAGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGc                                                                                                                                                                                                              >  1:176596/1‑251 (MQ=255)                       gggCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCg                                                                                                                                                                                                            >  1:220236/1‑251 (MQ=255)                           gTGCTGGGCGACTGGTCTCACCCGTACCTAACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCAACGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAg                                                                                                                                                                                                        <  2:237830/251‑1 (MQ=255)                                               cccGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCGGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAAtct                                                                                                                                                                                    <  1:45910/251‑1 (MQ=255)                                                   tACCTGACCATGGACTACAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAACGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGGATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGTCAGTCGATCAGGATGCCCTGCAAGCAAAATTTGCCGTAAGCACCGTTAAAGGCACAATCTCGCt                                                                                                                                                                                >  2:475199/1‑251 (MQ=255)                                                    ccctaaCCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCTAAATCATCGGCAACGGTCACCTGCACAAGGGGGCGAAGCC‑GTTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTg                                                                                                                                                                               <  2:217414/246‑1 (MQ=255)                                                    aCCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGGCGATCAGGTTGCACTGAAAGCAAAATTTGCCGTTAGCACCGTTAACGGCCCAATCTCGCTg                                                                                                                                                                               >  2:231411/1‑251 (MQ=255)                                                    aCCTAACCATGGACTCCAAAACTGAGGCAAACATCACCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTg                                                                                                                                                                               <  2:292514/251‑1 (MQ=255)                                                         aCCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGATGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGAGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGAAGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCCAAATTAGCCGTAAGCAACGTTAACGGCCCCATCTCGCTGGtgat                                                                                                                                                                          >  2:495622/1‑248 (MQ=255)                                                              ggACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAAGCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGa                                                                                                                                                                     <  1:61258/251‑1 (MQ=255)                                                                 cTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCAGCACCGGACACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCAAATCTCGCTGGTGATCTGGacc                                                                                                                                                                   >  2:97422/1‑250 (MQ=255)                                                                   tCAAAGCTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGaccac                                                                                                                                                                 <  1:518852/250‑1 (MQ=255)                                                                     aaaaCTGAAGCCAACATCATCCGCGCGCTGGGCATAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGTCTACCGTTCCGCGCTGGCGGAAGCGGAAGTTGAGTATTACGAAAAAACTTCTCCGTCCATCGCAGTTGCTTTTCAGGCGGTCGATCCGGATGCCCTGAAAACCAAATTTTCCATAAGCCAAGTTAACGGCCTCATCTCGCTgatggtctggcacactac                                                                                                                                                              >  2:116373/1‑234 (MQ=255)                                                                      aaaCTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGCCACCTGCACAAAGGCGCGAAGACGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGACGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTGTCAGGCGGTCGATCCGAATTCACTGAAAGCAAAATTTGCCGTAAGCACCGTTAACGGCCCAATCTCCCTGGTGATATGGACCACTACg                                                                                                                                                             >  2:617388/1‑251 (MQ=255)                                                                         cTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCCACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGGTGAGTATTACGACAACACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATATGGacaactaagccg                                                                                                                                                          >  2:242566/1‑241 (MQ=255)                                                                           gAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGCCGTTGCTTTTCAGGCGGTCGATCAGGATGc                                                                                                                                                                                                                                >  2:342838/1‑179 (MQ=255)                                                                           gAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGc                                                                                                                                                                                                                                <  1:342838/179‑1 (MQ=255)                                                                               ccAACATCATCCGCGCGCTGGGCAAAATAATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTATGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGATGATCTGGCCCACTACGCCGTGGAct                                                                                                                                                    >  1:125952/1‑251 (MQ=255)                                                                                    aacagccgctgtcTGGGCAAAATCATCGGCCACGGTCACCGTCACAAAGGCCCGAAGCCGGCA‑CACTGGTGCGTGGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAACCTGCTCCGGCCATCGACGTTGCTTTTCAGGCGGTCGAACAGGATGCACTGAAAGCAAAATTTGCCGTAAGCACCGTTAACGGCCCCATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTctgcc                                                                                                                                               <  1:29330/239‑1 (MQ=255)                                                                                      catcCGCGCGCGGCGCAAAATCATCGGCAACGGTCCCCTGCACAAAGGCGCGAAGCCAGTT‑CCCTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCGAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTctgcctg                                                                                                                                             <  2:280087/251‑1 (MQ=255)                                                                                          cgcgcgCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTT‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCa                                                                                                                                          >  1:295310/1‑250 (MQ=255)                                                                                             gcgcTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACcg                                                                                                                                      <  1:344042/251‑1 (MQ=255)                                                                                                      acaaTCATCGGCAACGGTCACCTGCACAAAGGCGCGACGCCGGAA‑CACTGGTGCGTTCACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACCCAACTTCACCGTCCATAGACGTTGCTTTTCAGGCGGTCGATCAGGATGTACTGAAAGCAAAACTTGCCGTAAGCAACGTTAACGGCCCACTCTCGCTGTCGATATGGACTACCACACCCTGGACTCTGCCTGCGCACCGCGCCAtctc                                                                                                                             >  2:287759/3‑251 (MQ=255)                                                                                                       aaaTCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGAAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATctct                                                                                                                            >  1:488134/1‑251 (MQ=255)                                                                                                        aaTCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTa                                                                                                                           <  1:97422/251‑1 (MQ=255)                                                                                                                        cACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCg                                                                                                            >  1:487567/1‑250 (MQ=255)                                                                                                                         aCCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGAc                                                                                                          <  1:53479/251‑1 (MQ=255)                                                                                                                          ccTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCACTAGCGGACGCGGAAGTTGAGTCTTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCCACCCGCGCCATCTCTATTGCACCTGCTTTCGACt                                                                                                         >  1:421319/1‑251 (MQ=255)                                                                                                                             gCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAACATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATg                                                                                                      >  1:188603/1‑251 (MQ=255)                                                                                                                             gCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATg                                                                                                      >  1:138441/1‑251 (MQ=255)                                                                                                                                caAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGATGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATgcgg                                                                                                   >  2:233262/1‑250 (MQ=255)                                                                                                                                  aaGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTg                                                                                                 >  1:612801/1‑251 (MQ=255)                                                                                                                                         gAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGCCGTTGAGTATTACGACCAAACTTCTCCGTCCATCGACGTTGCTTTTCCGGCGGTCGATCAGGATGCACTGCAAGCAAAATTTGCCGTAAGCAACGTTACCGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGAAATCTCTATTGCACCTGATTTCGACTAAGCGCTGGTGccgg                                                                                          >  2:283476/1‑248 (MQ=255)                                                                                                                                         gAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCACCGTTAACGGCCCCATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGc                                                                                             >  1:400745/1‑248 (MQ=255)                                                                                                                                          aaGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGAt                                                                                         <  1:600369/251‑1 (MQ=255)                                                                                                                                               cggtacACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGCGCTGATCGAc                                                                                     >  2:616229/6‑251 (MQ=255)                                                                                                                                                    cACTGGGGCGTTGCCGGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGAACAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTc                                                                                 <  2:295310/250‑1 (MQ=255)                                                                                                                                                         gggcgtggccTGCCGTGCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTGCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCg                                                                           <  2:3594/242‑1 (MQ=255)                                                                                                                                                            cGTTGACTGCCGCTCTGCGCCGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGGGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGa                                                                        <  1:193557/251‑1 (MQ=255)                                                                                                                                                                gtcTGCCGATCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACGTCTCCGACCAGCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTGCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTc                                                                     <  1:5346/248‑1 (MQ=255)                                                                                                                                                                  cTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTgg                                                                  >  1:449555/1‑251 (MQ=255)                                                                                                                                                                       gTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGATAACGGCCCAATCTCGCTGGTGATCTGGACCACTACCCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCCCCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTTTggtcaa                                                              >  2:333784/1‑246 (MQ=255)                                                                                                                                                                        ttCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGCACCGCGCAATCTCTATTGCCCCTGCTTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGaa                                                              >  1:408925/1‑249 (MQ=255)                                                                                                                                                                               cTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGt                                                     >  1:331542/1‑251 (MQ=255)                                                                                                                                                                                    ggaagcggaagTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAAGTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGaaa                                                <  1:475199/251‑1 (MQ=255)                                                                                                                                                                                      aagcggaagTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGc                                              >  1:104273/1‑251 (MQ=255)                                                                                                                                                                                          ggaagTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGCACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGGGCCGATCGACGGTCAGGCCGGGATTCTGGCGAAAGATCTg                                                       >  2:283003/1‑238 (MQ=255)                                                                                                                                                                                          ggaagTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTg                                                       <  1:283003/238‑1 (MQ=255)                                                                                                                                                                                              gTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGc                                       >  1:438627/1‑250 (MQ=255)                                                                                                                                                                                                tGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCCATCTCGCTGGGGATCTGGACCAATACGCCGTGGACTCTGCCTGCGAACCGCGCAATctct                                                                                                                            >  2:539503/1‑163 (MQ=255)                                                                                                                                                                                                tGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATctct                                                                                                                            <  1:539503/163‑1 (MQ=255)                                                                                                                                                                                                tGAGTATTACGAAAAAACTTCCCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTGCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAg                                     <  2:676222/250‑1 (MQ=255)                                                                                                                                                                                                   gTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGt                                  >  1:79123/1‑250 (MQ=255)                                                                                                                                                                                                      ttACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCg                              >  1:166640/1‑251 (MQ=255)                                                                                                                                                                                                          gACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGt                          >  1:253931/1‑251 (MQ=255)                                                                                                                                                                                                           aCAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGAGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCg                           >  1:430606/1‑249 (MQ=255)                                                                                                                                                                                                              aaaCTTCGCCGTCCACCGACTTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGTTGCAGAGCGACGGTCAGGCCGTGAGTCTGGCGAAAGAGCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAcc                      <  2:612801/251‑1 (MQ=255)                                                                                                                                                                                                                     tcCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTac                >  1:288808/1‑250 (MQ=255)                                                                                                                                                                                                                      ccGTCCATCGCCGTTGCTTTTCCGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAAGCTCTCTGGTGATCTGGACCACTACGCCGTGGACTCGGCCTGCGAACCGCGCAATCTCTATTGCACCTTATTTCGACTATGCGCTGGTGCAGATCGACGGGCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTacac              <  2:381234/251‑1 (MQ=255)                                                                                                                                                                                                                      ccGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTa                 <  1:242566/248‑1 (MQ=255)                                                                                                                                                                                                                       cGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACAcc             >  1:166281/1‑251 (MQ=255)                                                                                                                                                                                                                       cGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACAcc             >  1:482421/1‑251 (MQ=255)                                                                                                                                                                                                                        gTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGCTTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGCTCTGGTTTAAAGCGGAATGCAGCGTATCGGCgaggaccgttaaaccc            >  2:196756/1‑236 (MQ=255)                                                                                                                                                                                                                        gTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGTAGCGTATCGGCGTGAACGATTaaccc             >  2:541079/1‑246 (MQ=255)                                                                                                                                                                                                                          ccATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCAtt          <  1:236817/251‑1 (MQ=255)                                                                                                                                                                                                                           cATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGTCCGTGATTCTGGCGAAAGATCGGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtc         <  1:50283/251‑1 (MQ=255)                                                                                                                                                                                                                             ccGACGTTGCTTTTCAGGCCGTAGATCCGGAATCACTGAAATCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTTATCTGGACCACTACGCCGTTGACTCTGCCTGCGAACCGTGCAATCTCTATTGCACCTGATTTCGTCTATGCGCTGGTGCAGATCGACGTTCAGGTTGTGATTCTGGCGAAATATCTGGTTGAAAGCGTAATGCATCGTATCGGCTTGTCCGATTACTCCATtctc       <  2:468469/250‑1 (MQ=255)                                                                                                                                                                                                                                aCGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtctcagc    <  1:495622/251‑4 (MQ=255)                                                                                                                                                                                                                                 cGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtctcagc    >  1:59892/1‑247 (MQ=255)                                                                                                                                                                                                                                  gTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCAc  >  1:284427/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             AACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTT‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCAC  >  NC_000913/22800‑23274

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