Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,591,5130TG26.1% 303.5 / 57.2 309D86A (GAC→GCC) ECDH1_RS12815hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (103/90);  new base G (27/41);  total (151/158)
Fisher's exact test for biased strand distribution p-value = 6.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 310 total aligned reads displayed.

CCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAG  >  DH1/2591365‑2591662
                                                                                                                                                    |                                                                                                                                                     
    aaGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCga                                                                                                                                                 <  2:140286/151‑1 (MQ=255)
        gCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTCGAGCGTCAGCgag                                                                                                                                                >  7:1178131/1‑148 (MQ=255)
            aGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCgag                                                                                                                                                >  7:544982/1‑144 (MQ=255)
             gcgcgcAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGACTTGAGCTTCAGCCAGAGAGCGtt                                                                                                                                        >  3:2094230/1‑151 (MQ=255)
             cgtcttcgatgggtctAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGt                                                                                                                                         <  2:1045300/135‑1 (MQ=255)
              cgcgcAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGACCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGACTAGAGCGGCAGCGAGAGAGCGTTg                                                                                                                                       >  2:1271039/1‑151 (MQ=11)
                   aGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCa                                                                                                                                  <  1:1970236/151‑1 (MQ=255)
                      gcTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTGGAGCGGCAGCGAGAGAGCGTTGCCGCAGGc                                                                                                                               >  5:1386421/1‑151 (MQ=255)
                         gAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGCCTGGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAAc                                                                                                                            >  7:1598080/1‑151 (MQ=255)
                         gAGTAATCCTTCCCCCACCAACATCACATTCAAAACTACCTGAACAATCAAGCGAATCCGAAATCAAATATTCTGGCAATGCGGTTAAGGATGAGAAGCTTCGACTAAGGTTGGGGTCGAGCGTCAGCGAGAGAGCGTTGCCGCGGGCagc                                                                                                                            >  6:2156990/1‑149 (MQ=17)
                              ataattaccccACCACCATCACTTAAAAAAGTACCGGAATTCGCAAGCGAATCCCAAACGTAATATCCTCCCAAGGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGAccc                                                                                                                       <  1:773991/144‑1 (MQ=11)
                              aTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTGGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGAccc                                                                                                                       >  5:1739569/1‑151 (MQ=255)
                                ccTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCGCAAGCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGa                                                                                                                     <  1:410716/151‑1 (MQ=255)
                                   tCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAACCTAATCCGAAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCTACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAg                                                                                                                   <  8:255976/150‑1 (MQ=255)
                                   tCCCCCACCACCATCACTTTAAAAAGTCCCTTAACTATCAAGCTAATCCGAAATAAAATATTCTGCCCATGCGGGTAAGTATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGaa                                                                                                                    <  6:416701/149‑1 (MQ=14)
                                   cccccccccACCATCACTTTCAAAAGTCCCAGACCTCTCAAGCTAATCCGCAATCAATTATGGTCCCCATGGGGGGAAGGAGTAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAgg                                                                                                                  <  5:746812/150‑1 (MQ=14)
                                    cccccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCTTGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAgg                                                                                                                  <  4:404538/150‑1 (MQ=255)
                                     ccccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGACTAGAGCGCCAGCGAGAGAGCGTTGCCg                                                                                                                                    >  6:90480/1‑131 (MQ=255)
                                     ccccACCACCATCACTTTCAAAAGTCCCTGAACTATCAAGCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCaa                                                                                                                             <  2:1833744/138‑1 (MQ=255)
                                     ccccACCACCATCACTTTCAAAAGTCAGTGAAATCTCAAGCGAAACCGAAATCAAATATTCTGCCCAAGAGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTAGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGc                                                                                                                >  2:65317/1‑151 (MQ=11)
                                     acgcactacctttcttttCAAAAGTCCCTGAACTCACAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGc                                                                                                                <  8:1593677/136‑1 (MQ=14)
                                       ccaccaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCg                                                                                                               >  4:821942/1‑150 (MQ=255)
                                         accaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAg                                                                                                            <  8:1882815/151‑1 (MQ=255)
                                          ccaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTGGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGAcc                                                                                                                        >  2:810569/1‑138 (MQ=255)
                                           caccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcg                                                                                                          >  3:2257450/1‑151 (MQ=255)
                                            cccaTCACTTTCAAACGTCCCTCAACTCCCACCCCAACCCCCAATAAAATATTCGGCCCATGCGGGGACGGATGAGAAGCTTCGACTACGGTGCGGCGCGCGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgc                                                                                                         <  3:2173128/150‑1 (MQ=11)
                                            acgaTCTCTTTCAAAAGTCCCTGAACTTTCAATCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTATGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcg                                                                                                          <  2:1968516/147‑1 (MQ=11)
                                             ccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcg                                                                                                        >  4:1387942/1‑151 (MQ=255)
                                              caTCACTTTCAAAAGTCCCTGAACTCCCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                       <  3:1068721/151‑1 (MQ=255)
                                               aTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                       >  6:1046915/1‑150 (MQ=255)
                                               aTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                       >  3:1497680/1‑150 (MQ=255)
                                                tttatttCAAAAGTCCCTCAACTCTCATGCTAATCGGCAATCAAATTTTTTGCCAATCCGTGGAAGTAGGACAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCa                                                                                                      <  4:1078894/146‑1 (MQ=14)
                                                 aaCTTTCAAAAGTCCATTAACTTTCAATCGAATCCGCAATTAAATATTCTGCCCTTGCGGGGAAGGATGATATGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgc                                                                                                    <  8:2132999/150‑1 (MQ=11)
                                                  aCTTTTAAAATTCCCTTAACTTGCAAACGAATTCGCAATTAAATATTCTGCCCATGGGGGGTAGGATGTTAAGCTTCGACTTAGGTTCGGCTCGATCGTCAGCGATATTGCGTTGCCGCATGCAACGACCCGAAGGGCGAAGCGCGCAgcg                                                                                                   <  2:2188426/151‑1 (MQ=11)
                                                  aCTTTCAAATGTCCCTGAACTATCAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGATGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCTAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcg                                                                                                        <  2:275602/146‑1 (MQ=14)
                                                   cTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgcgc                                                                                                  >  7:1487948/1‑151 (MQ=255)
                                                   cTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgcgc                                                                                                  >  5:411209/1‑151 (MQ=255)
                                                   cTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTCGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgcgc                                                                                                  >  8:2199147/1‑151 (MQ=255)
                                                   atttAAAAAGTCCCTGAACTTTCAATCGAATCCGCAATAAAATATTCTGCCCATGCGGGTAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgcgc                                                                                                  <  8:665392/150‑1 (MQ=11)
                                                    tttCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgcgg                                                                                                  >  8:537083/1‑149 (MQ=255)
                                                     ttCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgcgc                                                                                                  <  4:1946084/149‑1 (MQ=255)
                                                      tCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTg                                                                                                >  3:289725/1‑150 (MQ=255)
                                                      tCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTg                                                                                                >  5:1481268/1‑150 (MQ=255)
                                                             tCCCTGAACTTTCCACCGANTCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCa                                                                                                                                                     <  6:326372/90‑1 (MQ=11)
                                                             tCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTAGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATc                                                                                         >  8:1111601/1‑150 (MQ=255)
                                                             tCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATcc                                                                                        >  6:1553721/1‑151 (MQ=255)
                                                              cccTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCt                                                                                       <  4:1766098/151‑1 (MQ=255)
                                                              cccTGAACTCCCCACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTGCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCt                                                                                       <  3:553998/151‑1 (MQ=255)
                                                                  gAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGACTTGAGCGGCAGCGAGAGAGCGTTGcc                                                                                                                                     >  4:1461466/1‑101 (MQ=11)
                                                                  gAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGACTAGAGCGCCAGCGAGAGAGCg                                                                                                                                          >  5:2105623/1‑96 (MQ=11)
                                                                  gAACTCTCAACCGAATCCGCATTCAAATATTCTGCCCATGCGGGTATGGATGATAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTcc                                                                                    <  2:433271/150‑1 (MQ=255)
                                                                    aCTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTccccc                                                                                 >  5:797366/1‑151 (MQ=255)
                                                                     ctctCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTcccccc                                                                                >  1:1002978/1‑150 (MQ=255)
                                                                      tctcAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTaagctcggaagcgc                                                                               >  7:615855/1‑139 (MQ=255)
                                                                      tctcAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccac                                                                               >  6:772434/1‑151 (MQ=255)
                                                                        ccaagccaccccccAACCAAATCTGCTCCCCAGGCGGGGAAGCAGGAGACCCTCCGACTAAGTTTCCGCGCGACCCGCCCCGAGACAGCGGGCCCGCCTGAAACGACCCGAACGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccacca                                                                             <  5:2157798/138‑1 (MQ=2)
                                                                         ccaacGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGAcc                                                                                                                        <  2:1582851/103‑1 (MQ=9)
                                                                         cAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccac                                                                            >  7:273920/1‑151 (MQ=40)
                                                                          aaacGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGc                                                                                                                                      <  4:158064/89‑1 (MQ=16)
                                                                           accgatcccGCAATAAAATATTCTGCCCATGCGGGGAAGCATGAGAAGCTTCGACTACGGTTCGCCGCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccac                                                                            <  3:388780/142‑1 (MQ=11)
                                                                           aacGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGc                                                                                                                                           <  4:1256450/84‑1 (MQ=19)
                                                                           aGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTccc                                                                                   >  3:1831979/1‑142 (MQ=38)
                                                                             cGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATc                                                                        <  5:1070341/151‑1 (MQ=9)
                                                                              gAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAACCGCGCAGCGCTGAGTAATCCttacaacggaacaaac                                                                        >  8:1031258/1‑136 (MQ=14)
                                                                               aaTCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATccc                                                                      >  8:1414051/1‑149 (MQ=38)
                                                                                aTCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCt                                                                     >  7:1769124/1‑151 (MQ=38)
                                                                                       ataaaaTATTCTGCCCTTGGGGGTAAGGTTGATAAGCTTCGACTACGGTTCGGCTCGAGGGTCTGTTATTGAGCGTGGCCGCTGGAAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCaaaa                                                              <  2:366933/148‑1 (MQ=11)
                                                                                         gtgctcttccgcgcttGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAg                                                             <  1:1641480/135‑1 (MQ=11)
                                                                                         cAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGc                                                            >  1:1481462/1‑151 (MQ=37)
                                                                                          aaaTATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCttccttaaat                                                            >  2:1655277/1‑142 (MQ=255)
                                                                                             tatTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttcc                                                                  >  5:1182487/1‑139 (MQ=25)
                                                                                              atTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCg                                                        <  7:2093392/150‑1 (MQ=2)
                                                                                                tCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGCCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGa                                                      <  8:205397/150‑1 (MQ=9)
                                                                                                   gCCCATGCGGGGAAGGATGTGAAGGTTCTACTATGGTTCGTCTCGAGCGTCTGCGATAGAGCGTTGCCGCAGGAAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAAtt                                                   <  8:396450/150‑1 (MQ=14)
                                                                                                   gCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGaa                                                     <  6:1931021/148‑1 (MQ=40)
                                                                                                   gCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGaa                                                     <  3:873868/148‑1 (MQ=38)
                                                                                                      cATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTGGAGCGGCAGCGAGAGAGc                                                                                                                                           >  1:1362752/1‑59 (MQ=255)
                                                                                                      cATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCt                                                                     >  6:2242206/1‑129 (MQ=37)
                                                                                                      cATGCGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGACTCGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATtctc                                                >  4:555146/1‑150 (MQ=9)
                                                                                                         gCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCGTCACAGACAAAAGGCCCGGAATTCTCAAAc                                            >  8:821034/1‑151 (MQ=255)
                                                                                                          ggggtAAGTATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACg                                           <  7:2168987/150‑1 (MQ=255)
                                                                                                          ggggtAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACg                                           <  6:1717880/150‑1 (MQ=255)
                                                                                                            gggAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGaa                                         <  1:928032/151‑1 (MQ=255)
                                                                                                             ggAAGGATGAGAAGCTTCGACTAAGGTGCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAAt                                        <  1:2158606/151‑1 (MQ=255)
                                                                                                                 ggATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGc                                    <  1:1059602/151‑1 (MQ=255)
                                                                                                                   aTGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGACTAATCCTTCCCCCCCCACCATCACATTCAAAAGCCCCGGAATTCTCAAACGAATCCTCa                                   >  5:2274423/1‑150 (MQ=255)
                                                                                                                     gcgatgtttCGACAAAGGTTCGCCGCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATc                                <  3:725369/144‑1 (MQ=11)
                                                                                                                     gagaAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAAt                                 <  2:2026635/150‑1 (MQ=255)
                                                                                                                      agaAGCTTCGACTAAGGTGCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCa                               <  5:868810/151‑1 (MQ=255)
                                                                                                                         aGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAat                            >  7:1292159/1‑151 (MQ=255)
                                                                                                                          gCTTCGTCTACGGTTGGGCGGGCGGGTCCGCGAGAGAGCGTTGCCGCCGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAata                           <  1:1888999/151‑1 (MQ=14)
                                                                                                                            ttCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCa                               <  4:2044850/145‑1 (MQ=255)
                                                                                                                            ttCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCa                               >  3:2045632/1‑145 (MQ=255)
                                                                                                                             tCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATACGCAATCAAATATTc                        >  1:1584684/1‑151 (MQ=255)
                                                                                                                               gACTATGGTTCGGCTCGCGCGCCAGCGAGTGAGCGTTGCGGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTg                      <  7:1409613/151‑1 (MQ=255)
                                                                                                                               gACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTg                      >  4:298557/1‑151 (MQ=255)
                                                                                                                               gACCAAGGTTCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCgagatctcaagatctg                      >  8:1311277/1‑136 (MQ=11)
                                                                                                                                accaaGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGc                     <  7:1885128/148‑1 (MQ=255)
                                                                                                                                   aaGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTTCTAAAAACACCGAAAGCTCAAACGCAGCCGCCATGACAAATTCTGcgc                   >  7:861231/1‑148 (MQ=11)
                                                                                                                                                    tCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTGCTTAAGGACGGCAATGCTGAAACGCCTCGGGAATGAAGGATTCTGGCCAAGCGGGGaacgcggcg  >  1:250809/1‑143 (MQ=2)
                                                                                                                                                    |                                                                                                                                                     
CCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAG  >  DH1/2591365‑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: ‑