Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,906,8370CG33.7% 115.6 / 27.9 90intergenic (+8/‑27)ECDH1_RS19315/ECDH1_RS19320tRNA‑Leu/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base C (45/10);  new base G (28/0);  total (77/13)
Fisher's exact test for biased strand distribution p-value = 1.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTCAAGTCCCCCCCCTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTTGATTT  >  DH1/3906690‑3906987
                                                                                                                                                   |                                                                                                                                                      
gggtcaacgccccccccTCGCACCAACAACCACGTTGATATTGCTCGCACTTTGCGAAGGGGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGat                                                                                                                                                     <  1:1130311/142‑1 (MQ=18)
         ccccccccTCGCACCAAAAAACAGGTTTAGATTGCTGGAAATGTGCGAAGGTTTAGGAAAGTGTAGACGATAAAGATTCAGGTGTTAGTGTAATTACGGACTTGGGGGTTAAAGTCGCCACCCACGCCCACAAAAAGGGGATATCagtaa                                                                                                                                             >  2:438241/1‑146 (MQ=255)
                cTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAAACGTGAAATCCAAAAAAGttac                                                                                                                                      >  4:447449/1‑147 (MQ=18)
                 tCGCACCAAAAACCACGTTGATATTGCTAGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATCATGGAGAGATGCAAAAAAGTataat                                                                                                                                    >  6:131795/1‑147 (MQ=18)
                  cGCACCAAAAACGACTTTGATATTGCTCGCAGTGGGCGAAGGTGGCGGAATTGCGAGCCGGGCTCCCTTGAGGTGTTCGGCTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATg                                                                                                                                   <  1:1301785/151‑1 (MQ=21)
                  cGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGAAGGGTACATCATAAAAAGTTAGTTg                                                                                                                                   >  5:392459/1‑151 (MQ=18)
                   gCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTGAAGTCCCCCCCCTCGCCCCAAAAGACTGTATATCAATTAAAGtttgtttc                                                                                                                                  >  1:1645763/1‑144 (MQ=18)
                   gCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTACTTACGGACGTGGGGCTTCAAGTCCCCCCCCTCGCACCAAAAAACGATAAATCAAAAAAAGttacatc                                                                                                                                   >  2:1202251/1‑144 (MQ=18)
                       caaccaaccggttgttATTTCTCGTGCTGGGGGAAGGTGGCGTAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTGCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATAAAAAAAAGTAAGATGACtgt                                                                                                                              <  5:473881/141‑1 (MQ=18)
                       cAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCGCACCAAATGACGGGAAATGAAAAAAAGTAAAATGACtgt                                                                                                                              >  1:1353295/1‑151 (MQ=18)
                                cGTTTATATTGGTGGGTCTTGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAgg                                                                                                                       <  7:634393/149‑1 (MQ=18)
                                  tttatatAGCTGGTACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTgg                                                                                                                    <  7:588275/147‑1 (MQ=255)
                                     atatTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCCCACGAGGCGATATCAAAAAAACTAAGAAGACTGTGCGAAGGTGGCggc                                                                                                                >  3:558850/1‑150 (MQ=255)
                                      aatatatgtgggtgtgCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTTTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGa                                                                                                                <  8:102864/135‑1 (MQ=21)
                                      aagacctgcccctaggCGAAGGTGGCGGAATTGGTATTCTTTCTATCTTCAGGTGTTAGTGTTCTTATGGACTTGGGGGTTCAAGTTCCCCCACTCGCACCATACTAGGATATATCATAAAAAGTAATATGACTCTGCGAAGGTGGCGGaa                                                                                                               <  2:1629147/137‑1 (MQ=12)
                                                 cTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGTGGGGATATCAAATAAAGTAAGATGACTGTTGTTTGGTGGCGGAATTgaaatacgt                                                                                                    >  3:765132/1‑143 (MQ=17)
                                                  tGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGTGGGTAAATCAAATAATGTTTGATGCGTGTTGTTTGGTGGCGGAATTGGTAgcaccc                                                                                                   >  1:471575/1‑146 (MQ=14)
                                                   cggCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTTTTTTTTTTATTATTGTGGTTGGCGTGCTAGTTCCACACCTCTAAACAATCTATGCGATATCTAATATATTACGATGAGTTTGCGAAGGTGGCGGAATTGGTAGAcgcg                                                                                                   <  6:194966/149‑1 (MQ=12)
                                                    tgCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCg                                                                                                                           <  3:101868/124‑1 (MQ=255)
                                                       gAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGAAGGGGATATCAAAAAATGTTTGTTGACTGTTGGTTGGTGGCGGAATTGCTataagttgtagc                                                                                              >  7:1429103/1‑140 (MQ=14)
                                                              gCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTGAAGTCCCCCCCCGCGCCCCAAACGAGGGGTTATCAAAAAAAGTTTGATGACTGTTGGTTGGTGGCGGAATTTGTTTACGCGGTTGCTTCAGgt                                                                                       >  7:1230158/1‑151 (MQ=18)
                                                               cGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGAGGGGAAATCATATAATGTTAGTTGACTGTGTGTAGGTGGCGGATTTGTTTGCCGGGCTATCTGCAGgtt                                                                                      >  1:1426947/1‑150 (MQ=14)
                                                               cGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGAGGCGATATCAAAAAAAGTTTGATGACTGTTCGAACGTGGCGGAATTGCTAGACTCACTATCATCAGgtt                                                                                      >  1:1434855/1‑150 (MQ=18)
                                                                  aaTTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTGAAGTCCCCCCCCGCGCCCCCAAAGAGGGGAAATCAAAAAAAGTTTGATGAATGTTGGTTGGTGGCGGAATTGTTATACGGGGTAGCATCAGGTGtt                                                                                    >  8:603422/1‑150 (MQ=17)
                                                                    tttttAGACGCTCTAGCGTCAGGTGGTACTGTGCTTACGGTCCTGGGGGTGCTACGCCCCCCCCTCGCAGCAAACGACTAGATATCAATAAAATTCATATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAg                                                                                  <  3:531354/146‑1 (MQ=16)
                                                                             gcgcTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTGAAGTCCCCCCCCACGCCCCAAAAGACGGGCTATACAAAAAAGTTAGATGACTGTTGGAATGTGGGGGCATTTGTTGCCGCGGTATCTTCAGGTTTTTGTGTCCTTACg                                                                        >  7:1486266/1‑151 (MQ=12)
                                                                              cgcTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTGAAGTCCCCCCCCCCGCCCCCAACGAGGGGATATAAAAAAAAGTTAGATGAATGTTGTTTGGTGGCGGGATTGGTATACGCGGTATCGTCAGGTGTTAGTGTACTTTCgg                                                                       >  7:18838/1‑151 (MQ=16)
                                                                              cgcTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGAGTGGATATCCAAAAATGTTTGATGCATGTGTTTTGGTGGCGGAATTTGTTTAAGGTGTAGCATCAGGTGTTAgtgtagttaggg                                                                       >  8:585403/1‑143 (MQ=9)
                                                                               gcTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCACGCACCAAAAGAGGCGAAATAAAAAAAAGTAAGATGAATGTGGGAACGTGGCGGAATTGGTATACGAGGTAGCTTCAGGTTTTAgtgtaaataagg                                                                       >  8:85048/1‑142 (MQ=21)
                                                                                  aGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCGCACCAAAAGAGGGGAAATGAAAAAAAGTAAGATGCATGTTGTTAGGTGGCGGAATTTGTATAATGGGTATCTTCAGGGGTTAgtgtaataaaggccgg                                                                   >  7:551370/1‑139 (MQ=12)
                                                                                     ttCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAAAGGTTAAATCAAATAAAGTTAGATGCCTGTGGGGAGGTGGAGGGATAGGTTGACGGGCGAGGTTGAGGTGTGAGTGTGCGTCCGGAcctttt                                                                >  2:504615/1‑146 (MQ=12)
                                                                                      tCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGAGGGGATATAAAAAAAAGTTTGTTGACTGTGGGTTGGTGGCGGCATTGGTAGACCCGCAAGCTTCACGAGTTTGTGTCGTTACGGACGTggtg                                                               >  1:1350881/1‑149 (MQ=21)
                                                                                      tCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAATAGGGTATATCAAATAAAGTTTGATGAATGATTTTTGGTGGCGGAATTAAAATAAGAGATAGGTACAGGTGGTTGTTTGCTTACGGAGGtttgg                                                               >  8:745628/1‑147 (MQ=2)
                                                                                      gcAGGTGGTCGTGCGCTTCCGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTACGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTgg                                                                 <  5:196298/148‑1 (MQ=25)
                                                                                       cAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGGGATATCAAAAAAAGTATGTTGACTGTGGGTTGGTGGCGGAATTGGTTTACGGGGTAGCTTCATGTGTGAGTGTACTTAAGGACGTggggg                                                              >  7:1329640/1‑151 (MQ=21)
                                                                                       cAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGGTAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTAATTACGGACGTGGGGGt                                                             >  1:849056/1‑151 (MQ=37)
                                                                                        accgttgaggttgcTTTCGGTCGTTGTGGTTCTAGTCCCCGCCCTCGCACTATCAGATTATATATCTATTATATTTAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGt                                                             <  2:1636641/138‑1 (MQ=9)
                                                                                           tgtTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGACGGTAGATCAAAAAAAGTTAGATGCCTGTGGGTTGGTGGCGGATTTGGTATAAGGGCTAGCTTCAGGTTTTACTGTAATTACGGCCGTGGGGGttg                                                           >  3:597606/1‑149 (MQ=16)
                                                                                            gtTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGAGGGTATTTCAAATAAAGTTAGATGAATGTTAGTAGGTGCCTGATTTGGTATACGAGGAAGCAGCATGTTGTAGTGAAAATAGGGCCGTGGGGGTTct                                                          >  7:517377/1‑149 (MQ=9)
                                                                                              tAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGGGATATCAAAAAATGTATGATGCATGTTGGATGGTGGCGGAATTGGTATAATCTGTAGCTTCAGGTGTTAGTGTAATTACGGAAGTGGGGGTTGAAGt                                                       >  5:811589/1‑151 (MQ=21)
                                                                                              tAGTGTTCTTACGGACGTGGGGGTTCAAGTACCCCCCCACGCACCAAAAGACGGGATATCAAAAAAAGTTAGATGACTGTGAGTAGGTGGCGGAATTAGTATAAGAAGTATCATCAGGAGTTAGAGTAATTAAGGAccggcgggtggatct                                                       >  8:934109/1‑137 (MQ=12)
                                                                                                gtgtTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCACGCACCAAACGAGGGGAGATGAAAAAATGTTAGATGACTGTGAGAATGTGGCGGAATTGGTATACGGGCTAGCTTCAGGTTTTTGTGTAGTTTCGGGCGTggggcgtgtaatcc                                                     >  7:1184735/1‑141 (MQ=19)
                                                                                                 tgtTCTTACGGACGTGGGGGTTAAAGTCCCCCCCCTCCCACCAAAAAACGGGAAATCAAAAAAAGTTAGATGAAAGTGGGAAAGTGCCGGAATTGATAAAAGAGATAGAAACAGGTGAAAGTGTACGTAAAGAAGTgggggatgcactcgc                                                    >  4:1660394/1‑141 (MQ=255)
                                                                                                  gtTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGGGATATCAAAAAAAGTTAGATGACTGTTATAAGCTGGCGGAATTGGTAGAATGGCTATCTTCATGTGTTAGTGTCATTACGGAAGTGGGGGTTGAAggaac                                                    >  7:1410145/1‑146 (MQ=25)
                                                                                                   ttCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGGGATATCAAAAAAAGTAAGATGAATGTTGGTAGGTGGCGGAATTGGAAGAAGCGCTAGCTTCAGGTGTTAGTGTACTAAAGGCCGTGGGGGTTGAAggcaaa                                                   >  7:215124/1‑145 (MQ=25)
                                                                                                    tCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGAGGGGAGATAAAAAAAAGTTAGATGCATGTTCGAAGGTGGCGGAATTGGTAGACGGGGGAGCTTCATGTGTTAGTGTAAATAAGGCAGTGGGGGttggaggcacccc                                                 >  8:143066/1‑140 (MQ=16)
                                                                                                      ttACGGACGTGGGGGTTCAAGTCCCCCCCCACGCCCCAAAAGAGGGGATTTCAAAAAAAGTTTGTTGCCTGTGGTAAGGTGGCGGCATTGGTATACGGGGTAGCTTCAGGGTTTAGTGTAATAAAGGCCCGGGGGGTGGAAGTcccccccc                                               >  5:1229854/1‑151 (MQ=255)
                                                                                                      ttACGGACGTGGGGGTTCAACTACCCCCCCTCGCACCAAACGAGGCGAAATCAAAAAAAGGAAGATGACTGTGCGAAAGGGGCGGAATTGGTAGACGAGCAAGCATCAGGTGTGAGTGTCCTTAAGGCCGTTGCGGTTACAGTccccaccc                                               >  7:1257209/1‑151 (MQ=25)
                                                                                                       tACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAAAAGAGGCGATATCAAAAAAATTAAGATGCATGTGGGTAGGTGGCGGAATTGGTATAAGCGGTAGCTTCAGGTGTTAGTGTCATTAAGGACCTGGGGGGTGAACTccccccccc                                              >  3:1355758/1‑150 (MQ=21)
                                                                                                        acgttcgtttgggTTCAATTCCCCCCCCTCGCACCAAACGAGGCGATATCAAATAAATTTAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTc                                             <  4:558372/141‑1 (MQ=21)
                                                                                                          ggACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGAGGGGAGATCAAAAAAAGTTAGATGAATGTGGGTAGGTGGCGGAATTGGTTGAAGGGGGAGCAGCAGGTTTGAGTGTACTTAGGCACGGGGGGGTTGAAggtaaagcgacggc                                           >  8:144993/1‑138 (MQ=16)
                                                                                                          ggACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGCGATATCAAAAAAAGTAAGATGACTGTGGGAAGGTGGCGGAATTGGTATACGCGCTAGCTTCAGGTTTTAGTGTCATTAAGGACGTGGGGGTTGAAGTccccccacgagc                                           >  3:1100876/1‑145 (MQ=25)
                                                                                                                             cccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGCAGTCCCCCCCCGCGCCCCAATTATGTTTACttgc                        >  7:22594/1‑149 (MQ=21)
                                                                                                                             cccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTcc                        >  3:695343/1‑151 (MQ=255)
                                                                                                                              ccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCt                       >  3:434486/1‑151 (MQ=255)
                                                                                                                              ccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGGCCCCCCCCCCGCCCCCATTATCTTTACTGccg                       >  7:1018170/1‑150 (MQ=255)
                                                                                                                              ccccccTCGCACCAAAAGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGTAACCCCCCACGCACCAATTAtaattgggtact                       >  4:819798/1‑140 (MQ=21)
                                                                                                                               cccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGGTAAGGTGGCGGATTTGGTAGACGCGATAGCTTCAGGTGTTAGTGTCATTCCGGACGTGGGGGTTCAAGTCCCCCCCATCGCACCAATTTTGTTTAATTCCtt                      >  5:1345508/1‑151 (MQ=21)
                                                                                                                                 cccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGTGCCCCCCCTGGCACCAATTATATTTTCTTCCgttg                    >  4:1105216/1‑150 (MQ=21)
                                                                                                                                  ccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAATTATCTTTACTTCCTTTCt                   >  6:634016/1‑151 (MQ=255)
                                                                                                                                  ccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCtttgt                   >  3:1589864/1‑149 (MQ=255)
                                                                                                                                  ccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCGCACCAATTATCTTTACGTCCTTTCt                   >  6:87872/1‑151 (MQ=255)
                                                                                                                                   cTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTTCTTCCTTTCtt                  >  1:604740/1‑151 (MQ=255)
                                                                                                                                   cTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCtt                  >  1:1452931/1‑151 (MQ=255)
                                                                                                                                    tCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCACGTCCCCCCCCCCGCACCCATTATATTTACATCCTTTCtgg                 >  2:440775/1‑149 (MQ=21)
                                                                                                                                    tCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCGTTCGTg                 >  6:457418/1‑151 (MQ=255)
                                                                                                                                    tCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCGCACCAATTATCTTTACATCCTTTCttt                 >  4:546803/1‑150 (MQ=255)
                                                                                                                                     cGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTgg                >  7:1619242/1‑150 (MQ=255)
                                                                                                                                     cGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGt                >  4:745486/1‑151 (MQ=255)
                                                                                                                                     cGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCGTTACTTCCTTTCTTGt                >  6:601075/1‑151 (MQ=255)
                                                                                                                                       cACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCCAGTGCCCCCCCGCGGACCAATTAAATTTACTTGCTTTGTTGttt              >  1:1643401/1‑151 (MQ=21)
                                                                                                                                        aCCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGtttg             >  6:1662388/1‑150 (MQ=255)
                                                                                                                                         ccAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAATATCTTTACTTCCTTTCTTGTttct            >  6:776141/1‑151 (MQ=255)
                                                                                                                                         ccAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCGCCCCCATTATGTTTACTTCCTTTATTGtttgt            >  4:445420/1‑149 (MQ=21)
                                                                                                                                         ccAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCCAGTGCCCCCCCGCGCCCCAATTATGTTTACTTCCTTTGGTGTttct            >  1:1108859/1‑151 (MQ=21)
                                                                                                                                         ccAAACGAGGCGATATCAAAAAAAGTAACATGACTGTGCGAAGGTGGCGGAATAGGAAGACGCGCTAGCTTCAGCTGTTAGTGTCATTACCGACCTCTGCGTTCAAGTCCCACCCCTCGCACCAATAATCTTTACGTCCTTTATTGttaat            >  2:301239/1‑148 (MQ=255)
                                                                                                                                           aaaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCGCCCCAATTATGTTTTAGTCGTTTCTTGtgtggtg          >  4:47395/1‑145 (MQ=18)
                                                                                                                                             aCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAATTATCTTTACTTCCTTTCTTGTTTGTTGCt        >  1:1032069/1‑151 (MQ=255)
                                                                                                                                              cGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGGCCCCCCCCGCGCACCAATTATGTTTACTGCCTTTGTTGTTTGTTccgt       >  8:757879/1‑149 (MQ=21)
                                                                                                                                              cGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATATTTAATTCCTTTATTTTTTATTACtt       >  6:366559/1‑151 (MQ=18)
                                                                                                                                              cGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGCGTTCAAGTCCCCCCCCTAGCACCAATTATCTTTACTTCCTATCTAGTTTCTTACt        >  6:815367/1‑150 (MQ=255)
                                                                                                                                               gAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTTg      >  3:435830/1‑151 (MQ=255)
                                                                                                                                               gAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTACCCCCCCAAGCACCAAATATCTTTACTTGCTTTATTGTTTATTACttt      >  4:367066/1‑150 (MQ=18)
                                                                                                                                               gAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGGCCCCCCCCGCGCACCAATTATCTTTACTGGCTTTCTTGTTTCTTCCTTg      >  3:1346616/1‑151 (MQ=255)
                                                                                                                                               gAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTGCAAGTCCCCCCCCTCGCCCCAATTATCTTTTCTTCCTTTCTTGTTTCTTCCTTg      >  4:123360/1‑151 (MQ=255)
                                                                                                                                                aGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGTGCCCCCCCGCGCCCCCATTTTGTTTTATTGCGTTATTGtttggtgctttt     >  8:915903/1‑142 (MQ=17)
                                                                                                                                                aGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAGTGCCCCCCCGCGCCCCAAATATATTTAAATCCGTTCTTGTTTGTTACTTg      >  8:82751/1‑150 (MQ=18)
                                                                                                                                                  gCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGTGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTTGAtt   >  4:787032/1‑151 (MQ=255)
                                                                                                                                                  gCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAATTTTCTTTACTTCCTTTCTTGTTTCTTCCTTGttt   >  4:680242/1‑151 (MQ=255)
                                                                                                                                                  gCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTTGAtt   >  2:728691/1‑151 (MQ=255)
                                                                                                                                                   cGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTTGAttt  >  1:230163/1‑151 (MQ=255)
                                                                                                                                                   cGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTGCCCCCCCGGGCACCAATTATATTTAATTCCTTTCTTGTTTCTTCCTTGAttt  >  7:1049913/1‑151 (MQ=255)
                                                                                                                                                   cGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGAGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTGGCACAAATTATCTTTACTTCCTTTATTGTTTCTTCATTGtttt  >  1:1203767/1‑151 (MQ=21)
                                                                                                                                                   |                                                                                                                                                      
GGTTCAAGTCCCCCCCCTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTTGATTT  >  DH1/3906690‑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: