Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	936,278	1	.	T	24.1%	138.4 / 65.7	54	intergenic (‑12/+22)	ECDH1_RS04620/ECDH1_RS04625	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base . (35/6);  new base T (12/1);  total (47/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-02

GCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCTTTTTTTTA  >  DH1/936141‑936428                                                                                                                                           |                                                                                                                                                       tcttcgcgaCCCCACTTCACCGTACTCCGTCGACTGCCTCATATTTCATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTtcatca                                                                                                                                            <  5:971152/142‑1 (MQ=14)                  cACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACAGTAATCC‑ATGACTTCATCACATTTTAATTCGTTGc                                                                                                                            >  4:775942/1‑149 (MQ=17)                                                        cGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCaaa                                                                                                                                          <  2:63857/97‑1 (MQ=35)                                                                           gacgacCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCCCGGACTGG‑GTTTCTTCTGCAATGTGGTATTGGTTGCGGCGGCCGGATTTGACCCGCCGGCCTTAGGGTTCTGATCGAGCCGAGCGAACAGgctgc                                                                 >  7:952632/1‑151 (MQ=2)                                                                                 cTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGGGGAGTGC‑TAGTCTTCTGCGAATTTTAATTGGTTGCGGGGTCGGGTTATGACCCGCCGCCCTTCGGGTTCTTAGCCCGCCGAGCGACCAGGCTGCAcccc                                                            >  7:1014158/1‑148 (MQ=2)                                                                                  ttCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACGTCATGAAATTTTa                                                                                                                                     >  5:909856/1‑76 (MQ=255)                                                                                                         tACCAGGCTGCTCCACCCCGCGTCAGGGTAGTGG‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCGGGATTTGAACCGACGACCTTAGGGTTATGATCCAGAAGAGCTACCAGGCTGCTCCCCCCCGAGTCAGTGGATTCGCACtat                                    >  3:1504848/1‑150 (MQ=17)                                                                                                                ctgctCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCTGATTTGAACCGACGCCCTTCGGGTTATGAGCCAGACGAGCTACCAGGCTGCTCCACCCCGCTTCGGTTGATGGGCACTATACTCCGCt                            >  1:461826/1‑151 (MQ=17)                                                                                                                ctgctCCACCCCGCGTCACCGTACTGC‑GTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGCCGAGCTACCAGGCTGCTCCCCCCCGCGTCCGTGGATGCGCACTATACTCCGcg                            >  7:1056501/1‑150 (MQ=255)                                                                                                                ctgctCCACCCCGCGTCACAGTAATCC‑ATTAATTCAGCAAAGTATAATTCATTCCGCAGGACGGATTTGAAACGACGACCTTCGGGTTATGAGCCAGACTAGATACCAGGCTGCTCCACCCAGCGTAAGAAGATGCGCACTATACTCCGCt                            >  6:1378821/1‑151 (MQ=12)                                                                                                                 tgctCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACATTAGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTc                           >  7:1166016/1‑151 (MQ=255)                                                                                                                    tCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTcgcg                        >  3:658780/1‑151 (MQ=255)                                                                                                                     ccACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGCGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGt                       >  5:61396/1‑151 (MQ=255)                                                                                                                      cACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGtt                      >  7:127704/1‑151 (MQ=255)                                                                                                                       aCCCCGCGTCACCGTACTGCTTTTACTTAATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCAGTGGAGTCGCAGTATACACCCCTCTCGttt                     >  3:1096078/1‑151 (MQ=11)                                                                                                                       aCCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGttt                     >  2:750206/1‑151 (MQ=255)                                                                                                                       aCCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGttt                     >  5:996259/1‑151 (MQ=255)                                                                                                                       aCCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATAATCCGCTCGCGttt                     >  6:575157/1‑151 (MQ=255)                                                                                                                         cccGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGATCGCGttttg                   >  1:480857/1‑150 (MQ=14)                                                                                                                         cccGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCCGGGAAGCCGACTATAcatccaaagagttttc                   >  7:500767/1‑136 (MQ=11)                                                                                                                         cccGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTc                   >  1:1270177/1‑151 (MQ=255)                                                                                                                          ccGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTc                   >  2:434263/1‑150 (MQ=255)                                                                                                                           cgcgTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCAt                 >  7:673001/1‑151 (MQ=255)                                                                                                                             cgTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCCTGGATGGGCAGTATACTACAATCGTGTTTTcgttg               >  5:719431/1‑147 (MQ=11)                                                                                                                             cgTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGc               >  3:936101/1‑151 (MQ=255)                                                                                                                               tCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTACGCTAGAGTTttcattca              >  8:1201263/1‑150 (MQ=14)                                                                                                                               tCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCa              >  6:1149148/1‑150 (MQ=255)                                                                                                                               tCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCa              >  5:1541545/1‑150 (MQ=255)                                                                                                                                cACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGTCGACCTTCGCGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAAc            <  5:1373289/151‑1 (MQ=255)                                                                                                                                 aCCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCAATATACACCAATCGCTTATTCATGCAAcc           >  3:1414503/1‑151 (MQ=11)                                                                                                                                 aCCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAAcc           >  7:946638/1‑151 (MQ=255)                                                                                                                                 aCCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAAcc           >  4:842146/1‑151 (MQ=255)                                                                                                                                  ccGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATTCGGAATATACTACGCTCGAGATTTCATGCAACCt          >  4:971276/1‑151 (MQ=11)                                                                                                                                  ccGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGCGCCCGACGAGCTACCAGCCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCt          >  2:1189175/1‑151 (MQ=255)                                                                                                                                  ccGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCt          >  6:281523/1‑151 (MQ=255)                                                                                                                                   cGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGTCCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGATCGCGTTTTCATGCACCCtt         >  1:301676/1‑151 (MQ=14)                                                                                                                                   cGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtt         >  2:810800/1‑151 (MQ=255)                                                                                                                                    gTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGCCGAGCTACCAGGCTGCTCCACCCCGCGTCCCGGGAGTGGGAGACTACTCCGCTCGCTTTTTcgttccgggggg        >  6:689598/1‑140 (MQ=11)                                                                                                                                    gTACTGC‑TTTACTTCCACAAATTTTTATTGGTTGCGGCGGCCGGATTTGAACCGACGACCTTCGGTTTATGAGCCCGACGAGCTAGCAGGCTTCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCttt        <  7:375072/151‑1 (MQ=21)                                                                                                                                     tACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCaa             <  1:1491998/145‑1 (MQ=255)                                                                                                                                     tACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTAGCGTTTTCATGCaaagatcc       >  2:1508779/1‑145 (MQ=14)                                                                                                                                     tACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGAGGATGCGCACTATACTCCCCTCGAGTTTTCATGCAACCtttt       >  2:876605/1‑151 (MQ=14)                                                                                                                                     tACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtttt       >  4:1475609/1‑151 (MQ=255)                                                                                                                                     tACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtttt       >  1:109663/1‑151 (MQ=255)                                                                                                                                     tACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtttt       >  7:1499944/1‑151 (MQ=255)                                                                                                                                      aCTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGGGGAAGCGGTCTATTCTCCGATCGAGttttgtttcacccgtttt      >  3:1040012/1‑137 (MQ=11)                                                                                                                                      aCTGC‑TTTACTTCATCAAATTTAAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCTGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCttttt      <  8:1062321/151‑1 (MQ=255)                                                                                                                                       cTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtttttt     >  2:916878/1‑151 (MQ=255)                                                                                                                                        tGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCttttttt    >  8:764472/1‑151 (MQ=255)                                                                                                                                        tGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCttttttt    >  7:124214/1‑151 (MQ=255)                                                                                                                                        tGC‑TTTACTTCATCAAATTTTAATTGGTTGCGCGGCCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGCTGCGCACTATACTCCGCTCGCGGTTTCATGCAACCtgttttt    >  6:97702/1‑151 (MQ=255)                                                                                                                                         gC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtttttttt   >  4:1380917/1‑151 (MQ=255)                                                                                                                                         gC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtttttttt   >  1:869463/1‑151 (MQ=255)                                                                                                                                         gC‑TTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCTTTTTTTTa  <  5:1306495/151‑1 (MQ=255)                                                                                                                                          c‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtttttttt   >  3:1025057/1‑150 (MQ=255)                                                                                                                                          c‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtttttttt   >  2:1137808/1‑150 (MQ=255)                                                                                                                                           |                                                                                                                                                       GCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCTTTTTTTTA  >  DH1/936141‑936428

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