Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,572,319	T→G	E395A (GAA→GCA)	glgB ←	1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,572,319	0	T	G	97.9%	134.6 / ‑3.6	48	E395A (GAA→GCA)	glgB	1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (22/24);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-01

CATATCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAATAACGGAAGTCGTCGC  >  NC_000913/3572087‑3572559                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 cATATCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTCCAATCCAGTAAAGCGCGt                                                                                                                                                                                                                                >  1:219819/1‑251 (MQ=255)    atCCTGCGGACGAGAAACGCCAGTGTAATCGGTATACTCCTCAGCCATTGTCCCCGCACCGGAAACCGGCTCACCAAGAATACGGTTGTTATTACGCAAGCATTCAATCGCTTCAAGATTCGCGCCCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATCCCAAAACGTGCAATCCAGTAAAGCGCGTTAc                                                                                                                                                                                                                             <  1:326249/251‑1 (MQ=255)       cTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACcg                                                                                                                                                                                                                           >  1:142384/1‑250 (MQ=255)           gtaCGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCCTCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACcgacg                                                                                                                                                                                                                        >  1:284385/3‑249 (MQ=255)           ggACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAGACGTGCAATCCAGTAAAGCGCGTTACCGACGAg                                                                                                                                                                                                                      <  1:73778/251‑1 (MQ=255)                     cGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAACGCGCGTTACCGACGAGGAAGTTACTg                                                                                                                                                                                                            >  1:375717/1‑251 (MQ=255)                            ggaaTCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCACGCCAAATTCGGTCGGGATCCACTCCCCCTCTTTACGGCTGGAGTCGCGATAAATCATTAAAGCCACCGCAACGCCACGAAGCGCATCAATACCAAAACGTGCACTCCACAAACGCGCGTTACCGAACAAGAAGTTACTgcctgtcc                                                                                                                                                                                                     >  2:164059/3‑244 (MQ=255)                            gAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCAc                                                                                                                                                                                                     <  1:254710/251‑1 (MQ=255)                                 aggTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCCAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACTCGCCGCGCATCTATCCCAAAACATGCAATACAGTTACGCGCGGTAACGACGAGGAAGTTACTGaatgccctg                                                                                                                                                                                                   >  2:110419/2‑240 (MQ=255)                                            ccaggcatttgcaccgcactgGAAATCTGCTCCCCAAGAATACGGTGGGTATTAAGTAAGAATTCAATCGCTTCAAGATTCCTGCTCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGAAGCGCAGCGCATCAATACCAAAACGGGCAATCCAGTAAAGCGCGTTAACGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTa                                                                                                                                                                                     <  2:375717/240‑1 (MQ=255)                                             caaccATTGTCACCGCACCGAAAACCTGCTCACAAAGAATACGTTTTGTATTACGCAAGAAGTCAATCGCTCCACGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCACATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAg                                                                                                                                                                                    <  2:454338/248‑1 (MQ=255)                                              aGCCATTGTCACCTCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGt                                                                                                                                                                                      >  1:258955/1‑248 (MQ=255)                                                           gcaccgGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCa                                                                                                                                                                      <  1:26884/251‑1 (MQ=255)                                                           caaccgGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAAGCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCa                                                                                                                                                                      <  2:142384/249‑1 (MQ=255)                                                              ccgGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTc                                                                                                                                                                   <  1:103614/251‑1 (MQ=255)                                                                             ccAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTc                                                                                                                                                    <  1:114119/251‑1 (MQ=255)                                                                                           ggTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCACGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGATCATAGTTGTAGAGCAGCGTGTTCCATGCCTGCTGATAGCCTTCAAGCGGATAGCtgg                                                                                                                                      >  2:226884/1‑250 (MQ=255)                                                                                                   gCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCCAAACGTGCAATCCCATAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATCGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATAc                                                                                                                              >  2:294592/1‑251 (MQ=255)                                                                                                     aaGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACa                                                                                                                             >  1:395670/1‑250 (MQ=255)                                                                                                     aaGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGAGCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCTCTGTGTTCATACaa                                                                                                                            <  1:353963/251‑1 (MQ=255)                                                                                                           gcAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCACCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCg                                                                                                                       <  1:22423/249‑1 (MQ=255)                                                                                                            cAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTg                                                                                                                     >  1:260072/1‑251 (MQ=255)                                                                                                                cGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCAt                                                                                                                 >  1:29481/1‑251 (MQ=255)                                                                                                                      aaGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGGCGCCACCGCATCGACGCGCAGCGCATCAGTACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGATGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGGGCAATCAAAtt                                                                                                           >  2:400111/1‑251 (MQ=255)                                                                                                                         aTTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCgg                                                                                                        <  1:117882/251‑1 (MQ=255)                                                                                                                                cgcccgccAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCAGCGACGCGCAGCGCAGCAATACCAAATCGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAATTTACTGACTTCACGGCGACCATAGTTGTAGAGCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGTGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAg                                                                                                    <  2:284385/248‑1 (MQ=255)                                                                                                                                       cAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGtca                                                                                          <  1:368081/251‑1 (MQ=255)                                                                                                                                                            cccccTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATTCAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTgg                                                                      <  1:300396/250‑1 (MQ=255)                                                                                                                                                            ccaactctTTACGGCTGTAGTCTCGATAAATCATTTACTCCACCGCATCGACGCGCAGCGACTCACTACCACAACGTGCAATCCAGTCAAGCGCGTTCCCGCCGCATAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCATTCCTGATTATAGCCTTCCCGCGGATCCCTGTGTTCATACAAGTTCTAGCCATCACATTCTGCAAACGCAGAGacatcttccggacattggccttcc                                                                     >  2:351797/5‑232 (MQ=255)                                                                                                                                                             ccccTCCTTATGGCTGTAGTCACGAGAAATCCTGGACGCCTCTGCATCTACGCGCAGCGCACCAATTCCAAAACTTGCACTCCAGGTGATCGCGTGACCGACGAGGAAGTTACTGCCGTCACGGCTGTCATAGTTGTAGCTTAGCGTGTTCCTGTCCGGATGATAGCCTTCACTCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAACGCCAAGTCATCAGTCGGTAAGCGGCCTGGCa                                                                    <  2:435117/251‑1 (MQ=255)                                                                                                                                                                               gTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGt                                                  <  1:51362/251‑1 (MQ=255)                                                                                                                                                                                 cgcgATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTAATACACGTTCGGGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGGGGCCGGGCACCCAGGCGAGAATCAAGTTc                                                >  2:311500/1‑251 (MQ=255)                                                                                                                                                                                     aTAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGAc                                            >  2:141630/1‑251 (MQ=255)                                                                                                                                                                                         aTCAATGACGCCATCGAAGCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCTCGTTACCGACGAGGAAGTTACTGTCTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTTATGATAGCCTTCACGCGGATCGCTGTGTGCATACAAGTTCTTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTTGTCTGGCACCCAGTCGAGAATCACGTTCATACCAGc                                        <  2:33040/251‑1 (MQ=255)                                                                                                                                                                                           cATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTg                                      <  1:433396/251‑1 (MQ=255)                                                                                                                                                                                              tGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCTTACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTAGAGCATCACGTTCAGACCAGCTGCgg                                   >  2:127335/1‑250 (MQ=255)                                                                                                                                                                                              tGACGCCACCGCAGCGAAGCGCAGCGCATCAATACCAAAACGTGCAAACCAGTAAAGCGCGTTAACTACGCGGAAGTTACTGACTTATCGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATTATAGCCTTCACGCGGAGCGCTTTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCg                                    <  2:367005/250‑1 (MQ=255)                                                                                                                                                                                                aCGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCgtgt                                 >  1:225424/1‑251 (MQ=255)                                                                                                                                                                                                  gCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAACTTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGc                               >  2:230465/1‑251 (MQ=255)                                                                                                                                                                                                    cACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCgg                             >  1:298147/1‑251 (MQ=255)                                                                                                                                                                                                          aTCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGATCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGACACCCGTTCAGACCAGCCGCGTGGGCGGCATCa                        >  2:323608/1‑250 (MQ=255)                                                                                                                                                                                                            cGACGCGCAGCGCATCAATTCCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATg                     <  1:348879/251‑1 (MQ=255)                                                                                                                                                                                                              aCGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCCGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGaa                   <  1:21084/251‑1 (MQ=255)                                                                                                                                                                                                                gcgaagcgcaTCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAAt                 <  1:164059/247‑1 (MQ=255)                                                                                                                                                                                                                       gcaTCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAATAACGGaa          <  1:127335/251‑1 (MQ=255)                                                                                                                                                                                                                        caTCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAATAACGGa           >  1:27842/1‑249 (MQ=255)                                                                                                                                                                                                                             aTACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGCAGTTACTTACGTCACGGCGACCATAGTTGTAGACCATAGTGTTCCAGTCCGGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAATTTCGTGCCATCGAATTCGGCAAGCTCAACGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGATCAGACCAGCGGCGTGTGCGGCATCATTCAAATAACGGAAttcgtc    <  2:410437/251‑1 (MQ=255)                                                                                                                                                                                                                               aCCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGAAAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAATAACGGAAGTCGTcgc  >  2:149212/1‑251 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 CATATCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAATAACGGAAGTCGTCGC  >  NC_000913/3572087‑3572559

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