Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,906,8370CG26.1% 140.1 / 18.4 92intergenic (+8/‑27)ECDH1_RS19315/ECDH1_RS19320tRNA‑Leu/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base C (44/21);  new base G (22/1);  total (70/22)
Fisher's exact test for biased strand distribution p-value = 9.93e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGTCCCCCCCCTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTTGATTT  >  DH1/3906694‑3906987
                                                                                                                                               |                                                                                                                                                       
cAAGTCACCCCCCTCGCTCCAAAAACCACGTTGATATTGGTCGCTCTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCa                                                                                                                                                  <  3:1562749/151‑1 (MQ=255)
actatcagcctccTCGCTGCAAGAAACTCGTTGTTATAGTTGACTCTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCa                                                                                                                                                  <  1:1140156/140‑1 (MQ=11)
 cacgacccccccTCGCACCAAACACCACGTTTATATTGCTCCCTCTCGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCaa                                                                                                                                                 <  8:244004/146‑1 (MQ=18)
 attgccccccccTCGCACCTAACACCTCGTTGATATTGCTCGCCCTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCaa                                                                                                                                                 <  2:320961/147‑1 (MQ=18)
     ccccccccTCGCACCAAAAACCTCGTTTATAGTCCTCCCAGTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCaaaaaa                                                                                                                                             <  7:1692084/151‑1 (MQ=21)
       ccccccTCGCACCAACCACATCGTTTAAATTGCTCTCTGTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTGCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAg                                                                                                                                           <  5:35954/151‑1 (MQ=18)
        cccccTCGCACGTAAAAACACGTTTATATTGCTCGCTCTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGt                                                                                                                                          <  1:3969/151‑1 (MQ=18)
         ccccTCGCACCAACAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTa                                                                                                                                         <  5:533060/151‑1 (MQ=255)
          cccTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCTATATCAAAAAAAGttt                                                                                                                                        >  5:520145/1‑149 (MQ=255)
           ccTCGCAGGAATAACGTCTTTGTTATTGCTCCGCCTCTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAg                                                                                                                                       <  3:233273/151‑1 (MQ=17)
           ccTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGTCGGTAATTCAAAAAAAGtttg                                                                                                                                       >  7:896734/1‑148 (MQ=18)
            cTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCCCACCAAAAAAATTGAAATCAAAAAAAGTatc                                                                                                                                       >  8:516863/1‑148 (MQ=18)
              cGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTACCCCCCCTCGCACCAAAAAACTAGAAATCAAAAAAAGtttgatc                                                                                                                                    >  7:1545752/1‑145 (MQ=18)
               gCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAATACGAGAATTCAAAAAAAGtttgttc                                                                                                                                    >  8:205802/1‑144 (MQ=18)
               gCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAAAGGGGATTTCAAAAAATGtttgttg                                                                                                                                    >  8:1580972/1‑144 (MQ=18)
               gCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTACTTACGGACGTGGGGGTTCAAGTCCCCCCCCAAGCACCAAACAACGCGATATAAAAAAAAAAAAGaaa                                                                                                                                    >  8:1246558/1‑148 (MQ=18)
                 aCCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGTGGAATTGGTAGACGCGCTAGCTTCAGGTCTTAGTGTACTTACGGACGAGGGGGTTCAAGTCCCCCCCCTCGCACCAAAATACGGGATATGAAATAAAGtttcttcagt                                                                                                                                 >  1:1184717/1‑142 (MQ=17)
                   cAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAATATGTCGATATCcaaaaaataaacactacagt                                                                                                                               >  3:1036913/1‑138 (MQ=18)
                    aaaaaaGACGTTTATATTGTTCCCTCTCTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACtgtg                                                                                                                              <  7:1365919/151‑1 (MQ=18)
                    aaaaaCCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGAGGCGAAATCAAATAAAGttacttgagtgtt                                                                                                                              >  3:1799495/1‑139 (MQ=18)
                      aaaCCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGAGGGGATATCAATTAATGTTAGTTGCCtgttgt                                                                                                                            >  7:283174/1‑148 (MQ=18)
                      aaaCCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGGTAGATCAAAAAAAGTatcatccgtgtgat                                                                                                                            >  4:441657/1‑138 (MQ=18)
                           aCGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAAACTAGGGGTGATCAAAAAAAGTTTGATGCCTGTGGGTAGGt                                                                                                                       >  3:403474/1‑151 (MQ=18)
                           aCGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCCGCCCCCAAAGAGGCGATATGAAAAAAAGTTAGATGACtgttggtaggt                                                                                                                       >  7:1379287/1‑143 (MQ=18)
                            ggTTGATATTGCTCGCTCTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTg                                                                                                                      <  7:761250/150‑1 (MQ=255)
                                tatagttctcgtccTGGGCGAATGTGGCGGAATTGGTAGTCGAGCTAGCTTCAGGGGTTATTGTTCTTACGGACGTGGGGGTTAAAGTCCCCCCCCTCGCACCAAACGAGGCTATATCAAAAAAGTAAGATGAGTGTGCGAAGGTGGCgg                                                                                                                  <  3:628482/137‑1 (MQ=18)
                                gATATTGCTCGCTCTTGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCg                                                                                                                   <  7:763555/150‑1 (MQ=255)
                                gATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAAAGGATAATTCAAAAAATGTAACATGAAtgatttttcatgcctg                                                                                                                  >  4:930648/1‑137 (MQ=14)
                                  tatttctctttgtttgCGAAGGTGTCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCACTCGCACCAAACGAGTCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGa                                                                                                                 <  4:264454/135‑1 (MQ=18)
                                   atTCCTCGATCTAGGCGACGGTGGCGGAATTGGTAGACGCGCTAGCGTCAGGTGGTAGTGTTCTTACGGACGTGGGGGTGCAAGTCCCCCCCCTCGCACCACCCGAGGCGATATCAAAAAAAGTAAGCTGACTGTGCGAAGGTGGCGGAAt                                                                                                               <  3:910187/151‑1 (MQ=18)
                                      gctgcgggggggCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTACTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTg                                                                                                             <  5:1100961/141‑1 (MQ=255)
                                      gctgcacgtgtgCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTTCTGTAATTACGGACGTCCCGGTTCAACCCCCCCCCCTCGCACCAACCGAGCCGATATCATATAAAGTAATATGACTCTGCGAAGGTGGCGGAATTg                                                                                                             <  1:1133084/139‑1 (MQ=17)
                                          caagttagCGAAGTGGTCGTAATTGGTACACGCGCTACCGGCAGGTGGTATTTTTCGTACGGACGATGCGGTTCAAGGCCCCCCCCTCGCACCACCCGAGTCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGa                                                                                                        <  3:873446/144‑1 (MQ=19)
                                            aCTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACTAGGGGAATTCAAATATAGTTAGATCCCTGTGGGTTGGTGGCGGATTTGCTTGAcg                                                                                                      >  1:969114/1‑151 (MQ=14)
                                                 gCGAAGGTGGCGGTTTTGTTAGAAGCTGTAGCTTCAGGTGGATGTGGTGTTATGTCACTGGTGTTGCGAGTCGCGCCCCTCGGACCTATCGTTCGGATTTCTTATAAAGTATTCTGTGTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTa                                                                                                 <  7:433952/151‑1 (MQ=255)
                                                   gAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAAATGAGGGGAGTTCCTAAAATGTTTGATGCGTGTTTTTTGGTGGCGGAATTGGTataagttgtatc                                                                                               >  3:870964/1‑140 (MQ=2)
                                                     aGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGGTATATCAAAAATAGTTAGTTGCCTGTGGGTTGGTGGCGGAATAggattcagttttagcat                                                                                             >  5:684742/1‑136 (MQ=17)
                                                            ggAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGGGATTTCCAAAAATGTTAGATGAATGAGGGTTGGTGGCGGAATTGGTTGACGAGGTAGCTTCAGgtgg                                                                                      >  7:451452/1‑150 (MQ=18)
                                                                   aatgtcgctCTAGCTTCTGCTGTTCCTCTTCTTACGGACTTGGGGGTTCAAGTCCCCCCCCTCGCACCAACCGACCCGATATCAAACATAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTc                                                                               <  1:1410628/146‑1 (MQ=25)
                                                                     aGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGGGAAATCAAAAAAAGTTAGATGAATGTTGTTTGGTGCCGGAATTGTAATAAGGGTTATCATCAGGTTTTAgtgaaa                                                                              >  7:667149/1‑147 (MQ=9)
                                                                          cgcTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAACTCCCCCCCCTCGCCCAAAACCACGCGATATAAAAAAAAATTATATTAAGGTGTTAAGCTGGCGGAATTCTTATACGGGGATGAATCAGGTTTTAgtgtaaaaaagg                                                                        >  7:1351866/1‑143 (MQ=9)
                                                                            cTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCACGCACCAAAAGACGCGAAATCAAAAAAAGTAAGATGCCTGTTGGAAGGTGGCGGAATTTGTAGAAGAGGAAGCTTCAGGTGTTAGTGTACtataggcc                                                                      >  2:615905/1‑144 (MQ=21)
                                                                               gCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAAAAGAGGGGAAATCAAAAAAAGTAAGATGCCTGTGGGATGGTGCCGGAATGGCTTTACGCGGTAGCTGCAGGTTTTAGTGTCATTATGGCCggg                                                                  >  7:540218/1‑149 (MQ=18)
                                                                                 ttCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAACAGGGGAAATCAAAAAAAGTAAGATGACTGTGGGAAGGTGGCGGAATTGTTATAAGCGCGAGCGGCATGTGTGAGTGTACTTACGGACGTggg                                                                >  7:1110307/1‑151 (MQ=21)
                                                                                   cAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTggggg                                                              >  3:626762/1‑151 (MQ=255)
                                                                                       tgtTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGGGATATCAAAAAAAGTTAGATGACTGTTCGATGGTGGCGGAATGGGTAGCAGGGGTAGCTTCATGTGTTACTGTACATACGGCCGTGGGGGttg                                                           >  5:248733/1‑149 (MQ=25)
                                                                                       tgtTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAATGTATGATGACTGTGGGTAGGTGGCGGAATTGGTTTCCGGGGTAGCTTCAGGTTGTAGTGTAGTTAGGGCCGTGGGGGttgc                                                          >  5:1520408/1‑149 (MQ=21)
                                                                                         ttAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGGTATATCAAAAAAAGTTTCATGAATGTGGTTAGGTGGCGGAATTGGTAGCAGAGGGAGCATCAGGTGTTAGGGTCATAACGGCCGTGCGGCttga                                                          >  1:1753637/1‑147 (MQ=12)
                                                                                         ttAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGCGATATCAAAAAAAGTATGATGCCTGTTTGTAGGTGGCGGAATTGGTTTACGGGGTATCTTCAGGTGTTAGTGTACTTACGGACGTGGGGGttgt                                                          >  7:1309278/1‑147 (MQ=25)
                                                                                         ttAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCTCGCCCCCAAAGATGGTATATGAAAAAAAGTAATATTACTTTTGTTATGTGGCGGAATTGGAATACGCTCTATCTACAGGTGTTACTGTAATTACGGACCTGCCCGTTCAAGt                                                       >  3:40632/1‑151 (MQ=255)
                                                                                           cctgtCCTTACGGACGTGGGGGGTCAAGTCCTTTTTTTATTTCTAATCGACGCGATATCATATAAAGTAAGATTACTCTTCTTTGGTTTCTGTACTGTTATACGATCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTc                                                      <  8:849892/149‑1 (MQ=2)
                                                                                             tgtTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAAAAGAGGGGATTTAAAATAAAGTAAAATGACTGTTATTAGGTGCCGGAATTGATATACTAAAAATCTTCAGGTATTAATGAAATTAAGGCCGTGGGGGtggaaggaa                                                     >  8:1639372/1‑142 (MQ=9)
                                                                                             tgtTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGCGATATCAAAAAAAGTTAGATGACTGTGGGTAGGTGGCGGAATTGGTATACGCGCGAGCTTCAGGTGGTACTGTACTTAAGGAAGTGGGGGTTCAAGtac                                                     >  5:921424/1‑148 (MQ=25)
                                                                                              gtTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGGGAGATCAAAAAAAGTTTGATGACTGTTTGAAGGTGGCGGAATTGGTTTACGAGGTAGCTTCAGGTGTTAGTGTACATTAGGCCGTCGGGGtggaaggcacc                                                   >  8:272956/1‑141 (MQ=19)
                                                                                              gtTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCGCACCAAAAGAGGGGATTTCAAAAAAAGTTTGATGAGTGTGGTTTGGTGGCGGATTTGGTTTAAGGTTTAGGTTCAGGTTTTCTTGTACTTAAGGCCggggggcgtggaggaca                                                    >  7:1112541/1‑139 (MQ=9)
                                                                                                      ggACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGCCTGTGCGATGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTGAAGGCCCCCCCCACGc                                           >  1:1263644/1‑151 (MQ=25)
                                                                                                      ggACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGAGATCAAAAAAAGTTAGATGCCTGTGGTATGGTGGCGGAATTCGTATAAGCGCGAGCTTCAGGTGTTAGTGTACTTACGGCCGTGGGGGTTCCAGGAAACCCCCTCGc                                           >  7:1475447/1‑151 (MQ=21)
                                                                                                      ggACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACCAGGGTATATCAAAAAATGTATGATGCCTGTGATAAGGTGGCGGAATTTGTTTCCTATGTAGCTTCATGTTGTAGTGTACAAAGCGACGTGGGGGTTGAAGtacaacggaacgc                                           >  7:1740802/1‑139 (MQ=12)
                                                                                                        aCGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATTTCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGGTAGCTTCAGGTGTTAGTGTACTTAAGGACGTGCGGGTTCAAGTCCCCCCCCTCGCAc                                         >  1:666089/1‑151 (MQ=255)
                                                                                                         cGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCGTACGGACGTGGGGGGTGACGTGCCCCCCCGCGcccc                                        >  7:890412/1‑151 (MQ=25)
                                                                                                          gTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGCCTGTGCGATGGTGGCGGAATTGGTATACGCGGTAGCTTCAGGTGTTAGTGTCCTTACGGCCGTGGGGGTTGAAGGCCCCCCCCGCGCAcc                                        >  5:874914/1‑150 (MQ=25)
                                                                                                                 ttCAAGTCCCCCCCCTCGCACCAAACGAGGCGAGATCAAAAAAAGTAAGTTGACTGTGAGATGGTGGCGGATTTGGTATACGCGGTAGCTTCAGGTGTTAGTGTCGTTACGGACGTGGGGGTTGAAGGCCCCGCCCGCGCCCCAAatatat                                >  8:270170/1‑149 (MQ=21)
                                                                                                                   cAAGTCCCCCCCCTCGCACCAAAAGAGGCGATATCAAAAAAAGTATGATGCGTGTTCTAAGGTGGCGGAATTGGTTGACGGGGTAGCTTCAGGTGTTAGTGTCGTTAGGGACGTGGGGGTTGAAGGGCACGCCCTAGGACCaaatgtattt                              >  2:322380/1‑143 (MQ=16)
                                                                                                                      gTCCCCCCCCTCGCACCAAAGGAGGGGAGATCATAAAAAGTATGATGAATGTGAGAAGGTGGCTGATTTGCTATACGAGATAGCTTCAGGAGGTAGTCTAATTAAGGACGTCGGGGTTCAAGGCACAACCCTGGCCCCAACTTTCTTTaat                           >  6:1299129/1‑149 (MQ=255)
                                                                                                                        ccccccccTCGCACCAAACGAGGGGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCGTTACGGACGTGGGGGTTGAAGGCCCCCCGCGGGCCCCGAGTATGTTttcttc                         >  4:995234/1‑151 (MQ=18)
                                                                                                                        ccccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTGGCCCCAATTATGTGttcttc                         >  1:1742167/1‑151 (MQ=255)
                                                                                                                         cccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTTTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGGCACCCCCCTCGCACCAATTATCTTTACTTcc                        >  5:1557098/1‑151 (MQ=255)
                                                                                                                         cccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAGTGCCCCCCCTCGCACCCATTATCTTTACTTcc                        >  5:1235150/1‑151 (MQ=255)
                                                                                                                         cccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGCACGTGGGGGTTCAAGTCCCCCCCCTCGGACCAATTATCTTTACTTcc                        >  3:864130/1‑151 (MQ=255)
                                                                                                                           cccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTccgt                      >  8:273277/1‑149 (MQ=255)
                                                                                                                           cccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGGCCCCCCCCGCGCCCCAATTATCTTTTCTTGCtt                      >  5:1665809/1‑151 (MQ=255)
                                                                                                                            ccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCttt                     >  8:1101629/1‑151 (MQ=255)
                                                                                                                            ccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAGTCCCCCCCCGCGCACCAATTATGTTTACTTCCtgt                     >  7:449952/1‑149 (MQ=255)
                                                                                                                             cccTCGCACCACACGCGGCGATATCATTATAACTAATATGCCTCTGCCCCTGTGGACCAATTTGTATACCAGCTTTCGTCACGTGTTCGCGCCTTCCAGCTCGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTc                    <  2:108942/151‑1 (MQ=255)
                                                                                                                             cccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTGAAGTGCCCCCCCCCGGCCCCAATATCTTTTCttacgttg                    >  7:1055727/1‑145 (MQ=21)
                                                                                                                              ccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCt                   >  2:84492/1‑151 (MQ=255)
                                                                                                                                tCGCCCCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGGGTTAGTGTCCTTACGGACGGGGGGGTTCCCCGCCCCCCCCGGGGCCGCATTATCGTTACTTGCTGTCtgg                 >  6:1655635/1‑149 (MQ=18)
                                                                                                                                tCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCttt                 >  5:233866/1‑150 (MQ=255)
                                                                                                                                 cGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGt                >  8:1129973/1‑151 (MQ=255)
                                                                                                                                  gCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCCAGTCCCCCCCCTCGCGCCAATTATCTTTAGTTGCTTTgttgtt               >  6:1551671/1‑151 (MQ=255)
                                                                                                                                  gCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAATTATCTTTTCTTCCGTTCTTGtt               >  3:1777999/1‑151 (MQ=255)
                                                                                                                                  gCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGtt               >  8:1591133/1‑151 (MQ=255)
                                                                                                                                    aCCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTGCCCCCCCGCGCACCAATTATCGTTAATTACTTTGTTGtttg             >  5:1629558/1‑150 (MQ=21)
                                                                                                                                       aaaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTttcttc          >  2:1640815/1‑151 (MQ=255)
                                                                                                                                          cGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCtt       >  1:1445289/1‑151 (MQ=255)
                                                                                                                                           gAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGACGGCCCCCCCCGCGCCCCAATTATCTTTACTTCCGTTCTTGGTTCTTCCTTg      >  1:469488/1‑151 (MQ=21)
                                                                                                                                           gAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTGCCCCCCCGCGCCCCAATTATGTTTACTTACTTTCTTGTTTCTTCCttt      >  5:700682/1‑150 (MQ=255)
                                                                                                                                           gAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCttt      >  2:456401/1‑150 (MQ=255)
                                                                                                                                           gAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCACGCACCAATTATCTTTACTTCCTTAGTTGGTTCTTCCTTg      >  1:496790/1‑151 (MQ=255)
                                                                                                                                            aGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTAATTACTTTCTTGTTTGTTCCTTg      >  6:92348/1‑150 (MQ=255)
                                                                                                                                            aGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTAATTACTTTCTTGTTTATTCCTTg      >  2:505534/1‑150 (MQ=255)
                                                                                                                                             ggCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCCACTGGCCCCCCTCGCCCCAATTATCTGTAGTTCCGTTCTTGTTTGTTCCTTGAt    >  5:372214/1‑151 (MQ=21)
                                                                                                                                               cGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTGCCCCCCCGCGCCCCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTTGAttt  >  5:1671783/1‑151 (MQ=255)
                                                                                                                                               cGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCCATTATCTTTAATTCCTTTCTTGTTTCTTCCTTGAttt  >  2:1782153/1‑151 (MQ=255)
                                                                                                                                               cGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCACGCACCAATTATCTTTACTTCCTTTCTTTGTTGTTCCTTGAttt  >  6:331865/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                                       
CAAGTCCCCCCCCTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTTGATTT  >  DH1/3906694‑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: