Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,780	T→C	T349T (ACT→ACC)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,780	0	T	C	95.4%	180.9 / ‑1.5	65	T349T (ACT→ACC)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base C (26/36);  total (26/39)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

GACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCT  >  NC_000913/4549‑5025                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       gACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACACAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGAGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCAGTGGATGATGAAACCACGCAACAGACCATGCGTGAGTTAAAAGAACTGGGCTTAACCTCGGAGCCGCACGCTGCCg                                                                                                                                                                                                                                     >  2:629754/1‑251 (MQ=255)   ccAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCCCCCAAAGCGCCTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCTCGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTa                                                                                                                                                                                                                                   <  1:74054/251‑1 (MQ=255)     aaCGAGAACGATACCGTGCCACGTTTCCTGAACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTAATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAgg                                                                                                                                                                                                                                 >  2:228428/1‑250 (MQ=255)            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)                                 tGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACg                                                                                                                                                                                                                                           <  1:100773/213‑1 (MQ=255)                                 tGCACGACGGACAGTGGTCACCAAAAGCGGCTCAGGCGACGTTATCCAACGCGATGG‑CCGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGCTGACACCACGCAACAGACAATGCGTGAGTTAAAAGAAATGGGCTACACCTCGGAGCCGCACg                                                                                                                                                                                                                                           >  2:100773/1‑213 (MQ=255)                                  cctcgacgGTCAGGGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATc                                                                                                                                                                                                    <  1:560991/248‑1 (MQ=255)                                     cgacgGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAg                                                                                                                                                                                                 <  1:374398/251‑1 (MQ=255)                                       acgGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCt                                                                                                                                                                                                                     <  1:206187/229‑1 (MQ=255)                                       acgGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCt                                                                                                                                                                                                                     >  2:206187/1‑229 (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)                                                          aGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTg                                                                                                                                                                                                                                        >  2:484051/1‑191 (MQ=255)                                                          aGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTg                                                                                                                                                                                                                                        <  1:484051/191‑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)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       GACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGG‑ACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAA‑AGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCT  >  NC_000913/4549‑5025

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