Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	DH1	3,906,837	0	C	G	27.0%	106.1 / 10.8	70	intergenic (+8/‑27)	ECDH1_RS19315/ECDH1_RS19320	tRNA‑Leu/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base C (35/11);  new base G (16/1);  total (57/13)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCAAGTCCCCCCCCTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTTGATTT  >  DH1/3906692‑3906987                                                                                                                                                  |                                                                                                                                                       gtaaaggacccccccTCGCAGGAATAACGTCGTTTATATTTCTCACAGTCGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGatat                                                                                                                                                   <  1:352455/143‑1 (MQ=17)      gTCCCCCCCCTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTGTTGCGGACGGGGGGGTTCAAGAGGCCCCCATCGGCCGCAGCTAGGCGATATCata                                                                                                                                               >  3:316753/1‑148 (MQ=21)          cccgccTCGCTCCAAAAACTACTTTTATAGTCGTCACTGTTTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTGCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAg                                                                                                                                          <  7:1164306/151‑1 (MQ=14)           cccccTCGCACCAAAAACCACGTTGATATTGCTAGGAGTGTGCGAAGGTGGAGGAATTGGTAGCAGCTATAGCTTCAGGTGTTAGTGTGCTTAAGGACGTGGGGGTTAAAGTCCCACCCCAAGCACAAAAAGAGACTATTTGAAATAAAGt                                                                                                                                         >  6:691853/1‑151 (MQ=255)            cgccTCGCACCTACAACCTCGTTTATATTTCTCGGTGTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTa                                                                                                                                        <  5:622831/149‑1 (MQ=18)              ccTCGCACCAACCACCACGTTGATATTGCTGGGTGTCGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAg                                                                                                                                      <  7:1129646/151‑1 (MQ=18)                tCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTGAACTCCCCCCCCGCGCCCCAAACGAGGGGAGTTCAAAAAAAGttaaat                                                                                                                                    >  7:59121/1‑146 (MQ=21)                  gCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAAAGCAGGATATATCaattaatgtttgttgc                                                                                                                                  >  5:1060977/1‑137 (MQ=18)                     ccAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGAAGTATAAATCaaaaaatattacttgactg                                                                                                                               >  2:946234/1‑138 (MQ=17)                      cAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAAAAGAGGTTATTTCAAAAAAAGTTACATgcttgt                                                                                                                              >  1:256100/1‑146 (MQ=18)                      cAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAATGTGGTTATATCAATAATTGTTAGATGACtgt                                                                                                                              >  1:1077727/1‑151 (MQ=18)                      cAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAAACGATAAATCAAAAAATGTAAGAtccatgt                                                                                                                              >  8:1291224/1‑145 (MQ=18)                       aaaaaCCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGGGAAATCAAATAAAGTTTGATGCCtgtt                                                                                                                             >  5:411934/1‑150 (MQ=18)                            ccACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGGGAGTTCCAATAAAGTTAGATGCGtgtgtgt                                                                                                                          >  3:597822/1‑146 (MQ=17)                            ccACGTTGATAATTCTCTCACTTGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAg                                                                                                                        <  1:997289/151‑1 (MQ=255)                            attcttttatatTGCTCGCCGTTGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAg                                                                                                                        <  3:1047855/143‑1 (MQ=21)                             cACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCGCACCAAAACAGGGGATATCAAATAAAGTTTGATgcatgttggtagt                                                                                                                       >  7:7696/1‑139 (MQ=17)                                      attcctctccctctgCGAAGGTGGCGGAATTGGTATTCGTTCTAGCTTCAGGTGGTAGTGTTCTTAAGGACTTGGGGGTTCATGTTCCCCCCCTCGCACCAAACGAGTAGATATCAATTAAAGTAAGCTGAGTTTGCGAAGGTGGCGGAAt                                                                                                              <  2:212422/137‑1 (MQ=12)                                      atTGCTGGCACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGaa                                                                                                               <  5:368839/150‑1 (MQ=255)                                           tCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACAAGGGGAGATCAATAAATGTATGATCACTGTGATATGCTGGCGGattctcat                                                                                                         >  3:76411/1‑144 (MQ=17)                                            ttatgtttgggAATGTGGCGGTATTTTTAGACTGTGTAGCTTTAGTTGTTTGTGTTCTTAGGGACGGGGCGGTGAAAGTCCCTCCCCTCGCCCTAAACGCGTGGATATCTATTAAAGTAATCTTTGTTTGCGAAGGTGGCGGAATTGGTa                                                                                                         <  2:759059/140‑1 (MQ=255)                                            gatggtatgCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGat                                                                                                                                                     <  1:999688/98‑1 (MQ=255)                                             gctcTGCGCGTAGGTGCCTGAATTGGTACAGGGGCATGCTCCAGGCGGTACTTTTACTACGGACGATGCGGTGCAAGGCACCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAATCTGCCTGTGCGAAGGTGGCGGAATTGGTAGa                                                                                                       <  5:14536/148‑1 (MQ=255)                                             gCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGAc                                                                                                                                 <  4:1209241/125‑1 (MQ=255)                                             gCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCACGCACCAAAAGAGGCGATATCAAAAAAAGTAAGATGAc                                                                                                                                 >  3:1211380/1‑125 (MQ=255)                                               ccTGGGCGAAGGGGGCGTATTTGGTAGACGCGTTCGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGagg                                                                                                     <  7:466047/150‑3 (MQ=25)                                                     cGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGGGATTTCCAAAAATGTTTGATGCCTGTTGTTTGGTGCCGGAATTGATatcctgactat                                                                                               >  1:1169156/1‑141 (MQ=11)                                                               ggAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGGTATATCAAAAAAAGTTTGTTGCCTGTGGGTTGGTGGCGGAATTGATATAATCTCCAGCATCAGgtgg                                                                                     >  3:909064/1‑150 (MQ=14)                                                                        aGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGGGAAGGTGGCGGAATTGCTAGACTCGCTAGCATCAGGAGTTAGTGTACt                                                                            >  5:387347/1‑151 (MQ=25)                                                                         gACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGGGATATCAAAAAAAGTAAGATGACTGTTGTAAGGTGCCGGATTATCTAGAATAGCAACCTTCAGGTGTTAgtcaacat                                                                           >  1:868824/1‑145 (MQ=17)                                                                           cgcgCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGTCGATATCAAAAAAAGTAACATGACTGTGAGAAGGTCGCGGAATGGGTATACGAGATATCATCAGGAGTAAgtgaaatta                                                                          >  6:378026/1‑144 (MQ=18)                                                                             cgcTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGCGAAATCAAAAAAAGTTAGATGAGTGTGGGTTGGTGGCTGAATTGGTTGCCGAGGGAGCTTCAGGTGGGCGTGTCATTAggg                                                                       >  2:938488/1‑151 (MQ=18)                                                                                    ttCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGGGATATCAAAAAAAGTAAGATGACTGTGAGATGGTGGCGGAATTTGTATACGAGATATCATCAGGTTTGAGTGTCAATAAGGACGTgcg                                                                >  3:222833/1‑149 (MQ=21)                                                                                    ttCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAAGAGAGATCAAAAAAAGTAACATGACTGTTGGAAGCTGGCGGAATACGTAGAAGCGATATCATCAGGTTTTACTGTACtataggacctgtg                                                                >  1:382534/1‑139 (MQ=12)                                                                                      cAGATGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACCAGGGGAGTTCAAAAAAAGTATGATGCCTGTTCGAACGTCGCGGAATTCCTATACGATGGAGCTGCAGGATTTACTGTCCTTCAGGCCGTgccgg                                                              >  1:26288/1‑147 (MQ=12)                                                                                       aGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTGCCCCCCCGCCCCCCAAAGGGGGCGAGATGAATAAATGTAACTTGACTGTGGTTAGTTGGCGGGTTTTGTATACGCTCGATCATCATGTGTTACTGAACATAAGGACGGGGGGGt                                                             >  1:811870/1‑151 (MQ=9)                                                                                          tgtTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTTAGATGCCTGTTGGTTGGTGGCGGAATTGGTTTACGGGGTATCTTCATGTTTTAGTGTAATTACGGCCGTGGGGGttgc                                                          >  1:465284/1‑149 (MQ=21)                                                                                          tgtTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCGCACCAAACGTGGGGAGATCAAAAAAAGTTTGTTGCCTGTGGGTTGGTGGCGGATTTGGTAGACTGGGTATCGTCAGGTTGTAGTGTAGTTACGGCCGTGGGGGttgt                                                          >  5:988595/1‑149 (MQ=16)                                                                                            ttAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTTAGATGCCTGTGCGATGGTGGCGGAATTGGTTGACGCGGTAGCTTCAGGTGGTAGTGTCCTTACGGCCGTGGGGGttga                                                          >  3:852427/1‑147 (MQ=37)                                                                                             tAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGGGAGATCAAATAAAGTTAGATGCCTGTGGGTAGGTGGCGGAATAGGTAGCCGTGGGAGCATCAGGTTTTAGTGTACTTAAGGCCGGGGGGGTTGAAgg                                                       >  3:1012101/1‑150 (MQ=19)                                                                                                    cTTACGGACGTGGGGGTTCAAGTCCCCCCCTTCGCACCAAAGGAGGCGATATCAAAAAATGTTAGATGACTGTGGGTAGGTGGCGGATTTGGTATACGAGGTAGCTGCAGGTTTTAGTGTCGTAACGGCCGTGGGGGTTGAAggcacaccc                                                >  7:1277619/1‑143 (MQ=21)                                                                                                     ttACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGTGGGGATATCAAAAAAAGTATGATGAATGTGGTTTGGTGGCGGATTTGGTAGACGGGGGAGCTGGAGGGGGGAGTGTACTTGGGGCCGTGGGGGtggtaggggcccggc                                               >  6:331955/1‑137 (MQ=12)                                                                                                       aCGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAAAAGAGGGTATATGAAAAAATGTTAGATGACTTTTGTTAGGTGGCGGAATTGGTATAAGGGTGAGCTTCATGTTGTAATGTAAAAAAGGACGTGGGGGTTGAAggaaacacacg                                              >  8:125122/1‑140 (MQ=12)                                                                                                            cGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGGTATTTCAAAAATTGTAAGTTGAATGTGGGATGCTGGCGGATTTGGTTGACGGGGGAGCTGCAGGTTTTAGTGTCCATACGGCAGTGGGGGTTGAAGTACACCCCCAGGCAcc                                        >  3:134367/1‑151 (MQ=255)                                                                                                            cGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCTATATCAAAAAAAGTTAGATGAGTGTGGGAAGGTGGCGGAATTGGAATACGCTCAAGCTTCAGGTGTTAGTGTAATTACGGACGTGGGGGTTGTAGTCAAACCACTCGCAcc                                        >  5:112200/1‑151 (MQ=21)                                                                                                               gggggTTCAAGTCCCCCCCCTCGCACCAAAAGAGGCGATATCAAAAAAAGTAAGATGACTGTGAGAAGGTGGCGGAATTGGTAGACGAGGTAGCTTCAGGTGTTAGTGTCCTTACGGCCGTGGGGGTTGAAGTACCACCCCACGccccca                                      >  4:635529/1‑148 (MQ=25)                                                                                                                     tCAAGTCCCCCCCCTCGCACCAAACGAGGGGATATCAAAAAAAGTAAGATGAATGTGCTAAGGTCGCGGATTTGGTAGACTAGATAGCTTCAGGTGTTAGTGTCGTTAGGGACGTGGGGGGTCAAGTCCCCCCCCTAGCCCCaaaattggt                               >  6:148788/1‑144 (MQ=21)                                                                                                                        aGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTGGAAGTCCCCCCCCTTGCCCCaactttttgttc                            >  7:909083/1‑141 (MQ=21)                                                                                                                              cccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAATTATCTTTACTTccgt                      >  1:604104/1‑149 (MQ=255)                                                                                                                              cccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGAATGTGCGAAGGTGGCGGAATTGGAAGACTCGATAGCTTCAGGTGTTAGTGAACTAACGCACGTGGGGGTTAAACACACCCCCCTCCCAACAAATaaaattaattaat                       >  8:235357/1‑138 (MQ=255)                                                                                                                               ccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTCCGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACGTCCttt                     >  4:57938/1‑151 (MQ=255)                                                                                                                               ccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGTCCCCCCCCTCGCACCCATTATCTTTACTTCCgtt                     >  8:1098078/1‑151 (MQ=255)                                                                                                                               ccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACGTCCttt                     >  6:389249/1‑151 (MQ=255)                                                                                                                                cccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTGAAGTGCCCCCCCGCGCACCAAATATATTTAcatacgttt                    >  2:1277663/1‑143 (MQ=21)                                                                                                                                   tCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTCGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCATTCTTg                 >  6:603936/1‑151 (MQ=255)                                                                                                                                     gCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTGCCCCCCCGCGCCCCAATTATCTTTAAGTGCTTTgttgtt               >  2:384217/1‑151 (MQ=21)                                                                                                                                     gCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAAGTATCTTTACTTCCTTTgttgtt               >  5:554051/1‑151 (MQ=255)                                                                                                                                     gCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGtt               >  1:541305/1‑151 (MQ=255)                                                                                                                                     gCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAGTCCCCCCCCGCGCCCCAATTATCttttggtgcgttgttgtt               >  8:839856/1‑136 (MQ=21)                                                                                                                                      cACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTACCATTCGTGttt              >  6:1287810/1‑151 (MQ=255)                                                                                                                                      cACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAACTGCCCCCCCGCGCCCCAATTATCTTTTCTTCCTTTCTTGttt              >  5:323436/1‑151 (MQ=255)                                                                                                                                         cAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCGCACCAATTATCTTTACTTCCTTTCTTGTttctt           >  6:393694/1‑151 (MQ=255)                                                                                                                                          aaaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCGCCCTCGCACCAATTATGTTTACTTCCTTTGTTGtttagtc          >  4:1254538/1‑147 (MQ=255)                                                                                                                                            aCGAGGCGATAGCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGTGCCCCCCCGCGCCCCAATTATGTTTACGTCCGTTCTTTTTTGTTGCt        >  5:924899/1‑151 (MQ=18)                                                                                                                                             cGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCtt       >  5:1176331/1‑151 (MQ=255)                                                                                                                                             cGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTGAAGTGCCCCCCCGCGCCCCAATTATCTTTTATTGCTTTGTTGTTTGTTGCtt       >  7:1120334/1‑151 (MQ=18)                                                                                                                                              gAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCttt      >  3:275711/1‑150 (MQ=255)                                                                                                                                              gAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCttt      >  1:116115/1‑150 (MQ=255)                                                                                                                                               aGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTGAAGTGCCCCCCCGCGCCCCCATTATCTTTTATTCCTTTGTtgtttgttggttgt     >  7:38287/1‑145 (MQ=18)                                                                                                                                                ggCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACAAATTAACTTTACTTCCTTTCTTGTTTATTCCTTGAt    >  2:1100036/1‑151 (MQ=255)                                                                                                                                                  cGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGGTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACGAATATGGTTTACTTCCTTGGTTGTTGGTACGTTGAttt  >  6:1079170/1‑151 (MQ=17)                                                                                                                                                  |                                                                                                                                                       TTCAAGTCCCCCCCCTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTTGATTT  >  DH1/3906692‑3906987

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