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	95.6%	207.0 / ‑0.3	68	D363D (GAT→GAC)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (26/39);  total (28/40)
Fisher's exact test for biased strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

CCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGC  >  NC_000913/4579‑5054                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          ccTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTg                                                                                                                                                                                                                                     >  1:273863/1‑250 (MQ=255)            gTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTCTCCAACGCGATGGACGGGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGaa                                                                                                                                                                                                                         <  1:398352/251‑1 (MQ=255)            gTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGTTGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGa                                                                                                                                                                                                                          >  2:116093/1‑250 (MQ=255)                 tGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGCAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATgg                                                                                                                                                                                                                    >  2:418795/1‑251 (MQ=255)                     cACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACACCTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCtt                                                                                                                                                                                                                 <  1:34561/250‑1 (MQ=255)                      ccccAAAGCGACTCAGGCGTCGTTATCCAACGCGATGGCCGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCCCGCTGCCGTAGCGTATCGTGCGTTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGt                                                                                                                                                                                                               <  1:297593/250‑1 (MQ=255)                        ccAAAGCGACTCAGGCGACGTTATCCACCGCGATGGACGTGAGCCGGCCGAACAACTGTCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTc                                                                                                                                                                                                             <  1:414253/251‑1 (MQ=255)                           aaGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGGGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTc                                                                                                                                                                                                          <  1:463741/251‑1 (MQ=255)                            ggtgtctagggCGACGTTACACAACGTGAGGGACGGAAGTCAGCCGAACAACTGGCCGCGGGGGGAAGAGGTGTACCGCCGCAAAAACTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTTCACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGAGAATATGGCTTGTTCCTCg                                                                                                                                                                                                         <  2:284611/242‑1 (MQ=255)                              cGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGGGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGc                                                                                                                                                                                                       >  2:368653/1‑251 (MQ=255)                                ccTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGGGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCAc                                                                                                                                                                                                     <  1:116093/250‑1 (MQ=255)                                  tCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACcg                                                                                                                                                                                                   >  1:272649/1‑251 (MQ=255)                                       cGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCAt                                                                                                                                                                                              >  1:293130/1‑251 (MQ=255)                                        gACGTTACCCAAGGCGTGGGAGGTGAGTCAGCCGAACAACTGGCCGGGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATc                                                                                                                                                                                             <  2:272649/251‑1 (MQ=255)                                               tCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGa                                                                                                                                                                                       >  1:312740/1‑250 (MQ=255)                                                      cgATGGACGTGCGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTTTTCCCCGGCCCCGCCCATCCGGCGAAATTTaa                                                                                                                                                                                >  2:57865/1‑250 (MQ=255)                                                       gATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGt                                                                                                                                                                                                               >  1:300640/1‑218 (MQ=255)                                                          ggACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCAGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGACAttttagg                                                                                                                                                                              >  1:479725/1‑244 (MQ=255)                                                           gACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAagagag                                                                                                                                                                          <  1:478783/251‑1 (MQ=255)                                                             cGTGAGTCAGCCGAAAAGCTGGCCGCGTGGGGAAGAGTGATGCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCg                                                                                                                                                                        <  2:127693/251‑1 (MQ=255)                                                                    cAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAg                                                                                                                                                                  <  1:172918/250‑1 (MQ=255)                                                                         gAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATtc                                                                                                                                                            <  1:394235/251‑1 (MQ=255)                                                                           acaacTGGCCGCGTGTGGAAAAGTGGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATtctc                                                                                                                                                          <  2:293130/251‑1 (MQ=255)                                                                             agcgggCCGGGTGTGGAATGGTTGTTCCGCCGCAAAAGCTGGCAACTGAAAGGGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCGCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCgg                                                                                                                                                        <  2:241495/247‑1 (MQ=255)                                                                             agagggtggggtgtGGAAGAGTTGTCCGGCCGCAAAATCGGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCgg                                                                                                                                                        <  2:229784/241‑1 (MQ=255)                                                                             aacTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCgg                                                                                                                                                        <  1:416985/251‑1 (MQ=255)                                                                                gggCCGCTTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGa                                                                                                                                                     <  1:234080/250‑1 (MQ=255)                                                                                 ggCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGa                                                                                                                                                     <  1:80541/250‑1 (MQ=255)                                                                                      ggtgtGGACGAGGTGTCCCGCCGCAAAATCTGGCAACTGAAAGAGCTTGGTTACGCAGCCGTGGATGATGACACCACGCACCAGACAATGCGTGAGTTAACAGCCCTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTCAACCCAGGCGAATATGGCTTGTTCCTCGGCACCGCTCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGtt                                                                                                                                               <  2:107271/250‑1 (MQ=255)                                                                                         gtgGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTgg                                                                                                                                             <  2:216008/250‑1 (MQ=255)                                                                                          gggaggagTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAt                                                                                                                                           <  2:72559/246‑1 (MQ=255)                                                                                           ggAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAt                                                                                                                                           <  1:73750/250‑1 (MQ=255)                                                                                            gAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTTCTCGGCACCGCGCCTCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTggtact                                                                                                                                         >  2:95219/1‑247 (MQ=255)                                                                                             aaGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCCTCCGGCGAAATTTAAAGAAAGCGTGGAAGCGCTTCTCGGTGAAACGTTGGATCTg                                                                                                                                        >  1:156937/1‑251 (MQ=255)                                                                                             aaGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTg                                                                                                                                        <  1:451975/251‑1 (MQ=255)                                                                                                               aaaTCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGGGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGaa                                                                                                                      <  2:273863/251‑1 (MQ=255)                                                                                                                     ggCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGc                                                                                                                 <  1:474321/250‑1 (MQ=255)                                                                                                                        aaCTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑At                                                                                                            >  1:58289/1‑251 (MQ=255)                                                                                                                               agagCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCt                                                                                                     <  2:312740/251‑1 (MQ=255)                                                                                                                                gagCTGGGTTAGGCAAGCGTGGATGATGAAACCACGCAAAAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCtt                                                                                                    <  2:478511/251‑1 (MQ=255)                                                                                                                                       gctATGCAGCCGGGGATGATGAAACCACGCAACAGACAATGCGTGCGTTAAGAGTACGGGGCTACACCTCGGGCCCGCACGCTGCCGTAGCGTATCGGGGGCTGCGTGACCAATTGACTCCAGCCGAGTATAGCTTTTTTCAGGGCACCGCCCAAAGCGCGACATTTTAAGAGAGCGTGGAGCCCGTTCTCGGGGCAAAGTTTGCACGTCCAAAAGAGATGGCAGAACGAGTCG‑ATTTACCCTTGTTTTc                                                                                              >  2:35770/3‑250 (MQ=255)                                                                                                                                       gTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTca                                                                                             <  2:442370/251‑1 (MQ=255)                                                                                                                                         tATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTcac                                                                                            <  1:95219/250‑1 (MQ=255)                                                                                                                                                ccGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTg                                                                                    <  1:256979/251‑1 (MQ=255)                                                                                                                                                  gTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCCCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAACTTGAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGcc                                                                                  >  1:156666/1‑251 (MQ=255)                                                                                                                                                    ggATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGTCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCccg                                                                                <  2:25999/251‑1 (MQ=255)                                                                                                                                                         atgaAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGAt                                                                           >  1:362427/1‑251 (MQ=255)                                                                                                                                                             aaCCACGCAACAGACAATGCGTGAGTTCAAAGAACTGGGCTACACCTCGGAACCGCACGCTGCCGTAGCGTATCGAGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCCCGGCACCGCGCCTCCGGCGAAATTTAAAGAGGGCGTGGAAGCGATTCTCGGAGAAACGGTGGATATGCCAAAAGAGGTGGCAGAACGAGATG‑ATTTTCCCTTGTTTTCACCTATACTGCCCCCCGAttttt                                                                       >  2:294296/1‑250 (MQ=255)                                                                                                                                                                  cgcaacagacaagggggtatttAAAAGAACTGGTCTACCCCCCGGAGCCGCACGCGGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATAGGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACTTTGGATCTGCAAAAAGAGCTGGCAGATCGTGCTG‑ATTTCCCCTTGCTTTAACATAATCTGCCCGCCGATTTTGCtgcg                                                                  <  2:362427/231‑1 (MQ=255)                                                                                                                                                                    cAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATAGGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGGGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCtgcgtt                                                                <  1:418795/251‑1 (MQ=255)                                                                                                                                                                      aCAGACAATGAGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGAGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTTAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCtgcgttgc                                                              >  1:83719/1‑251 (MQ=255)                                                                                                                                                                             aggCGTGAGTTAAAAGAACGGGGCTACACCCCTGAGCCGCACGCGCCCGGAGCGTTTCGTGCGCTGCGTGACCATTTGAATCCGGGCGAATATGGTTTGTTCCCCGCCACCGCGCATCCTGCGAAATTTAAAGAGAGCGGGGAAGCAATTCTCGGTGAAACGTTTGATCTGCCAAAAGAGCTGGCAGAACGTG‑TGTATTTCCCCTTGCTTTCACATATCCTGCCCGCCGATTTTGCTGCGTTGCGTAAAtt                                                       <  1:60698/249‑1 (MQ=255)                                                                                                                                                                             aTGCGTGAGTTACAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCg                                                                                                                                                                                                         >  2:478645/1‑106 (MQ=255)                                                                                                                                                                             aTGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGGAGCGGATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCg                                                                                                                                                                                                         <  1:478645/106‑1 (MQ=255)                                                                                                                                                                                 gTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATtgatg                                                   >  1:320584/1‑251 (MQ=255)                                                                                                                                                                                         aaaGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGGGTATCGTGCTCTGCGTGACCAGTTGAATCCAGGAGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCCCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAatca                                           <  1:23503/251‑1 (MQ=255)                                                                                                                                                                                               cTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTACTCGGCACCGCGCATCCGGCGAAAATTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAAAGTTGGATCTGCCAAAAAAGCTGGCAGAACGGGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATAATCAGt                                      >  2:49592/1‑250 (MQ=255)                                                                                                                                                                                                       cacctCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGGTCTGCCAAAAGAGCTGGCAGAACGGGCAG‑ATTTACCCTTGCTCTCACATAATCTGCCCGCCGATTTTGTTGCGGTGCGTAAATTTATAGTGAATAATCAGTAAAATCTAt                             >  2:407594/5‑251 (MQ=255)                                                                                                                                                                                                             gcaGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCCGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCCCCGCGCATCCGGCGAAATTTAACGAGACCGTGGAAGCGATTCTCGGCGCAACGTTGGATCTGCCAAAAGAGCTGGCAGAACCCGCTG‑CTTTACCCCTGCCTTCCCATCATCCGCCCGCCCCCTTTGCTGCGTTGCCTACATTGCTGCTGaaccacccaccaaacctactccctc                       >  2:308322/3‑228 (MQ=255)                                                                                                                                                                                                              gAGCCGCACGCTGCCGTAGCGTAGCGTGCGCTGCGTGACCAGTTGAACCCAGGCGAATATGGCTTATTCCTCGGCCCCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTAt                      <  1:244796/251‑1 (MQ=255)                                                                                                                                                                                                                    cACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGGGAATATGGCTTGTTCCTCGGCCCCGCGCAGCCGGCGAAATTTAAAGAGAGCGGGGAAGCGTTTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTGACCCTTGCTTGCACAAAATCGGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAAt                <  1:57865/251‑1 (MQ=255)                                                                                                                                                                                                                            cGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGACACGTTGGATCTTCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCTTCAGTAAAATCTATTCATTATCTCAATCAGGCCgg        >  1:245834/1‑251 (MQ=255)                                                                                                                                                                                                                               agcgtATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGc           >  2:25561/5‑245 (MQ=255)                                                                                                                                                                                                                               agcgtATCGGGCGCTGCGTGACCAGTTGAATCCAGGCGAATAGGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGc           <  1:25561/241‑1 (MQ=255)                                                                                                                                                                                                                                gCTTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAACCGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTCCCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGttt    <  1:147435/251‑1 (MQ=255)                                                                                                                                                                                                                                 cgtATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGtttt   >  2:88303/3‑250 (MQ=255)                                                                                                                                                                                                                                  gtATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGc  >  2:72357/2‑251 (MQ=255)                                                                                                                                                                                                                                    aTCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGCGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTCTTCATTATCTCAATCAgg            >  1:46239/1‑239 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          CCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGC  >  NC_000913/4579‑5054

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