Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,647,428	C→T	A377V (GCG→GTG)	gor →	glutathione reductase (NADPH)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,647,428	0	C	T	100.0%	196.0 / NA	57	A377V (GCG→GTG)	gor	glutathione reductase (NADPH)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/29);  total (28/29)

AAGGTATTTACGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGACA‑CCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGA  >  NC_000913/3647185‑3647668                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 aaGGTATTTACGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGACA‑CCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGtat                                                                                                                                                                                                                                           >  1:349890/1‑251 (MQ=255)  aGGTATTTACGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGACA‑CCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGtat                                                                                                                                                                                                                                           >  1:104141/1‑250 (MQ=255)   ggTATTTACGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGACA‑CCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGCTGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGtatcc                                                                                                                                                                                                                                         <  1:363428/251‑3 (MQ=255)    gTATTTACGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGACA‑CCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCCCCGTGATGTATAcc                                                                                                                                                                                                                                        >  1:17195/1‑251 (MQ=255)         tACGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGCCA‑CCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAACGTGTATAAATCCTCTTTCACCGTGATGTATAccgccgc                                                                                                                                                                                                                                   >  1:309937/1‑250 (MQ=255)         tACGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGACA‑CCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGt                                                                                                                                                                                                                                   >  2:72752/1‑251 (MQ=255)           cgcgGTGGGCGATAACACGGGTGCAGTGGAGCTGACA‑CCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGAATACCGCCGTc                                                                                                                                                                                                                                  >  2:193080/1‑250 (MQ=255)               gTGGGCGATAACACGGGTGCAGTGGAGCTGACA‑CCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTcaccac                                                                                                                                                                                                                             >  1:11110/1‑251 (MQ=255)                 gggCGATAACACGGGTGCAGTGGAGCTGACA‑CCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTc                                                                                                                                                                                                                           >  1:344252/1‑251 (MQ=255)                     gATAACACGGGTGCAGTGGAGCTGACA‑CCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCg                                                                                                                                                                                                                       >  2:113271/1‑251 (MQ=255)                              ggTGCAGTGGGGCTGACA‑CCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGcgtaaatac                                                                                                                                                                                                              <  2:17195/251‑9 (MQ=255)                              ggTGCAGTGGAGCTG‑CAGCGGGTGGAAGTTGCAGGGGGTCGCCTTCTCTCTGAACGCCTTTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCACCATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATCCCGCCGTCGCCACTCACCGCCAGCCGTg                                                                                                                                                                                                              <  2:247634/251‑1 (MQ=255)                                     tGGAGCTGACA‑CCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATg                                                                                                                                                                                                       >  1:446043/1‑251 (MQ=255)                                            gACA‑CCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTg                                                                                                                                                                                                 <  1:577143/250‑1 (MQ=255)                                                      ggCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCCCCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTg                                                                                                                                                                                       >  1:309269/1‑251 (MQ=255)                                                      ggCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGtt                                                                                                                                                                                        <  1:456132/250‑1 (MQ=255)                                                       gCAGTTGCGGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTgg                                                                                                                                                                                      <  2:434680/251‑1 (MQ=255)                                                       gCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGGGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTgg                                                                                                                                                                                      <  2:45590/251‑1 (MQ=255)                                                       gCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTgg                                                                                                                                                                                      <  2:344252/251‑1 (MQ=255)                                                              cggcggggcgcgttctctcTAAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTCAAGCAACATTCCGCACGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaa                                                                                                                                                                               <  2:152208/238‑1 (MQ=255)                                                                gCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaaga                                                                                                                                                                             <  2:130161/251‑1 (MQ=255)                                                                      cGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCAACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaagaga                                                                                                                                                                           >  1:202508/1‑247 (MQ=255)                                                                              ctctGAGCGCCCGTTTAAAAACAGGCCGGATGAGCATCTGGATAACAGCAACATTCCGCCCGTGGTCTTCACCCATCCGCCGATTGGTACGGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCTCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCgg                                                                                                                                                                 <  1:318407/249‑1 (MQ=255)                                                                                      gCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGc                                                                                                                                                       >  1:421031/1‑251 (MQ=255)                                                                                              aataaCAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCtt                                                                                                                                               <  1:117765/251‑1 (MQ=255)                                                                                               aaaccAAGCCGGAGGAGCATCTGGATTACAGCAACATTCCGACGGTGGTCTGCACCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCttt                                                                                                                                              <  2:24153/247‑1 (MQ=255)                                                                                                     aGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATg                                                                                                                                        >  2:326448/1‑251 (MQ=255)                                                                                                      gCCGGCTGAGCATCTGGATTACAGCAACATTCCGCCCGAGGTCTTCAGCCATCCGCCGCTTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGCAGATTGTCGGTATTAACGGCATTGGCTTTGGTCTgg                                                                                                                                       >  2:249975/1‑251 (MQ=255)                                                                                                         ggATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACg                                                                                                                                    <  1:225843/251‑1 (MQ=255)                                                                                                            tGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATACGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGaa                                                                                                                                  >  1:380858/1‑250 (MQ=255)                                                                                                                aaTCTGGATTAAAGAAAAATTCCGACCGGGGTCTTCAGCCATCCGCCGATGGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGGATGGCGACGATCAGGGGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGtt                                                                                                                             <  2:446043/250‑1 (MQ=255)                                                                                                                             gCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTcgcg                                                                                                                >  2:36210/1‑251 (MQ=255)                                                                                                                                  aTTCCGACCGGGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGgc                                                                                                           <  2:345507/251‑1 (MQ=255)                                                                                                                                     ccGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCt                                                                                                        <  1:364826/251‑1 (MQ=255)                                                                                                                                        aCCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGa                                                                                                      >  2:42420/1‑250 (MQ=255)                                                                                                                                           gTGGTCTTCAGCCACCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCCCCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGAt                                                                                                  <  2:306187/251‑1 (MQ=255)                                                                                                                                             ggTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATgg                                                                                                <  1:206313/251‑1 (MQ=255)                                                                                                                                                  tCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGc                                                                                            <  1:26028/250‑1 (MQ=255)                                                                                                                                                     gCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAAcc                                                                                        >  2:69758/1‑251 (MQ=255)                                                                                                                                                       cATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCa                                                                                       >  1:417311/1‑250 (MQ=255)                                                                                                                                                                         gTTGGTTTAACGGAACCGGAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCAGGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTGTGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGAGGTGGGCAACCAAAAAAGACTTCGACAAt                                                                      <  2:554264/249‑1 (MQ=255)                                                                                                                                                                                    ggAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTc                                                         >  1:79025/1‑251 (MQ=255)                                                                                                                                                                                      aaCCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCAc                                                       >  2:363957/1‑251 (MQ=255)                                                                                                                                                                                          gCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCa                                                    >  1:70820/1‑250 (MQ=255)                                                                                                                                                                                                        ggtatggcaacactCAGGTGAAATTGTATAAATCCCCTTTCACCGTGATGTATTCCGCCGTCACCAATCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGATATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACAGCGGCAGAagag                                     <  2:397712/238‑1 (MQ=255)                                                                                                                                                                                                         gTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGt                                    >  1:513100/1‑251 (MQ=255)                                                                                                                                                                                                           aTGGCGACGATCAGGTGAAAGTTTATAAATCCTCTTTCCCCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTc                                  <  2:349890/251‑1 (MQ=255)                                                                                                                                                                                                             ggCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGt                                >  1:201403/1‑251 (MQ=255)                                                                                                                                                                                                              gCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGGTGCAGGGCTTCGCGGTGGCGCTGAATATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTg                               <  1:343916/251‑1 (MQ=255)                                                                                                                                                                                                               cgacgtgcAGGTGAAAGGGTATAAATCCTCTTTCGCCGTGATTTATCCCGCCTGCACCACTCCCCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTGCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGGCTTCGACAATACCGTCGCCATTCACCCGACGGCGGCAGAAGAGTTCGTGa                              <  2:404046/244‑1 (MQ=255)                                                                                                                                                                                                                gacgaTCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGAc                             <  1:144686/251‑1 (MQ=255)                                                                                                                                                                                                                      cAGGTGAAAGTGTATAAATCCTCTTTCCCCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCTCCCAACGGCGGCAGAAGAGTTCGTGACAATGc                        <  1:350870/250‑1 (MQ=255)                                                                                                                                                                                                                          tGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTa                    >  2:529895/1‑250 (MQ=255)                                                                                                                                                                                                                          tGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCAAACGGCGGCAGAAGAGTTCGTGACAATGCGTTa                    >  1:41425/1‑250 (MQ=255)                                                                                                                                                                                                                            aaaGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAAt                 >  2:113001/1‑251 (MQ=255)                                                                                                                                                                                                                                      aaTCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGc       <  1:190139/251‑1 (MQ=255)                                                                                                                                                                                                                                              ttCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAaga  <  1:72752/248‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 AAGGTATTTACGCGGTGGGCGATAACACGGGTGCAGTGGAGCTGACA‑CCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGA  >  NC_000913/3647185‑3647668

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