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	100.0%	189.3 / NA	57	E395A (GAA→GCA)	glgB	1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (32/25);  total (32/25)

GAAGCCCAGACCGCCCATATCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAATAACGGAAGTCGTC  >  NC_000913/3572072‑3572557                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               gaatcctaGACCGCCCCTATCCTGCGGACGAGAACCGCAAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATGCAATCGCTTCAAGATTCTCGCGCCCGCCAAAGTCGTTCGGGATCCACGCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCaa                                                                                                                                                                                                                                             <  2:103529/244‑1 (MQ=255)     cccAGACCGCCCATATCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATcc                                                                                                                                                                                                                                          >  2:94139/1‑250 (MQ=255)          aCCGCCCATATCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTaa                                                                                                                                                                                                                                     >  1:278086/1‑250 (MQ=255)                   atCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTAc                                                                                                                                                                                                                           >  1:217368/1‑251 (MQ=255)                     ccTGCGGACGAGAAACGCCAGGGAAATCCGAAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCTCTCCCCCTCTTTACGGCTGTATTCGCGTTAAATCATTGACGCCACCGCATCGACGCGCAGCGCCTCAATACCAACACGTGCAATCCAGTAAAGCGCGTTAccc                                                                                                                                                                                                                         >  2:116294/1‑250 (MQ=255)                              gagaAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGa                                                                                                                                                                                                                  >  1:186366/1‑249 (MQ=255)                                  aaCGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACAGCACAGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTAGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGAAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGAGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGATAc                                                                                                                                                                                                            >  1:236482/1‑251 (MQ=255)                                       cAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTGCTGACt                                                                                                                                                                                                       <  1:94139/251‑1 (MQ=255)                                        aGGGAAATAGGTAGACTCCTCAGACATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTACCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCGTCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGaatt                                                                                                                                                                                                      >  2:207371/1‑248 (MQ=255)                                             aaTCGGTAGACGCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCTCTGCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCATCGCATCAATACCAAAACTTGCAATCCAGTAAGGCCCGTTACCGACGAGGAAGTTACTGACTTCACgg                                                                                                                                                                                                 <  1:182202/251‑1 (MQ=255)                                                 ggTAGACTCCTCAGCCATTGGCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGAc                                                                                                                                                                                             <  1:188649/251‑1 (MQ=255)                                                    aGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCGGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCAt                                                                                                                                                                                          <  1:79821/251‑1 (MQ=255)                                                     gACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCAt                                                                                                                                                                                          >  2:234829/1‑250 (MQ=255)                                                                   tGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGc                                                                                                                                                                            <  1:138838/250‑1 (MQ=255)                                                                       accgcaccgGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAgc                                                                                                                                                                                                                                  <  1:12202/192‑1 (MQ=255)                                                                       accgcaccgGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAgc                                                                                                                                                                                                                                  >  2:12202/1‑192 (MQ=255)                                                                       accgcaccgGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGtt                                                                                                                                                                       <  1:185135/251‑1 (MQ=255)                                                                           caccgGAAGCATGCTCACCAAGAATACGGTTGGTATTACGCACGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCAACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATACAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCCCGGCGACCATAGTTGTAGATCAGCGTGTTcc                                                                                                                                                                     >  2:342377/1‑249 (MQ=255)                                                                             ccgGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTc                                                                                                                                                                 <  1:321747/251‑1 (MQ=255)                                                                                 aaaCCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtg                                                                                                                                                              >  1:335892/1‑250 (MQ=255)                                                                                  aaCCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtgat                                                                                                                                                            >  1:316593/1‑251 (MQ=255)                                                                                   aCCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtgatg                                                                                                                                                           >  1:167194/1‑251 (MQ=255)                                                                                         tCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGcc                                                                                                                                                     >  1:251422/1‑251 (MQ=255)                                                                                          cACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCTCCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTGACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGc                                                                                                                                                      <  1:243810/249‑1 (MQ=255)                                                                                           aCCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTGCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCtt                                                                                                                                                   <  2:121589/251‑1 (MQ=255)                                                                                             cAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTc                                                                                                                                                  >  1:118954/1‑250 (MQ=255)                                                                                               aGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCAGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCaac                                                                                                                                               >  2:59425/1‑249 (MQ=255)                                                                                                 aaTACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTTCGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCAcgcg                                                                                                                                             <  2:278086/251‑1 (MQ=255)                                                                                                   tACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCCTCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGCGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGGGGTCCAGTCCTGATGATAGCCTTCAAGCgg                                                                                                                                            >  2:271526/1‑250 (MQ=255)                                                                                                        ttGGTATTACGCAAGAATTCAATCGCATCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATAAGCGTGTTACAGTCCTGATGATAGCCTTAACGCGGATCGCt                                                                                                                                      >  1:237115/1‑251 (MQ=255)                                                                                                          ggTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCACGACAACTTCGTTCGGGATCAACTCCCCCTCTTTCCGGCTGTAGTCGCGAAAAATCATTGCCGCCACCGCATCGACGCGCCGCGCCTCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTACCGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATAGCtgt                                                                                                                                    >  2:353495/1‑251 (MQ=255)                                                                                                           gTATTACTCAAGAATTCAATCGCTTCAAGATTCGCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCTTCAATACCAAAACGTGGAATCCAGTAAAGCGCTTTACCGACGAGGAAGTTACTGACCTCACGGCGACCATAGTTGTAGATCAGCTTGTGCCAGTCCTGATGATAGCCTTCACGCGGATCGCtgtg                                                                                                                                   <  1:353495/251‑1 (MQ=255)                                                                                                            tATTACGCAAGAATTCAATCGCTTCAAGCTTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCTCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCCCGGCGACCATAGTTGTAAATCCGCTTGTTCCAGTCCTGATGATCGCCTTCACCCGGATCGCtgttt                                                                                                                                  >  2:186215/1‑249 (MQ=255)                                                                                                             aTTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTgtgt                                                                                                                                  <  1:148902/250‑1 (MQ=255)                                                                                                                 cGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGGAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATAACAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCAt                                                                                                                              >  2:98490/1‑250 (MQ=255)                                                                                                                    aaGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACa                                                                                                                           >  1:186837/1‑250 (MQ=255)                                                                                                                           cAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTg                                                                                                                   <  1:305551/251‑1 (MQ=255)                                                                                                                              ccGCTTCAAGACTCTCGGGCCCTCCAAATTCGTTCGAGATGCACTCCCGCTCCTGACGGCTGTAGGCGCAATAAATCAGTGACGCCACCGCATCGACGCGCAGCTCATCAATTCCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTGCACGGCGACCATATTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGGTCGTGCCa                                                                                                                <  2:105093/250‑1 (MQ=255)                                                                                                                                gCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCCCTCGCCCTCTTTACGGATGTAGTCGCGACAAATCATTGACCACAACGCAACGACGCGCAAAGCATCACTACCAAAAAGTGGGATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGAAGATCAGCGTGCACCAGGCCTAATGATAACCTTCACGCGGATCGGTGTGCTCATACAAGTTCGTGCCATc                                                                                                              >  1:156001/1‑251 (MQ=255)                                                                                                                                   tCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCaa                                                                                                            <  2:89378/250‑1 (MQ=255)                                                                                                                                                      cAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGtca                                                                                        <  1:245527/251‑1 (MQ=255)                                                                                                                                                           tcgttcgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAg                                                                                    <  2:339998/250‑1 (MQ=255)                                                                                                                                                                     tCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTgg                                                                         >  1:294257/1‑251 (MQ=255)                                                                                                                                                                                    tACGGCTGTATTCGCGATAAATCATTGACGCCCCCCCTTCGACGCGCAACTCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGGTCGCTGTGTTCATACAAGTTCGTGCCATCAAAGTCTGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTTGCACCCAGTc                                                           <  2:313759/250‑1 (MQ=255)                                                                                                                                                                                    tACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCATCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGGTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTc                                                           <  2:167194/250‑1 (MQ=255)                                                                                                                                                                                        gCTGTAGTCGCGATAAATCATTGCCGCCACCGAATCGCAGCGCAGCGCATCAATACCAAAACGTGCAATCAAGTAAAGCGCGATACCGACGCGGAAGTTACTGACTTAACGGCGACCATAGTTGTAGATCAGAGTGTTCCAGTCCTGATGATAGCCTTAACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCCAATTCGGAAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCAACCAGTAgaga                                                       >  1:325715/1‑250 (MQ=255)                                                                                                                                                                                                  cgATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTc                                              >  2:74616/1‑249 (MQ=255)                                                                                                                                                                                                               gcGCCACCGCAGCGACGCGCAGCGCATCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCTACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGGCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCgtgt                               <  2:150679/250‑1 (MQ=255)                                                                                                                                                                                                                                cgcagcgcaTCCATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGCCGAGGAAGTTACTGACTTCACGGCGCCCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGACATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAAt               >  2:52501/1‑250 (MQ=255)                                                                                                                                                                                                                                cgcagcgcaTCAATACCAACACGTGCAATCCAGTACAGCGCGTTACCGAAGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGCTCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGCGAATCACGATCAGACCAGCTGCGTGTGCGGCATCAATGTAAt               >  2:42595/1‑250 (MQ=255)                                                                                                                                                                                                                                 gcagcgcaTCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCAt                                                                                                               >  2:267507/1‑153 (MQ=255)                                                                                                                                                                                                                                 gcagcgcaTCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCGGTGTTCATACAAGTTCGTGCCAt                                                                                                               <  1:267507/153‑1 (MQ=255)                                                                                                                                                                                                                                   agcgcaTCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCCCGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAATAACg           >  2:331884/1‑251 (MQ=255)                                                                                                                                                                                                                                         tCAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAATAACGGAAgtc     >  1:191705/1‑251 (MQ=255)                                                                                                                                                                                                                                          cAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGCAGTTACTGCCTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACACGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAATAACGGAAgtcg    >  2:172243/1‑251 (MQ=255)                                                                                                                                                                                                                                          cAATACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCAGCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAATAACGGAAgtc     <  1:30375/250‑1 (MQ=255)                                                                                                                                                                                                                                            aTACCAAAACGTGCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAATAACGGAAgtcgtc  <  2:294257/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               GAAGCCCAGACCGCCCATATCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACGCGGATCGCTGTGTTCATACAAGTTCGTGCCATCAAATTCGGCAAGCGCAAAGTCATCAGTCGGGAAGTGGCCTGGCACCCAGTCGAGAATCACGTTCAGACCAGCTGCGTGTGCGGCATCAATGAAATAACGGAAGTCGTC  >  NC_000913/3572072‑3572557

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