Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH1936,2781.T24.1% 142.1 / 68.7 58intergenic (‑12/+22)ECDH1_RS04620/ECDH1_RS04625tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base . (39/5);  new base T (11/3);  total (50/8)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-02

TGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCTTTTTTTT  >  DH1/936190‑936427
                                                                                          |                                                                                                                                                      
tGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTAc                                                                                                                                                  <  1:209999/94‑1 (MQ=31)
tGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTAc                                                                                                                                                  >  2:210159/1‑94 (MQ=31)
  gTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCCCGGAGTGCGTGAATTCTGCACTATATAATTGGTTGCGGGGGGAGTTTTTGACCCGCCGCCATTCGGGttctt                                                                                       >  5:619944/1‑148 (MQ=2)
                  ttGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCCAGTAATCGATGTATTCTTCGATGTATTTTTCATTGCGCGGGAGGGATTAGCACCGACGCCCTTAGGGTTATTAGTCCGACGAGCTaa                                                                        >  8:1205133/1‑149 (MQ=2)
                              ccTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCAGCTATGTTTAATTGGTTGCGGCGGCCGGTTTTGAACCGACGCCCTTCGGGTTATGATCCCGCCGAGAGACCAGGCTGCTcccc                                                           >  7:81622/1‑149 (MQ=11)
                                                gacgaGCTACCAGGCTGCTCCACCCCGCGTCCCGGGAGTCCGTGACTTCTGCTATGTTTATTTGGTTGCGCGGGGGGGATTTGACCCGCCCCGCTTGGGGTTATGATCCCGACGAGCTACAAGGCTGCTCCCCCGCGCGTcgggcgatg                                           >  1:111190/1‑141 (MQ=2)
                                                       tACCAGGCTGCTCCACCCCGCGTCACCGTAGTCCATAAATTCATCAAATTTTAATTGGTTGCGGCGTACTGATATGAACCGACGACCTTCGGCTTATGATCCCGCAGAGCAACAAGGCTGCTCCAAAAATAGTCAGTGGATGCTCACtat                                   >  4:1355039/1‑150 (MQ=9)
                                                         ccAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACt                                >  3:432076/1‑151 (MQ=255)
                                                             gctgctCCACCCCGCGTCACCGTAATGCGTTACTTCATCAAATTTTAATTGGTTGCGGGGGCAGGATTTGACCCGACGACATTCGGGTTATGAGCCCGAAGAGCTACCAGGCTGCTCCACCCCGCGTCAGTGGATTCGCACTATACTCCGc                            >  4:281375/1‑151 (MQ=17)
                                                               tgctCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACCAGCTACCAGCCTGCTCCACACCGAGTCCGTGGATGCGCACTATACTCCGCTc                          >  4:979352/1‑151 (MQ=255)
                                                                gctCCACCCCGCGTCACCGTAGTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCAGGATTTGAACCGACCACCTTCGGGTTATTATACCGACGAGCTACCAGGCTGCTCCCCACGTAGTCGGTTGATGAGCACTATACTCCGCt                           >  4:1162486/1‑149 (MQ=17)
                                                                 ctCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTcgc                        >  5:892979/1‑151 (MQ=255)
                                                                  tCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGACGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACAAGGCTGATCCACCCCGCGTCAGCGGATGCGCACAATTCTGCCCTcg                         >  1:875796/1‑149 (MQ=11)
                                                                  tCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGCCCTTCGGGTTATGAGACCGACGCGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTcgcg                       >  2:1285667/1‑151 (MQ=255)
                                                                   ccACCCCGCGTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGt                      >  2:416946/1‑151 (MQ=255)
                                                                    cACCCCGCGTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGAAGAGCTACCAGGCTGCTCCACACCGCGTCAGTGGATGCGCACTATACTCCGCTCGCGtt                     >  1:613916/1‑151 (MQ=255)
                                                                      ccccctctCACCGTAGTGCTTTACCTCCTGAAAATTTTATTGGTTGCGGGGGCAGGAGTTGAACCGACGACCTTCGGGTTATGAGCACGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTTGATGCTCACTATACTCCGCTGGCGtt                     <  3:1341430/142‑1 (MQ=17)
                                                                      ccccGCGTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACAGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGtttt                   >  6:389940/1‑151 (MQ=255)
                                                                       cccGCGTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTc                  >  6:1289240/1‑151 (MQ=255)
                                                                         cgcgTCCCCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGACCCGCCGCCCTCCGGGTTAGGAGCCCGCCGCGCTCCCAGGCGGAGCCCCCCCGCGGCCGGGGATGCGCCCTATCCTCCGCTCGCGtttgc                  >  1:517811/1‑147 (MQ=18)
                                                                         cgcgTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGAGGAATCGCAGTATACTACTCACGCGTTTTCAt                >  1:17268/1‑151 (MQ=1)
                                                                         cgcgTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGCCGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCAt                >  2:1249923/1‑151 (MQ=255)
                                                                         cgcgTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCAt                >  7:1257989/1‑151 (MQ=255)
                                                                          gcttCACGGTACTGCTTTTACTTCATAAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATg               <  2:17201/148‑1 (MQ=14)
                                                                          gcgTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATTAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACAGGAATCCCAGTATACTccaaaagatatttaatt               >  6:567751/1‑136 (MQ=14)
                                                                          gcgTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCCCt                           >  1:49065/1‑139 (MQ=11)
                                                                          gcgGCACCGTACTGCTTTAAATGCATCAAATTTTATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCt                           <  2:48900/139‑1 (MQ=14)
                                                                            gTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCCTGGAGGCGCACAATTCTACCCTAGAGTTTTCAttc              >  6:363678/1‑148 (MQ=11)
                                                                            gTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGc              >  4:136888/1‑150 (MQ=255)
                                                                             tCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACtat                                   >  3:560528/1‑128 (MQ=11)
                                                                             tCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGAAGCGGACTATACAACGCACGAGTTttcattca             >  4:1211902/1‑150 (MQ=11)
                                                                             tCACCGTACTGCTTTTACTTCATAAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACtat                                   <  4:560530/128‑1 (MQ=11)
                                                                             tCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCa             >  1:688640/1‑150 (MQ=255)
                                                                             tCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCa             >  8:363850/1‑150 (MQ=255)
                                                                             tCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCa             >  5:989183/1‑150 (MQ=255)
                                                                              cACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCCTGGAAGCGAAATATACAACGATCGATATTTCATTCAAc           >  6:1203914/1‑151 (MQ=11)
                                                                              cACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAAc           >  8:890561/1‑151 (MQ=255)
                                                                                ccGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGCGCTACCAGGCTGCTCCACCCCGCGTCCGTGGAGGCGTACTATACACCGCTGGCGTTTTccttgcggcg         >  8:1080697/1‑142 (MQ=11)
                                                                                ccGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCt         >  2:1360667/1‑151 (MQ=255)
                                                                                ccGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCt         >  5:1355973/1‑151 (MQ=255)
                                                                                ccGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCt         >  4:587960/1‑151 (MQ=255)
                                                                                 cGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGGGGATTCGGACTATACTCCGATAGATTTTTCAt                >  7:1236732/1‑143 (MQ=1)
                                                                                 cGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGGGGATGCGCACTATACTCCGCTCGCGTTTTCATggcgccgt        >  5:433251/1‑144 (MQ=14)
                                                                                 cGTACTGCTTTTACTTAATAAATTTTTATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCAt                <  8:1236059/143‑1 (MQ=17)
                                                                                 cGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtt        >  6:915438/1‑151 (MQ=255)
                                                                                 cGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtt        >  7:630878/1‑151 (MQ=255)
                                                                                  gTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCttt       >  5:1344208/1‑151 (MQ=255)
                                                                                   tACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtttt      >  6:742809/1‑151 (MQ=255)
                                                                                   tACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtttt      >  1:1070271/1‑151 (MQ=255)
                                                                                   tACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGCATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtttt      >  5:1111601/1‑151 (MQ=255)
                                                                                    aCTGCTTTTAATAATAAAATTTTTATTGGTTGCGGTGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCttttt     <  5:568273/151‑1 (MQ=21)
                                                                                     cTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtttttt    >  6:1351297/1‑151 (MQ=255)
                                                                                     cTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtttttt    >  8:403161/1‑151 (MQ=255)
                                                                                      tGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTAATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCCTCCGTGGATGCGCACTATACTCCGCTCGCCTTTTCATGCAACCtttttgt   >  4:601768/1‑149 (MQ=255)
                                                                                      tGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGAACTTCGGGTTATGAGCCCGACCAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGGTTTCATGCAACCttttttt   >  8:662269/1‑151 (MQ=255)
                                                                                        gCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtttttttt  >  8:1329139/1‑151 (MQ=255)
                                                                                        gCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtttttttt  >  3:45453/1‑151 (MQ=255)
                                                                                        gCTTTACCTAATTAAATTTTTATTGGTTGCGGGGGCCGGATTTGAACCGACGACGTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCttttttt   <  3:155854/150‑1 (MQ=255)
                                                                                         cTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCtttttttt  >  7:496599/1‑150 (MQ=255)
                                                                                          |                                                                                                                                                      
TGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGTGGATGCGCACTATACTCCGCTCGCGTTTTCATGCAACCTTTTTTTT  >  DH1/936190‑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: