Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,184,5850CA27.9% 159.3 / 24.3 81intergenic (+8/‑30)ECDH1_RS15795/ECDH1_RS15800tRNA‑Gln/tRNA‑TERM
Reads supporting (aligned to +/- strand):  ref base C (28/29);  new base A (20/2);  total (49/32)
Fisher's exact test for biased strand distribution p-value = 6.46e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATCCCGTCGTAGCCACCAAATTCTGAATGTATCGAATATGTTCGGCAAATTCAAAACCAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGCGAGCTTTTTGCTTTTCTGCGTTC  >  DH1/3184436‑3184734
                                                                                                                                                     |                                                                                                                                                     
gAATCCCGTCGTAGCCACCAAATTCTGAATGTATCGAATATGTTCGGCAAATTCAAAACCAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCa                                                                                                                                                      <  8:221631/151‑1 (MQ=255)
  aTCCCGTCGTAGCCACCAAATTCTGAATGTATCGAATATGTTCGGCAAATTCAAAACCAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAg                                                                                                                                                    <  6:944359/151‑1 (MQ=255)
    cccGTCGTAGCCACCAAATTCTGAATGTATCGAATATGTTCGGCAAATTCAAAACCAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGa                                                                                                                                                   <  6:735000/150‑1 (MQ=255)
    cccGTCGTAGCCACCAAATTCTGAATGTATCGAATATGTTCGGCAAATTCAAAACCAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGAc                                                                                                                                                  <  8:14804/151‑1 (MQ=255)
         cgtAGCCACCAAATTCTGAATGTATCGAATATGTTCGGCAAATTCAAAACCAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTATAAAAGatccctt                                                                                                                                             >  8:558305/1‑145 (MQ=21)
         cgtAGCCACCAAATTCTGAATGTATCGAATATGTTCGGCAAATTCAAAACCAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTCATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTAATAAAGaaccat                                                                                                                                              >  8:56196/1‑145 (MQ=255)
           tAGCCACCAAATTCTGAATGTATCGAATATGTTCGGCAAATTCAAAACCAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAAATTATTCCAGAAGATTAc                                                                                                                                            >  7:23965/1‑150 (MQ=255)
              ccaccaCATGATGAATGTATCCAATAGGAGCGGCACATTCAACCCCAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTg                                                                                                                                        <  3:933121/151‑1 (MQ=21)
                accaAATTCTGAATGTATCGAATATGTTCGGCAAATTCAAAACCAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTa                                                                                                                                      <  6:429021/151‑1 (MQ=255)
                 acaAATTCTGAATGTATCGAATATGTTCGGCAAATTCAAAACCAATTTGTTGGTGTATCCCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTaa                                                                                                                                     <  3:942659/150‑1 (MQ=255)
                  caAATTCTGAATGTATCGAATATGTTCGGCAAATTCAAAACCAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAg                                                                                                                                    <  3:879698/151‑1 (MQ=255)
                  caAATTCTGAATGTATCGAATATGTTCGGCAAATTCAAAACCAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAg                                                                                                                                    <  3:874461/151‑1 (MQ=255)
                   aaaTTCTGAATGTATCTAATATGTTCGGCAAATTCAAAACCAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTAGCGCAGGAGAGTTGTTGAAGAAGCTTCACTTGTaggg                                                                                                                                   >  8:944066/1‑148 (MQ=21)
                      ttCTGAATGTATCGAATATGTTCGGCAAATTCAAAACCAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAAATTATAAAAGATGGCTACGTTGtcatttat                                                                                                                                >  1:328008/1‑144 (MQ=18)
                          gAATGTATCGAATATGTTCGGCAAATTCAAAACCAATTTGTTGGGGTATCGCCAAGCGGTAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCAccc                                                                                                                             <  7:798579/149‑1 (MQ=255)
                           aaTCTATCCAATATGTTCTGCAAATTCAAACCCAATTTGTTTGTGTATCTCCAAGCCGTAATGCACCTGATTCTGATTCCGGCATTCCTAGTTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAg                                                                                                                           <  2:57325/151‑1 (MQ=25)
                                    aaTATGTTCGGCAAATTCAAAACCAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTa                                                                                                                   <  6:66831/150‑1 (MQ=255)
                                     atatTTTCGGAAAATTAATATCTTCTTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATc                                                                                                                 <  4:173052/151‑1 (MQ=21)
                                            cGGCAATTTCAAATCGAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAg                                                                                                           <  6:1056179/150‑1 (MQ=255)
                                               aaaaTTCAATACAAATTTGTTGGGTTATCGCAAAGCTGTAAGGCAACGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTTTAAGTGCACCGAGTTGGGGTATCGCCAAGCGGt                                                                                                       <  6:92542/150‑1 (MQ=25)
                                                ataTTCAAAACCACTTTGTTTGTGTATCTCCAACCTTTAAGGAACATGATTCTGTTTCAGGCTTTACGAGGTTCGAATCCTCGTAACCCAGACAATTTATTAAAGACGCTTAAATTTTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTa                                                                                                      <  4:616323/149‑1 (MQ=255)
                                                  aTTCAAAACCAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCACCGTACCTGGTAAGTGTACCGATTGTGGGTACCTGCAAGGGGTa                                                                                                      >  8:594459/1‑149 (MQ=11)
                                                    tCAAAACCAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATAAAAGAGGCGTACCTTGTTAGTGGAGCTATTTGGCGTTTCGCCAAGCGGTAAGGc                                                                                                  >  7:929/1‑151 (MQ=18)
                                                     cAAAACCAGTTTGTTTGGGTCTTTCAACGCTGTAAGGCACCTGATGCTGTTTCCGGCACGCATCGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTAACTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGc                                                                                                  <  5:722068/150‑1 (MQ=255)
                                                         aCCAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTAAAGAAGCTTACCTTGTAAGTGCACCAAGTTTGGGTATCGCCAAGAGGTAAGGCAACgg                                                                                             >  3:177133/1‑151 (MQ=25)
                                                          cccgtttATTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCt                                                                                                                                          <  4:814400/101‑1 (MQ=255)
                                                          cccggttattgGGGTATCGCCAAGCGGTAAGGCACCGGTTTTTGATACCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAAGTTATTAAAGACGGTTACCt                                                                                                                                          >  3:814749/9‑105 (MQ=11)
                                                             aTTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTAAAGAGGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCTGTAAGGCACCGGATTc                                                                                         >  4:1070540/1‑151 (MQ=255)
                                                                 attgGGGTATCGCCAAGCGGTAAGGCACCGGTTTGTGATACCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTg                                                                                                                        >  2:1067031/2‑116 (MQ=255)
                                                                 attgGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTg                                                                                                                        <  1:1067424/115‑1 (MQ=255)
                                                                    ggggTATCGCCAAGCGGTAAGGCACCGGTTTTTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCg                                                                                                                >  8:339805/1‑121 (MQ=255)
                                                                    ggggTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACTTTATTCAAGCCGCTTACCTTGTGAGTGCTCCCTGTTGGGGTATCGGCAAGCGGGAAGGCTCCggtgtctgcttgc                                                                                  >  7:548598/1‑140 (MQ=18)
                                                                    ggggTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCg                                                                                                                <  7:340160/121‑1 (MQ=255)
                                                                         aTCGCCAAGCGGTAAGGCACCGGTTTGTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGt                                                                                                                                   >  2:103779/1‑97 (MQ=255)
                                                                         aTCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGGTACGGAGTGCTTTATCCGAGTTGGCGTATCGCCAAGCGGTAAGGCGCCGGagacaatttactgcat                                                                             >  5:907462/1‑136 (MQ=9)
                                                                         aTCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGt                                                                                                                                   <  1:103841/97‑1 (MQ=255)
                                                                            gCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTAAAGCCCCTTACCTTGTAAGTGAACCGAGTTGGGGTATCGCCAAGCGCTAAGGCACCGGAGTCTAATTCCGGCCTTcc                                                                          >  1:186456/1‑151 (MQ=25)
                                                                               aaGCGGTAAGGCACCGGTTTTTGATACCGGCATTCCCTGGTTCGAATCCTCGTACCCCAGCCAATTTATTc                                                                                                                                                       >  3:27013/1‑71 (MQ=1)
                                                                                   cgcaacccACCGGATGCTGATACCGGCATTCCGAGGGTCGAATCCTCGTACGCCAGACAATTTATTCAAGCCGCTTACCTTGTAACTGCACGGAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTc                                                                   <  1:554553/144‑1 (MQ=25)
                                                                                    tgaagcctccGGACTCTGTTACCGGCATTCCGCGCTTCGAATCCTCCTACACCATCATATTTATTCAATACTCTTAACTTGTAAGTGCACCGAGTTGGGGTTTCGCCAAGCGGTAAGGCACCGGATTCTTATTCCGGCATTCCGAgg                                                                      <  6:603099/139‑1 (MQ=18)
                                                                                     cacgcccccGGATGCGGATACCGGCACTCCGCGCTTCGAATCCTCGTACCCCAGACAATTTATTCAAGCCGCTTACCTTGTAATTGCAACCTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGa                                                                 <  5:659613/144‑1 (MQ=21)
                                                                                         cctaaggATTCTGTGTCCGGCATTAATCGGGTCGAATCCTCGTAACCCTGCATATTTCTTATATTCTCGTCCCATGTAATTCCACTTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATcc                                                             <  3:1032659/146‑1 (MQ=2)
                                                                                          cACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTAATCAAGATCCGTACCTTGAAAGTGCACCGAGTTTGGGTATCTACAATGTCTAATGCTACGGATACTAATTGCGGCATTcaatgttgcgattccc                                                            >  7:460227/1‑136 (MQ=12)
                                                                                          cACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAAATTATAAAAGAAGCATACGTAGTAAGTTCACCCAGTGGGGGTATCGCGAATCGCTAAGGCTCAGGATTCTAATACGAGCATTCAGAGGttgcaaacc                                                             >  1:515030/1‑143 (MQ=16)
                                                                                          aacaggATTGTGACTCAGGCATTACTAGCTTCGAAGCCTCGAAAGCCTGACACTTTATTCTATCCGCTTACATTGTAATTGCAATCTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCt                                                            <  8:391440/147‑1 (MQ=9)
                                                                                           aaccgATCCTGACACCGGCACTCCCACCTTCGAATCCTCGTACACCTTAAAATTTATTCTAGCCGCTTACCTTTTAACTGCACCGAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTc                                                           <  6:559947/147‑1 (MQ=17)
                                                                                            ccGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTAAAGAGGCTTACCTTGTAAGTGAACCGAGTTGGGGTATAGACAAGCTGTAAGGCGCCGGTTTCTAGTTCCGTCATTCCGAGGTTCgtctccggg                                                          >  3:196238/1‑143 (MQ=25)
                                                                                                tttttGATACCGGCATTCCCAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCaa                                                                                                                                                     <  2:670387/52‑1 (MQ=17)
                                                                                                   cTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTAAAAAAGCAGCTTACGTGGTAAGTGCACCCATTTGGGGTATCGCCAAGCGGTAAGGCACGGGAGTCTAATTCCGGCATTCAGAGTATCGAAACCGCGTAccgcc                                                   >  7:460903/1‑148 (MQ=21)
                                                                                                     gATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAAATTAAAAAAGAACCGTACCTGGTGAGTGGACCCAGTTGGGGTAGCGGCAAGGGGTAATGCCACGTAGTCTACTTGCCTCATTCCGAGTATCGAGACCACGTCccccctc                                                 >  7:880091/1‑148 (MQ=9)
                                                                                                         ccGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTACGGCATTCCGAGTTTCGAATCCTCGTACCCCAGCcac                                              >  5:1039812/1‑150 (MQ=255)
                                                                                                                 ccGAGGTTCGAATCCTCGTACCCCAGCCAAATTATAAAAGCAGCTTACCTTGTTAGTGATCCGAGTAGGGGTAAAGCCAATCTCTAAGGCACAGTTTTCTACTTAATTCATTCAGAGTCTCGACTCAAGGTCCCCCCGCAACATTaaaaa                                      >  5:720375/1‑150 (MQ=255)
                                                                                                                  cGAGGTTCGAATCCTCGTACCCCAGCCAATTTATAAAAGAAGCTTTCCTTGTAAGTGCAGCGATTTGGGGTAGAGACAATCTCTAATGCGCGGGTTACTAATTGCGTCATTCAGAGTGTCGACTCCGCGTACCGCAGCCGCATTaaaaaat                                    >  5:713859/1‑150 (MQ=255)
                                                                                                                  cGAGGTTCGAATCCTCGTACCCCAGCCAAATTAATAAAGAAGCATACCTAGTAAGTGAACCAATTAGGGGTATATCCAAGCGCTAATGCAACGTATACTACTAAAATCATTCAGAGTATCCAATCCTCCTACCCCAGCGAAAATaaaaaac                                    >  8:892711/1‑150 (MQ=255)
                                                                                                                   gAGGTTCGAATCCTCGTACCCCAGCCAATTTATTAAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATAGACAAGCGGTAATGCACCGGATTCTACTTGCTGCATTCCGAGGATCGAATCCTCGTACCCCAGCCACATTAAAAAAGc                                   >  1:802592/1‑151 (MQ=255)
                                                                                                                    aGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGTTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCt                                  >  3:122344/1‑151 (MQ=255)
                                                                                                                         cGAATCCTCGTACCCCAGCCAATTTATAAAAGACCCATACCTGGTACGTGCAGCCAGTTGGGGTATCGCCAATCGGTAAGGCACAGGAGTCTAATTCAGGCATTCCCAGTATCGAAGCCTCGTGCCCCAGCGACATTAAAAAACCTCGCtt                             >  5:712786/1‑151 (MQ=255)
                                                                                                                         cGAATCCTCGTACCCCAGCCAATTTAATAAAGAGGCTTACCTAGTAAGTTAACCAAGTTTGCGTATCTACAATCTGTAATGCACAGGATACTAATTAAATCATTCAGTGTATCAAATACTACTAAACCAGCAACATTAAAAAAGCTCTCtt                             >  8:155659/1‑151 (MQ=255)
                                                                                                                          gAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCATCGTACCCCAGCCACATAAAAAAACCTCGCTTc                            >  4:188028/1‑151 (MQ=255)
                                                                                                                          gAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTc                            >  2:994192/1‑151 (MQ=255)
                                                                                                                          gAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGACAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGATCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTc                            >  3:95762/1‑151 (MQ=255)
                                                                                                                           aagccccctACCCCCGCCAATTTATTTTATTATCTTAACTTGTAACAGCTAGCTGTTGTGGTATCTCCAATCTGTTAGGCAACTGATTCTTATTCCGTCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCg                           <  2:101182/143‑1 (MQ=9)
                                                                                                                           aaTCCTCGTACCCCAGCCAATTTATACAAGCACCTTACCTTGTAAGTTAACCCAGTAGGCGTAACGACAAGCTGGAATGCACAGGATACTCATTCGGGCATTCAAAGGATCGAAGCAACGTGCCCCAGCGACAATAAAAAAGCTAGCTACg                           >  6:682269/1‑151 (MQ=255)
                                                                                                                              ccTCGTACCGCAGCCAATATATTCAAGCCGCTTCCCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGc                         <  3:429372/150‑1 (MQ=255)
                                                                                                                                  gTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGGTTCGGCGAGCt                    >  7:586410/1‑151 (MQ=255)
                                                                                                                                   tACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGCGAGCtt                   >  5:866935/1‑151 (MQ=255)
                                                                                                                                   tACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGCGAGCtt                   >  3:588962/1‑151 (MQ=255)
                                                                                                                                   tACACCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGAGGTAAGGCACCGGATTCTAATTCCGTCATTCCTAGGTTCGAATCCTCTAACCCCAGCCACATTAAAAAAGCTCGCTTCGGCTAGCtt                   >  3:1055824/1‑151 (MQ=255)
                                                                                                                                    aCCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGCGAGCttt                  >  7:720731/1‑151 (MQ=255)
                                                                                                                                    aCCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGCGAGCttt                  >  5:450165/1‑151 (MQ=255)
                                                                                                                                       ccAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGCGAGCTTTTTg               >  2:449305/1‑151 (MQ=255)
                                                                                                                                       ccAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGCGAGCTTTTTg               >  5:609459/1‑151 (MQ=255)
                                                                                                                                        cAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGCGAGCTTTTTGc              >  6:180180/1‑151 (MQ=255)
                                                                                                                                         aGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGCGAGCTTTTTGCt             >  8:392693/1‑151 (MQ=255)
                                                                                                                                           ccAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGCGAGCTTTTTGCttt           >  3:462605/1‑151 (MQ=255)
                                                                                                                                              ctatcttCAAGCCGCGTCCCATGCAACAGCAACGAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGCGAGCTTTTTGCtttt          <  3:125220/144‑1 (MQ=18)
                                                                                                                                               tttATTCAAGACGCTTACCTTGTAAGTGCACACAGTTGGGGTATCGACAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGCGAGCTTTTTGCATTTCTg       >  4:134960/1‑151 (MQ=255)
                                                                                                                                                ttATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGCGAGCTTTTTGCTTTTCTGc      >  7:298912/1‑151 (MQ=255)
                                                                                                                                                 tATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGCGAGCTTTTTGCTTTTCTGCg     >  6:241204/1‑151 (MQ=255)
                                                                                                                                                 tATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGCGAGCTTTTTGCTTTTCTGCg     >  1:136996/1‑151 (MQ=255)
                                                                                                                                                    tCAAGTCGGTTACCTTGTAATTGCACTCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGCGAGCTTTTTGCTTTTCTGCGtt   <  4:177118/150‑1 (MQ=255)
                                                                                                                                                    tCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGCGAGCTTTTTGCTTTTCTGCGttc  >  6:462688/1‑151 (MQ=255)
                                                                                                                                                     cAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGCGAGCTTTTTGCTTTTCTGCGttc  >  3:456161/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAATCCCGTCGTAGCCACCAAATTCTGAATGTATCGAATATGTTCGGCAAATTCAAAACCAATTTGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTGATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCAATTTATTCAAGACGCTTACCTTGTAAGTGCACCCAGTTGGGGTATCGCCAAGCGGTAAGGCACCGGATTCTAATTCCGGCATTCCGAGGTTCGAATCCTCGTACCCCAGCCACATTAAAAAAGCTCGCTTCGGCGAGCTTTTTGCTTTTCTGCGTTC  >  DH1/3184436‑3184734

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: