Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	936,278	1	.	T	25.0%	108.3 / 53.5	44	intergenic (‑12/+22)	ECDH1_RS04620/ECDH1_RS04625	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base . (29/4);  new base T (10/1);  total (39/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-02

CTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCTTTTTTTT  >  DH1/936165‑936427                                                                                                                   |                                                                                                                                                      cTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCCAGTAATGC‑GATAATTCATCAAATTATAtttccttgacgggtactt                                                                                                                  >  8:422928/1‑136 (MQ=14)                                                         cTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTTCTTGATGAAATTTTTATTGGTTg                                                                                                                            >  5:1010950/1‑85 (MQ=14)                                                         cTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTg                                                                                                                            <  6:1009960/85‑1 (MQ=255)                                                                 tATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCCGGTACTCC‑TATACTTCTACAAAATTTTATTCGTTGCGGCGGAGTGATTAGCACCTACGGCCTTCGGGTTATGAGCCCGCAGAGCAAACAGGCTGCTCCCCCCCGAGTcc                                                  >  7:690770/1‑151 (MQ=2)                                                                            cgaGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATgcgcc                                       >  7:427371/1‑150 (MQ=255)                                                                                 tACCAGGCTGCTCCACCCCGCGTCCCGGTAGGGG‑GTGACTTCAGCAATGTTTTATTGGTTGCGGGGGGCTGTTATGAACCGCCGCCCTTCGGCTTGTGATCCCGCAGAGCTACCACGCTGATCCACCCCGCGTCCGTCGATggtccctatc                                  >  3:598812/1‑142 (MQ=9)                                                                                 tACCAGGCTGCTCCACCCCGCGTCACAGTAGTCG‑TTTACTTCATCAAATTTTAATTGGTTGCACGGGCAGGATTAGAACCTCCGACCTTAGGGTTATGATCCAGACGAGCTACCAGGCTCCTCCAGACATAGTCAGTGGATGAGCCCtat                                   >  2:724794/1‑150 (MQ=9)                                                                                        ctgctCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTCTACTCCGcg                           >  6:761506/1‑150 (MQ=255)                                                                                          gctCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTcg                         >  6:651746/1‑151 (MQ=255)                                                                                            tCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTcgcg                       >  5:614614/1‑151 (MQ=255)                                                                                            tCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTcgcg                       >  3:991399/1‑151 (MQ=255)                                                                                            tCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCACGCTGCACCACCACGAGTCAGTGGATGCGCACTATACTCCGCTcgcg                       >  4:1042523/1‑151 (MQ=255)                                                                                            tCCACCCCGCGTCACAGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGACCGACGAGATACCAGGCTGCTCCACCCCGAGTCCGTGGATGCGCACTATACTCCGCTcgcg                       >  4:1037496/1‑151 (MQ=255)                                                                                             ccACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGt                      >  1:508438/1‑151 (MQ=255)                                                                                              cACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGtt                     >  5:550609/1‑151 (MQ=255)                                                                                              cACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGtt                     >  6:297769/1‑151 (MQ=255)                                                                                              cACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCAGCTCGCGtt                     >  2:905689/1‑151 (MQ=255)                                                                                               cccccGCGTCACCGTACTGC‑TTTACCTCAACAAAATTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGttt                    <  3:758078/150‑1 (MQ=255)                                                                                               aCCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGttt                    >  6:29935/1‑151 (MQ=255)                                                                                                ccccGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGGGGATGCGCACTATACTCCg                             >  3:372518/1‑141 (MQ=11)                                                                                                ccccGCGTCACCGTACTGC‑TTTTATTTATTAAATTTTTATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCg                             <  4:372380/141‑1 (MQ=255)                                                                                                     gctCACCGTACTGC‑TTTTCCTTAATAAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGc              <  3:916668/149‑1 (MQ=17)                                                                                                     cgTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGAAGCGCACTATACTCCGCTAGCGTTTTccttc              >  4:397668/1‑147 (MQ=14)                                                                                                     cgTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGCCCTTCGGGTTATGATCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGc              >  3:932506/1‑151 (MQ=255)                                                                                                      gTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGTAATCGAAGAATACTACGCTAGCGTTTTCATgg              >  6:715815/1‑149 (MQ=11)                                                                                                      gTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGc              >  2:887402/1‑150 (MQ=255)                                                                                                      gTCACCGTACTGC‑TTTACTTAATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACAGCCGACCTTCGGGTTATTAGCCCGACGAGCTACAAGGCTGCTCCACCCCCCGTCCCAGGATGCGCACTATACTCCGCTCGCGTTTTCATGc              >  7:362056/1‑150 (MQ=18)                                                                                                        cACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTACGCTAGCTTTTTCAttc              >  4:916022/1‑146 (MQ=14)                                                                                                          ccGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCt         >  5:467008/1‑151 (MQ=255)                                                                                                          ccGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCt         >  2:789058/1‑151 (MQ=255)                                                                                                           cGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCCGGGAAGCCAAGAATACTACGCTAGAGTTTTCATTCAACCtt        >  6:860821/1‑151 (MQ=11)                                                                                                           cGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtt        >  7:601429/1‑151 (MQ=255)                                                                                                            gTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGGGGAGGCGGACTATACTCCGCTCGATTTTTCATGCACCCgtt       >  7:451568/1‑151 (MQ=11)                                                                                                            gTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCCTGGATGCGCACTATAc                                 >  1:964467/1‑125 (MQ=11)                                                                                                            gTACTGCTTTTAATTAATAAAATTTT‑ATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATAc                                 <  2:963916/125‑1 (MQ=14)                                                                                                            gTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCttt       >  5:860341/1‑151 (MQ=255)                                                                                                            gTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCttt       >  5:1064386/1‑151 (MQ=255)                                                                                                             tACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGAGGATGCGAACTATACAACGCTCGCGTATTCATTCAACCttt       >  6:336975/1‑150 (MQ=14)                                                                                                             tACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACAGGATGCGGAATATACTGCGATAGATTTTTCATTCCACCtttt      >  6:570750/1‑151 (MQ=11)                                                                                                             tACTGCTTTTACTTCATCAAA‑TTTGATTCGTTGCGGCGCACGCATTTGAACCGACGACCTTCGGGTTATGACCCCGACGAGCTACCAGCCTGCTCCACCCCGCGTCCCAGGAAGCCAAGTATACTACACTATAGAATTCATGCAAccgtac      >  7:127707/1‑147 (MQ=9)                                                                                                              aCTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGACCCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCACCCttttt     >  6:9334/1‑151 (MQ=255)                                                                                                               cTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtttttt    >  5:696954/1‑151 (MQ=255)                                                                                                                 gC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtttttttt  >  6:941764/1‑151 (MQ=255)                                                                                                                 gC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCACGACGACCTACAAGGATTCTCCACCCAGCGTACGTGGATGCGCACTATACTCAGCTCGCGTTTTCATGCAACCtttttttt  >  8:899749/1‑151 (MQ=21)                                                                                                                   |                                                                                                                                                      CTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCTTTTTTTT  >  DH1/936165‑936427

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