Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,646,970	G→T	Q224H (CAG→CAT)	gor →	glutathione reductase (NADPH)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,646,970	0	G	T	94.4%	104.1 / NA	36	Q224H (CAG→CAT)	gor	glutathione reductase (NADPH)
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (21/13);  total (22/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

ACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCT‑GTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAGCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGG  >  NC_000913/3646732‑3647202                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         ccccGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCT‑GTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAAc                                                                                                                                                                                                                              <  1:101591/250‑1 (MQ=255)         atTCCGGGCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCT‑GTTTGTGCGTAAACATGCTCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATcc                                                                                                                                                                                                                       <  1:13204/250‑1 (MQ=255)                  gTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGGTTAAAGGCCCCGGCGCGGAAACGCATCT‑GTTTGTGCGTAAAAATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAAAGCTGGGCGGAGTGGTGAACGCCGAAGGCCCGCCTCTGCACCCCAACGCCATCCCGAAAGCggg                                                                                                                                                                                                             >  1:67688/1‑250 (MQ=255)                     gAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGGGGGATTAAAGGCCTCGGCGCGAAAAAGCATCTTGTTTG‑GGGGAAAAATGCGCCGCTGCGCAGCCTCGCCCCGGTGATTTCCGAAACGCTGGTCGAAGGGATGAACGCCCAAGGCACGCATCTGCACACAAACACCCTCCCGGAAGCgggagg                                                                                                                                                                                                          >  1:369301/1‑247 (MQ=255)                           ggTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGGGTGATTAACGGCCTCGGCGCGAAAACGCATCT‑GTTTGTGCGTAAAAATGCACCGCTGCGCAGCTTCGACCCGGTGGTTTACGAAACGCTGGGCGAAGTGGTGAAAGCCCAAAGCCCGCCTCCGCACACCAACACCATTCCCAAAACGGGAGTGgaaa                                                                                                                                                                                                     >  1:8053/1‑250 (MQ=255)                              aTTGATTCTGATGGCTTCTTCGCCCTTCCTGCGTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCT‑GTTTGTGCGGAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATa                                                                                                                                                                                                  <  1:162443/250‑1 (MQ=255)                                                               ttGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGGGATTAAATGTCTCAGCGCGAAAAGCCATCT‑GTTTGTGCGTAAAAATGAGCCGCTGCGGCGCTTTGACACGATGATTTCCCAAAAACTGGTCGAAGGGGGGAACGCCGAAGGGCCGCATATGGCCACCAAAGCCATCCCGAAAGCGGGAGTGAAAAATACCGATGGGAGCCTGACGCTGGGGCTGGAAGATgg                                                                                                                                                                >  1:36894/1‑251 (MQ=255)                                                                       gcgcGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCT‑GTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTg                                                                                                                                                        <  1:11489/251‑1 (MQ=255)                                                                                      tgGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGAGCGAAAACGCATCT‑GTTAGTGAATAAACATGCGCCGCTGCGCCGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGGGATGAACGCCGGAGGCCCGCATCTGCACACCAACGCCATCCCCCAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGGCGCAGCGACACGGTGGATTGcc                                                                                                                                         >  1:252891/1‑251 (MQ=255)                                                                                       tgctctgtttaatgtcccGTGGAGCTGGCGGGCGTTATTAACGGCCGCGGCTCGAAAACGCATCT‑GTTTGTGCGGAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATCTACGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTTACGCTGGAGATGGAATATGGTCGCAGTGAAACGGTGGATTGcc                                                                                                                                         <  2:44752/234‑1 (MQ=255)                                                                                         cgcgGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCT‑GTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGa                                                                                                                                      <  1:293263/251‑1 (MQ=255)                                                                                                 aCATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCT‑GTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGGGGGTTGCCTGATTTGGGCg                                                                                                                              >  1:13653/1‑251 (MQ=255)                                                                                                  cATCGCCGTTGAGCTGGCGTGCGTGATTAACGGCCTCGGCGCGAAAACGCATCT‑GTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCg                                                                                                                              >  1:274607/1‑250 (MQ=255)                                                                                                    gcGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCT‑GTTTGTGCGTAAACAGGCGCCGCTGCGCAGCTTCGACCCGATGAGTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAGCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGAtg                                                                                                                           <  2:171285/250‑2 (MQ=255)                                                                                                          ttGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCT‑GTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGGGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTcgc                                                                                                                     >  1:2114/1‑251 (MQ=255)                                                                                                                     ggcgTGATTAACGGCCTCGGCGCGAAAACGCATCT‑GTTTGTGCGGAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGACACGATGGTCGAAGTGATGAACGCCAAAGGCCCGCATCGGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGAAGCCTAACGCTGGACATGGCAGATGGTCGCAGTGAAACGGAGGATTGCCTGATTTGGGCGCTGGGTCGCGAggcggcca                                                                                                           >  2:295683/1‑243 (MQ=255)                                                                                                                                       ggggcgAAAACGCTGTT‑GTTGGTGCGTAAACAGGCGCCGTTGCGCAGCTTCGACCCGATGATTTCCGAAACGCGGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCACCGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGGGCCTGCCAATGACAACATCAACCTg                                                                                          <  2:67688/246‑1 (MQ=255)                                                                                                                                               aaCGCATCT‑GTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCt                                                                                 <  1:207904/250‑1 (MQ=255)                                                                                                                                               aaCGCATCT‑GTTTGTGCGAAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGCAGTGATGAACGCCGACGCCCCGCATCTGCACACCAACGCAATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGGGGATTGCCGGATTTGGGCGATTGGTCGGAAGCCTGCAAATGACCACATCACCCTAGAAGCCGCt                                                                                 >  2:207904/1‑250 (MQ=255)                                                                                                                                                aCGCATCT‑GTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGGCCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCACCCTGGAAGCCGCTgg                                                                               <  1:194097/251‑1 (MQ=255)                                                                                                                                                  cgcatcgttttGTGCTTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGGAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGc                                                                              <  1:11329/243‑1 (MQ=255)                                                                                                                                                       t‑GTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCATGACGCTGGAGCTGGATGATGGCCGCAGAGAAACGGTGGATTGCCTGATTTGGGAGATTGGTCGCGAGCCTGCCAATGACAACATCAACCCGGAAGCCGCTGGAGTTaa                                                                         >  1:379988/1‑250 (MQ=255)                                                                                                                                                                  tAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGaa                                                               >  1:374964/1‑250 (MQ=255)                                                                                                                                                                      cATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACAGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGGCAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAAAGAAAAAgg                                                          >  1:1378/1‑251 (MQ=255)                                                                                                                                                                          cgccgcTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCtat                                                      >  1:155965/1‑251 (MQ=255)                                                                                                                                                                          cgccgcTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCACGCATCTGCACACCAACGCCATCCCGACAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGGGGATTGCCTGCTTTGGGCGATTGGTCGCGAGCCTGCCAAGGACAACATCAACCTGGAAGCCGCAGGCGTTAAACCTAACGAAAAAGGCtat                                                      >  2:185123/1‑251 (MQ=255)                                                                                                                                                                                           gACCCGCTGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTATTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCTAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAATGCTATATCGTCGTCGATAAata                                     <  2:155965/251‑1 (MQ=255)                                                                                                                                                                                           gACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAat                                      >  1:379783/1‑250 (MQ=255)                                                                                                                                                                                               cGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTAGCGAGCCTGCCAATGACAACATCAACCTGGAAGACGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCa                                  >  1:144118/1‑250 (MQ=255)                                                                                                                                                                                                    atTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAacac                            >  1:259043/1‑251 (MQ=255)                                                                                                                                                                                                     tttCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAGCTGCACACCAACGCCATCCCGAGAGCGGTAGTGAAAAATACCGCTGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACAcc                           >  2:160927/1‑251 (MQ=255)                                                                                                                                                                                                                  ggTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGt               >  1:315717/1‑250 (MQ=255)                                                                                                                                                                                                                      gAAGTGATGAACGCCGAAGGTCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTa          <  1:296269/251‑1 (MQ=255)                                                                                                                                                                                                                              gaACGACGAAGGACCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGCTAAATATCAAAACACCTATATTGAAGGTATTTACGCGGTgg  >  1:77217/1‑251 (MQ=255)                                                                                                                                                                                                                                  gCCGAAGGCCCGCATCTGCACACCCACGCCATCCCGAAAGCGGTAGTGCAAAATCCCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGACACGGTGGATTGCCTGCATTGGGCGATTGGTCGCGAGCCTGCCACTGACAACATCACCCTGGAAGCCGCTGGCGTTAAAACTACCGAAAAAGGCTATATCGTCGTCGATAAATATCAAAaca                             >  2:175259/1‑220 (MQ=255)                                                                                                                                                                                                                                  gCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAaca                             <  1:175259/220‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         ACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCT‑GTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAGCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGG  >  NC_000913/3646732‑3647202

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