Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,906,8370CG23.4% 203.5 / 17.0 122intergenic (+8/‑27)ECDH1_RS19315/ECDH1_RS19320tRNA‑Leu/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base C (62/20);  new base G (25/0);  total (97/25)
Fisher's exact test for biased strand distribution p-value = 3.23e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

GTTCAAGTCCCCCCCCTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTTGATTT  >  DH1/3906691‑3906987
                                                                                                                                                  |                                                                                                                                                      
tgtccacgactcccccTCGCTCGAATAAACTCTTTGATATTCCTGTCTGTCTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGat                                                                                                                                                     <  8:480979/139‑1 (MQ=14)
cggcaatgcacccgccTCGGTCGAACAACCTCTTTTAAATTGCTGCCTCTCGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGata                                                                                                                                                    <  3:291935/141‑1 (MQ=14)
         cccccccTCGGACCAACAACCTCGTTTTTATTCTTGCCTGTCTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCaaaaaaa                                                                                                                                           <  7:1013740/151‑1 (MQ=17)
         cccccccTCGCTCGAAAAACAACTTTTATATTCCTGCCCGTCGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCaaaaaaa                                                                                                                                           <  1:1211789/151‑1 (MQ=17)
           cccccTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGACGGGAAATCAAAAAAAGt                                                                                                                                         >  3:525821/1‑151 (MQ=255)
            ccccTCGCACCAAAAAACACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTAAATTAACACCCCACGCACCAAAAGAGGGTAGATCAAATAAAGtt                                                                                                                                        >  4:318911/1‑150 (MQ=18)
             ccctccctCGAAACACGACTTTTAAATTGCTGTCCCTCGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTaa                                                                                                                                       <  7:1015162/151‑1 (MQ=17)
             cccTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGAAGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGAGGTTATATCATTAAAAGttt                                                                                                                                       >  5:829105/1‑149 (MQ=21)
                 cGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTTGTAGACGAGCTACCTTCACTTGTAAGTGTTCTTAATGACGAGGGGGTTCAAGTACCCCCCCTGGCAAAAAAAGAGGCGCTATCACAAAAAGTAAGCTg                                                                                                                                   >  3:907913/1‑151 (MQ=255)
                  ggagcAATAACTACTTTGATATTGGTGAGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGa                                                                                                                                  <  3:754997/147‑1 (MQ=17)
                  gCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCTATATCAAAAAATGTTTGATg                                                                                                                                   >  7:287897/1‑150 (MQ=255)
                   cACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGNGGTTCAAGTCCCCCCCCNCGCANCAANAAANNTNANNTCanaanaannttgttc                                                                                                                                   >  3:2017415/1‑141 (MQ=18)
                   cACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTAAAGTCCCCCCCCTCGCACCAAAAAAGGGGATATCAAAAAAAGTTAGATGAc                                                                                                                                 >  7:409784/1‑151 (MQ=255)
                   cACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTACTTACGGACGTGGGGGTTCAACTCCCCCCCCACCCACCAAAAAACGTGATTTCAAAAAAATTataatac                                                                                                                                  >  1:711850/1‑144 (MQ=18)
                    aCCAAAAACCACGTTGATATTGCTCGAACTGGCCGAAGGTCGCGGAATTGCTAGACGCGCTACCTTCAGGTGTTAGTGTACTTACGGACGTGGGGGTTCAAGTCCCCCCACTCGCACCAAATAAGGAGAATTCaaaaattgtatcatcact                                                                                                                                >  5:1441021/1‑138 (MQ=17)
                      cAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAATAAAGCGATGTAAAAAAAAGGAAggacaggc                                                                                                                               >  1:1276136/1‑143 (MQ=18)
                         aaaCCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGTGGGGAATTCAATTATTGTTAGATGAGtgtggt                                                                                                                           >  1:1433942/1‑149 (MQ=17)
                          aaCCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACATGTAGAAATCAAAAAATGTAAGTTCCCtgtgggt                                                                                                                          >  3:1748000/1‑148 (MQ=17)
                          aaCCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTGCTTACGGACGTGGGGGTTCAAGTCCCCCCCCACCCACCAAAAGTGGATATATCAATAAAAGTTAGATgcattagat                                                                                                                           >  7:534037/1‑142 (MQ=17)
                            ccACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGGGAAATCAAATAATGTATCATGAGtgagttatt                                                                                                                        >  6:2136169/1‑144 (MQ=17)
                             cACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGAAGGATAAATCCAAAAAAGTTACATGACtgatgtttgg                                                                                                                       >  2:1006435/1‑143 (MQ=17)
                                 ttGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTACTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGACGGGAGATCAAAAAAAGTTAGATGACtgtgataacgtccc                                                                                                                   >  2:551/1‑141 (MQ=18)
                                   gcttttgctgtcccTGGGCGAAGGTGGCGGATTTGGTAGACGCGCTAGCTTAAGGTGTTAGTTTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCg                                                                                                                  <  7:1268219/137‑1 (MQ=21)
                                   gATATTGCTCGCAGTTGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCg                                                                                                                  <  3:759301/150‑1 (MQ=255)
                                   gATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCACCCACCAAAAAACGAGAAATAAAAAAAAGTAACATGAAtgatataaaatcacgc                                                                                                                 >  8:821167/1‑137 (MQ=17)
                                   gATATTGCTCCCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCg                                                                                                                  <  7:310270/150‑1 (MQ=255)
                                        ttttttatttatgCGAAGGTGGCGGAATTGGAATACGCGCTAGCTTCAGTTTTTATTGTCCTTCCGGAGGTGTGGTTTCAAGTCCCTCCCCTCGCACTAACTGCGCCGATATCAAAAAAAGTAATATGCCTGTGCGAAGGTGGCGGAAtt                                                                                                             <  6:1741072/138‑1 (MQ=12)
                                            gtcagttttgtaaGGAGGCGGAATTGGTAGACGCGCTAGCTTCATCTGGTCGTGTGCTTACGGACGTGGTCGTGCTAGTCCCCCCCCTCGCACGAACCGCGTCGATATCAAATAAAGTAATATGACTGTGCGAAGGTGGCGGAATTGGTa                                                                                                         <  3:255269/139‑1 (MQ=16)
                                             gCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCaaaaaaa                                                                                                                                           <  5:1822015/115‑1 (MQ=255)
                                             gCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCaaaaaaa                                                                                                                                           >  6:1818650/1‑115 (MQ=255)
                                             atcgtttgCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTa                                                                                                                                        <  2:418698/111‑1 (MQ=255)
                                                 tGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCGCCCCAAAACAGGAGATTTGAAAAAATGTATCATGCATGTTGTATGGTGGCGGAAtttggttcattg                                                                                                   >  5:1234851/1‑141 (MQ=9)
                                                    gCGTAGGTTTCGCAAGTGGTAGACGCGCTAGCTGCTGCTGATAGGGTTCTTCCGGACGTGGGGGTGCTATTCCCCCCCCTCGCACCAAACGACGAGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTa                                                                                                <  3:908397/151‑1 (MQ=25)
                                                     cGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGGGAGATCAAAAATAGTTAGATGCCTGTTGTATGGTGCCGGAATTAGtttcataggg                                                                                                 >  4:664432/1‑140 (MQ=14)
                                                      gAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCAAGCACCCAACAAGGGGATATACAAAAATGTTTGATGAATTTGTGAATTTGCCGGATTTGTTataatatcgat                                                                                               >  7:687090/1‑140 (MQ=2)
                                                       aaGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTACCCCCCCACCCACCAAAACAGAAGATAAGAAAATAAGTAAAAAGAGTATTAGAAAGTGTAGGAATAGGTTGTCGCCCTAga                                                                                              >  2:1455028/1‑149 (MQ=2)
                                                              cGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAAACTGAAAAACGAAAAGAGTAAGATAAGGGTAGGAAGGGGTCGGCATTGGTTCACGCGCtgggggcaggtt                                                                                      >  2:1920489/1‑140 (MQ=2)
                                                               ggAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGGTTTATCAAAAAAAGTTAGTTGGCTTTGGGTTGGTGGCGGAATTGGTAGACGCGGTAGCTTGAGgtgt                                                                                     >  7:1562636/1‑151 (MQ=18)
                                                                aaaTTCGAAGACTCTCTAGCTGAAGCTGTTCCTCTTCTTACGGAGTTGGCGGTTGTAGTTCACCGCCTCGCACTATAAGATTATATATCAAATATAGTAAGCTGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGtt                                                                                    <  6:908937/150‑1 (MQ=9)
                                                                  aTTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGCTGGCGGAATTGCAAGAAGATCTAGCTTCATGAGTTa                                                                                   >  1:384209/1‑150 (MQ=21)
                                                                            gcgcTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTAAAGTCCCCCCCCTCGCACCAAAAGAGGGGATATCAAAAAAAGTTAGATGCCTGTGGGTAGGTGGCGGATTTCATATACTCTCGAGCTTCAGGTGTTAGTGTCATTACg                                                                        >  7:924963/1‑151 (MQ=18)
                                                                               cTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGACGCTATATCAAAAAAAGTAAGATGACTGTGCTTAGGTGGCGGAATTGCAATAAGCGCAAGATTCATGTTTTAgtgaaattaaggac                                                                     >  3:651884/1‑140 (MQ=18)
                                                                               cTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCCCCCCCCAAAAAGGGGATATCAAAAAAAGTATGTTGGCTGTGGGTAGGTGCCGGAATTTGTTAACGCGGTAGCGTCATGGGTTAgagtaaatacggcc                                                                     >  3:2261314/1‑141 (MQ=12)
                                                                                    ttCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAAAGGAGGGGATATGAAAAAAAGTTAGATGACTGTTTTTAGGTGGCGGATTTTGTTTACGGGGTAGCAACATGGGTTTGTGTCCATTCGGACCTggg                                                                >  7:959529/1‑151 (MQ=12)
                                                                                          tgtTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATGGAGGGGATTTCATAAAAAGTTACTTGAGTGTTGGTTGGTGGCTGATTGGGTTTCCGGGCGAGCATCAGGTGGGAGTGTAGATAGGGCCGTGGCGGttgg                                                          >  1:188042/1‑149 (MQ=2)
                                                                                              aGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGCGATATCAAAAAAAGTTAGATGACTGTGGGAAGGTGGCGGAATTGGTAGACGCGGTAGCTTCAGGTGTTAGTGTAATTAAGGCCGTGGGGGTTCTAggc                                                      >  5:684088/1‑149 (MQ=25)
                                                                                               gtgtTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGa                                                                      <  5:63789/134‑1 (MQ=255)
                                                                                               gtgtTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCCCACCAAACGAGGCGATATCAAAAAAAGTAACATGACAGTGGGAAGGTGGCGGAATTGGTAGACGCGCTAGCGTCAGGTGTTAGTGTCCTTACGGa                                                                      >  6:63675/1‑134 (MQ=37)
                                                                                                 gtTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGAGGGGATATCAAATAAAGTTTGATGAATGTTGTAAGGTGGCGGAATTGGTATACTGGGGAGCATCATGTTTTAGTGTACTAAAGGCAGTGGGGGTTGAAggaca                                                    >  7:1729318/1‑146 (MQ=16)
                                                                                                 gtTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAAGACGCGAGATCAAAAAATGTAAGTTGACTGTTGGAAGGTGGCGGAATTGCTTGACGGGGTATCTTCAGGTGTTACTCTCTATTCGGACGTGGGGGTTGAAggcac                                                    >  5:2275958/1‑146 (MQ=19)
                                                                                                    cTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCGCACCAAAAGAGGGTATTTCAAATAAAGTTTGATGCATGTTGTTTGGTGGCGGAATTTAAATAAGGGTAAGCAGCAGGTGGTAATGTCCTTAGGGACGtggtgggtggagtggcccgg                                                >  6:96381/1‑137 (MQ=2)
                                                                                                     ttACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAAAAGAGGGGATATCAAAAAAAGTTTGATGCCTGTTGTAAGGTGCCGGAATTGGTATACGAGGTATCTTCATGTTTTAATGTAAATAAGGCCGGGGGGGTTGaatcaaaaaac                                                >  1:737930/1‑141 (MQ=12)
                                                                                                     ttACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGGGATATCAAAAAAAGTTTGATGACTGTTGTATGGTGGCGGAATTGGTATAAGGGGTATCATCATGTTTTAGTGTAATAACGGACGTGGGGGttgtagtaaaaccg                                                >  3:1007731/1‑138 (MQ=16)
                                                                                                      tACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTAAGGACGTGGGGGTTGAAGTccccccccc                                              >  3:361036/1‑150 (MQ=255)
                                                                                                           aCGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAATAGTAAGATGACTGTGGTAAGGTGGCGGAATTGGTATACTCGGGAGCTTCAGGTGTTACTGTCCTTCCGGAAGTGGGTGGTCACCTCCCCCCCCTccccc                                         >  3:1836920/1‑147 (MQ=25)
                                                                                                             gTGGGGGTTCAAGTCCCCCCCCTCGCACCAAAGGAGGCGAGATCAATAAAAGTAAGATGAGTGTTAGAAGGTGGCTGATTAGGATGCCGCAGTATCTTCATGTGTAAGTGTATAAGAGGACGTGGGGGTTGAAGTCCCAGCACTAGCAcc                                        >  4:812882/1‑150 (MQ=255)
                                                                                                             gTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGGGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCGTTACGGCCGTGGGGGTTGAAGGCCCCCCCCTCGcccc                                        >  1:427703/1‑150 (MQ=25)
                                                                                                              tGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGCCTGTTGGTAGGTGGCGGAATTGGTAGACGGGGTAGCTTCAGGTGTTAGTGTCGTTACGGACGTGGGGGTTGAAGGcccccgactggcccca                                       >  1:1792995/1‑139 (MQ=21)
                                                                                                               gggcgtgcccgcccccccccccGCACCAAAATTTCACATATAAAAAAAATTACGCTCACGGGGGTAATGTGGGGGTATTTGTAGACGGGGGATCTTCAGTTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAt                                     <  7:2035647/139‑1 (MQ=2)
                                                                                                                gcgggtCACGTCCCCCCCCTCGCACCACACCAGCCGAAATCCAACATCCTCATATGACTGTGCCCAGGTGCCGTCAGTGGTAGAGGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAtt                                    <  7:1234645/146‑1 (MQ=9)
                                                                                                                       ccgcccccccccccccACCCACCCACCCCATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTCGCAGACGCGCTAGCTTCAGCTCTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTa                             <  3:357487/147‑1 (MQ=17)
                                                                                                                       aaGTCCCCCCCCTCGCACCAAACGAGGCGATATCAATAAAAGTAAGTTGACTGTGCGAAGGTGGCGGAATTGGTAGACTCGCGAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGTCCCCCCCCTCGCACCAATTATGtttg                             >  5:1347393/1‑150 (MQ=21)
                                                                                                                       aaGTCCCCCCCCTCGCACCAAAAGATTATATTTTTTAATTATTCTGATCACTGTGAGAAGTTGTCTGAATTTGTATATGATTTAGTTTAAGTTGTTAGTGTGATTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTa                             <  8:851517/151‑1 (MQ=9)
                                                                                                                        aGTCCCCCCCCTCGCACCAAAGGAGGGGATATCAAAAAAAGTAAGATCCGTGTGGGATGGTGGCGGAATTGGTTTACTCTGTAGCTTCAGGTGTTAGTGTCGTTTAGGACGTGGGGGTTGAAGTGCCCCCCATAGCCCCAAGTATATTTa                             >  4:167122/1‑150 (MQ=16)
                                                                                                                           ccccccccTCGCACCAAACGAGGGGATATCAAAAAAAGTAAGATGCCTGTGCGAAGGTGGCGGAATTGGTAGAAGCGCGAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCGAGGCCCCCCCCTCGCCCCAAGTCTGTTTACTTc                         >  4:927871/1‑151 (MQ=21)
                                                                                                                           ccccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGGCCCCCCCCGCGCCCCAATTTTTTTTTCttg                         >  1:753396/1‑150 (MQ=21)
                                                                                                                            cccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAATTATCTTTACTTcc                        >  1:1638986/1‑151 (MQ=255)
                                                                                                                            cccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTcc                        >  3:121987/1‑151 (MQ=255)
                                                                                                                            cccccccTCGCACCAAACGAGGCGAGATCAAAAAAAGTATGATGACTGTGCGAACCTGGCGGAATTGTTAGACGCGATAGCTTCAGGTGTTAGTGTCCTTAAGCACGTGCGGGTTCAAGACCCCGCACTGGCCCCAATTAAGTTTACTGcc                        >  6:152575/1‑151 (MQ=21)
                                                                                                                             ccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAATTATCTTTACTTCCt                       >  3:1303551/1‑151 (MQ=255)
                                                                                                                             ccccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCt                       >  2:2139038/1‑151 (MQ=255)
                                                                                                                              cccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCCCCAATTATCTTTACTTCCtt                      >  1:1323874/1‑151 (MQ=255)
                                                                                                                               ccccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCGCCCCAATTATCTTTTCTTGCttt                     >  7:568846/1‑151 (MQ=255)
                                                                                                                                cccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGCAGGAATTGGTAGACGCGCTAGCATCAGCTCTTAGTGTCCTTACGGACGTGGGGGTTCAATTACCCCCCCACGCACCAATTATCTTTAATTACttt                     >  2:2039867/1‑150 (MQ=21)
                                                                                                                                 ccTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGGAGTCCCCCCCCCCGCCCCAATTATATTTAAttgcgtttt                   >  7:836529/1‑144 (MQ=21)
                                                                                                                                   tCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTGAAGTCCCCCCCCGCGGACCAATTATGTTTACTTGCTTTCTTg                 >  8:299416/1‑151 (MQ=21)
                                                                                                                                   tCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACCGGCTAGCTTCAGGAGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCATTACTTCCTTTATTg                 >  6:694768/1‑151 (MQ=255)
                                                                                                                                     gCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTgttgtt               >  2:1053866/1‑151 (MQ=255)
                                                                                                                                      cACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCGCGCCCCAATTATCTTTACTTCCTTTCTTGttt              >  5:487464/1‑151 (MQ=255)
                                                                                                                                      cACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTGAAGGCCCCCCCCGCGCACCCAATATATTTAAGTACGTTGTTGttt              >  7:2035717/1‑151 (MQ=18)
                                                                                                                                      cACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTCGCGGAATTGGTAGACCCGCTAGCTTCAGGTCTTAGTCTCCTTACGGACTTCGGGCTTCAAGACACCCCCCACGCAAAAATTATATTTaaataaataaatgatt              >  2:1973223/1‑136 (MQ=255)
                                                                                                                                       aCCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTGCAGTGCCCCCCCGCGCCCCAATTTTATTTAATTGCGTTGttttttt             >  7:504128/1‑150 (MQ=18)
                                                                                                                                        ccAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTGCCCCCCCTCCCCCCAATTATCTTTTCTTGCTTTCTTGTttct            >  3:1858028/1‑151 (MQ=255)
                                                                                                                                         cAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTGCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTttctt           >  7:1535385/1‑151 (MQ=255)
                                                                                                                                         cAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTtgtttgtt           >  6:1179507/1‑151 (MQ=255)
                                                                                                                                         cAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTttctt           >  7:1535387/1‑151 (MQ=255)
                                                                                                                                         cAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTttctt           >  3:1598602/1‑151 (MQ=255)
                                                                                                                                         cAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAGTCCCCCCCCGCGCCCCAATTATCTTTACTTCCTTTGTtgtttgtt           >  8:1190572/1‑151 (MQ=255)
                                                                                                                                          aaaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGTCCCCCCCCGCGCCCCAATTATCTTTACTTCCTTTGTTGTttattt          >  1:1112014/1‑150 (MQ=255)
                                                                                                                                          aaaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTTTttcttg          >  4:206072/1‑150 (MQ=255)
                                                                                                                                          aaaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTttattc          >  4:869922/1‑151 (MQ=255)
                                                                                                                                           aaCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTGTTGTTTCTTcc         >  5:632492/1‑151 (MQ=255)
                                                                                                                                            aCGAGGCGATATGATAAATAGTAAGTTGACTGGTCCACGGTGGCCTTCCTCGGTGACGCGCTAGCTTCAGGTGTTAGTGTGCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCt        <  7:197324/151‑1 (MQ=18)
                                                                                                                                            aCGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCt        >  7:2131253/1‑151 (MQ=255)
                                                                                                                                              gAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTGCCCCCCCTCGCCCCAATTATCTTTACTTCCTTTCTTGTTTGTTCCTTg      >  2:1183154/1‑151 (MQ=255)
                                                                                                                                               aGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTGAAGTGCCCCCCCTCGCCCCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTTgt     >  1:2161102/1‑150 (MQ=255)
                                                                                                                                               aGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGAGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGGACCAATTAAGTTTACTTCCATTTGTGTTTATTCCttc      >  1:498988/1‑149 (MQ=18)
                                                                                                                                                ggCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTGCCCCCCCTCGCACCAATTATCTTTACTTGCTTTCTTGTTTGTTCCTTGAt    >  5:1763681/1‑151 (MQ=255)
                                                                                                                                                  cGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGTTGTTAGTGTCCTTACGGACGTGGGGGTTAAACTCCCCCCCCTCGCACCAATAATATTTAATTCCTTTATTGTTTATTAGTTGAttt  >  1:1807506/1‑151 (MQ=18)
                                                                                                                                                  cGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTTGAttt  >  1:2095923/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GTTCAAGTCCCCCCCCTCGCACCAAAAACCACGTTGATATTGCTCGCACTGGGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTTCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAAACGAGGCGATATCAAAAAAAGTAAGATGACTGTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCAATTATCTTTACTTCCTTTCTTGTTTCTTCCTTGATTT  >  DH1/3906691‑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: ‑