Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,591,5130TG27.7% 245.6 / 49.3 249D86A (GAC→GCC) ECDH1_RS12815hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (93/53);  new base G (21/35);  total (137/112)
Fisher's exact test for biased strand distribution p-value = 8.99e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

CCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAG  >  DH1/2591365‑2591662
                                                                                                                                                    |                                                                                                                                                      
  cGAACGGGGAATCGCGTACCGATCAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGc                                                                                                                                                    <  7:1222646/151‑1 (MQ=17)
        gCGAAGCCCGGAGCCCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCgagaga                                                                                                                                              <  3:1383452/151‑1 (MQ=255)
         cGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTAGAGCGCCAGCgagag                                                                                                                                               >  2:1219843/1‑149 (MQ=255)
                     aaggttcgaaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGc                                                                                                                                            <  7:1626738/132‑1 (MQ=255)
                      gcTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGCCTGGAGCGGCAGCGAGAGAGCGTTGCCGCAGGc                                                                                                                                >  8:541522/1‑151 (MQ=255)
                       ggttctAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTATCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCa                                                                                                                               <  4:447735/146‑1 (MQ=255)
                       gggtcgaaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCa                                                                                                                               <  5:1087042/145‑1 (MQ=255)
                            tAATCCTTCCCCCACCAGCATCTATTTCAAAAGTCCCTGAACTGTCAAGCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGa                                                                                                                           <  2:930762/150‑1 (MQ=255)
                              aTCCTTCCCCCACAACCATCACTTTCAAAAGTCCCTGAAATCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTGGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGAccc                                                                                                                        >  4:58441/1‑151 (MQ=255)
                                ccTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTATCAACCTAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCg                                                                                                                       <  2:1178601/150‑1 (MQ=255)
                                   tCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGaacc                                                                                                                   >  1:1535189/1‑149 (MQ=255)
                                   tCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGACTAGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAgg                                                                                                                   >  2:866840/1‑151 (MQ=255)
                                    cccccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCATCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAggg                                                                                                                  >  7:1063641/1‑151 (MQ=255)
                                    cccccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAggg                                                                                                                  <  2:674520/151‑1 (MQ=255)
                                    cccccACCACCATCACTTTCAAAAGTCCCTGAACTATCAATCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAgg                                                                                                                   <  6:560188/150‑1 (MQ=255)
                                       ccaccaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCg                                                                                                                >  8:1014738/1‑150 (MQ=255)
                                       ccaccaccaTCAATTTCAAAAGTCCATGAACTTACAATCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGTATTAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGATAGAGCGTTGCCGCAGGAAACGACCCGAAGGGCGa                                                                                                               <  8:1659496/151‑1 (MQ=11)
                                           cagattCAGTTTCAAACGGCCGTGACCTCCCAAGCGAAGCGCCATGCAATTATTGTGCCCCTGCGGGGATGGCTTAGAAGCTGCGACAACGGTTCGCCTCGCGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcg                                                                                                           <  3:786932/146‑1 (MQ=2)
                                              caTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                        >  7:1578817/1‑151 (MQ=255)
                                              caTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                        >  3:673297/1‑151 (MQ=255)
                                              caTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                        >  6:312221/1‑151 (MQ=255)
                                                tCACTTTCAAAATTCGATTAACTATCAAGCTAATCCGCAATCAATTATTCTGCCCTTGGGGGGAAGGTTTAGAAGATTGGACTATGGTTCGGCTCGAGGGTCAGCGAGAGTGCGTTGCGGCTGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                        <  8:479258/149‑1 (MQ=14)
                                                  aCTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTGGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgc                                                                                                     >  4:1182482/1‑150 (MQ=255)
                                                   cTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgcgc                                                                                                   >  5:452497/1‑151 (MQ=255)
                                                    tttCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCCCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCt                                                                                                  >  8:92090/1‑151 (MQ=255)
                                                      tCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTg                                                                                                 >  1:680082/1‑150 (MQ=255)
                                                      tCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGACTTGAGCGCCa                                                                                                                                                      >  2:628305/1‑97 (MQ=11)
                                                        aaatttctcTGAACTTACAACCTAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGAAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAg                                                                                               <  8:126763/142‑1 (MQ=11)
                                                         aaaGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGGCTAGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgc                                                                                                            >  2:1588483/1‑136 (MQ=255)
                                                          aaGTCCCTTAACTCTCAACCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTaa                                                                                            <  8:838696/151‑1 (MQ=255)
                                                           gttctcTGAACTCCCAACCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAAt                                                                                           <  6:1435158/146‑1 (MQ=11)
                                                           cgTCCCTGAACTCCCAACCGAAGCCCCAATAAAATATTCTGCCCATGCGGGGAAGGAGGAGACGCTGCGACTAAGGTGCGGCTCGAGCGTCAGCGAGAGAGCGTGGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAAt                                                                                           <  3:1033248/150‑1 (MQ=14)
                                                           aGTCCCTGATCACTCACGCCAATCCGCAATCAAATATTCTGCCCATCCGGGTAAGGATGAGAAGCTTCTACTAAGCTTCGCCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAAt                                                                                           <  6:503570/151‑1 (MQ=11)
                                                             tCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATcc                                                                                         >  5:1405130/1‑151 (MQ=255)
                                                             gcccTGAACTCGCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATcc                                                                                         <  5:194239/150‑1 (MQ=255)
                                                              cccTGAACTCTCAAGCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCt                                                                                        <  4:302059/151‑1 (MQ=255)
                                                                    aCTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTccccc                                                                                  >  5:784791/1‑151 (MQ=255)
                                                                    aCTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTccccc                                                                                  >  3:947591/1‑151 (MQ=255)
                                                                     cggagcaccGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTGCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCcca                                                                                 <  1:320207/143‑1 (MQ=11)
                                                                     atttcAATCGAATCGTCAATCAAATATTCTGCCCATGGGGGGATGTTTTTTAATCTTCTACTATGGTTCGGGTGTAGCGTCAGCTAGAGTGCGTTGCGTCAGGCAAGTACCCGAAGGGCGAAGCGCGCAGCTCTGAGTAATCCTTCCCcca                                                                                 <  8:1584288/148‑1 (MQ=9)
                                                                        tcAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccacc                                                                               >  5:738624/1‑150 (MQ=255)
                                                                        tcAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccacc                                                                               >  1:97465/1‑150 (MQ=255)
                                                                         cAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccac                                                                             >  2:1532802/1‑151 (MQ=39)
                                                                          aaGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccacc                                                                            >  4:235185/1‑151 (MQ=40)
                                                                          aaGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccacc                                                                            >  6:1622597/1‑151 (MQ=255)
                                                                            acGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCgag                                                                                                                                                 <  4:786968/79‑1 (MQ=255)
                                                                             gcctttctcAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGGGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCAt                                                                          <  4:138541/142‑1 (MQ=11)
                                                                             cGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATc                                                                         >  1:1561836/1‑151 (MQ=38)
                                                                              gAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATcc                                                                        >  2:588725/1‑150 (MQ=38)
                                                                               tctCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTTGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATcg                                                                        <  6:1029888/148‑2 (MQ=16)
                                                                               aaTCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGACTAGAGCGTCAGCGAGAGAGCGTTGACGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATcgc                                                                       >  8:1351592/1‑149 (MQ=12)
                                                                                aTCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCt                                                                      >  7:406866/1‑151 (MQ=38)
                                                                                    gCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttcct                                                                  >  8:303461/1‑148 (MQ=25)
                                                                                     cAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTCCa                                                                 >  2:1669751/1‑150 (MQ=25)
                                                                                      aaTCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttgctt                                                                 >  3:938802/1‑146 (MQ=255)
                                                                                       aTCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATACTTCCCACACAACCATCAATTACaaa                                                               >  8:191256/1‑150 (MQ=255)
                                                                                         ataaTTTTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGc                                                            <  1:1224452/149‑1 (MQ=17)
                                                                                          aaaTATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttacttaaat                                                             >  6:1394382/1‑142 (MQ=255)
                                                                                          aaaTATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTGCATAAga                                                            >  3:516435/1‑149 (MQ=255)
                                                                                             tatTCTGCCCATGCGGAGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCttccttaaacatcg                                                         >  8:251568/1‑139 (MQ=11)
                                                                                              atTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGc                                                                                                                                            <  8:1606416/67‑1 (MQ=255)
                                                                                              atTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttccttaaagaaac                                                         >  5:1589645/1‑138 (MQ=255)
                                                                                              atTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCg                                                        <  4:969003/150‑1 (MQ=37)
                                                                                               ttCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttccttaaagactga                                                        >  3:896744/1‑137 (MQ=11)
                                                                                               ttCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTCCataaataaaga                                                       >  1:1003922/1‑143 (MQ=11)
                                                                                                 cagCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTCCataaagacacaaa                                                     >  1:229033/3‑141 (MQ=11)
                                                                                                  tGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTACagaaataaacaaaa                                                    >  6:1256523/1‑140 (MQ=11)
                                                                                                   gCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttcctta                                                                >  1:275011/1‑133 (MQ=255)
                                                                                                   gCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGaa                                                     <  3:443668/148‑1 (MQ=16)
                                                                                                     ccATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTACaaaaaaaatcaaaagctc                                                >  8:1123038/1‑137 (MQ=11)
                                                                                                     ccATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATtctc                                                >  2:1029037/1‑151 (MQ=255)
                                                                                                       aTGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACttacttaaataat                                                          >  6:1182179/1‑129 (MQ=255)
                                                                                                       aTGCGGGGAAGGATGAGAAGCTGCGACTAAGGTTCGGCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCaa                                              <  5:436289/151‑1 (MQ=255)
                                                                                                        tGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCttccttaaataatc                                                         >  5:1064841/1‑128 (MQ=11)
                                                                                                           ggggAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTTCTAAAGAACGGGAAATCTGATACgc                                          >  3:212578/1‑150 (MQ=11)
                                                                                                               aaGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATcc                                      <  6:698592/151‑1 (MQ=255)
                                                                                                                aGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTCCTAAAAACCACGGaagagcctcgggtacgg                                      >  6:928541/1‑136 (MQ=11)
                                                                                                                  gATGAGAAGCTTCGACCAAGGTTCGCCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCa                                   <  2:480168/151‑1 (MQ=16)
                                                                                                                   aTGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCaa                                  <  4:265989/151‑1 (MQ=255)
                                                                                                                    tGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGAAATGATTACCCAACAAGACTGAAATACAAAAGACCCGGAATTATCAAACGAAGCGGCAAt                                 >  6:710040/1‑151 (MQ=14)
                                                                                                                     gagaAGTTTCGACCAAGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATc                                <  7:1359192/151‑1 (MQ=255)
                                                                                                                       gaAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCa                               >  1:471010/1‑150 (MQ=255)
                                                                                                                         accTTCGACTATGGTTCGGCTCGTGGTTCCGGGATTGAGCGTTGCCGCAGCAAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCCATCAAat                            <  7:191351/149‑1 (MQ=14)
                                                                                                                         aGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCaaa                             >  4:1257814/1‑150 (MQ=255)
                                                                                                                         aGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAat                            >  2:274490/1‑151 (MQ=255)
                                                                                                                            ttCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGCCCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATAtt                         >  3:288578/1‑151 (MQ=255)
                                                                                                                             tCGACTAAGTTTCGACTCGAGCGTCAGCTATAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAatat                          <  8:1453485/149‑1 (MQ=255)
                                                                                                                             tCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAatat                          <  1:1029870/149‑1 (MQ=255)
                                                                                                                             tCGACAAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAatat                          <  1:1415189/149‑1 (MQ=255)
                                                                                                                             tCGACAAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATAtt                         <  1:1625340/150‑1 (MQ=255)
                                                                                                                              cGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCt                       >  3:833916/1‑151 (MQ=255)
                                                                                                                              accctcttgggatgttCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCt                       <  4:1146704/136‑1 (MQ=11)
                                                                                                                               gcccaaGGTGCGCCGCGCGCGCCAGCGAGAGAGCGTGGCCCCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTg                      <  2:1377454/147‑1 (MQ=11)
                                                                                                                                 caaaGGTTCGTCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGAAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGcc                    <  4:1098758/149‑1 (MQ=255)
                                                                                                                                   aaGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGAGCTGCGCTGAGTAAGCGTTACCCCAACACGCTCAGGTGCAGAAGACCCGGAAGTGGGACCCGGATCCGGAATCAAATATTCTGccgg                  >  4:887219/1‑149 (MQ=2)
                                                                                                                                    aGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCaa                 <  7:686316/151‑1 (MQ=255)
                                                                                                                                     ggTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCCTCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCa                  >  8:7161/1‑149 (MQ=255)
                                                                                                                                        tCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCg              <  3:1259082/150‑1 (MQ=255)
                                                                                                                                         ggtcTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCggg            <  6:1011613/148‑1 (MQ=255)
                                                                                                                                               gAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGa      <  5:1315266/151‑1 (MQ=38)
                                                                                                                                                    |                                                                                                                                                      
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: