Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH1936,1690T.40.6% 148.3 / 107.7 70intergenic (‑13/+21)ECDH1_RS04615/ECDH1_RS04620tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base T (2/39);  new base . (3/25);  total (6/64)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

GAATGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCT  >  DH1/936019‑936253
                                                                                                                                                      |                                                                                     
gAATGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCTTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCt                                                                                       <  5:661649/151‑1 (MQ=255)
 aaTGGCGCTTTTTTATGGGAGACCAGAACAAAAAAAACACCGGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtt                                                                                      <  3:517861/151‑1 (MQ=255)
 aaTGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGTGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtt                                                                                      <  8:1306361/151‑1 (MQ=255)
 aaTGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtt                                                                                      <  4:1297845/151‑1 (MQ=255)
 aaTGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtt                                                                                      <  2:1138384/151‑1 (MQ=255)
 aaTGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCACGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtt                                                                                      >  3:780471/1‑151 (MQ=255)
  aTGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttt                                                                                     <  3:665552/151‑1 (MQ=255)
  aTGGCGCTTTTTATGGTAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtttt                                                                                    <  4:664226/151‑1 (MQ=255)
    ggggattttttATGGTAGACATGAAAAAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttta                                                                                    <  2:708957/145‑2 (MQ=18)
    ggCGCTTTTTTATGGGAGACCAGAAACAAATAAACACCCGTTATGTTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttta                                                                                    <  2:112243/150‑2 (MQ=255)
    ggCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCGGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTa                                                                                   <  2:462368/151‑1 (MQ=255)
    ggCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtttt                                                                                    <  7:844049/150‑1 (MQ=255)
    ggCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttta                                                                                    <  4:956202/150‑2 (MQ=255)
    ggCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTa                                                                                   <  2:965566/151‑1 (MQ=255)
    ggCGCCTTTTTATGGGAGACCAGAAACAAAAAACCAGCGGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGGTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtttt                                                                                    <  4:670915/150‑1 (MQ=255)
     gcgggtttttATGGGAGACCAGAAACTAAAATACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGTCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACACCGCGTCACCGTACTGCTTTTAc                                                                                  <  6:671554/146‑1 (MQ=18)
      cgcgTTTTTATGGGAGACCAGCAAAAAAAAAACACCGGATAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtttactt                                                                                 <  3:262561/151‑5 (MQ=17)
      cgcTTTTTTATGGGAGACCAGAAACAAAAAAACACCGGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACt                                                                                 <  6:1119179/151‑1 (MQ=255)
       gcTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACtt                                                                                <  6:1207506/151‑1 (MQ=255)
       gcTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACtt                                                                                <  8:739321/151‑1 (MQ=255)
       gcTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACtt                                                                                <  7:907666/151‑1 (MQ=255)
       gcTTTTTTATGGCAGACCTGAACCAAAAAAAGCGCCTTGAGGTTGATCATTTATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCtttactt                                                                                 <  8:64483/150‑5 (MQ=14)
        cTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTtc                                                                               >  2:987056/1‑150 (MQ=255)
         ttttttATGGGAGACAAGAAACAAAAAAACAGCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtc                                                                               <  6:712878/150‑1 (MQ=255)
          tttttATGTGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCTCATTTGAACCGACGACCTTCGTGTTATGAGCCCTACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtc                                                                               <  5:1012666/149‑1 (MQ=255)
          tttttATGGGACACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtca                                                                              <  8:654803/150‑1 (MQ=255)
          tttttAGGGGAGACCTGATAAAACAAATCAGGCCATAGGGTGTTAATTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTtca                                                                              <  7:551782/149‑1 (MQ=11)
          tttttAGGGGAGACCAGAACCCCCAAAACCCCCGTTAGGGTGTTCAATCATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTtcat                                                                             <  5:548772/150‑1 (MQ=17)
           ttttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCGATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtcat                                                                             <  3:630916/150‑1 (MQ=255)
           ttttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtcat                                                                             <  3:256338/150‑1 (MQ=255)
            tttATGTGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtcatc                                                                            <  4:470700/150‑1 (MQ=255)
            tttATGGGAGACCTGAAACATATAATCATTCGTTAGGTTGTTCAATAATTGGTTGCGGGGTGCGGTTTTGAACCGAGGAGGTTGGGGTTATTTGCCCGACGAGCTAGCATGCTGCTCGTCCCCTCGTGACGGTACTGCTTTTACTtcatc                                                                            <  6:468935/150‑1 (MQ=255)
            tttATGGGAGACCAGAAACAATAAAAGACTCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTtcatca                                                                           <  6:996210/150‑1 (MQ=255)
            tttATGGGAGACCAGAAACAAAAAACCACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCGGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtcatc                                                                            <  6:478409/150‑1 (MQ=255)
             ttATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTtcatca                                                                           <  3:1132139/150‑1 (MQ=255)
              tATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCaa                                                                          <  7:649985/150‑1 (MQ=255)
              tATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCaa                                                                          <  7:118376/150‑1 (MQ=255)
              tATGGGAGACCACAAACAACAACACCCCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCaaa                                                                         <  3:982436/150‑1 (MQ=17)
               aTGGGAGCCCAGAAACAAAAAAACACGCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAtt                                                                      <  5:1018613/151‑1 (MQ=255)
               aTGGGAGACCAGAAACAACAAAACACCCGTCAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAt                                                                       <  1:5453/150‑1 (MQ=255)
               aTGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCaaa                                                                         <  7:1322232/150‑1 (MQ=255)
               aTGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCaaa                                                                         <  6:128899/150‑1 (MQ=255)
                ttggAGCCCAGAAACAAAAAAACACCTGTTAGGTTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAtt                                                                      <  5:920026/148‑1 (MQ=17)
                tGGGAGACCTGAAACAATAAAAGCGCTGTTTGGTTGTTCATTTATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAtt                                                                      <  8:271360/150‑1 (MQ=11)
                 gggAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAAttt                                                                     <  5:690887/151‑1 (MQ=255)
                 ggccccccAGAAACAAAAAAACACTGTTTATGGTTTACAATAATTGGGTGGGGGGGCCGGTATTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAtttt                                                                    <  3:1290430/145‑2 (MQ=14)
                  ggagccctcaaaCAAATAAACAATCTTTTGGTTGTTAAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAtttta                                                                   <  6:536072/142‑3 (MQ=11)
                  ggAGACCCGAACCACCAACACCCGCGAGAGGGTGTTCAATCATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAtttta                                                                   <  7:787507/151‑3 (MQ=11)
                  ggAGACCAGCAACAAAAAAACACCCGTTAGGGTGTTCAATAATGGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTg                                                                    <  2:531105/151‑1 (MQ=255)
                  ggAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTg                                                                    <  2:105309/151‑1 (MQ=255)
                   gagaCCAGAAACAAAAAATGACCTTTTATGGTGTTAAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAttttaa                                                                  <  2:523080/151‑4 (MQ=14)
                   gagaCCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGa                                                                   <  1:582997/151‑1 (MQ=255)
                    agaCCAGAAACACAAACACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAttttaat                                                                 <  1:987032/151‑5 (MQ=255)
                    agaCCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTACACAAAttttaat                                                                 >  3:167895/1‑147 (MQ=17)
                     gaCCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGAtt                                                                 <  5:658822/151‑1 (MQ=255)
                     gaCCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGAtt                                                                 <  6:1339915/151‑1 (MQ=255)
                     cacctcaaccAAAAAAACACCCGTTTGGTTGTTAATTTATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAttttaatt                                                                <  4:168148/142‑6 (MQ=11)
                     caCCGGAAACAAAAAAACAGTGGTTAGCGTGTTCAATTATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAttttaatt                                                                <  6:590624/150‑6 (MQ=11)
                     caCCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAttttaatt                                                                <  5:596835/150‑6 (MQ=255)
                     caCCAGAAAAAAAAAATCCCCCGTGATGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAttttaatt                                                                <  8:450596/150‑6 (MQ=14)
                      aCCTGAAAAAATAAAACAGTTTTTTGGTTGTTATATTATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAttttaattg                                                               <  4:780214/151‑7 (MQ=11)
                      aCCAGAAACAATAAAAGACCTTATAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAttttaattg                                                               <  1:772194/151‑7 (MQ=11)
                      aCCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGCTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGAtn                                                                 <  4:129766/150‑2 (MQ=255)
                       ccAGAAACAAATAAACAGCCTTTAGGTTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAttttaattgg                                                              <  6:1000057/151‑8 (MQ=14)
                         aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGtt                                                             <  7:328057/151‑1 (MQ=255)
                         aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGtt                                                             <  1:313241/151‑1 (MQ=255)
                         aGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGtt                                                             <  4:1207776/151‑1 (MQ=255)
                          gattcAAATAAACAGCTTTTTGGTTGTTAATTTATTGGTTGCGGGGGCGGGATTTGTCCCGACGCCCTTGGGTTTATTAGCCCGACCAGCTTCCAGGCTTCTCCAGCCCGCGTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTg                                                            <  4:11034/147‑1 (MQ=9)
                          gAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTgg                                                               <  5:475162/148‑1 (MQ=255)
                          gAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTgg                                                               >  6:474766/1‑148 (MQ=255)
                          gAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTg                                                            <  2:1022209/151‑1 (MQ=255)
                          gAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTg                                                            <  1:1161268/151‑1 (MQ=255)
                          gAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGCTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTg                                                            <  3:1291029/151‑1 (MQ=255)
                          gAAAAAAAAAAACACCCTTTTTGTTGTTAAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTg                                                            <  2:1093174/151‑1 (MQ=11)
                           aaaaaaatactcacCCTTTTTGGTGTTAAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACTAGCTACCAGGCTTCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTg                                                            <  8:792761/139‑1 (MQ=11)
                           aaaCAAAAAACCACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACGTTCGGGTTATGAGCGCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGc                                                           <  3:88501/151‑1 (MQ=255)
                           aaaCAAAAAAACAGCCGTTTGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAATTTTAATTGGTTGc                                                           <  7:836800/151‑1 (MQ=14)
                           aaaCAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGc                                                           <  7:85990/151‑1 (MQ=255)
                           aaaCAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAttaaatgttcggg                                                           >  4:406434/1‑139 (MQ=17)
                              cAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCggg                                                        <  7:113361/151‑1 (MQ=255)
                               aaaaaaaCACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCgggg                                                       <  8:123811/151‑1 (MQ=255)
                                       acCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCt                                                                                       <  1:849942/112‑1 (MQ=255)
                                       acCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCt                                                                                       >  2:850123/1‑112 (MQ=255)
                                             tttgttgtTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTACTTCATCAAAttttaattggt                                                             <  8:962577/125‑9 (MQ=11)
                                                ggTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttta                                                                                    >  3:706282/1‑105 (MQ=40)
                                                ggTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCttta                                                                                    <  4:705876/106‑2 (MQ=40)
                                                        ataatTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCCGGGTACTCGGTGTAGTGACCGAAATTTGATTCGTTGCGGGGGCGTGATTTGAAGCGCCGCCCtt                                >  7:541820/1‑149 (MQ=11)
                                                                                    ccGACGACGTTCGGTTTATGAGCCCGTCTAGCTACCATTGTTCTCCACCCCGCTTCACAGTACTGCTTTTACTTCTTCAAATTTTATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGgct  <  7:312128/151‑1 (MQ=11)
                                                                                                    tATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTAATTAATCAAATTTAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTa                         >  2:711453/1‑112 (MQ=255)
                                                                                                                                                      |                                                                                     
GAATGGCGCTTTTTTATGGGAGACCAGAAACAAAAAAACACCCGTTAGGGTGTTCAATAATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGTCACCGTACTGCTTTTACTTCATCAAATTTGATTGGTTGCGGGGGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCT  >  DH1/936019‑936253

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: