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	97.7%	130.1 / NA	44	Q224H (CAG→CAT)	gor	glutathione reductase (NADPH)
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (21/22);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAGCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATAA  >  NC_000913/3646723‑3647209                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                gtcCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGGGATGAACGCCGAAGGCCCGCATCTg                                                                                                                                                                                                                                              >  1:90136/1‑251 (MQ=255)                atatTCCGGGCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCTGTTGAGCTGGCGGGCGGGATTAAAGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGGAAACCTGCGCCCCTGCGGCGCTTCCACCCGGTGATTTCCGAAAACCTGGTCGGAGGGATGAACGCCCAAGGCCCGCCTCTGCACACCAAcccccac                                                                                                                                                                                                                               >  2:122298/1‑245 (MQ=255)                           gTGGAATACGGTAGTGAGTCTGATGGCTTCTTCGCCCTTCCTGCTGGGCCAGAGCGCGGGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATGAACGGCGTCGGCGCGAAAACGCATCTGGTTGTGCGTAAACAAGCGCCGCGGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGGGATGAACGCCGACGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGt                                                                                                                                                                                                                    <  2:160190/251‑1 (MQ=255)                               aaTACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAAGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGGAGTg                                                                                                                                                                                                                >  1:95013/1‑251 (MQ=255)                                                   ggCTTCTCCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAATGCGGTAGTGAAAAATACCGATGGTAGcc                                                                                                                                                                                             <  1:216394/250‑1 (MQ=255)                                                         ttcGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGGTTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGGAAAAATGCGCCGCCGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGGTGAACGCCGAAGGCCCGCCGCTGCACCCCCACGCCCTCCCGAAAGCGGTAGGGGAAAATAccagagaacggaaacgcgg                                                                                                                                                                                      >  2:147012/1‑234 (MQ=255)                                                         ttcGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGa                                                                                                                                                                                                     <  1:147012/236‑1 (MQ=255)                                                         ttcGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGGGATTAAAGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGACGGGATGAACGCCGAAAGCCCGCATCTTCACACCAACGCCATCCCGAAAGCGGGAGGGAAAAATCCCGATGGTAGCCTTACgct                                                                                                                                                                                      >  1:177783/1‑251 (MQ=255)                                                                 tCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTCCCTCGCCGTTGAGCTGGCGGGCGTGATTAAAGGCCTCGGCGCGAAAACGCCTCTGTTTGTGCGTAAACAAGCCCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGGGATGAACGCCGAAGGCCCGCATCTGCCCACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCCGACgctggggctgg                                                                                                                                                                              >  1:78460/1‑251 (MQ=255)                                                                              gagCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGGTGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGGAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGgcgc                                                                                                                                                                   >  1:216470/1‑249 (MQ=255)                                                                                    gTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGTTGGAAGATGGTCGCAGTGAAAc                                                                                                                                                           <  1:163939/251‑1 (MQ=255)                                                                                             gttgGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATGCCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTg                                                                                                                                                  <  2:93659/251‑1 (MQ=255)                                                                                               tgGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGc                                                                                                                                                 <  1:102899/250‑1 (MQ=255)                                                                                                     gggTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAAAGGTGGATTGCCTGAttt                                                                                                                                          >  1:72095/1‑251 (MQ=255)                                                                                                            aTCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGGGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCCTCCCTAAAGCGGTAGGGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCCGTGACACGGTGGATTGCCTGGTTTGGGCGAt                                                                                                                                   >  2:233431/1‑251 (MQ=255)                                                                                                                 cGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTc                                                                                                                              <  1:203854/251‑1 (MQ=255)                                                                                                                                      tAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGacaac                                                                                                          >  1:42102/1‑250 (MQ=255)                                                                                                                                       aaCGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACTCCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTTGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACAt                                                                                                        <  2:177783/251‑1 (MQ=255)                                                                                                                                        aCGGCCTAGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATc                                                                                                       <  1:122298/251‑1 (MQ=255)                                                                                                                                          ggCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGGGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAACAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCa                                                                                                      >  1:215896/1‑250 (MQ=255)                                                                                                                                             cgcGGCGCGAAAACGCATCTGTTTTGGCGTAAACATGCGCCGCGGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAATATGGTCGCAGTGAAACGGTGGATTGCCGGATTTGGGCGATTGGGCGCGAGCCTGCCAATGACAACATCAACCt                                                                                                  <  2:39454/249‑1 (MQ=255)                                                                                                                                                 gcgcGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTgg                                                                                                                                        >  1:120380/1‑209 (MQ=255)                                                                                                                                                 gcgcGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTgg                                                                                                                                        <  2:120380/209‑1 (MQ=255)                                                                                                                                                        aaCGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTg                                                                                       <  2:216470/251‑1 (MQ=255)                                                                                                                                                        aaCGCATCTGTGTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGATAACATCACCCTGGAAGCCGCTg                                                                                       <  1:212569/251‑1 (MQ=255)                                                                                                                                                                 gTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTaaa                                                                               <  2:72095/250‑1 (MQ=255)                                                                                                                                                                         gTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGCGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACCTCAACCTGGAAGCCGCTGGCGTTACAACTAACGa                                                                       >  1:222645/1‑250 (MQ=255)                                                                                                                                                                             aCATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGaaaa                                                                    >  1:83178/1‑249 (MQ=255)                                                                                                                                                                             aCATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAg                                                                  <  1:65561/251‑1 (MQ=255)                                                                                                                                                                                  cgccgcTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGGGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCtat                                                             <  2:69130/251‑1 (MQ=255)                                                                                                                                                                                   gccgcTGCGCAGATTCGACCCGATGATTTCCGCAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCtat                                                             >  1:69864/1‑250 (MQ=255)                                                                                                                                                                                    ccgcTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTatat                                                           <  2:90136/251‑1 (MQ=255)                                                                                                                                                                                       cTGCGCCGCTTCGCCCCGCTGATTTCAGAAACGCTGGTCGAAGTGCTGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCACACTGGAAGCCGCTGGCGTTAAAACTAACGCAAAATGCTATCtcgt                                                        >  2:41995/1‑251 (MQ=255)                                                                                                                                                                                            cAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCTCACCAACGCCATCCCTAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATAtcgtcgtcg                                                   <  2:75057/251‑1 (MQ=255)                                                                                                                                                                                             aGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAACGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATAtcgt                                                        >  1:245958/1‑245 (MQ=255)                                                                                                                                                                                             aGCTTCGACCCGATGATTGCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAAGACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATGGCCTGATTTGGGCGATTGGTCGCTATCCTGCCAATGACAACATCAACCTGGAAGCCGCTTGCGTTAAAACTAACGAAAAAGGCTATAtcgt                                                        <  2:245958/245‑1 (MQ=255)                                                                                                                                                                                                         atgatTTCCGACAAGCCGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAAAGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCaaa                                       >  2:158771/1‑250 (MQ=255)                                                                                                                                                                                                           gatTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAaca                                    <  2:95013/251‑1 (MQ=255)                                                                                                                                                                                                            atTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCAGCACACCAACGCCATCCCGAAAGCGGTGGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAacac                                   >  1:171067/1‑251 (MQ=255)                                                                                                                                                                                                                ccGAAACGCTGGTAGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGACATCCCGAAAGCGGTAGTGAAAAATACAGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCACTGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAACCACAAat                               >  2:177714/1‑251 (MQ=255)                                                                                                                                                                                                                               aagggatgaACGCAGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCGGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTAc                <  2:78460/247‑1 (MQ=255)                                                                                                                                                                                                                                       aaCGCCGAAGGCACGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTGAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGAggg        >  1:203312/1‑251 (MQ=255)                                                                                                                                                                                                                                            cgaagggtcGCATCTGCACACCAACGCCATCCCGAAAGCGGTCGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCAGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGAATTTACGCGGTGGGCGAt    >  2:204515/9‑250 (MQ=255)                                                                                                                                                                                                                                             gAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCATTGAAACGGTGGATTGCCTGATGTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACGAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATaa  <  1:41995/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                GTCCGAGCCACCCGGATATTCCGGGCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCGACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAGCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATAA  >  NC_000913/3646723‑3647209

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: ‑