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

AATGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCA  >  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)
         ttttttTGGGAGACCAGTAACAATAAATGCGGTGTTTTGTTGTTCTATTATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTtcat                                                                         <  8:760717/150‑1 (MQ=11)
         tttttATTGTAGACCAGAAACAAAAAAACACGCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtca                                                                          <  8:561974/150‑1 (MQ=255)
         tttttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGTTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTtcat                                                                         <  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)
          ttttATGGGAGACCTGATACAAAAACTCACGTGTTAGGGTGTTCATTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTtcatc                                                                        <  4:103706/150‑1 (MQ=14)
          ttttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtcat                                                                         <  5:672204/150‑1 (MQ=255)
           tttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtcatc                                                                        <  5:231674/150‑1 (MQ=255)
           tttATGCGAGCCCAGAACCAAAAACACACCCGATATGATGATCTATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTtcatca                                                                       <  3:1481686/150‑1 (MQ=14)
           tctATGGGAGACCTGATACAAAAAAACCCGTGTTATGTTGTTCAATTATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTtcatca                                                                       <  8:1418100/148‑1 (MQ=14)
            ttATGGGAGACCAGCTACAATTAATGAGCCTTTATGTTGTTCAATTATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCaa                                                                      <  2:1403486/150‑1 (MQ=11)
            ttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttttcgtgctgca                                                                       >  1:565449/1‑141 (MQ=255)
            ttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtcatca                                                                       <  8:654393/150‑1 (MQ=255)
             tATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATGaaa                                                                     >  5:1213502/1‑150 (MQ=255)
              aTGGGAGACCAGAAACAAAAAAACACCTTTTAGGTTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAt                                                                   <  6:1211873/150‑1 (MQ=255)
              aTGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttttcgtgatgaacat                                                                    >  7:1419711/1‑139 (MQ=255)
              aTGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCaaa                                                                     <  6:75863/150‑1 (MQ=255)
              aTGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAAt                                                                   <  4:1519559/151‑1 (MQ=255)
              aTGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACCTCATCAAAt                                                                   <  7:1098413/150‑1 (MQ=255)
               tggcacaCCAGAACAAAAAAAACACCCTTTAGGTTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAtt                                                                  <  2:1477290/144‑1 (MQ=14)
                gggAGACCAGAAACAAAAAACCACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAAtt                                                                  <  6:1483790/150‑1 (MQ=255)
                gggAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGTGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAAtt                                                                  <  4:1078343/150‑1 (MQ=255)
                gggAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAAttt                                                                 <  6:1460227/151‑1 (MQ=255)
                gggAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAAttt                                                                 <  6:1452218/151‑1 (MQ=255)
                gggAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAAtt                                                                  <  4:658530/150‑1 (MQ=255)
                gggAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAAtt                                                                  <  8:1524390/150‑1 (MQ=255)
                gggAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAAtt                                                                  <  6:1443997/150‑1 (MQ=255)
                  gagaCCAGAAACAAGAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGa                                                               <  4:1022313/151‑1 (MQ=255)
                   agaCCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGAt                                                              <  4:588411/151‑1 (MQ=255)
                   accccAGAAACAATAACAGACCGGTTAGGGTGATCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAttttaat                                                             <  3:534569/148‑5 (MQ=14)
                    gaCCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGAtt                                                             <  3:159832/151‑1 (MQ=255)
                    caCCAGAAACAAAAAAACACCCGTTAGGTTGTTCATTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAttttaatt                                                            <  6:1207723/150‑6 (MQ=14)
                    caCCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGAtt                                                             <  7:1394659/150‑1 (MQ=255)
                     cccAGAAACAAAAAATCACCTGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAttttaattg                                                           <  4:332934/150‑7 (MQ=14)
                     acccgatagaaaaaaaCCCACGATAGGCTGTACAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAttttaattg                                                           <  1:593392/142‑7 (MQ=11)
                     aCCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTg                                                            <  5:519833/151‑1 (MQ=255)
                     aCCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTg                                                            <  2:387904/151‑1 (MQ=255)
                      ccAGCAACAAAAAAACCCCTTTTATGGTGTTAAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAttttaattgg                                                          <  4:1154667/151‑8 (MQ=11)
                      ccAGAAACAACAACACACCCGTGAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAttttaattgg                                                          <  7:1570277/151‑8 (MQ=14)
                      ccAGAAACAAAAAAACACCCGTTAGGTTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAttttaattgg                                                          <  3:1584313/151‑8 (MQ=255)
                      ccAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTgg                                                           <  7:674821/151‑1 (MQ=255)
                      ccAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTgg                                                           <  6:397089/151‑1 (MQ=255)
                       cAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGATGCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCAGCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGt                                                          <  2:844315/151‑1 (MQ=255)
                        cgttaCAAAAACACATCTGTTTGCTTGTTAAATTATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGtt                                                         <  6:897415/147‑1 (MQ=11)
                        cgAACCAACAAAACCCCCTTCAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGtt                                                         <  1:409433/150‑1 (MQ=11)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGGCCGACGAGCGACCAGTCTGCTCCCCCCCCCGTCGCGCGACTGCGGTACTTTCCTAAAATTTCATTGGtt                                                         >  7:660080/1‑151 (MQ=17)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGtt                                                         <  1:110411/151‑1 (MQ=255)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGtt                                                         <  8:628882/151‑1 (MQ=255)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGtt                                                         <  7:317006/151‑1 (MQ=255)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGtt                                                         <  6:375597/151‑1 (MQ=255)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGtt                                                         <  4:694034/151‑1 (MQ=255)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGtt                                                         <  1:870801/151‑1 (MQ=255)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAttt                                                                 <  3:1343848/142‑1 (MQ=255)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATAAAAttt                                                                 >  4:1342466/1‑142 (MQ=255)
                        aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGGTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGtt                                                         <  4:1496175/151‑1 (MQ=255)
                         gAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTg                                                        <  3:560281/151‑1 (MQ=255)
                          aaaCACAAAAACACCCGTTAGGTTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGc                                                       <  1:1199743/151‑1 (MQ=14)
                          aaaCAAAAAAACACCCGTTAGGGTGTTCAATAATTTGTTGCGGGGGACGCATTTGAACAGACGACCTTCGTGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGc                                                       <  4:309683/151‑1 (MQ=255)
                          aaaCAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTCGTGCTTCAAATTTGATTGGTTGc                                                       >  7:1230996/1‑151 (MQ=255)
                          aaaCAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGc                                                       <  8:810975/151‑1 (MQ=255)
                          aaaCAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGc                                                       <  8:759117/151‑1 (MQ=255)
                           aaCAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCg                                                      <  8:1662868/151‑1 (MQ=255)
                           aaCAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCg                                                      <  8:1287158/151‑1 (MQ=255)
                            aCAAAAAAACACCCGTTAGGGTGTTCAATACTTGGGTGCGGGGGCCGGAGTTGAACCGACGACCTTGGGGTGATGAGCCCGACGAGCTACCAGGCGGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCgg                                                     <  1:1068997/151‑1 (MQ=18)
                              aaaaaatCACCCTTTAGGGTGTTAAATAATTGGTTGCGTGGGCCGGATTTTAACGGACTACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACACCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCgggg                                                   <  4:492677/151‑1 (MQ=17)
                                aaaaaCACCCGTTAGGGTGTTCAATAATTGGTTGCGGTGGCCGGATTTGAACCGACTACCTTCGGTTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGc                                                 <  4:1578309/151‑1 (MQ=255)
                                 aaaaaaGCCGTTAGGGTGTTAAATAATTGGTTGCGGGGGGCGGATTCGACCCGACGACCTTCGGGTTATGAGCCCGACGAGCTCCCAGTCTCCTCCACCCCTCGTCACCGTACTGCTTTTACTTCATCAAATTTTATTGGTTGCGGGGGcc                                                <  4:1407454/151‑1 (MQ=17)
                                 aaaaCACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGcc                                                <  3:542770/151‑1 (MQ=255)
                                 aaaaCACCCGTTAGGGTGTTCAATAATTGGTGGCGGGGGCCGGATTTGACCCGACGTCCTTCGGGTTATGAGCCCGACGAGCTACCAGTCTGCTCGACCCCGCGTCACCGTACTGCTTTTACTTCATGAAATTTGATTGGTTGCGGGGGcc                                                <  6:1495055/151‑1 (MQ=18)
                                  aaaCACCCGTTAGGGTGTTCAATAATTGTTTGCGGGGGCCGGACTTGAACCGCCGAGCTTCGTGTTATTAGCCCGACGAGCTACCAGGGTTCTCGAGCCCGCGTCACCGTACTGCTTTTACTTCATGAAATTTGATTGGTTGCGGGGGCCg                                               <  7:630970/151‑1 (MQ=255)
                                    acacCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtt                                                                                  <  3:879160/115‑1 (MQ=255)
                                    acacCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtt                                                                                  >  4:877873/1‑115 (MQ=255)
                                    acacCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGa                                             <  8:1228113/151‑1 (MQ=255)
                                    acacCATTTAGGGTTTTAAATAATTGGTTGCGGGGGCCGGATTCGACCCGACGCCCTTCGGGTTATCAGCCCGACGAGCTACCAGGCTGCTCCACCCCCCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGa                                             <  7:180646/151‑1 (MQ=17)
                                         cGTGAGGGAGTTCAATAATTGGTTGCGGGGGCCGGAGTTGAACCGTCGACCTTCGGGTTATCAGCCCGTATAGCTACGTTTCTTCTCCACCCCTCGTCAACGTAATGCGTTTACATCTTGATATTTGATTGGTTGCGGGGGCCGGATTTGa                                        <  5:1624451/151‑1 (MQ=9)
                                                          atTGGTTGCGGGGGCCGGAGTTGACCCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTTCTCCACCCCGCGTCACCGTACTGCTTTTACTTCCTCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGt                       <  7:96553/151‑1 (MQ=255)
                                                                                gagcggacgacCTTCGGGTTAAGAGCCCGAAGAGCTACCAGGCTCCTCGACCCCGCTTAACCGTACTGCTTTTACTTCATCTTATATGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCa  <  7:90304/145‑1 (MQ=14)
                                                                                                                                                     |                                                                                 
AATGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCA  >  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: