Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	22,945	T→A	V185V (GTT→GTA)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	22,945	0	T	A	92.2%	135.3 / 3.4	52	V185V (GTT→GTA)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base A (22/25);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

GAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTT‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAG  >  NC_000913/22702‑23162                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                           gAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGt                                                                                                                                                                                                                    <  1:235062/251‑1 (MQ=255)    aGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCg                                                                                                                                                                                                                 <  1:86760/251‑1 (MQ=255)                aTACGGTAAACCGGGAGAGAACTTCACCGCCGCCGAATTCCGCGCCAAGTACCGCGAATACGCGGCGACACAGGTTGACGGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAAATGAAGCCAACAACATCCACGCGCTGGGAAAAATCATCGGCACCGGTCACCTGCACAAAGCCGCGAAGCCGGTC‑CCCTGTTGCGTTGACTGCCGtt                                                                                                                                                                                                     >  2:380803/1‑251 (MQ=255)                  aCGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCt                                                                                                                                                                                                   >  1:270039/1‑251 (MQ=255)                     gTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTgcg                                                                                                                                                                                                <  1:287690/251‑1 (MQ=255)                       aaGCCGGGTGAGAAATTCACCGCCGCCGAGGTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACAATGGACTTCCAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCTTTCTGCGCt                                                                                                                                                                                              >  2:103580/1‑251 (MQ=255)                        aGCCGGGTGAGAAATTCACCGCCACCGAGTTCCGCACCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTAGGCGACTGGGCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGTCAACATCATCCGCGCGCGGGGCAACATCATCGGCAACGGACACCTGCACAAAGACGCGACGCCGGTA‑CACTGACGAGATGACTGCCGTACTGCGCTg                                                                                                                                                                                             >  2:160371/1‑251 (MQ=255)                         cccggtggagaACTTCCCCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CTCCCGTACCTGACCATGGACTTGAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTgg                                                                                                                                                                                            <  2:300123/244‑1 (MQ=255)                                gagaAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGgcggaag                                                                                                                                                                                      <  1:180126/250‑1 (MQ=255)                                   ataTTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCAGCGACCCAGGTTGACGGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACATCAAAACTGAAGCCAACATCATCAGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGACCAAAGGCGCGAAGCCGGTA‑CGCTGGTACGATGACTGCCGTTCTACGATGgcggaagt                                                                                                                                                                                     >  2:126546/3‑247 (MQ=255)                                   aaaTTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGgcggaagcgga                                                                                                                                                                                  <  1:84819/251‑1 (MQ=255)                                    actccACCCCCGCCGAGTTTCTCGACAAGTGCCGCGAATACGCGGGGACCCAGGTTGCCGGGCAACGCAAAG‑ACTTTATCCGGCTGGGCGTGCTGGGCGACTGGTCGACGCCGGTAC‑TGACCATGGACTTCAAAAATGAAGCCAACATCATCCGGGCGCTGGGCAAAATCAGCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGgcggaagcggaa                                                                                                                                                                                 <  2:377527/247‑1 (MQ=255)                                       ccccagccgcctAGTGCCGCGCCAAGAGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACACAAAG‑ACTTTATCCGTCTGGGCGTGATGTGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCAGCGGCAACGGTCCCCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGtt                                                                                                                                                                              <  2:94353/246‑1 (MQ=255)                                           cgccgccgAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAAAGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAgg                                                                                                                                                                          >  1:160808/1‑250 (MQ=255)                                            gccgccgAGTTCCGCGCCAAGTGCCGCGAATACGCGGTGACCCAGGTTGACGGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCAGTACCTGACCCTGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATAGGCAACGGTAACCTGCACAAAGGCGCGAAGCCGGGA‑CGATGGTGCCTTGACTGACGTTCTGCGCTGacggacgatgaagtggcgt                                                                                                                                                                          >  2:242473/1‑236 (MQ=255)                                            cccgccggtttCCGCGCCAAGTACCGCGAATACGCGGTGACCCAGGTTGACGGTCAACGCAAAG‑CCTTTATCCGTCTGGGCGTGCTGGGCGTCTGGTCG‑CACCCGTCCCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGTTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTa                                                                                                                                                                         <  2:127970/242‑1 (MQ=255)                                             ccgccgAGTTCCGGGCCAAGTGCCGCGAATATGCGGCGACCCAGGTTGACGGTCAACGCAAAG‑ACTTTATCCGGCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCAGCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGgcggaagcgttatcctttgat                                                                                                                                                                        <  2:321962/251‑13 (MQ=255)                                                          cgcCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACtt                                                                                                                                                           >  1:66356/1‑251 (MQ=255)                                                           gcCAAGTGCCGCGAATACGTGGCGACCCAGGTTTACGGTTAACGCAAAG‑ACTTTATCCGTCTGTGCGTGTTGGTCGACGTGTCG‑CACCCGTACCTGTCCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGAAAAATCATCGTCAATGGTCACCTGCACAAAGGCGCGAAGACGGTA‑CACTGTTGCGTTGACTGCAGTTCTGCGCTGGCTGAAGCGGAAGTTGAGTATTACGACAAAACTtc                                                                                                                                                          <  1:776/251‑1 (MQ=255)                                                            ccAAGTGCCGCGAATACGCGGCGACCCAGGTTGACTGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGCCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTT‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTtct                                                                                                                                                         >  1:80197/1‑251 (MQ=255)                                                                    cgcgAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCAt                                                                                                                                                 >  1:283388/1‑251 (MQ=255)                                                                        aaTACGCGGCGACCCAGGTTGACGGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTCCCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCCCCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCTTTTACTGCCGTTCTGCGCTTGCGGAAGCGGAAGTTGAGTAGTACGACAAAACTTCTCCGTCCATCGAc                                                                                                                                             <  1:337960/251‑1 (MQ=255)                                                                             gcggcgACCCAGGGTGACGGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGc                                                                                                                                        >  1:91679/1‑251 (MQ=255)                                                                               ggcgACCCAGGTTGACGGTCAACGCCAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGAAACGGACACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCtt                                                                                                                                      >  2:113525/1‑251 (MQ=255)                                                                                       aGGTTGACGGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTT‑CACTGGTGCGTTGACTGCCGTTCTGCTCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGc                                                                                                                               >  1:81508/1‑250 (MQ=255)                                                                                                  cAACGCAAAGTACTT‑ATCCGGCTGGGCTTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGAAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAgg                                                                                                                    <  2:145386/250‑1 (MQ=255)                                                                                                    aCGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGgcggaagcgg                                                                                                                                                                                   <  1:202786/185‑1 (MQ=255)                                                                                                        aaaG‑AATTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACt                                                                                                             >  1:300165/1‑251 (MQ=255)                                                                                                         aaagtccttATCCGTCTGGGCGTGCTGGGCGCCTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACt                                                                                                             <  2:66356/244‑1 (MQ=255)                                                                                                           g‑ACTTTATCCGTCTGGGCGTGCTGGGCGACGGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCATGCGGTCGATCAGGATGCACTGaa                                                                                                          <  1:41255/251‑1 (MQ=255)                                                                                                                     cGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAAtt                                                                                                 >  1:168669/1‑251 (MQ=255)                                                                                                                      gTCTGGGGGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACCGCAGCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGTACAAAGGCGCGAAGCCGGTA‑CCCTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTACATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACGGAAAGCAAAAtgt                                                                                                <  2:86392/251‑3 (MQ=255)                                                                                                                      gTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGAATTACGCCAAAACTTCTCCGTCCATCGAAGTTGCTTTTCAGGCGGTAGCTCAGGATGCACTGCAAGCAAAAttt                                                                                                >  2:319312/1‑251 (MQ=255)                                                                                                                                   gggCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACGGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCaa                                                                                     <  1:232271/249‑1 (MQ=255)                                                                                                                                     gCGACTGGTCG‑CACCCGTATCTTACCATGGACTTAAAAACTGAAGCCAACATCATACGCGCGCTTGGAAAAATCATCGGCAACGGTCACCGGCACAAAGGCTCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGTAGGCACTGAAAGCAAAATTTGCCGTAAGCAACggt                                                                                 <  2:91329/251‑3 (MQ=255)                                                                                                                                           ggTCT‑CACCCGTACCTTACCATGGACTTAAAACCTGAAGCCACCATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTTCGCTGGCGGAAGCGGAAGTGGAGTATTACGACAAAACTTCTCCGTACATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACg                                                                             <  2:308637/249‑1 (MQ=255)                                                                                                                                              cG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCGGGGCAAAAGCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTGGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTGCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGGAAGCAACGTTAACGGCCCa                                                                        <  1:39459/251‑1 (MQ=255)                                                                                                                                                           tGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGtg                                                            >  1:378096/1‑250 (MQ=255)                                                                                                                                                           tGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGtg                                                            >  2:73616/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)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                           GAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAG‑ACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCG‑CACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTT‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAG  >  NC_000913/22702‑23162

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