Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH1936,1690T.28.9% 192.9 / 80.8 78intergenic (‑13/+21)ECDH1_RS04615/ECDH1_RS04620tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base T (1/53);  new base . (2/20);  total (4/74)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

AATGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCA  >  DH1/936020‑936249
                                                                                                                                                     |                                                                                 
attGGCTCTTTTTTTTGGGAGACCACTACCAATAAAACAGTTGTTTGGTTGTTAAATTATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCt                                                                                   <  2:1033272/148‑1 (MQ=14)
aaTGGCGCTTTTTTATGGGAGACCCGTTACAATAAAACACCGTTTAGGGTGTTCAATTATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtt                                                                                  <  4:174696/151‑1 (MQ=18)
aaTGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtt                                                                                  >  4:534156/1‑151 (MQ=255)
aaTGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtt                                                                                  >  5:376273/1‑151 (MQ=255)
aaTGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtt                                                                                  <  3:756351/151‑1 (MQ=255)
aaTGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtt                                                                                  <  2:33905/151‑1 (MQ=255)
aaTGCCGCGTTTTTATGGGAGACCAGAAAAAAAAAAACCCCCTTCATGGTGTTCAATCATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtt                                                                                  <  7:1650269/151‑1 (MQ=17)
 aTGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACGGTTTAGGTTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttt                                                                                 <  2:654550/151‑1 (MQ=255)
 aTGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttt                                                                                 <  5:327060/151‑1 (MQ=255)
 aTGGCGCTTTTTTATGGGAGACCACAAACAAAAAAACACCCTTTAGGTTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttt                                                                                 <  5:640189/151‑1 (MQ=255)
  tGGCGCTTTTTTATGGGAGCCCAGAAACAAAAAAACACCCGTTAGGGTGTTAAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttt                                                                                 <  3:172061/150‑1 (MQ=255)
   ggatcTTTTTTATGGGAGACCAGTAACAAATAAAGATTTTATTTGTTGTTCAATTATTGGTTTCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtttac                                                                               <  6:1155931/147‑3 (MQ=14)
   ggCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtttt                                                                                <  8:637386/150‑1 (MQ=255)
   ggCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTa                                                                               <  6:1582072/151‑1 (MQ=255)
   ggCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTa                                                                               <  3:1162843/151‑1 (MQ=255)
    gtacgatgttATGCGACACCACTACCAACAAAACATCCTACAGCTTGTTAAATCATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtttact                                                                              <  3:528651/143‑4 (MQ=11)
    gcgcTTTTTTATGGGAGACCAGAAACAAAAAAACCCCCTTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtttact                                                                              <  5:665642/151‑4 (MQ=255)
    gcgcTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTAc                                                                              <  8:823949/151‑1 (MQ=255)
      tcTTTTTTTTGGGAGACCAGAAACAATAAATCAGTTTTTATGTTGTTATATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtttactt                                                                             <  2:696731/149‑5 (MQ=14)
      gcTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACtt                                                                            <  7:1173992/151‑1 (MQ=255)
       gttttttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtc                                                                           <  1:846516/150‑1 (MQ=255)
         ttttttTGGGAGACCAGTAACAATAAATGCGGTGTTTTGTTGTTCTATTATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTtcat                                                                         <  8:760717/150‑1 (MQ=11)
         tttttATTGTAGACCAGAAACAAAAAAACACGCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtca                                                                          <  8:561974/150‑1 (MQ=255)
         tttttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGTTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTtcat                                                                         <  4:481636/150‑1 (MQ=255)
         tttttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtttacttgat                                                                          >  2:1198664/1‑143 (MQ=255)
         tttttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtca                                                                          <  7:530669/150‑1 (MQ=255)
         tttttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtca                                                                          <  7:530672/150‑1 (MQ=255)
          ttttATGGGAGACCTGATACAAAAACTCACGTGTTAGGGTGTTCATTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTtcatc                                                                        <  4:103706/150‑1 (MQ=14)
          ttttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtcat                                                                         <  5:672204/150‑1 (MQ=255)
           tttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtcatc                                                                        <  5:231674/150‑1 (MQ=255)
           tttATGCGAGCCCAGAACCAAAAACACACCCGATATGATGATCTATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTtcatca                                                                       <  3:1481686/150‑1 (MQ=14)
           tctATGGGAGACCTGATACAAAAAAACCCGTGTTATGTTGTTCAATTATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTtcatca                                                                       <  8:1418100/148‑1 (MQ=14)
            ttATGGGAGACCAGCTACAATTAATGAGCCTTTATGTTGTTCAATTATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCaa                                                                      <  2:1403486/150‑1 (MQ=11)
            ttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttttcgtgctgca                                                                       >  1:565449/1‑141 (MQ=255)
            ttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtcatca                                                                       <  8:654393/150‑1 (MQ=255)
             tATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATGaaa                                                                     >  5:1213502/1‑150 (MQ=255)
              aTGGGAGACCAGAAACAAAAAAACACCTTTTAGGTTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAA‑t                                                                   <  6:1211873/150‑1 (MQ=255)
              aTGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttttcgtgatgaacat                                                                    >  7:1419711/1‑139 (MQ=255)
              aTGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCaaa                                                                     <  6:75863/150‑1 (MQ=255)
              aTGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑t                                                                   <  4:1519559/151‑1 (MQ=255)
              aTGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACCTCATCAAA‑t                                                                   <  7:1098413/150‑1 (MQ=255)
               tggcacaCCAGAACAAAAAAAACACCCTTTAGGTTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAA‑tt                                                                  <  2:1477290/144‑1 (MQ=14)
                gggAGACCAGAAACAAAAAACCACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑tt                                                                  <  6:1483790/150‑1 (MQ=255)
                gggAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGTGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑tt                                                                  <  4:1078343/150‑1 (MQ=255)
                gggAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑ttt                                                                 <  6:1460227/151‑1 (MQ=255)
                gggAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑ttt                                                                 <  6:1452218/151‑1 (MQ=255)
                gggAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑tt                                                                  <  4:658530/150‑1 (MQ=255)
                gggAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑tt                                                                  <  8:1524390/150‑1 (MQ=255)
                gggAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑tt                                                                  <  6:1443997/150‑1 (MQ=255)
                  gagaCCAGAAACAAGAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGa                                                               <  4:1022313/151‑1 (MQ=255)
                   agaCCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGAt                                                              <  4:588411/151‑1 (MQ=255)
                   accccAGAAACAATAACAGACCGGTTAGGGTGATCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAA‑ttttaat                                                             <  3:534569/148‑5 (MQ=14)
                    gaCCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGAtt                                                             <  3:159832/151‑1 (MQ=255)
                    caCCAGAAACAAAAAAACACCCGTTAGGTTGTTCATTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAA‑ttttaatt                                                            <  6:1207723/150‑6 (MQ=14)
                    caCCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGAtt                                                             <  7:1394659/150‑1 (MQ=255)
                     cccAGAAACAAAAAATCACCTGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAA‑ttttaattg                                                           <  4:332934/150‑7 (MQ=14)
                     acccgatagaaaaaaaCCCACGATAGGCTGTACAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAA‑ttttaattg                                                           <  1:593392/142‑7 (MQ=11)
                     aCCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTg                                                            <  5:519833/151‑1 (MQ=255)
                     aCCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTg                                                            <  2:387904/151‑1 (MQ=255)
                      ccAGCAACAAAAAAACCCCTTTTATGGTGTTAAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAA‑ttttaattgg                                                          <  4:1154667/151‑8 (MQ=11)
                      ccAGAAACAACAACACACCCGTGAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAA‑ttttaattgg                                                          <  7:1570277/151‑8 (MQ=14)
                      ccAGAAACAAAAAAACACCCGTTAGGTTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAA‑ttttaattgg                                                          <  3:1584313/151‑8 (MQ=255)
                      ccAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTgg                                                           <  7:674821/151‑1 (MQ=255)
                      ccAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTgg                                                           <  6:397089/151‑1 (MQ=255)
                       cAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGATGCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCAGCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGt                                                          <  2:844315/151‑1 (MQ=255)
                        cgttaCAAAAACACATCTGTTTGCTTGTTAAATTATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGtt                                                         <  6:897415/147‑1 (MQ=11)
                        cgAACCAACAAAACCCCCTTCAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGtt                                                         <  1:409433/150‑1 (MQ=11)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGGCCGACGAGCGACCAGTCTGCTCCCCCCCCCGTCGCGCGACTGCGGTACTTTCCTAAAA‑TTTCATTGGtt                                                         >  7:660080/1‑151 (MQ=17)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGtt                                                         <  1:110411/151‑1 (MQ=255)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGtt                                                         <  8:628882/151‑1 (MQ=255)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGtt                                                         <  7:317006/151‑1 (MQ=255)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGtt                                                         <  6:375597/151‑1 (MQ=255)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGtt                                                         <  4:694034/151‑1 (MQ=255)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGtt                                                         <  1:870801/151‑1 (MQ=255)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAA‑ttt                                                                 <  3:1343848/142‑1 (MQ=255)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATAAAA‑ttt                                                                 >  4:1342466/1‑142 (MQ=255)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGGTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGtt                                                         <  4:1496175/151‑1 (MQ=255)
                         gAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTg                                                        <  3:560281/151‑1 (MQ=255)
                          aaaCACAAAAACACCCGTTAGGTTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGC‑TTTACTTCATCAAATTTTAATTGGTTGc                                                       <  1:1199743/151‑1 (MQ=14)
                          aaaCAAAAAAACACCCGTTAGGGTGTTCAATAATTTGTTGCGGGGGACGCATTTGAACAGACGACCTTCGTGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGc                                                       <  4:309683/151‑1 (MQ=255)
                          aaaCAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTCGTGCTTCAAA‑TTTGATTGGTTGc                                                       >  7:1230996/1‑151 (MQ=255)
                          aaaCAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGc                                                       <  8:810975/151‑1 (MQ=255)
                          aaaCAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGc                                                       <  8:759117/151‑1 (MQ=255)
                           aaCAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCg                                                      <  8:1662868/151‑1 (MQ=255)
                           aaCAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCg                                                      <  8:1287158/151‑1 (MQ=255)
                            aCAAAAAAACACCCGTTAGGGTGTTCAATACTTGGGTGCGGGGGCCGGAGTTGAACCGACGACCTTGGGGTGATGAGCCCGACGAGCTACCAGGCGGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCgg                                                     <  1:1068997/151‑1 (MQ=18)
                              aaaaaatCACCCTTTAGGGTGTTAAATAATTGGTTGCGTGGGCCGGATTTTAACGGACTACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACACCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCgggg                                                   <  4:492677/151‑1 (MQ=17)
                                aaaaaCACCCGTTAGGGTGTTCAATAATTGGTTGCGGTGGCCGGATTTGAACCGACTACCTTCGGTTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGc                                                 <  4:1578309/151‑1 (MQ=255)
                                 aaaaaaGCCGTTAGGGTGTTAAATAATTGGTTGCGGGGGGCGGATTCGACCCGACGACCTTCGGGTTATGAGCCCGACGAGCTCCCAGTCTCCTCCACCCCTCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTTATTGGTTGCGGGGGcc                                                <  4:1407454/151‑1 (MQ=17)
                                 aaaaCACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGcc                                                <  3:542770/151‑1 (MQ=255)
                                 aaaaCACCCGTTAGGGTGTTCAATAATTGGTGGCGGGGGCCGGATTTGACCCGACGTCCTTCGGGTTATGAGCCCGACGAGCTACCAGTCTGCTCGACCCCGCGTCACCGTACTGCTTTTACTTCATGAAA‑TTTGATTGGTTGCGGGGGcc                                                <  6:1495055/151‑1 (MQ=18)
                                  aaaCACCCGTTAGGGTGTTCAATAATTGTTTGCGGGGGCCGGACTTGAACCGCCGAGCTTCGTGTTATTAGCCCGACGAGCTACCAGGGTTCTCGAGCCCGCGTCACCGTACTGCTTTTACTTCATGAAA‑TTTGATTGGTTGCGGGGGCCg                                               <  7:630970/151‑1 (MQ=255)
                                    acacCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtt                                                                                  <  3:879160/115‑1 (MQ=255)
                                    acacCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtt                                                                                  >  4:877873/1‑115 (MQ=255)
                                    acacCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGa                                             <  8:1228113/151‑1 (MQ=255)
                                    acacCATTTAGGGTTTTAAATAATTGGTTGCGGGGGCCGGATTCGACCCGACGCCCTTCGGGTTATCAGCCCGACGAGCTACCAGGCTGCTCCACCCCCCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGa                                             <  7:180646/151‑1 (MQ=17)
                                         cGTGAGGGAGTTCAATAATTGGTTGCGGGGGCCGGAGTTGAACCGTCGACCTTCGGGTTATCAGCCCGTATAGCTACGTTTCTTCTCCACCCCTCGTCAACGTAATGCGTTTACATCTTGATA‑TTTGATTGGTTGCGGGGGCCGGATTTGa                                        <  5:1624451/151‑1 (MQ=9)
                                                          atTGGTTGCGGGGGCCGGAGTTGACCCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTTCTCCACCCCGCGTCACCGTACTGCTTTTACTTCCTCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGt                       <  7:96553/151‑1 (MQ=255)
                                                                                gagcggacgacCTTCGGGTTAAGAGCCCGAAGAGCTACCAGGCTCCTCGACCCCGCTTAACCGTACTGCTTTTACTTCATCTTA‑TATGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCa  <  7:90304/145‑1 (MQ=14)
                                                                                                                                                     |                                                                                 
AATGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAA‑TTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCA  >  DH1/936020‑936249

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