Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,591,5130TG27.5% 124.3 / 33.7 170D86A (GAC→GCC) ECDH1_RS12815hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (55/48);  new base G (14/25);  total (78/92)
Fisher's exact test for biased strand distribution p-value = 8.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

CCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAG  >  DH1/2591366‑2591662
                                                                                                                                                   |                                                                                                                                                      
    aGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGACTTGAGCGCCaagctc                                                                                                                                                 >  1:927943/1‑146 (MQ=11)
             cgctcAGCTCTTAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTg                                                                                                                                        <  7:605503/151‑1 (MQ=255)
               cgccgcGCGGACGAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGcc                                                                                                                                      <  3:335621/151‑1 (MQ=255)
                      ggtacgaaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCa                                                                                                                               <  5:515566/145‑1 (MQ=255)
                       gttcgaaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCaa                                                                                                                              <  3:690822/146‑1 (MQ=255)
                          cgaaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCGCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACg                                                                                                                            <  7:237882/148‑1 (MQ=255)
                            aaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGAcc                                                                                                                         >  7:110647/1‑151 (MQ=255)
                               ccTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGa                                                                                                                      <  4:789221/151‑1 (MQ=255)
                               ccTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCCCAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATTAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGa                                                                                                                      <  7:530027/151‑1 (MQ=255)
                                cTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGa                                                                                                                      >  7:752591/1‑150 (MQ=255)
                                  tCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAg                                                                                                                    <  2:438856/150‑1 (MQ=255)
                                  tCCCCCACCACCATCACTTTCAAAAGTAACAGAAAGCTCTAGCGAATCCGCAATGAAATATTATGCCCAAGCGGTTAAGGATGAGAAGATGCGAAGAGGGTTCGGCTGGAGCGCCAGCCAGACAGCGTTGCCGCAGGCCCCGACCCGAAgg                                                                                                                   >  5:738581/1‑151 (MQ=2)
                                   cccccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGACTAGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAggg                                                                                                                  >  4:170565/1‑151 (MQ=255)
                                   cccccACCACCATCACTTTCAAAAGTCAAGGAACACTCAAGCGCCGCCGAAATCAAATAATCTGCCAATGCGGGCAAGGATGAGAAGCTTCTACGAAGATTCGTCTATAGCGGCAGCGAGCGAGAGTTGAAGCAGACAACGACCCGAAggg                                                                                                                  >  6:1151154/1‑151 (MQ=2)
                                       caccaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGa                                                                                                               >  7:169125/1‑150 (MQ=255)
                                       caccaccaTCACTTTCAAAAGGCCCTGAACGCACAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGaa                                                                                                              <  3:541517/151‑1 (MQ=14)
                                       aaccaccaTCACTTTCAAAAGTCCCTGAACTTTCAATCGAATCCGCTATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGaa                                                                                                              <  2:881335/150‑1 (MQ=255)
                                         ccaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgc                                                                                                            >  3:1194156/1‑151 (MQ=255)
                                          cccctgccgatacAAAAGGCCCTGAACTCCCCACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTACGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcg                                                                                                           <  1:815153/139‑1 (MQ=11)
                                          caccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcg                                                                                                           >  5:1030460/1‑151 (MQ=255)
                                           accaTCACTTTCAAAAGCCCCTCACCTCCCACCCCAATCCGCAATAAAATATTCTGCCCATGCGGGTAAGTATGAGAAGCTTCGACTAAGGTGCGGCTCCCGCGCCCGCGAGAGAGCGTGGCCGCAGGCAACGACCCGAAGGGCGAAgcgc                                                                                                          <  7:154488/151‑1 (MQ=11)
                                             caTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGAAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                        <  2:817808/151‑1 (MQ=255)
                                             caTCACTTTCAAAAGTCCCTGAACTCCCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTGCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                        <  5:1229133/151‑1 (MQ=255)
                                              aTCACTTTCAAAAGGCCCTGAACTCGCCAGCGAAGCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCa                                                                                                       <  5:107876/151‑1 (MQ=255)
                                                cACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgc                                                                                                     >  8:580206/1‑151 (MQ=255)
                                                cACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgc                                                                                                     >  6:150601/1‑151 (MQ=255)
                                                 aCTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGGCTGGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGgagacgaaccgaccc                                                                                                    >  5:345604/1‑140 (MQ=14)
                                                   tttCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCt                                                                                                  >  1:1017532/1‑151 (MQ=255)
                                                   tttCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGACAGCGAGAGAGCGTTGACGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCt                                                                                                  >  8:241037/1‑151 (MQ=255)
                                                    ttCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTg                                                                                                 >  1:286852/1‑151 (MQ=255)
                                                    tgcAAACGACCCTCAACTCGCAACCCAAGCCGCAATAAAATATGCGCCCCATGCGGGGAAGGATGAGACGCTTCGACTAAGGGGCGGCGCCAGCGGCCGCGAGAGCGCGTGGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCt                                                                                                  <  5:699871/148‑1 (MQ=9)
                                                     tCAAAAGTCCCTGAACTTTCAACCGAATCCTCATTTAAATTTTCTGACTTTGCGGGGAAGTATGTGAAGCTTGGAATATGGTTCGGCTGGAGCGTCAGCGTTATAGCGTTGCCGCATGCAACGACCCGAAGGGCGAAGCTCGCAGCGCTg                                                                                                 <  6:33478/150‑1 (MQ=2)
                                                      cAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGATCCGAAGGGCGAAGCGCGCAGCGCTg                                                                                                 >  3:653766/1‑149 (MQ=255)
                                                      cAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTg                                                                                                 >  6:227691/1‑149 (MQ=255)
                                                      aaaaaGTCCCTGAATTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGTCTTGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAg                                                                                               <  2:286295/150‑1 (MQ=255)
                                                       aaaaGTCCCTGTACTTTCAAGCGAATCGGCAATCAATTATTCTGCCGATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAg                                                                                               <  6:622508/150‑1 (MQ=255)
                                                       aaaaGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGt                                                                                              >  6:270435/1‑151 (MQ=255)
                                                       aaaaGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGt                                                                                              >  2:619118/1‑151 (MQ=255)
                                                       aaaaGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGt                                                                                              >  8:105366/1‑151 (MQ=255)
                                                       aaaaGTCACTGAAATTTCAATCTAATCATAAATCAAATTTTCTGCACATGCGGGGAAGGATGAGAAGCTTCGACTAAGTTTGGTCTCGAGCTTCAGATAGAGAGCGTTGCCGCAGGAAACGACCCGAAGGTCGAAGCGCGCAGCGCTGAGt                                                                                              <  4:677695/151‑1 (MQ=14)
                                                        aaaGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGt                                                                                              >  6:612975/1‑150 (MQ=255)
                                                        aaaGTCCCTGAACTCTCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTa                                                                                             <  4:748124/151‑1 (MQ=255)
                                                         aaGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTa                                                                                             >  7:333420/1‑150 (MQ=255)
                                                          aGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAAt                                                                                           >  5:401601/1‑151 (MQ=255)
                                                              ccTGAACTCCCCACCCAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAACCTTCGACTACGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCtt                                                                                       <  1:836196/151‑1 (MQ=11)
                                                                 acacttccAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTccc                                                                                    <  1:190610/144‑1 (MQ=11)
                                                                  aaCTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTcccc                                                                                   >  7:329637/1‑151 (MQ=255)
                                                                    atttcaatataaTCCTAAATTAAATATTCTTACCATGCTGGTAATTATTATTAGATTCTACTTAGGTTCGTATCTAGCTTCATGTATAGAGCTTTGCCTCAGGAAACGACCCGAAGTGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCcca                                                                                 <  2:1095626/141‑1 (MQ=2)
                                                                        caacctaaTCCCCAATCAAATATTCTGCCCATGCGGGTAAGGATGAGAAGCTTCGACAAAGGTTCGTGTTGTGCGTCAGCGAGAGAGCGTTGCCGCAGGAAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccacc                                                                               <  8:751361/143‑1 (MQ=11)
                                                                         aaGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccacc                                                                            >  4:572306/1‑151 (MQ=40)
                                                                          aGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccacc                                                                            >  1:183974/1‑150 (MQ=40)
                                                                           gCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccacca                                                                           >  2:147821/1‑150 (MQ=255)
                                                                           ccGAATCCCCAATCAAATATTCTGCCCATGCGGGGAAGGAGGACAAGCTGCGACAAAGGTGCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCAt                                                                          <  1:662408/150‑1 (MQ=11)
                                                                            cGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATc                                                                         <  8:869621/151‑1 (MQ=14)
                                                                             gAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTGGAGCGCCAGCg                                                                                                                                                   >  6:404135/1‑76 (MQ=11)
                                                                              aaTCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATcgc                                                                       >  3:753971/1‑149 (MQ=38)
                                                                                 ccGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCtt                                                                     >  8:354292/1‑150 (MQ=38)
                                                                                  cGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCCCCACCATCGCttgc                                                                   >  6:343394/1‑149 (MQ=255)
                                                                                    cAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTGGAGCGTCACCGAGAGCGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCCCGACCATcccatccat                                                                 >  3:880527/1‑143 (MQ=11)
                                                                                     aaTCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttcctt                                                                 >  3:709513/1‑146 (MQ=25)
                                                                                      aggatcttctgttctcATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCa‑aaa                                                              <  4:991050/136‑1 (MQ=11)
                                                                                      aTCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCa‑aaa                                                              <  3:843045/151‑1 (MQ=25)
                                                                                           atatTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTCCt‑aaaaaaa                                                          >  4:1110955/1‑146 (MQ=11)
                                                                                            tatTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCACCGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTCCATAAAGcacc                                                         >  3:1073683/1‑151 (MQ=11)
                                                                                              ttCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCCCCATCCCTTTcttaactccgcc                                                        >  5:519955/1‑139 (MQ=11)
                                                                                              ttCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttgcttaaagaagcc                                                        >  7:1266057/1‑137 (MQ=11)
                                                                                                  gCCCATGCGGGTAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATtc                                                  <  8:4633/151‑1 (MQ=255)
                                                                                                  gCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCTAGATAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATtc                                                  <  4:262743/151‑1 (MQ=255)
                                                                                                    ccATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTACa‑taaacaatccaaaactc                                                >  5:560186/1‑137 (MQ=11)
                                                                                                    ccATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATtctc                                                <  6:128258/151‑1 (MQ=18)
                                                                                                      aTGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATtctc                                                <  2:734811/149‑1 (MQ=255)
                                                                                                      aTGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCTCGCAGCGCTGAGTAATCCTTCCCCAACCACCATCGCTTACa‑aaaaaaatcaaaagctat                                               >  2:489777/1‑135 (MQ=11)
                                                                                                      aTGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTACt‑aaaaaaatcaaaaactga                                               >  2:966375/1‑135 (MQ=11)
                                                                                                      aTGCGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCaa                                              <  4:373060/151‑1 (MQ=25)
                                                                                                        gCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCaaa                                             <  7:667566/150‑1 (MQ=255)
                                                                                                          gggtAAGTATTAGAAGCTTCGACTATGGTTCGCCTCGAGCGTCAGCGAGAGAGGGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGa                                          <  2:1264196/151‑1 (MQ=14)
                                                                                                          ggggAAGGATGATAAGCTTCGACAAAGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGa                                          <  8:515681/151‑1 (MQ=255)
                                                                                                              aaGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATc                                       <  4:1088184/150‑1 (MQ=255)
                                                                                                                 gATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATCCGCa                                   <  5:169027/151‑1 (MQ=255)
                                                                                                                 gATGAGAAGCTTCGACCAAGGTTCGTCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATCCGCa                                   <  6:805347/151‑1 (MQ=16)
                                                                                                                   tGAGAAGCTTCGACGAAGGTTCGCCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATCCGCaa                                  <  5:781248/150‑1 (MQ=255)
                                                                                                                    gagaAGCTGCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATCCGCAAt                                 <  1:627797/150‑1 (MQ=255)
                                                                                                                     agaAGCTTCGACCAAGGTTCGTCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATCCGCAATCa                               <  2:1251307/151‑1 (MQ=16)
                                                                                                                     agaAGCTTCGACAACGGTTCGCCTCGAGCGCCCGCGTGAGAGGGTTGCCGCCGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATCCGCAATCa                               <  6:1128303/151‑1 (MQ=11)
                                                                                                                      gaAGCTTCGACTAAGGTTCGGCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATCCGCAATCaa                              <  1:415966/151‑1 (MQ=255)
                                                                                                                         gCTTCGACTAAGGTTCGTCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAata                           <  2:1089049/151‑1 (MQ=255)
                                                                                                                           cgcGACAAAGGGGCGCCGCGCGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATAtt                         <  1:1027209/149‑1 (MQ=14)
                                                                                                                            tCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTACA‑AAATAAAACAAAAACTCAAACGACTCCGCAAGCAAatata                         >  3:937620/1‑149 (MQ=14)
                                                                                                                              gACCAAGGTTCGCCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTg                      <  1:1181146/151‑1 (MQ=16)
                                                                                                                              gACAAAGGACCGGCTCGAGCGTCACCGAGAGACCGGTCCCGCCGGAAACGACCCGAACGGCGAAGCGCGGAGCGCTTAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTg                      <  3:527411/151‑1 (MQ=14)
                                                                                                                               accaaGGTGCGCCGCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCTAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGc                     <  7:1084868/148‑1 (MQ=11)
                                                                                                                                 tAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACAAACACTTTAA‑AAAGCCACGAAATTCTCAAAAGAATGAGCAATCAAATATACAGccc                   >  6:71373/1‑151 (MQ=14)
                                                                                                                                   aGGTTCGTCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCaa                 <  6:146403/151‑1 (MQ=255)
                                                                                                                                     gTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGc               <  3:1171355/151‑1 (MQ=255)
                                                                                                                                      ttCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCCCCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCg              >  5:84219/1‑151 (MQ=255)
                                                                                                                                       tCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCgg             >  5:350542/1‑151 (MQ=255)
                                                                                                                                             cGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGaa         <  4:103829/149‑1 (MQ=38)
                                                                                                                                              gAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑TAAGCACCGGAATTCTCAAACGAATCCGGAATGAAATATTCTGCCCACCCGGGGAAgc       >  7:59153/1‑149 (MQ=14)
                                                                                                                                               aGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCTGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGa          <  2:213132/146‑1 (MQ=38)
                                                                                                                                                   tCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑TAAGCCCCGGAATTCTCAAACGAATCCGCAATCAATGAGTCTGCCCACGCGGGGAAGGATgcg  >  1:438435/1‑148 (MQ=14)
                                                                                                                                                   |                                                                                                                                                      
CCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCA‑AAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAG  >  DH1/2591366‑2591662

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