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	89.1%	132.5 / 7.2	55	A216A (GCA→GCG)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (1/5);  new base G (25/24);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01

CAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTT‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGT  >  NC_000913/22799‑23284                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       cAACGCAAAGTACTT‑ATCCGGCTGGGCTTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGAAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAgg                                                                                                                                                                                                                                              <  2:145386/250‑1 (MQ=255)       aaaG‑AATTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACt                                                                                                                                                                                                                                       >  1:300165/1‑251 (MQ=255)        aaagtccttATCCGTCTGGGCGTGCTGGGCGCCTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACt                                                                                                                                                                                                                                       <  2:66356/244‑1 (MQ=255)          g‑ACTTTATCCGTCTGGGCGTGCTGGGCGACGGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCATGCGGTCGATCAGGATGCACTGaa                                                                                                                                                                                                                                    <  1:41255/251‑1 (MQ=255)                    cGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAAtt                                                                                                                                                                                                                           >  1:168669/1‑251 (MQ=255)                     gTCTGGGGGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACCGCAGCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGTACAAAGGCGCGAAGCCGGTA‑CCCTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTACATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACGGAAAGCAAAAtgt                                                                                                                                                                                                                          <  2:86392/251‑3 (MQ=255)                     gTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGAATTACGCCAAAACTTCTCCGTCCATCGAAGTTGCTTTTCAGGCGGTAGCTCAGGATGCACTGCAAGCAAAAttt                                                                                                                                                                                                                          >  2:319312/1‑251 (MQ=255)                                  gggCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACGGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCaa                                                                                                                                                                                                               <  1:232271/249‑1 (MQ=255)                                    gCGACTGGTCGCACCCGTATCTTACCATGGACTTAAAAACTGAAGCCAACATCATACGCGCGCTTGGAAAAATCATCGGCAACGGTCACCGGCACAAAGGCTCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGTAGGCACTGAAAGCAAAATTTGCCGTAAGCAACggt                                                                                                                                                                                                           <  2:91329/251‑3 (MQ=255)                                          ggTCTCACCCGTACCTTACCATGGACTTAAAACCTGAAGCCACCATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTTCGCTGGCGGAAGCGGAAGTGGAGTATTACGACAAAACTTCTCCGTACATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACg                                                                                                                                                                                                       <  2:308637/249‑1 (MQ=255)                                             cGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCGGGGCAAAAGCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTGGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTGCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGGAAGCAACGTTAACGGCCCa                                                                                                                                                                                                  <  1:39459/251‑1 (MQ=255)                                                         tGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGtg                                                                                                                                                                                      >  2:73616/1‑250 (MQ=255)                                                         tGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGtg                                                                                                                                                                                      >  1:378096/1‑250 (MQ=255)                                                           aCCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGtgat                                                                                                                                                                                    >  1:310785/1‑248 (MQ=255)                                                                      cAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGaccact                                                                                                                                                                          <  1:113525/250‑2 (MQ=255)                                                                         aaCTGAAGCCCACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCAAATCTCGCTGGTGATCTGGACCACTACGc                                                                                                                                                                      >  1:94221/1‑251 (MQ=255)                                                                                  cAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGActc                                                                                                                                                             <  1:10627/251‑1 (MQ=255)                                                                                     catcatCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGCCTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTctgc                                                                                                                                                          >  1:367100/1‑251 (MQ=255)                                                                                        catcCGTGCGCAGGGAAAAATCATCGGCAACAGGCTCCTGCACAAGTGCGCGAAGCC‑GTTACACAGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTCTTACGCAAAAACTTCTCCGGCCATCGACGTTGCTTTTCAGGCGGTCGATTAGGATGCACTGAAAGCAAAATTTGTCGTAATCAACGTTAACGGCCCAATCTCGCTGGTTATCTGTACCACTACGCCGTGGACTctgcatg                                                                                                                                                       <  2:369210/251‑4 (MQ=255)                                                                                               gcgcTGGGCAAAAGCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACcg                                                                                                                                                <  1:319312/251‑1 (MQ=255)                                                                                                cgcTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTAGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACcgc                                                                                                                                               >  1:292481/1‑251 (MQ=255)                                                                                                   ttggCAAAATCATCGGCATCGGTCACCTGCACCAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTAGTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGACGATCGATCAGGCTGAACTGAAAGCAAAATTTGCCGTAAGCAACGCTACCGGCCCAATCACGCTGGTTATATGGAACACTACGCCGTGGACTCGCCCTACGAACCgcgc                                                                                                                                             >  2:288086/3‑250 (MQ=255)                                                                                                      gCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAAtc                                                                                                                                         >  1:84339/1‑251 (MQ=255)                                                                                                           atcatcGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTAt                                                                                                                                    <  1:380803/251‑1 (MQ=255)                                                                                                                  gCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTT‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGACTGCCAACCGAGCAATCTCTATTGCAca                                                                                                                              >  2:78572/1‑249 (MQ=255)                                                                                                                    aaCGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCAcctg                                                                                                                            >  1:295584/1‑248 (MQ=255)                                                                                                                                                      cACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCCCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCTCTGGTGCAGATCGACGGTCa                                                                                          <  1:176412/251‑1 (MQ=255)                                                                                                                                                      cACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTc                                                                                           <  1:381459/250‑1 (MQ=255)                                                                                                                                                       aCTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAg                                                                                         <  1:242473/251‑1 (MQ=255)                                                                                                                                                            tGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGt                                                                                    >  1:123931/1‑251 (MQ=255)                                                                                                                                                             gcggtGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTg                                                                                   <  1:28281/247‑1 (MQ=255)                                                                                                                                                              cGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGACCGTGa                                                                                  <  2:80197/251‑1 (MQ=255)                                                                                                                                                                ttGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCGGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTCTTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCCGGCCGTGAtt                                                                                >  1:318516/1‑251 (MQ=255)                                                                                                                                                                ttGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGAtt                                                                                >  1:283002/1‑251 (MQ=255)                                                                                                                                                                    ctgacGTTCTGCGCTGGCGGAAGCGGCAGTTGAGTATTACGACAAAACGTCTCCGTCCATCGACGGTGCTTTTCAGGCGGTCGATCCGGATGCACTGACAGCAACATTTGCCGTAAGCAACGTTAACGGCCCTATCTCGCTGGTGATCTGGACAACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTTTTGCACCATATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGcagctcttct                                                                              >  2:270170/5‑240 (MQ=255)                                                                                                                                                                      gCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTAAGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCACGCAATCTCTATTTCACCTGATTTAGACTATGCGCTGGTGCAGAACGACGGGCAGGCCGTGGTTCTGGCg                                                                          >  1:193842/1‑251 (MQ=255)                                                                                                                                                                              gcgcTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGACTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCACTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGa                                                                     <  1:78572/248‑1 (MQ=255)                                                                                                                                                                                ccTGGAGGAAGCGGAATTTGAGTTTTACGACAAAACTTCTCCTGCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTTGCGAAAGAGCTgg                                                                <  2:300165/250‑1 (MQ=255)                                                                                                                                                                                     cggaagcggcagTTGAGTATTACGACAAAACTTCTCCGTCCATAGACGTTGCTTTTCAGGCGGTCGATCAGGATGCAATGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGAACCTGATTTCGACTATGCGCTGGTGAAGATCGACGTCCAGGCCGTGATTCTGGCGAAAGATCTGGTTGa                                                            >  2:362690/1‑250 (MQ=255)                                                                                                                                                                                         agcggaagTTGAGTATTACGACAAAACTTCTCCGTGCATCTACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCCCGCCGTTGACTCTGCCTGCCAACCGCGCAATCTCTATTGCAGCAGATTTCGACTATGCGCTGGTGCAGAGCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCg                                                       <  2:81508/251‑1 (MQ=255)                                                                                                                                                                                            ggaagTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGTTTCTGGCGAAAGATCTGGTTGAAAGCGTa                                                     >  2:193650/1‑250 (MQ=255)                                                                                                                                                                                               agTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGc                                                 <  1:144319/251‑1 (MQ=255)                                                                                                                                                                                                gTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCa                                                <  1:113558/251‑1 (MQ=255)                                                                                                                                                                                                  tGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGAGGTCGATCAGGCTGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGAGATCTGGACAACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATATCTATTGCAACTGATTTCGACTATGCGCTGGTGCAGATCGCCGGACAGCCCGTGATTCTGGCGAGAGATCTGGTTGGAAGCGGAGTGCAGc                                              >  2:850/1‑251 (MQ=255)                                                                                                                                                                                                                aaaCTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAcc                                >  1:230837/1‑251 (MQ=255)                                                                                                                                                                                                                aaaCTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAcc                                <  1:70087/251‑1 (MQ=255)                                                                                                                                                                                                                         cGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACAcc                       <  1:73616/251‑1 (MQ=255)                                                                                                                                                                                                                         cGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACAcc                       >  1:137541/1‑251 (MQ=255)                                                                                                                                                                                                                                    gTTGATTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGAAACGTTAACGGCCCACTCTCGCTGGGGATCTGGACCACTACGACGTGGCCTCTGCCTGCGAACCTCGCAATCTCTATTGCACCTTATTTCGACTATGCGCTGGCGCAGAGCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTTAAAGAGTAATGCAGCTTCTCGGCGTAGCCGATTaccacaatttcgacac            >  2:25393/1‑237 (MQ=255)                                                                                                                                                                                                                                      tGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTg                                   <  1:293802/226‑1 (MQ=255)                                                                                                                                                                                                                                      tGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGagatcggaagagcgacgggtcggg           >  2:293802/1‑227 (MQ=255)                                                                                                                                                                                                                                          tttcAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTaa       <  1:211986/247‑1 (MQ=255)                                                                                                                                                                                                                                            tcAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTa                           >  1:207447/2‑228 (MQ=255)                                                                                                                                                                                                                                            tcAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCGCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTa                           <  2:207447/227‑1 (MQ=255)                                                                                                                                                                                                                                              aGGCGGTCGATCAGGGTGCACTGAAAGCAAAATTTGCCGTACGCAACGTTAACCGCCCAAGCTCGCTGGTGCTCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTGTGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTaaaatgt  <  1:346754/251‑4 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       CAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTT‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGT  >  NC_000913/22799‑23284

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