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	95.5%	122.1 / ‑3.9	44	Q224H (CAG→CAT)	gor	glutathione reductase (NADPH)
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (17/25);  total (17/27)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

GGGCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAGCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATAACACGGGT  >  NC_000913/3646745‑3647216                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       gggCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGGTTAACGGCCTCGGGGCGAAAACGCATCTGTTTTTGCGTAAACATGCGCCGCTGCGCAGCCTCG‑ACCCGATGGTTTCCGGAAAGCTTGTCGGAGTGGTGAACGCCGAAGGCCCGCATCTGCACCACAACGCCATCCCgcaa                                                                                                                                                                                                                                >  1:424359/1‑247 (MQ=255)                 aTTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGGGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACAGGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAATGCGGTAGTGAAAAATAc                                                                                                                                                                                                               <  2:318213/251‑1 (MQ=255)                  gtgagtCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTGGAGCTGGCGGTCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGGAAACATGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCTGTAGTGAAAAATAc                                                                                                                                                                                                               <  2:321002/245‑1 (MQ=255)                            tGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGCGCTGGCGGGCGGGATTAACGGCCCCGGGGCGCAAACGCATCTGTTTGTGCGTAAACCTGGGCCGCTGCGCAGCCTCGGCCCCGTTATTTTCCGAAA‑GCTGGGCGGAGGGATGGAAGCCGAAGGCCCGCCTCTTCGCACAACCCCCATACCGCAAGCGGGAGTGGAAAAAACCGATGGtgcgc                                                                                                                                                                                                    >  2:33375/1‑247 (MQ=255)                             ggCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGGAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGGAGTGATGGACGCCGAAGGCCCGCCTCTGCACACCAACGCCATCCCGGAAGCGGTAGTGAAAAATAACGATGGGAGCCt                                                                                                                                                                                                   >  1:49892/1‑251 (MQ=255)                             ggCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGcc                                                                                                                                                                                                    <  1:381963/250‑1 (MQ=255)                                           tCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTAAATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACgctggagctg                                                                                                                                                                                      <  2:49892/250‑1 (MQ=255)                                            ccTGCTTGGCCAGAGCGCGTGGCGGTTGGTGGCGCGGGTTACATCGCCGGTGAGCTGCCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAGGCGGTAGTGAAAAATACCGATGGTAGCCTGACgctggagctgg                                                                                                                                                                                     <  2:97494/250‑1 (MQ=255)                                                           cgcgTGGCGGTGGTTGGCGCTGGTTACATCGCAGTTGAGCTGGCGGGAGTGATGAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGa                                                                                                                                                                     <  2:125296/251‑1 (MQ=255)                                                           ccgtgtacggtggtgtgcgcGGGGTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGTCGCGAAAACGCATCTTTTTGTGCGTAAACAGGCGCCGCTGCGCAGCTTCG‑ACCCTATGATTTCCGAAACGCGGGTCGAAGTGATGAACTCTGAAGGCACGCATCTGCACACCAATGCCAGCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACTATGGTTTTCGAAGATGGTCGCAGTGa                                                                                                                                                                     <  2:213490/235‑1 (MQ=255)                                                                  gggtggttgGGGCGGGGTTCATCGCCGTGGAGCTGGCGGGCGTGATTAACGGCGTCGGAGCGAAAGCGCAGCTGTTTGTGCGTAAACATGCTCCGCTGCGCAGCGTCG‑ACCCGATGATTGCTGAAACGCTGGTCGAAGTGATGAACGTCGAAGGCCCGCAGATGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGACTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTg                                                                                                                                                              <  2:146905/246‑1 (MQ=255)                                                                       gttggggcgggttcacgcGCCGGTGAGCGGGCGGGCGTGATTAACGGCATCGGCGCGAAAACGCATCTGTTTGTGCGTAAACAGGCGCCGCGGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTg                                                                                                                                                         <  2:254593/234‑1 (MQ=255)                                                                         tgGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGAGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGcc                                                                                                                                                       >  1:74647/1‑251 (MQ=255)                                                                                     cATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCgg                                                                                                                                           >  1:503004/1‑250 (MQ=255)                                                                                        cGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCG‑CCCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAAAGCCATCCCGAAAGCGGTAGTGAAAAATACCGCTGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTg                                                                                                                                        >  1:216989/1‑251 (MQ=255)                                                                                          ccGTTTAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCTGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGt                                                                                                                                      <  1:413279/251‑1 (MQ=255)                                                                                                                tAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGCCGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGCTTTGGGCGATTGGTCGCGAGCCTGCCAATGacaac                                                                                                                 >  2:208446/1‑250 (MQ=255)                                                                                                                 aaCGGGCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCTGTTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCTCGAAAGCGGTAGTGAAAAATACCGATGGTAGCTTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCATGCCAATGACAACAt                                                                                                               <  1:336133/251‑1 (MQ=255)                                                                                                                     gCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAAc                                                                                                           <  1:464593/251‑1 (MQ=255)                                                                                                                      cccgggCGCGAAAACGCATTTGTTTGGGCTTAAACATGCCCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGGGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAAcc                                                                                                          <  2:365264/247‑1 (MQ=255)                                                                                                                           gcgcGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGaa                                                                                                     <  1:128245/251‑1 (MQ=255)                                                                                                                              cgAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGcc                                                                                                  <  1:85696/251‑1 (MQ=255)                                                                                                                               gAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCCGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCg                                                                                                 >  2:484134/1‑251 (MQ=255)                                                                                                                                  aaCGGATCGGTTGGTGCGTAAACAGTCGCCGCTGCGCAGCTTCG‑ACCCGAGGATTTCCGAAACGCAGGACGAAGTGATGAACGATGAAGGCCCGCATTCGAACGCCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCGTGACGCTGGAGCTGGAAGATGGTAGCAGTGAAACGGGGGATTGCCTGATGGGGGCGATTGTTCACGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTg                                                                                              <  2:142883/251‑1 (MQ=255)                                                                                                                                   aCGCATCTGTTTGGGGGTAAACATGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCGGGTCGAAGTGATGACCGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGTAAGATGGTCGCAGGGAAACGGTGGATTGCCTGATTGGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTgg                                                                                             <  2:52792/251‑1 (MQ=255)                                                                                                                                        tttttgtgtgCTTAAACATGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACTCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGCCGCAGTGAAACGGTGGATTGCCGGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGtt                                                                                         <  2:371704/244‑1 (MQ=255)                                                                                                                                                    tAAACATGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGaa                                                                             >  1:92768/1‑250 (MQ=255)                                                                                                                                                      aaCATGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAGCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGaaaaa                                                                          <  1:229581/251‑1 (MQ=255)                                                                                                                                                                         cTTCG‑ACCCGATGGTTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGACAGAGGTAGCGAAAAATCCCGATGGTAGCCTGACGCTGGAGCTGGCAGATGGACGCAGTGAAACGGGGGCTTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCCACCTGGAAGCCGCTGGCGTTAAAACTAACGCAAACGGCTATATAGTCGTCGAt                                                        >  2:503360/1‑250 (MQ=255)                                                                                                                                                                             g‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACACCATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTCTATCGTCGTCGATAAatt                                                   >  2:132766/1‑250 (MQ=255)                                                                                                                                                                               aCCCGATGAGTTCCGAAACGCTGGGCGAAGGGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGGGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATAGCGTCGTCGATAAatat                                                  <  2:503004/251‑1 (MQ=255)                                                                                                                                                                                             gAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAatat                                    <  1:5935/251‑1 (MQ=255)                                                                                                                                                                                             gAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAAGACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGTCTATATCGTCGTCGATAAATATCAAAACACCAatat                                    <  1:161294/251‑1 (MQ=255)                                                                                                                                                                                                 cGCTGGGCGAGGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCGGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGCTACAACTAACGAAAAAGGCTATATCGGCGTCCATAACTATAAAAACACCactcttgc                                 >  2:211591/1‑243 (MQ=255)                                                                                                                                                                                                     ggTCGAAGTGAGGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTa                            <  1:397865/251‑1 (MQ=255)                                                                                                                                                                                                      gTCGAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTAt                           <  1:33375/251‑1 (MQ=255)                                                                                                                                                                                                        ggAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTAttt                         <  2:148762/250‑1 (MQ=255)                                                                                                                                                                                                         gAAGTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTAttt                         >  1:148762/1‑250 (MQ=255)                                                                                                                                                                                                            gTGATGAACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGAGGTAGTGAAACATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCATAACACCAATATTGAAGGTATTTaagc                     >  2:505587/1‑248 (MQ=255)                                                                                                                                                                                                             tgatgaACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAGGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTAcgcg                    >  1:131013/1‑251 (MQ=255)                                                                                                                                                                                                               atgaACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGt                  <  1:484134/251‑1 (MQ=255)                                                                                                                                                                                                                 gaACGCCGAAGGCCCGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTgg                <  1:81276/251‑1 (MQ=255)                                                                                                                                                                                                                        gAAGGCACGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGCTAAATATCAAAACACCAATATTAAAGGTATTTACGCGGTGGGCGATa          >  2:219102/1‑250 (MQ=255)                                                                                                                                                                                                                               cGCATCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAACCTAACGAAAAAGGCTATATCGTCGTCGCTAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATAACACgggg  >  2:337454/1‑250 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       GGGCGTGGAATACGGTATTGATTCTGATGGCTTCTTCGCCCTTCCTGCTTTGCCAGAGCGCGTGGCGGTTGTTGGCGCGGGTTACATCGCCGTTGAGCTGGCGGGCGTGATTAACGGCCTCGGCGCGAAAACGCATCTGTTTGTGCGTAAACATGCGCCGCTGCGCAGCTTCG‑ACCCGATGATTTCCGAAACGCTGGTCGAAGTGATGAACGCCGAAGGCCCGCAGCTGCACACCAACGCCATCCCGAAAGCGGTAGTGAAAAATACCGATGGTAGCCTGACGCTGGAGCTGGAAGATGGTCGCAGTGAAACGGTGGATTGCCTGATTTGGGCGATTGGTCGCGAGCCTGCCAATGACAACATCAACCTGGAAGCCGCTGGCGTTAAAACTAACGAAAAAGGCTATATCGTCGTCGATAAATATCAAAACACCAATATTGAAGGTATTTACGCGGTGGGCGATAACACGGGT  >  NC_000913/3646745‑3647216

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