Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,822	T→C	D363D (GAT→GAC)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,822	0	T	C	98.3%	185.0 / ‑4.7	59	D363D (GAT→GAC)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (24/33);  total (24/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.34e-01

TTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGC  >  NC_000913/4576‑5062                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   tttCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGt                                                                                                                                                                                                                                                >  1:504130/1‑251 (MQ=255)       cctcgacgGTCAGGGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATc                                                                                                                                                                                                                                          <  1:560991/248‑1 (MQ=255)          cgacgGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAg                                                                                                                                                                                                                                       <  1:374398/251‑1 (MQ=255)                      ggcACCAAAAGCGACTAAGGGTACGTAATCCAACGCGAGGG‑AGGTGGGTCCGCCAAACAACTGGCCGTGTGTGGAATAGGTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAAGGCGTGAGTTAAAAGAACTGGGCTCCACCTCGGAGCCGCCCGCTGCCGTAGCGTAGCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑Ct                                                                                                                                                                                                                          <  2:61563/249‑1 (MQ=255)                            cAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTcc                                                                                                                                                                                                                    <  1:364255/251‑1 (MQ=255)                               aGCGCCTCAGGCGACGTTAGCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGGTGAAGCCAGGCGAATATGG‑CTTGTTCCTCg                                                                                                                                                                                                                 <  1:629754/251‑1 (MQ=255)                                gCGACTCAGGCGATGTTATCAAACGC‑ATGGGCCGGGAGTCATCTGACCACCTGTCCGCGTGGGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAATCCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACCCCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATTG‑CTTGTTCCTCgg                                                                                                                                                                                                                <  2:658559/251‑1 (MQ=255)                                           gACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGACACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGAGTATCGTGCGCTGCGTGACCAGTTGATTCCAGGCGAATATGG‑CTTGTTCCTCGGCCCCGCACATc                                                                                                                                                                                                     >  2:44062/1‑251 (MQ=255)                                              gttctCCAACGCGATGG‑ATGTGAGCCAGCCGAACACCTGGCCGCGGGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCgg                                                                                                                                                                                                  <  2:393817/247‑1 (MQ=255)                                              gTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCgg                                                                                                                                                                                                  <  1:14968/251‑1 (MQ=255)                                                 aTCCACCGCGAGGG‑ACGTGGGTCAGCCGAACAACTGGCAGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCg                                                                                                                                                                                                <  1:467124/250‑1 (MQ=255)                                                         gcgtgggacgtgttgcAGCCGAACAATGGGCCGGTTGTGGAAGAGTTGGGCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGAGGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTTTCGTGCTCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTaa                                                                                                                                                                                        <  2:587991/236‑1 (MQ=255)                                                                        tggccgaaccaggggCCGTGTGTGGAAGAGTTGTCCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAg                                                                                                                                                                          <  2:368343/237‑1 (MQ=255)                                                                        cAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTGCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAAGGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGc                                                                                                                                                                         <  1:44062/251‑1 (MQ=255)                                                                            cGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCACGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCCCGGGCCCGCGCATCCGGCGACATTTAAAGAAGGCGTggtaacgagt                                                                                                                                                                     >  2:577042/1‑243 (MQ=255)                                                                               acaacTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAACTCTGGCAACTGCAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGCCAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCACATCCGGCGAAATTTTAAGAGAGCGTGGTAGCGATtctc                                                                                                                                                                  >  2:80458/1‑251 (MQ=255)                                                                                cagcggccggggtggggAAGATTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCg                                                                                                                                                                 <  2:393309/236‑1 (MQ=255)                                                                                     ggCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTg                                                                                                                                                              <  1:18396/249‑1 (MQ=255)                                                                                     ggCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGa                                                                                                                                                             >  1:476729/1‑250 (MQ=255)                                                                                       ccGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAAc                                                                                                                                                          >  1:48505/1‑251 (MQ=255)                                                                                          ggtgtgtaggagTTGCTCAGCCGCAAAATCTGGCCACTGAAAGAGCTGGGTTAGGCAGCCGTGGATGATGATACCACGCAACAGACAATGCGTGAGGTAAAAGAACTGGGCGACACCTCGGAGCCGCACGCTGCCGTAGCGTAGCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACggt                                                                                                                                                       <  2:232509/242‑3 (MQ=255)                                                                                          ggtgtGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGGGTGAGTTAAAAGAACTGGGCTACACTGCGGATCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGAGCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCGCGGCGCCGCGCAGCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGtt                                                                                                                                                       <  1:45971/250‑1 (MQ=255)                                                                                           gtgttgaaggttgtgtCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTAGCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTg                                                                                                                                                      <  2:319905/236‑1 (MQ=255)                                                                                                        ttcccgccgCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCaaa‑a                                                                                                                                         <  2:48505/247‑1 (MQ=255)                                                                                                          tCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTTGATCTGCCAAA‑aga                                                                                                                                       <  1:213553/250‑1 (MQ=255)                                                                                                              ccgcAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTgg                                                                                                                                  <  1:593325/251‑1 (MQ=255)                                                                                                               cgcAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGc                                                                                                                                 >  1:591809/1‑251 (MQ=255)                                                                                                                         ggCAACTGAAATAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACCGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGt                                                                                                                          <  1:370668/248‑1 (MQ=255)                                                                                                                         ggCAACTGAAAGAGCTGGTTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGc                                                                                                                        <  1:503454/250‑1 (MQ=255)                                                                                                                           cAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTg                                                                                                                      >  1:162144/1‑250 (MQ=255)                                                                                                                              cgtgaaGAGCTGGGTTATTCAGCCGTGGATTATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCCCGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTTCCAAA‑AGAGCTGGCAGAACGTGCTGAttt                                                                                                                  <  2:104366/247‑1 (MQ=255)                                                                                                                                    gagCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACACtt                                                                                                            >  2:185040/1‑251 (MQ=255)                                                                                                                                          ggTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGGTTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTc                                                                                                      >  2:524145/1‑251 (MQ=255)                                                                                                                                             tATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTcac                                                                                                    <  1:228428/250‑1 (MQ=255)                                                                                                                                                    ccGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTg                                                                                            >  1:453980/1‑251 (MQ=255)                                                                                                                                                    ccGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑ATTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGAGAAAAGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTg                                                                                            >  1:158657/1‑251 (MQ=255)                                                                                                                                                     cGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATcagc                                                                                           >  2:303781/1‑248 (MQ=255)                                                                                                                                                     cGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGc                                                                                           <  1:174490/251‑1 (MQ=255)                                                                                                                                                        ggATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCccg                                                                                        <  1:274851/251‑1 (MQ=255)                                                                                                                                                                    cccGCAACAGACAATGCGTGAGTTAAAAGAACTGTGCTACACCTCGGAGCCGCACGATGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTTTCGGTGAAACGTGGGATCTGCCAAA‑AGGCCTGCCTGAACGGCCTGATTTACCCTGGTTTTCACATAATCTTCCCGCCGTTTTtgtgg                                                                            >  2:54891/3‑248 (MQ=255)                                                                                                                                                                      cGCAACAGACAGTGCGTGTGTTAAAAGAACTGGGCTTCACCTCGGCGCCGCACGCTGCCGTGGGGTATCGTGCGCTGCGTGACGATTGGAACCCAGGAGAATTTGG‑CTTGTTCCCCGGCACCGCGCACCCGGCGAAATTTAAAGAGTGCGTGGAAGCGTTTCTCGGTTAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGACCGTGCGGATTGCCCCTTGCTTTCAATTAATCTGCCCGCCGATTTTGCtgcg                                                                          <  2:658717/251‑1 (MQ=255)                                                                                                                                                                         taCAGAAAATGCGTGTGTTAAAAGAACGGGGCTACACCTCGGGGCCGCTCGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGGACGGAGCAACCGGCGAAATTTAAAGAGAGCGTGGAAAGGCTTCTCGGTGAAACGTTGGGTATACCAAACAG‑GCTGGCAGAAAGGTCTGATTTAACCTTTCTGTCACATAATTTGCCCACCGAATTTGCtgtgttt                                                                       >  2:72285/2‑250 (MQ=255)                                                                                                                                                                          aCAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCg                                                                                                                                                                        <  1:417640/154‑1 (MQ=255)                                                                                                                                                                                aaTGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAAt                                                                >  1:188015/1‑251 (MQ=255)                                                                                                                                                                                     gTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGAGATTCTCGGTGAAACGTTGGCTCTTACAAA‑AGAGCTGGCAGAACGTGCTGTTTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATtgatg                                                           >  2:205781/1‑251 (MQ=255)                                                                                                                                                                                            aaaaGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGCTTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAatc                                                    >  1:240447/1‑251 (MQ=255)                                                                                                                                                                                                gAACTGGGCTACACCTCGGGGCCGCACGCTGCCGTAGCGTATCGGGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCCCCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAACCGTTGGATCTGCAAAA‑AGAGCTGGCAGAACGTGCTGATTTCCCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAAtcatca                                                <  1:339977/251‑1 (MQ=255)                                                                                                                                                                                                 accTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTCTCGTGCGCTGAGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTTAAGAGAGCGTGGAAGCGCTTCTCGGTGAAACGTTGGATCTGCCAAA‑AAAGCTGGAAGAAAGCTCTGATTTACCCTTTCTTTCACATAATCTGCCCGCCGGTTTTTCTCGGTTTGGTAAATTGATGATGAAtcatcac                                               >  2:158798/3‑250 (MQ=255)                                                                                                                                                                                                    ggggttACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGGGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTaa                                            <  1:164088/246‑1 (MQ=255)                                                                                                                                                                                                         tACACCTCGGAGCCGCACGCTGCCGTAGCGTAACGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGAATCTCGGTGAAACGTTGGATCTGCCAAC‑AGAGCTGGCAGAACGTGCTGCTTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTTGGTTGCGTAAATTGGTGATGAATCCTCAGTAAAATCt                                       >  2:509513/1‑251 (MQ=255)                                                                                                                                                                                                            ccctCGGAGCCGCACGCGGCCGGAGCGTATCGTGCGATGCGTGACCAGTTGAATCCAGGTGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACTTGGGTTCTGCCAAA‑AGAGCTGGCAGAACGTGTTGATTTACCCTTGCTTTCACAGAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTAtt                                    <  1:581461/248‑1 (MQ=255)                                                                                                                                                                                                            acctCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTAtt                                    >  1:182973/4‑251 (MQ=255)                                                                                                                                                                                                                        cACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCt                                                                                                                                                 <  1:5011/131‑1 (MQ=255)                                                                                                                                                                                                                        cACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCt                                                                                                                                                 >  2:5011/1‑131 (MQ=255)                                                                                                                                                                                                                               ccGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTTCTGCGTTCCGTAAATTGATGATGAATCATCAGTAAAATATATTTATTATATCACTCATGCCg                 >  2:311030/1‑251 (MQ=255)                                                                                                                                                                                                                                  tggcgtagcgggcgcTGCGTGACCAGTTGAATCCAGGCAAAT‑TGGTCTTGTTCCTCGGCACCGCGGAGCCGGGGAAATTGAAAGAGAGCGTGGAAGCGATTCCCGGTGAAAGTTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTGGCTTTCAAATAACCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCGCAATCAGTCCggg               <  1:215517/239‑1 (MQ=255)                                                                                                                                                                                                                                      gtATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGc          <  1:185040/250‑1 (MQ=255)                                                                                                                                                                                                                                            ggcgcTGCGGGACCAGTTGAATCCAGGCGAATATGG‑TTTGTTCCTCGGACCCGCGCATCCGGCGAAATTTAAAGAGAGCGGGGAAGCGAGTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTAt    <  1:478084/250‑1 (MQ=255)                                                                                                                                                                                                                                              cgcTGCGTGACCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGGAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGc  >  1:185524/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   TTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGG‑CTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGC  >  NC_000913/4576‑5062

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