Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,804	T→G	A357A (GCT→GCG)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,804	0	T	G	95.2%	185.5 / ‑1.1	63	A357A (GCT→GCG)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base G (25/35);  total (25/38)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

ACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCG  >  NC_000913/4560‑5047                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     ccGATACCGTGCCACGTTTCCTGCCCCACGGTCAGTCGTCACCCAAATGGACTCAGGCGACTTAATCCAACGAGATGG‑ACGTAAGGCAGCCGAACAACTGGCCGCGTGTGAAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGACGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGt                                                                                                                                                                                                                                                <  1:268457/250‑1 (MQ=255)  cGATACCGTGCCACTTTTCCTGCACGACGGTCAGTGGCCACCCAAAGCGACTCAGGCGACGTTACCCAACGCGATGG‑ACGTGAGTCAGCCTAACAGCTGGCCGCTTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTg                                                                                                                                                                                                                                               <  1:129698/251‑1 (MQ=255)           gCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTg                                                                                                                                                                                                                                      >  1:369344/1‑251 (MQ=255)            ccactgtcccggCACGACGTTCGGTGGTCACCCAAAGCGACTCAGGCGGCGTTACCCACCGCGGGGG‑ACGTGATACAGCCGAACAACTGGCCGCTTGTGGAAGAGTTGTGCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTAGCGTGCGCTGCGTGa                                                                                                                                                                                                                                     <  2:428663/240‑1 (MQ=255)              ccGTTTCCTGCACGACGGTCAGTGGTCCCCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGacc                                                                                                                                                                                                                                   <  1:372779/250‑3 (MQ=255)                 tttCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGt                                                                                                                                                                                                                                >  1:504130/1‑251 (MQ=255)                 gttCCTGAACGATGTTCGGTGCTCACCAAAAGCGACTCAGGCGCCGTTATCAAAAGCGGTGG‑ACGGGAGTCAGCCGAACAACTGGCCGCGTGGGGAATAGTTGCTCCGCCGCAAAATCTGGCAACTGAAAGATCGGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTTTCGTGCGCGGCGTGacca                                                                                                                                                                                                                                  <  2:532088/248‑4 (MQ=255)                       cctcgacgGTCAGGGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATc                                                                                                                                                                                                                          <  1:560991/248‑1 (MQ=255)                          cgacgGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAg                                                                                                                                                                                                                       <  1:374398/251‑1 (MQ=255)                            acgGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCt                                                                                                                                                                                                                                           >  2:206187/1‑229 (MQ=255)                            acgGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCt                                                                                                                                                                                                                                           <  1:206187/229‑1 (MQ=255)                                      ggcACCAAAAGCGACTAAGGGTACGTAATCCAACGCGAGGG‑AGGTGGGTCCGCCAAACAACTGGCCGTGTGTGGAATAGGTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAAGGCGTGAGTTAAAAGAACTGGGCTCCACCTCGGAGCCGCCCGCTGCCGTAGCGTAGCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCt                                                                                                                                                                                                           <  2:61563/249‑1 (MQ=255)                                            cAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTcc                                                                                                                                                                                                     <  1:364255/251‑1 (MQ=255)                                               aGCGCCTCAGGCGACGTTAGCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGGTGAAGCCAGGCGAATATGGCTTGTTCCTCg                                                                                                                                                                                                  <  1:629754/251‑1 (MQ=255)                                                gCGACTCAGGCGATGTTATCAAACGC‑ATGGGCCGGGAGTCATCTGACCACCTGTCCGCGTGGGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAATCCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACCCCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATTGCTTGTTCCTCgg                                                                                                                                                                                                 <  2:658559/251‑1 (MQ=255)                                                           gACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGACACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGAGTATCGTGCGCTGCGTGACCAGTTGATTCCAGGCGAATATGGCTTGTTCCTCGGCCCCGCACATc                                                                                                                                                                                      >  2:44062/1‑251 (MQ=255)                                                              gttctCCAACGCGATGG‑ATGTGAGCCAGCCGAACACCTGGCCGCGGGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCgg                                                                                                                                                                                   <  2:393817/247‑1 (MQ=255)                                                              gTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCgg                                                                                                                                                                                   <  1:14968/251‑1 (MQ=255)                                                                 aTCCACCGCGAGGG‑ACGTGGGTCAGCCGAACAACTGGCAGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCg                                                                                                                                                                                 <  1:467124/250‑1 (MQ=255)                                                                         gcgtgggacgtgttgcAGCCGAACAATGGGCCGGTTGTGGAAGAGTTGGGCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGAGGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTTTCGTGCTCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTaa                                                                                                                                                                         <  2:587991/236‑1 (MQ=255)                                                                                        tggccgaaccaggggCCGTGTGTGGAAGAGTTGTCCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAg                                                                                                                                                           <  2:368343/237‑1 (MQ=255)                                                                                        cAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTGCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAAGGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGc                                                                                                                                                          <  1:44062/251‑1 (MQ=255)                                                                                            cGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCACGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCCCGGGCCCGCGCATCCGGCGACATTTAAAGAAGGCGTggtaacgagt                                                                                                                                                      >  2:577042/1‑243 (MQ=255)                                                                                               acaacTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAACTCTGGCAACTGCAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGCCAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCACATCCGGCGAAATTTTAAGAGAGCGTGGTAGCGATtctc                                                                                                                                                   >  2:80458/1‑251 (MQ=255)                                                                                                cagcggccggggtggggAAGATTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCg                                                                                                                                                  <  2:393309/236‑1 (MQ=255)                                                                                                     ggCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTg                                                                                                                                               <  1:18396/249‑1 (MQ=255)                                                                                                     ggCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGa                                                                                                                                              >  1:476729/1‑250 (MQ=255)                                                                                                       ccGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAAc                                                                                                                                           >  1:48505/1‑251 (MQ=255)                                                                                                          ggtgtgtaggagTTGCTCAGCCGCAAAATCTGGCCACTGAAAGAGCTGGGTTAGGCAGCCGTGGATGATGATACCACGCAACAGACAATGCGTGAGGTAAAAGAACTGGGCGACACCTCGGAGCCGCACGCTGCCGTAGCGTAGCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACggt                                                                                                                                        <  2:232509/242‑3 (MQ=255)                                                                                                          ggtgtGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGGGTGAGTTAAAAGAACTGGGCTACACTGCGGATCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGAGCAGTTGAATCCAGGCGAATATGGCTTGTTCCGCGGCGCCGCGCAGCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGtt                                                                                                                                        <  1:45971/250‑1 (MQ=255)                                                                                                           gtgttgaaggttgtgtCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTAGCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTg                                                                                                                                       <  2:319905/236‑1 (MQ=255)                                                                                                                        ttcccgccgCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCaaa‑a                                                                                                                          <  2:48505/247‑1 (MQ=255)                                                                                                                          tCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTTGATCTGCCAAA‑aga                                                                                                                        <  1:213553/250‑1 (MQ=255)                                                                                                                              ccgcAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTgg                                                                                                                   <  1:593325/251‑1 (MQ=255)                                                                                                                               cgcAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGc                                                                                                                  >  1:591809/1‑251 (MQ=255)                                                                                                                                         ggCAACTGAAATAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACCGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGt                                                                                                           <  1:370668/248‑1 (MQ=255)                                                                                                                                         ggCAACTGAAAGAGCTGGTTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGc                                                                                                         <  1:503454/250‑1 (MQ=255)                                                                                                                                           cAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTg                                                                                                       >  1:162144/1‑250 (MQ=255)                                                                                                                                              cgtgaaGAGCTGGGTTATTCAGCCGTGGATTATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCCCGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTTCCAAA‑AGAGCTGGCAGAACGTGCTGAttt                                                                                                   <  2:104366/247‑1 (MQ=255)                                                                                                                                                    gagCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACACtt                                                                                             >  2:185040/1‑251 (MQ=255)                                                                                                                                                          ggTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGGTTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTc                                                                                       >  2:524145/1‑251 (MQ=255)                                                                                                                                                             tATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTcac                                                                                     <  1:228428/250‑1 (MQ=255)                                                                                                                                                                    ccGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTg                                                                             >  1:453980/1‑251 (MQ=255)                                                                                                                                                                    ccGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGATTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGAGAAAAGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTg                                                                             >  1:158657/1‑251 (MQ=255)                                                                                                                                                                     cGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATcagc                                                                            >  2:303781/1‑248 (MQ=255)                                                                                                                                                                     cGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGc                                                                            <  1:174490/251‑1 (MQ=255)                                                                                                                                                                        ggATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCccg                                                                         <  1:274851/251‑1 (MQ=255)                                                                                                                                                                                    cccGCAACAGACAATGCGTGAGTTAAAAGAACTGTGCTACACCTCGGAGCCGCACGATGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTTTCGGTGAAACGTGGGATCTGCCAAA‑AGGCCTGCCTGAACGGCCTGATTTACCCTGGTTTTCACATAATCTTCCCGCCGTTTTtgtgg                                                             >  2:54891/3‑248 (MQ=255)                                                                                                                                                                                      cGCAACAGACAGTGCGTGTGTTAAAAGAACTGGGCTTCACCTCGGCGCCGCACGCTGCCGTGGGGTATCGTGCGCTGCGTGACGATTGGAACCCAGGAGAATTTGGCTTGTTCCCCGGCACCGCGCACCCGGCGAAATTTAAAGAGTGCGTGGAAGCGTTTCTCGGTTAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGACCGTGCGGATTGCCCCTTGCTTTCAATTAATCTGCCCGCCGATTTTGCtgcg                                                           <  2:658717/251‑1 (MQ=255)                                                                                                                                                                                         taCAGAAAATGCGTGTGTTAAAAGAACGGGGCTACACCTCGGGGCCGCTCGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGGACGGAGCAACCGGCGAAATTTAAAGAGAGCGTGGAAAGGCTTCTCGGTGAAACGTTGGGTATACCAAACAG‑GCTGGCAGAAAGGTCTGATTTAACCTTTCTGTCACATAATTTGCCCACCGAATTTGCtgtgttt                                                        >  2:72285/2‑250 (MQ=255)                                                                                                                                                                                          aCAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCg                                                                                                                                                         <  1:417640/154‑1 (MQ=255)                                                                                                                                                                                                aaTGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAAt                                                 >  1:188015/1‑251 (MQ=255)                                                                                                                                                                                                     gTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGAGATTCTCGGTGAAACGTTGGCTCTTACAAA‑AGAGCTGGCAGAACGTGCTGTTTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATtgatg                                            >  2:205781/1‑251 (MQ=255)                                                                                                                                                                                                            aaaaGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGCTTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAatc                                     >  1:240447/1‑251 (MQ=255)                                                                                                                                                                                                                gAACTGGGCTACACCTCGGGGCCGCACGCTGCCGTAGCGTATCGGGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCCCCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAACCGTTGGATCTGCAAAA‑AGAGCTGGCAGAACGTGCTGATTTCCCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAAtcatca                                 <  1:339977/251‑1 (MQ=255)                                                                                                                                                                                                                 accTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTCTCGTGCGCTGAGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTTAAGAGAGCGTGGAAGCGCTTCTCGGTGAAACGTTGGATCTGCCAAA‑AAAGCTGGAAGAAAGCTCTGATTTACCCTTTCTTTCACATAATCTGCCCGCCGGTTTTTCTCGGTTTGGTAAATTGATGATGAAtcatcac                                >  2:158798/3‑250 (MQ=255)                                                                                                                                                                                                                    ggggttACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGGGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTaa                             <  1:164088/246‑1 (MQ=255)                                                                                                                                                                                                                         tACACCTCGGAGCCGCACGCTGCCGTAGCGTAACGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGAATCTCGGTGAAACGTTGGATCTGCCAAC‑AGAGCTGGCAGAACGTGCTGCTTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTTGGTTGCGTAAATTGGTGATGAATCCTCAGTAAAATCt                        >  2:509513/1‑251 (MQ=255)                                                                                                                                                                                                                            ccctCGGAGCCGCACGCGGCCGGAGCGTATCGTGCGATGCGTGACCAGTTGAATCCAGGTGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACTTGGGTTCTGCCAAA‑AGAGCTGGCAGAACGTGTTGATTTACCCTTGCTTTCACAGAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTAtt                     <  1:581461/248‑1 (MQ=255)                                                                                                                                                                                                                            acctCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTAtt                     >  1:182973/4‑251 (MQ=255)                                                                                                                                                                                                                                        cACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCt                                                                                                                                  >  2:5011/1‑131 (MQ=255)                                                                                                                                                                                                                                        cACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCt                                                                                                                                  <  1:5011/131‑1 (MQ=255)                                                                                                                                                                                                                                               ccGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTTCTGCGTTCCGTAAATTGATGATGAATCATCAGTAAAATATATTTATTATATCACTCATGCCg  >  2:311030/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     ACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCG  >  NC_000913/4560‑5047

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