Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,906,8370CG33.3% 104.2 / 19.6 76intergenic (+8/‑27)ECDH1_RS19315/ECDH1_RS19320tRNA‑Leu/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base C (30/16);  new base G (23/0);  total (59/17)
Fisher's exact test for biased strand distribution p-value = 6.43e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGGTTCAAGTCCCCCCCCTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTTGATTT  >  DH1/3906687‑3906987
                                                                                                                                                       |                                                                                                                                                      
ggccccgccaGTACCCCCCCTCCGACCTACAACCTCTTTTATATTCCTGTCAGTCGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGc                                                                                                                                                        <  3:1298135/142‑1 (MQ=14)
 gggcggcctcgacgcccccTCGGACCAATAACCACGTTGATATTGCTCTCCCTTGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCg                                                                                                                                                       <  7:334946/137‑1 (MQ=18)
  gtgTTCAACTCACCCGCCTCGGACCTAAAAACACTTTTATATTTCTCTCACTCGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGa                                                                                                                                                      <  8:1400770/149‑1 (MQ=17)
     cccaactccggcgccTCGCAGCAACCACGACTTTGATAGTTCTCCCCGTCGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGatat                                                                                                                                                   <  4:695692/138‑1 (MQ=14)
       caactcctgccccTCGCACGAATAAAGACGTTTATATTTCTGTCTCTTTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGatat                                                                                                                                                   <  8:1315493/140‑1 (MQ=14)
                  ccTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCCCACCAAATAACGGGACATCcaaaaatgattc                                                                                                                                      >  4:903164/1‑145 (MQ=18)
                    gcccACCAACAAACTCGTTTATATTCCTGTGTCTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGAt                                                                                                                                    <  7:1325149/148‑1 (MQ=17)
                     cGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTACTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCACGCACCAAAAGACGATAAATCAAAAAAAGTAAAAt                                                                                                                                    >  4:770864/1‑150 (MQ=18)
                      gCACCAAAAACCCCTTTTATATTGCTCGCCCTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGa                                                                                                                                  <  7:859678/151‑1 (MQ=255)
                      gCACCAAAAAACACGTTGATATTTCTCGCACTGGGCGAAGGTGGCGGTATTTGTAGACGCGCTAGCGTCAGGTGTTAGTTTTCTTACGGACGTGGGGGGTGAAAT‑GCCCCCCCCCGCCCCCAAAGAGGGGAAATCaaaaaaatttttttgg                                                                                                                                  >  4:1519511/1‑142 (MQ=255)
                       cACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCCCCCCAAAAGAGTGTAGATCCAATAAAGTatcatgc                                                                                                                                  >  4:1533768/1‑144 (MQ=18)
                       cACCAAAAACCACGTTGATATTGCTCGCACTCGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTACTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCACGCACCAAAACACGCGAGATCAAAAAAACTAAAAtca                                                                                                                                  >  4:1273028/1‑148 (MQ=18)
                        aCCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCACGCACCAAAAAAAGTTAAATAAAAAAAAGTTAAATGaat                                                                                                                                >  8:825688/1‑149 (MQ=17)
                        aCCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGGTGAAGT‑ACCCCCCCGCGCACCAAAACAGGGGATTTACAATAAAGTTTGATGACt                                                                                                                                >  7:70937/1‑151 (MQ=18)
                          cAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTAAAGT‑ACCCCCCCAACCACCAAAAAAAGATAAAAAAAATAAAGTAAGATGaattt                                                                                                                              >  2:499780/1‑147 (MQ=17)
                          cAAAAACCACGTTGAAATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCCCCAAAAAACGGTATATCAAATAAAGTTAGATGAgtgt                                                                                                                              >  8:1000255/1‑151 (MQ=18)
                           aaaaaaCACTTTGATATTGATCGGTCTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACtgtg                                                                                                                             <  3:512471/151‑1 (MQ=255)
                             aaaCCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAAAACGGGAAATCCATAAAAGTTAgttccatgtttt                                                                                                                           >  8:504880/1‑140 (MQ=17)
                               aCCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTAAAGT‑CCCCCCCCGCGCCCCAAACTAGGGTATATCAAAAAAAGTTATTTGACTTTGTGat                                                                                                                         >  7:706247/1‑150 (MQ=18)
                                  aCGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCGCGCACCAAAAGAGGGGATTTCAAAAAAAGTTAGTTGACGTTGAGAAGGt                                                                                                                      >  7:76273/1‑151 (MQ=18)
                                      tGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCCCCAAAAGTCGTGAAATCAATAAAAGTTTGTTGACtgtttgttgttggcg                                                                                                                  >  7:1322751/1‑139 (MQ=17)
                                        atatTGCTCGCAGTTTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGTGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGa                                                                                                                <  1:538662/151‑1 (MQ=255)
                                         tattattgttccTGGGCGAAGGTGGCTGAATTGGTAGACGCGCTTGCTTCAGGTGTTAGTGGTCTTACGGAGGTGGGGGTGCAAGT‑CCCCCCCCTCGCACCAACCGCGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGa                                                                                                                <  5:868941/139‑1 (MQ=21)
                                         atgatttgtatgtttgCGAAGGTGTGGTAATTTTTAGACGACATAGCTTTATGTGTAATTGTTCTTACGGACGTGGTGGTGCAATT‑CACACCCCTCGGACCAACCGATCCGATATCAAAAAAAGTAATATGACTGTGCGAAGGTGGCGGa                                                                                                                <  8:1298842/135‑1 (MQ=255)
                                          atTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGGTAAAGT‑CCCCCCCCGCGCCCCAAAACACGGGATATCAAAAATAGTTAGGTGCCTGTTGTTAGGTGGCGGAAt                                                                                                              >  1:1538775/1‑151 (MQ=17)
                                            acctccaagttagCGAAGGGCGCGGAATTGGTAGACCCTCTAGCTTCAGGTGTAAGTGTTCTTACGGACGTGTGGGTGCAAGT‑CCCTCCCCTCGCACCAAACGCGGCGAAATCAAAAAAAGTAAGCAGAGTGTGCGAAGGTGGCGGAAtt                                                                                                             <  8:1209670/138‑1 (MQ=255)
                                               tCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTAAAGTCCCCCCCCCACCCACCAAAACCGGGGCTAACAAAAAAAGTTATATCACTGTTGGAATGTTGCGCAGTTgtt                                                                                                          >  7:1325024/1‑149 (MQ=9)
                                                  atgtactcGAAGGTGGAGGAATTGGTAGACGCGGTAGCTTTTGGTTTTTTTGTAATTATTGACGTTTGGGTGCAAGT‑CCCTCGCCTGGGACTAACCGCTTCGATATCTTTTAAAGTAATTTGTGTTTGCGAAGGTGGCGGAATTGGTAGa                                                                                                       <  6:1130919/143‑1 (MQ=255)
                                                   aCTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCACGCACCAAAAAAGTAAAAAACAAAAAGAGTAATTTAACTGAAAGAAGGTCCCCGAATTGGTataa                                                                                                      >  6:19002/1‑147 (MQ=11)
                                                    cTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAACACGGGATATCAAAAAAAGTTTGATGACTGTGGGTAGGTGGCGGAATTCGTataatc                                                                                                    >  7:1097330/1‑146 (MQ=17)
                                                      gggCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCCCCAAAAAACGGGAAATAAAAAAAAGTATGATGAATGTTATATCATGGCGGAATTggattacaaa                                                                                                   >  8:1124945/1‑142 (MQ=14)
                                                          gAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAg                                                                                                                                      <  2:835283/111‑1 (MQ=255)
                                                           aaGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCACGCACCAAAACAGGCGATATCAAAAAAAGTTTGATGCATGTGGTTTGGTGGCTGAATTGGTTGAAGGGGTAGCt                                                                                             >  7:244611/1‑151 (MQ=14)
                                                            aggggccGGAAGTGGGAGACGGTGTAGCTTCAGGTGATAGTGTTCTTCCGCACGGGGTGGTTGACTT‑CCCCGCCCTCGCACCAACCTCCGATATATCAAAAAAAGTAATCTGACTCTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCtt                                                                                            <  7:353080/145‑1 (MQ=255)
                                                            aGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAAAAGGATATATCAAAAAAACGAAGATTACTGTAAGAAGGTGGCGGAATTGGTAGaacgccttggtt                                                                                            >  4:1247010/1‑140 (MQ=17)
                                                                     aaatgctacacccgcTAGCGGCAGTAGTTAGTGTTCATCGGGCCGTGGGGGTGGCAGT‑CCCCCCCCTCCCCCGACCCCACCCGATCTCAATCAAATTCATATCACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTa                                                                                   <  7:384799/139‑1 (MQ=9)
                                                                                 cgcTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAAGAGGGTATATCAAAAAAAGTTTGATGCATGTTATATGGTGCCGGAATTTATATAAGAGGAATCTTCATGTGTTActgtaaatacgc                                                                       >  5:973063/1‑143 (MQ=9)
                                                                                  gcTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCACGCCCCAAACGAGGGGATATCATAAAAAGTAAGTTGCATGTTCTAAGGTGGCGGCATTGGTAGAAGCACTAGCTGCAGGTGTGAgtgtacatatggc                                                                      >  3:375006/1‑142 (MQ=18)
                                                                                    tAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTGAAGT‑CCCCCCCCGCGCCCCAAAAGATGGGAAATCAAAAAAAGTTAGATGACAGTTGGTAGGTGGCGGAATTGCTAGACAGGCTAGCATCATGTGTTAgtgtaagaaaggacc                                                                    >  7:56849/1‑140 (MQ=255)
                                                                                    tAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGCGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTg                                                                                                                     >  3:299526/1‑102 (MQ=255)
                                                                                    tAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTg                                                                                                                     <  4:301743/102‑1 (MQ=255)
                                                                                            ggTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAGGAGGGGAGTTCAAAAAATGTTTGATGAGTGTGGGATGGTGGCGGAATTGGTTGAAGGTGTAGCTTCATGTGTTAGTGTAATTATGGCCGTGGGGGtg                                                            >  3:121610/1‑150 (MQ=16)
                                                                                                ttAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGGGATATCAAAAAAAGTTAGATGCCTGTGCTTAGGTGGCGGAATTGGTATACTGGGGATCTTCAGGTGGTAGTGTCCTTAAGGCCGTGGGGGttga                                                          >  7:490993/1‑147 (MQ=25)
                                                                                                    tgtTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGGGATTTCAAAAAAAGTAAGATGACTGTGGGAAGGTGGGGGTATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTAGGGACGTGGGGGTTGAAGTcc                                                     >  5:716741/1‑150 (MQ=37)
                                                                                                    tgtTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGGGATATCAATAATAGTTTGATGCCTGTGTTTTGGTGGCGGATTTGGTATACGGGGTAGCTTCAGGTTTTACTGTACTAAGGGCCGTGGGGGttgtaggca                                                     >  5:771163/1‑143 (MQ=16)
                                                                                                        cTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCCCCAAAGGACGGGATATCAAAAAAAGTTTGATGCGTGTGGGTAGGTGGCGGATTTGGATGCCGGGGGAGCGTCAGGTGTTCGTGTAATTAGGGCCGGGGGGGtggtaggccgcccc                                                >  2:467832/1‑138 (MQ=12)
                                                                                                         ttACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAAGAGGGGAAATCAAATAAAGTTAGATGCCTGAGAGAAGGTGGCGGAATTGGTAAACGCGGGATGGTCCGGTTTTAGTGTCCATAAGGACGTGGGGGGGGAAGTcccccccc                                               >  4:514952/1‑151 (MQ=19)
                                                                                                             ggACGTGGGGGTTCAAGT‑CCCCCCCCTCGCCCCAAACTAGGCGATATCAAAAAAAGTAAGATGACTGTTGGAAGGTGCCGGAATAGGTAGACTAGCTAGCTGCAGGTGTTAGTGTCAATACGGACGTTGGGGTGGAAGTccaccccaagcc                                           >  1:1070240/1‑146 (MQ=25)
                                                                                                              gACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAAGGAGGCGAGATCAAAAAAAGTAAGATGACTGTGGGAAGGTGGCGGAATTGGTAGAAGGGCTAGCATCAGGTGTTAGTGTCGTTACGGACGTGGGGGTTGTAGTcaacgccatggcc                                          >  3:289762/1‑139 (MQ=25)
                                                                                                                        ttCAAGT‑CCCCCCCCTCGCACCAAAGGAGGGGATATCAAAAAATGTATGATGAGTGTTGGATGGTGGCGGAATGGGTTGCAGGGGGATCTGCAGGTGTTAGTGTCGTTAGGGCCGGGGGGCTGGTAGGACCCCCCCTGGCCCCAATTAttt                                >  3:938294/1‑149 (MQ=255)
                                                                                                                         tCAAGT‑CCCCCCCCTCGCACCAAACGAGGGGATATCAAAAAAAGTAAGATGACTGTGGGTAGGTGGCTGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTTCGGACGTGGGGGTTGAAGTCCCCCCCCAGGCCCCaagtgtgtt                               >  7:1004008/1‑144 (MQ=21)
                                                                                                                                ccccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAATTATCTTttcttc                         >  7:376525/1‑151 (MQ=255)
                                                                                                                                 cccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCACCCCCAAGCCCCAATTATCTTTACTTcc                        >  1:605814/1‑151 (MQ=255)
                                                                                                                                  ccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGCTAGACGCGCTATCTTCAGGTGTTAGTGTCCTTACGGACGTGGAGGTTAAAGGACCCCCCCTAGCAACAAATATCttaaagtac                        >  2:1490655/1‑143 (MQ=21)
                                                                                                                                   cccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCtt                      >  3:726008/1‑151 (MQ=255)
                                                                                                                                    ccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTGAAGGCCCCCCCCGCGCCCCAATTATGTTTACTTCCgtt                     >  5:62700/1‑151 (MQ=21)
                                                                                                                                      ccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCt                   >  4:907825/1‑151 (MQ=255)
                                                                                                                                      ccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTAAAGTCCCCCCCCGCGCCCCAATTATGTTTAATTCCattat                   >  6:28130/1‑149 (MQ=21)
                                                                                                                                        tCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTg                 >  3:710201/1‑151 (MQ=255)
                                                                                                                                        tCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTg                 >  4:135052/1‑151 (MQ=255)
                                                                                                                                        tCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTACAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTg                 >  2:1189070/1‑151 (MQ=255)
                                                                                                                                          gCACCAAACGAGTCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTAAGGACGTGGGGGTTAAACTCCCCCCCCGCGCACCAAATATCATTACTTACTTTCTTGtt               >  8:895827/1‑151 (MQ=21)
                                                                                                                                          gCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGtt               >  5:500994/1‑151 (MQ=255)
                                                                                                                                          gCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAGTCCCCCCCCGCGCCCCCATTATCTTTACTTCCTTTCTtttt               >  8:23296/1‑151 (MQ=255)
                                                                                                                                              cAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCGCCCCCATTATCTTTACGTGCTTTCTTGTttctt           >  5:53512/1‑151 (MQ=255)
                                                                                                                                               aaaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTGCTTTCTtgtttgttg          >  6:141319/1‑150 (MQ=255)
                                                                                                                                               aaaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGGCCCCCCCCTCGCACCAATTATATTTACTTCCTTTCTTGTttattc          >  4:1487802/1‑151 (MQ=255)
                                                                                                                                               aaaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGCTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCACCCCCTCGCACCAAATATCTTTACTTCCTTTCTTGTttcttc          >  6:1260209/1‑151 (MQ=255)
                                                                                                                                                aaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGTCCCCCCCCGCGCCCCAATTATGTTTTATTACTTTGTTTTTTTTTcc         >  3:1541622/1‑151 (MQ=18)
                                                                                                                                                aaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCGCCCCCGCGCACCAATTATATTTTCATCCTTTGGTGtttggtc          >  8:552086/1‑146 (MQ=18)
                                                                                                                                                 aCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGTCCCCCCCCTCGCCCCAATTATCTTTACTTCCTTTCTTGTTTCTTGCt        >  6:1334965/1‑151 (MQ=255)
                                                                                                                                                  cGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAATTATGTTTACTTCCTTTCTTGTTTGTTCGtt       >  4:905813/1‑151 (MQ=255)
                                                                                                                                                  cGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGGTTCTTCCtt       >  4:320082/1‑151 (MQ=255)
                                                                                                                                                   gAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCCAGTCCCCCCCCGCGGCCCAATTATCTTTTCTTCCTTTCTTGTTTCTTGCttt      >  5:882345/1‑150 (MQ=21)
                                                                                                                                                   gAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTAtggtt                              >  4:244229/1‑123 (MQ=255)
                                                                                                                                                   gAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCttt                              <  3:242181/127‑1 (MQ=255)
                                                                                                                                                   gAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTGTTGTTTcttgcttg      >  6:1161481/1‑151 (MQ=255)
                                                                                                                                                       cGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCACGTCACCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTTGAttt  >  2:1102769/1‑151 (MQ=255)
                                                                                                                                                       cGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAATTATGTTTACTTCCTTTCTTGTTTGGTCCTTGAttt  >  1:988973/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
GGGGGTTCAAGTCCCCCCCCTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTTGATTT  >  DH1/3906687‑3906987

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑