Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,906,8320CA36.1% 108.7 / 18.3 92intergenic (+3/‑32)ECDH1_RS19315/ECDH1_RS19320tRNA‑Leu/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base C (37/16);  new base A (30/0);  total (76/16)
Fisher's exact test for biased strand distribution p-value = 3.45e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTTCAAGTCCCCCCCCTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTT  >  DH1/3906689‑3906982
                                                                                                                                                |                                                                                                                                                      
gtgTTCAAGTCCCCCCCCTCGCACCAAAAACCACTTTGATATTGCTCTCTCTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGa                                                                                                                                                 <  8:54036/149‑1 (MQ=255)
     cAAGTCCCGCGCCTCGGTCCAAACACCTCGTTGATATTGCTGTCTGTCGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATATCa                                                                                                                                            <  4:2050165/151‑1 (MQ=17)
     cAAGTCCCGCCCCTCGCAGCAATAACGTCGTTGTAATTTGTCTCAGTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATATCa                                                                                                                                            <  7:385767/151‑1 (MQ=17)
         gcaccccccTCGCAGCTACAACCTCGTTTTTATTGCTGCGAGTTGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATATCaaaaa                                                                                                                                        <  1:1044235/148‑1 (MQ=17)
           cccccccTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGATGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACCGAAGTGGGGGTTCAACT‑CCCCCCCCACGCCAAAAACGCCGCGATATGaaaaaa                                                                                                                                       >  3:1018952/1‑150 (MQ=21)
            cgcgccTCGCAGCAAACACCTCGTTTATATTGCTCCCTGTTTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAg                                                                                                                                     <  5:852275/147‑1 (MQ=17)
              ccccTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACTAGGCTAATTCAAATAAAGtt                                                                                                                                   >  5:1956755/1‑150 (MQ=255)
               cccTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAACT‑CCCCCCCCACGCCCCAAAGCAATAGATATAAAAAAAAGttt                                                                                                                                  >  1:1008533/1‑149 (MQ=18)
                ccTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCGTCATGTGGTAGTGTTCTTACGGCAGTGGGGGTTCAAGT‑CCCCCCCCTGGCAAAAAACGAGGCGATATCAAAAAAAGTa                                                                                                                                   >  1:1056573/1‑149 (MQ=25)
                   cGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAGAACGATAGATCAATAATAGTTAGATg                                                                                                                              >  4:1225143/1‑151 (MQ=18)
                    gCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAAGACGTTAAATCAAAAAATGTTAGAtcc                                                                                                                             >  1:102988/1‑149 (MQ=18)
                     cACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTGAAGT‑CCCCCCCCTCGCCCCAAACCACGGGATATCAAATAAAGTTTGATgcc                                                                                                                            >  1:1689348/1‑149 (MQ=18)
                      acgaataaagtcGTTTATATTGCTGTCTCTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACt                                                                                                                           <  3:1375148/140‑1 (MQ=18)
                      aCCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCCCGCACCAAAATAGTTAATATCaaaatttgtaattagtt                                                                                                                            >  3:1909560/1‑138 (MQ=18)
                      aCCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAGCACGTGATATCCATAATTGTTAGATGACt                                                                                                                           >  5:199504/1‑151 (MQ=18)
                      aCCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAAAAGTCGAAATCAAAAAAAGTAAGATgcat                                                                                                                           >  7:1414268/1‑148 (MQ=255)
                             acctcttttatatTGCTGGGTGTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGaa                                                                                                                    <  1:4569/142‑1 (MQ=18)
                             aCCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAATATAGGGGAATTCAATTAAAGTTTGATGACtttttttt                                                                                                                    >  5:36972/1‑144 (MQ=17)
                               cACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCCCCAAAAAACTGTATATCATATAAAGTAAGTTGCCtgttatttgc                                                                                                                  >  3:291578/1‑144 (MQ=17)
                               cACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAAAAGTAGAAATCAATTAATGTAAGAtcaatgatttt                                                                                                                     >  6:1691395/1‑138 (MQ=17)
                                     gATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAGTAGGGTATTTCAAATAATGTTTGATGACtgttttttgctcgcgg                                                                                                            >  3:1342986/1‑138 (MQ=17)
                                       tattcatggaagggggCGCAGGTGGGGGCTTTGGGCGTCGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGa                                                                                                           <  5:1412893/137‑1 (MQ=21)
                                       tatTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAAAAGGGTAAATCAAATAAAGTTAGATGCATGTTGTAAGGTGCCGGa                                                                                                           >  6:1827650/1‑150 (MQ=17)
                                        attattgtttgtttgCGAATGTTTCGGAATTGTTAGACGCTCTAGCTTCAGGTGTTATTTTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACGAAACGAGGCGATATCAAAAAAAGTAAGATGACTTTGCGAAGGTGGCGGaa                                                                                                          <  6:111886/136‑1 (MQ=18)
                                        atatatgatagttttcGAAGGTTGCGGAATTGGTAGACCTGCTATCTTGAGGTGTTAGTTTTCTTACTGACGTTGTGGTGAAAGT‑CCCCCGACTCGCACCAAACGAGCCGATATCAAAAAAAGTAAGCTGACTGTGCGAAGGTGGCGGAAt                                                                                                         <  6:1622670/136‑1 (MQ=255)
                                        atTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCACGCACCAAAATAGGGGAAATCCAAAAATGTTTGATGCATGTGTGTTGGTGGCGGatt                                                                                                         >  5:448941/1‑149 (MQ=14)
                                        atTGCGCGCACTGGGCGAAGTTGGCGTAATTTGTAGACGAGCTAGCTTCAGTAGTTATTGTTCTTACGGACGTCGCTCCTCAAGT‑CCCCCCCCTCGCACCAACCGACTATATATCAAAAAAATTAAGATGACTGTGCGAAGGTGGCGGAAt                                                                                                         <  2:1057122/151‑1 (MQ=18)
                                         gtgctccaccTGAGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTGCAGGTCGTAGTGTTCCTACGGACGTGGGGGTGCAAGG‑CCCCACCCTCGCACCAACCGCGGCGATATCAAAAAAAGTAAGCTGACTGTGCGAAGGTGGCGGAAtt                                                                                                        <  5:832617/142‑1 (MQ=18)
                                          tGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAAGAGGCGATATCAAAAAAAGTAAGATGACTTTTGGTAGGTGGCGGAAttt                                                                                                       >  8:108204/1‑150 (MQ=255)
                                          tGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGGTGAAGT‑CCCCCCCCGCGCCCCCAAAAAGGGGTTTTCCAATAAATTTAGATGACTGTGGGTAAGTGGCGGAAttt                                                                                                       >  3:2025388/1‑150 (MQ=17)
                                           gCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTGAAGT‑CCCCCCCCACCCACCAAAAAAGTCGAAATAAAAAAAAATATGATGACTGTTATAAGGGGGCAGAAATgg                                                                                                      >  8:815948/1‑151 (MQ=17)
                                           gCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAAGAGGTGAGATCAAATATTGTTTGTTGACTGTGTGTTGTTGGCGGATTTgc                                                                                                      >  1:954657/1‑150 (MQ=14)
                                              gatggtatgCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCACGCACCAAACGCGCCGATATCaa                                                                                                                                           >  6:1337888/9‑111 (MQ=11)
                                                cACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTGAAGT‑CCCCCCCCACGCACCAAAACAAGAGATATAAAAAAAAGTTAGTTGAATGTTTTTTGGTGGCGGAATTTGTTGa                                                                                                  >  5:1853897/1‑150 (MQ=14)
                                                       cGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAAAAGGGGATATCAAAAAATGTTAGATGCATGTTGTTTGGTGCCGGAATTggaataagttaaat                                                                                          >  6:535824/1‑139 (MQ=14)
                                                       cGAAGCTGTCGCTATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTCCGGACGTGGTTGTTCTAGT‑CCCACCCCTCGCACCATACGAGTCGATATCAATAAAATTTAGATGACTGTGCGAAGGTGTCGGAATTGGTAGACGCGCTAg                                                                                          <  4:861610/151‑1 (MQ=21)
                                                        gAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAATAAGGGGAATTCAAATATTGTTTGATGCCTGTTTGTTGGTGGCGGAttttgaatacgtattatc                                                                                         >  3:1998138/1‑136 (MQ=2)
                                                        gAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAAAACGGGAAATCAAAAAAAGTTAGATGAGTGTTTGTAGGTGCCGGAATTGGTATActaggaagc                                                                                         >  5:1490614/1‑143 (MQ=14)
                                                         aaGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAGGTGGGTAGTTCAAAAAAAGTTTGTTGAGTGTGTTTTGGTGGCGGAAttactttaattgttatc                                                                                         >  3:1468604/1‑135 (MQ=11)
                                                             tGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTGAAGT‑CCCCCCCCGCGCACCAAAATAGGGGATATCCAAAAAAGTTTGTTGCCTGTTGTTTGGTGGCGGAATTGGTTTACGCTGTAGCTTc                                                                                      >  7:587737/1‑149 (MQ=14)
                                                               gCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAAGAGGGGATATCAAAAAAAGTTAGATGCCTGTTGTATGGTTGCTGAATTGGTATCCGATCTAGCGTCATgt                                                                                  >  1:855582/1‑151 (MQ=17)
                                                                   aattgtttttgtatcTATCTTCAGCTGTTCGTGTTCTTACGTACTTTTTGGTGGATTT‑CCCTCCACTCGGATCATACGATGTGATATCTATTTAAGTATTTTTTGTTTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTa                                                                              <  8:806538/137‑1 (MQ=0)
                                                                        gTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGGGAGATCAAAAAAAGTAAGTTGCCTGTTGTAAGGTGGCGGATTTGGTAGACTCGAGATCTTCAGGTTTTAgtgt                                                                          >  5:2035043/1‑150 (MQ=18)
                                                                               cgcTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTGAAGT‑CCCCCCCCGCGCCCCAAAAGATGGGATTTCAAAAAAATTTATATGACTGTGGGTAAGTGGCGGAATTTGTTTACGCGCGAGCGTCAGGTGTTAgtgatcattcgc                                                                  >  5:807713/1‑142 (MQ=16)
                                                                               cgcTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAAGACGCGATATCAAAAAAAGTAAGATGAATGTGGGAAGGTGGCGGAATTGGTAGACGAGATAGCATCAGGTGTTAGTGTACTTAAgg                                                                  >  6:583784/1‑151 (MQ=25)
                                                                                  tAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAAAAGGAGAAATCAAAAAAAGTTAGATGACTGTTATAAGGTGGCGGAATTGGTATAATAGATAGCTTCATGTTTTAgtgaaataaaggccg                                                               >  5:1475456/1‑139 (MQ=17)
                                                                                   aGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGGTAGTTCAAAAATAGTTAGATGCATGTGGGTTGGTGGCGGAATTGCTTTACGCTCGAGCTGAACGTTTTAGTGTAATTACggcccg                                                              >  5:923692/1‑147 (MQ=12)
                                                                                     cTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAATAGTGTAAATCAAAAAATGTTAGATGACTGTGTTAAGGTGGCGGAATTGATTTAATCTAAAGCATCAGGTGTTACTGTACATACGGACGTgg                                                            >  5:1042028/1‑151 (MQ=12)
                                                                                     cTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAACT‑CCCCCCCCCCGCCCCAAAACAGGCGATTTCAAAAAAAGTAAGATGCCTGTGGGAAGGGGGCGGAATTGGTTGACGGGCGAGCGTCAGGTGTTTGTGTACTAACGGCCgggg                                                            >  1:2053828/1‑151 (MQ=255)
                                                                                         cggtgtcggttgcgctACGGACGTGGGCGTGCACGT‑CCCCCCCCGCGCACCAACCGCCCCGATATCTATTATAGTACGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGt                                                        <  1:463287/136‑1 (MQ=17)
                                                                                          ggTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAGGAGGCGATATCAAAAAAAGTAAGATGACTGTTCGAAGGTGGCGGAATTGGTATACTCGCTAGCTTCAGGTGTTACTGTACTTACGGACGTGGGGGtt                                                       >  5:1285747/1‑151 (MQ=37)
                                                                                            tgtTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAAGAGGGGAAATAAAAAAAAGTTAGATGACTGTGTGAAGGTGGCGGAATTGGTATAAGAGGTATCATCATGTTTTAATGTAATTAAGGACGgggggggtgg                                                     >  8:831294/1‑147 (MQ=16)
                                                                                                aGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAAGAGGGGAGTTCAAAAAAAGTTAGATGAATGTTGGTAGGTGGCGGATTTAGTATAATGTCTATCATCATGTGTTATTGTAATAACGGAAGTGCGGGtggtaggc                                                 >  3:1795233/1‑144 (MQ=12)
                                                                                                 gtgtTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCCCCAAAAGAGGGGATATCAAAAAAAGTTAGATGACTGTTGGAAGGTGGCGGAATTGGTTGCCGCGGGAGCTTCAGGTGTTAGTGTCCTTACGGCCGGGGGGGtggaaggcc                                                >  3:433277/1‑143 (MQ=25)
                                                                                                  tgtTCTTACGGACGTGGGGGTTCAACT‑CCCCCCCCGCCCCCCAAAAGACGGGATATCAATAAATGTTAGATGCCTGTTGTAAGGTGGCGGAATGGGTAGGCGGGCGACCGTCAGGGGGTAGGGGCATTTCGGCCGGGCGGGtggccgtgc                                                >  2:1634994/1‑142 (MQ=255)
                                                                                                       ttACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAAGAGGAGATATCAAAAAAAGTTAGATGAATGTTAGTAGGTGGCGGAATTGGTAGAAGAGCTATCTTCAGGTTTTAGTGTACTTACGGACGTGGGGGGTTAAGTaccccccc                                          >  5:1998763/1‑151 (MQ=21)
                                                                                                       ttACGGACGTGGGGGTTCAACT‑CCCCCCCCTCGCACCAAAAGAGGAGAAATCAAAAAAACTAACATGAATGAGAGAAGCTGGCGGAATTGGTATACTGACAAGCTTCAGGAGGTACTGTCCTAACGGACGTggggggtgaacgctcaccg                                           >  2:444874/1‑136 (MQ=16)
                                                                                                         aCGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAAGAGGGTATATCAAAAAAAGTTAGTTGACTGTTGGAAGGTGGCGGAATTTTTATAAGTGGTATCTTCATGTTTTAGTGTAATTAAGGACGTGGGGGGTGAAGtaaaaaccaa                                         >  1:1638239/1‑141 (MQ=16)
                                                                                                          cGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTTAGATGCCTGTTGGATGGTGGCGGAATTGGTAGACGGGCGAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGTAGGCCCCCCCCACg                                       >  1:314218/1‑151 (MQ=25)
                                                                                                          cGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACAAAACGAGGCTATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCg                                       <  5:247703/151‑1 (MQ=255)
                                                                                                            gACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATTTCAAAAAAAGTAAGATGACTGTGCGATGGTGGCGGAATTGGTTGACGAGGTAGCTGCAGGTGTGAGTGTCGTTACGGCCGTGGGGGTGGAAGGcccagccccggcc                                     >  5:1793544/1‑141 (MQ=25)
                                                                                                               gTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAATAGGATATATCAAAAAATGTAAGATGAGTGTGATAAGGTCGCGGAATTGGTATACGGTGTATCTTCAGGTGTTAGTGTAGTTAAGCACGTGGGGGTGGAAGtaaaagcaatagcccc                                   >  4:624569/1‑135 (MQ=255)
                                                                                                                      ttCAAGT‑CCCCCCCCTCGCACCAAAGGAGGCGATATCAAAAAAAGTAAGATGACTGTTCTAAGGTGGCGGATTTGGTATACGCGCTAGCTTCAGGTGGTAGTGTCATTACGGACGTGGGGGTTGAACGCCCCCCCCACggccccagacttt                           >  3:2025511/1‑139 (MQ=21)
                                                                                                                      ttCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGATTTGGTAGACGCGGTAGCTTCAGGTGTGAGTGTCCTTTCGGCCGTGGGGGTGGGAGGCCACCCCCTCGCCCCaagtttgt                           >  6:569757/1‑145 (MQ=25)
                                                                                                                              ccccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAATTTTCTTTACTTc                    >  5:1073488/1‑151 (MQ=255)
                                                                                                                              ccccccccTCGCACCAAACGACGCGATATCAAAAAAAGTAAGATGACTGTGCGAACGTGGCGGAATAGGTAGACGCGCGAGCTTCAGGAGTTAGTGTCCTAACGGACGTCAGGGTTCAAGTGCACGCCCAGGCCCCAACTATCTGTACTTc                    >  2:1498919/1‑151 (MQ=255)
                                                                                                                               cccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGGGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCACGCCCCAATTATCTTTACTTcc                   >  1:462330/1‑151 (MQ=255)
                                                                                                                               cccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAACGTGGCGGAATTGGTAGACGAGCTAGCTTCAGGTGTAAGTGTCCTTACGCACGTGCGGGTTCAAGTCCCACCCCACGCACCAAATATCTTTACTTc                    >  6:1977204/1‑150 (MQ=255)
                                                                                                                                ccccccTCGCACCAAAGGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTCGTAGACGAGCTAGCTTCAGGTGTTAGTGTCCTAACGGACGTGCGGGTTCAACTCCCCCCCAACGCACCAAAAAAGATTACTTCCt                  >  2:1902874/1‑151 (MQ=18)
                                                                                                                                ccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCCATTATCTTTACTTCCt                  >  4:1815261/1‑151 (MQ=255)
                                                                                                                                 cccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCGCAATTATCTTTACTTCCtt                 >  6:853658/1‑151 (MQ=255)
                                                                                                                                 cccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCACTCGCACCAATTATCTTTACTTCCtt                 >  2:1051000/1‑151 (MQ=255)
                                                                                                                                  ccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCttt                >  6:1357148/1‑151 (MQ=255)
                                                                                                                                   cccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAATTATCTTTACTTCCTTTc               >  1:1755107/1‑151 (MQ=255)
                                                                                                                                   cccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGGCCCCCCCCTCGCCCCAATTATCTTTACTTCCTTTc               >  5:277131/1‑151 (MQ=255)
                                                                                                                                    ccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGTGCCCCCCCTGGCCCCAATTATGTTTACTTGCTTTCt              >  3:1649128/1‑151 (MQ=21)
                                                                                                                                    ccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCGCCCCCCGCGCACCAATTATCTTTACTTCCTTTCt              >  4:1397312/1‑151 (MQ=255)
                                                                                                                                       cGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGt           >  6:260620/1‑151 (MQ=255)
                                                                                                                                        gCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCGTATCTTGtt          >  4:653267/1‑151 (MQ=255)
                                                                                                                                        gCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCGCCCCAATTATCTTTACTTCCTTTCTTGtt          >  1:1234002/1‑151 (MQ=255)
                                                                                                                                        gCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTGAAGTCCCCCCGCGCGAACCAATAAAGTTTAGTTGCttgattttt          >  4:1207361/1‑144 (MQ=18)
                                                                                                                                         cACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCtttttt         >  2:212493/1‑151 (MQ=255)
                                                                                                                                           ccAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTTTttct       >  7:1564027/1‑151 (MQ=255)
                                                                                                                                            cAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCGCCAAGTTTCTTTTCTTCCTTTCTTGtttatg      >  6:1401760/1‑148 (MQ=255)
                                                                                                                                            cAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAATTATCTTTACTTCCTTTCTTGTttctt      >  3:285521/1‑151 (MQ=255)
                                                                                                                                            cAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTttctt      >  1:799993/1‑151 (MQ=255)
                                                                                                                                             aaaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTGGCACCAATTATCTTTACTTCCTTTCTTGTttctt      >  6:1306157/1‑150 (MQ=255)
                                                                                                                                             aaaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTttcttc     >  2:83817/1‑151 (MQ=255)
                                                                                                                                             aaaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTGGAGTCCCCCCCCACGCCCCAATTATCTTTTCTTGCGTTGTTGgttgtt      >  8:770897/1‑150 (MQ=18)
                                                                                                                                              aaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGTGCCCCCCCTCGCCCCAATTATCTTTACTTCCGTTGTTGTttcttac    >  6:47229/1‑149 (MQ=21)
                                                                                                                                              aaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTcc    >  8:1939902/1‑151 (MQ=255)
                                                                                                                                              aaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAGTCCCCCCCCTCGCACCAATTATCATTAATACCATAGTTGTttcttcg    >  7:1002228/1‑150 (MQ=18)
                                                                                                                                                cGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAGTTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAAATATCTTTACTTCCTTTCTTGTTTCTTCCtt  >  8:1251693/1‑151 (MQ=255)
                                                                                                                                                cGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGTGCCCCCCCGCGCACCAATTATGTTTACTTGCTTTCTTGTTTCTTGCtt  >  4:657220/1‑151 (MQ=21)
                                                                                                                                                cGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTGCCCCCCCTCGCACCAATTATATTTACTTCGTTTCTTGTTTCTTGCtt  >  1:1985361/1‑151 (MQ=255)
                                                                                                                                                cGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCGCTCGCACCAAGTATGTTTACTTCCTTGGGTGTTTGTTCGtt  >  4:1250336/1‑151 (MQ=18)
                                                                                                                                                |                                                                                                                                                      
GGGTTCAAGTCCCCCCCCTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTT  >  DH1/3906689‑3906982

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