Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,591,5130TG24.2% 286.9 / 48.6 297D86A (GAC→GCC) ECDH1_RS12815hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (89/96);  new base G (23/36);  total (130/167)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 298 total aligned reads displayed.

CCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATG  >  DH1/2591365‑2591660
                                                                                                                                                     |                                                                                                                                                    
catgAAGGGTGATGCGCGTAGCTGTTACTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCa                                                                                                                                                    <  8:344082/148‑1 (MQ=17)
          gAAGCGCGCAGCGCTGAGTAATCCTTCGCCCACGACGATCAGTTTCTAAAGTCCCGGAAGTCTCAAGCGCATCCGCAATCAAATATTATGCCCATGCGGGGCAGGATGCGCAGCTTCGCCTAAGGTTCGGCTGGCGCGGCAGCgcgcgcgc                                                                                                                                          >  1:100300/1‑144 (MQ=14)
          caacgtctgagaggtGAGGAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGc                                                                                                                                          <  3:531257/137‑1 (MQ=14)
             cgttgtcgatgtttgTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTTTCAATCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGtt                                                                                                                                       <  8:110067/137‑1 (MQ=14)
                gttagaaggttgTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTATCAAGCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGcc                                                                                                                                    <  6:835663/140‑1 (MQ=14)
                cgaaGCGCGGACTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGcc                                                                                                                                    <  1:352755/148‑1 (MQ=255)
                 gcAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTGGAGCGGCAGCGAGAGAGCGTTGCCg                                                                                                                                   >  2:798858/1‑151 (MQ=255)
                  gagacgggtcgaaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCGCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGcc                                                                                                                                    <  5:166311/138‑1 (MQ=14)
                    gcacgtacgaaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAg                                                                                                                                <  1:1714867/142‑1 (MQ=255)
                       ggttctAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAACGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCa                                                                                                                             <  8:2193486/146‑1 (MQ=255)
                       ccgcgcaacccatCCCCCACCACCAACACTTGCACCCCACCCGCCACGCTCAAGAGACACAGCACTCAACTATTCTGCCCATCCGCTGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCa                                                                                                                             <  3:1996703/139‑1 (MQ=9)
                        tGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGACCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGACTTGAGCGCCAGCGAGAGAGctaggaagagcgtcg                                                                                                                            >  2:35328/1‑137 (MQ=11)
                         taGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTTTCAAGCGAATCCGCAATGAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCaa                                                                                                                            <  8:1582173/149‑1 (MQ=255)
                         gacgaaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCCAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAAc                                                                                                                           <  5:1337394/147‑1 (MQ=255)
                           gTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGCCTGGAGCGGCAGCga                                                                                                                                                >  5:2049814/1‑128 (MQ=255)
                           cgaaTCCATCCCCCACCACCATCACTTTCAAAACTCCCTGAACTCACAACCCAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGa                                                                                                                         <  5:75316/149‑1 (MQ=14)
                            tAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGAc                                                                                                                        >  1:1915446/1‑151 (MQ=255)
                             aaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGa                                                                                                                         <  3:1419979/149‑1 (MQ=255)
                             aaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCCCAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGAcc                                                                                                                       <  3:1805165/151‑1 (MQ=255)
                             aaTCCTTCCCCCACCACCATCACTATCAAAAGTCCCTGAACTCACAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGGTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGa                                                                                                                         <  5:64519/149‑1 (MQ=255)
                                 cTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTAGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGa                                                                                                                    >  8:1854069/1‑150 (MQ=255)
                                   tCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAACCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGaa                                                                                                                   <  6:1428491/149‑1 (MQ=255)
                                    cccccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAggg                                                                                                                <  5:1459699/151‑1 (MQ=255)
                                    cccccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTAGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAAAGAACCGAAgcg                                                                                                                >  1:909074/1‑149 (MQ=14)
                                    ccaccaccaccaGCACTTTCAAAAGGCCCTTAAGTCACAGGCGATGCGCCATGCCATTTTGCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAgg                                                                                                                 <  7:1840305/147‑1 (MQ=14)
                                       ccaccaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCg                                                                                                                                                 <  5:1108817/115‑1 (MQ=255)
                                          ccaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGGCTAGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgc                                                                                                          >  6:2193167/1‑151 (MQ=255)
                                             ccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGcacatcgc                                                                                                        >  1:1900202/1‑143 (MQ=255)
                                              caTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                      >  7:1346928/1‑151 (MQ=255)
                                               ggcccTTTCCAAAGCCCCTCAACTCCCAACCCAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGACGCTTCGACTAAGGTGCGGCGCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCa                                                                                                     <  7:647707/147‑1 (MQ=11)
                                               aTCCCTTTCAAAAGTCCCTCAACTCTCACCCGAACCCCCAATCAAATATTCTGCCCATGCGGGCAAGGATGAGAAGCTGCGACTAAGGTGCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCa                                                                                                     <  1:742498/151‑1 (MQ=14)
                                               aTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                      >  3:1298377/1‑150 (MQ=255)
                                                 cACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgc                                                                                                   >  2:54723/1‑151 (MQ=255)
                                                 cACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgc                                                                                                   >  8:439319/1‑151 (MQ=255)
                                                 cACTTTCAAAAGTCCCTGAACTCCCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgc                                                                                                   <  5:1303452/151‑1 (MQ=255)
                                                  aCTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgcg                                                                                                  >  1:991540/1‑151 (MQ=255)
                                                  aCTTTCAAAAGTCACTGAACTCTCAAGCCAATCCGCAATCAAATATTCTCCCCAAGCGGGGTAGGATGAGAAGCTTCGACTAAGCTTCGGCTCGAGCGACAGCGAGAGATCGTTGCCGCAGGCAACGACCCCAACGGCGAAGCCCGCAgcg                                                                                                  >  6:1110757/1‑151 (MQ=255)
                                                     ttCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCACAtcgga                                                                                                >  8:242408/1‑148 (MQ=255)
                                                     ttCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGACTTGAGCGCCAGc                                                                                                                                                  <  4:1829493/100‑1 (MQ=11)
                                                      tCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTg                                                                                               >  6:173689/1‑150 (MQ=255)
                                                      tCAAAAGTCCCTGAACTCGCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTg                                                                                               <  3:969237/150‑1 (MQ=255)
                                                      tCAAAAGTCCCTGAACTCACAACCGAATCCGAAATAAAATATTCTGCCCATGGGGGTAAGGATGAGAAGCTTCGACTACGGTTCGGCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGAAACGACCCGAAGGGCGAAGCGCGCAGCGCTg                                                                                               <  2:1901218/150‑1 (MQ=11)
                                                        aaaaGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTGGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCAGAAGGGCGAAGCGCGCAGCGCTGAGt                                                                                            >  5:186830/1‑151 (MQ=255)
                                                         aaaGTCCCTGAACTTCCAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTa                                                                                           <  6:291259/151‑1 (MQ=255)
                                                         aaaGTCCCTGAACTCCCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTGCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTa                                                                                           <  3:1021806/151‑1 (MQ=255)
                                                           aGTCTCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAAt                                                                                         >  3:829570/1‑151 (MQ=255)
                                                             tCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATcc                                                                                       >  2:170271/1‑151 (MQ=255)
                                                               ccTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCtt                                                                                     >  3:1504799/1‑151 (MQ=255)
                                                                cTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTc                                                                                    >  4:324301/1‑151 (MQ=255)
                                                                  gAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTccc                                                                                  >  1:1175768/1‑151 (MQ=255)
                                                                  caaCTCTCAACCGAATCCGCAATCAAATATTCTCCCCATGCGGGGAAGTATTAGAAGCTTCGCCTAAGGTTCGGCTCGACCGTCAGCGAGAGAGCGTTGCCGCAGGAAACGACCCGAAGGTCGAAGCGCGCAGCGCTGAGTAATCCTTccc                                                                                  <  6:1384327/150‑1 (MQ=14)
                                                                   aaCTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTcccc                                                                                 >  1:131002/1‑151 (MQ=255)
                                                                       cccAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTGCGACTAAGGTTCGGCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCcca                                                                               <  1:634999/147‑1 (MQ=255)
                                                                         cAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccac                                                                           >  7:1614456/1‑151 (MQ=40)
                                                                            gCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTTGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCAt                                                                        <  4:827563/151‑1 (MQ=2)
                                                                             ggAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATc                                                                       <  2:1618131/150‑1 (MQ=14)
                                                                              gAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATc                                                                       >  7:973136/1‑150 (MQ=38)
                                                                                aTCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGACTAGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGc                                                                                                               >  6:1938560/1‑108 (MQ=255)
                                                                                 tCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCtt                                                                   >  6:1603360/1‑151 (MQ=38)
                                                                                  ccGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACtt                                                                   <  2:2102629/150‑1 (MQ=14)
                                                                                  ccGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACttt                                                                  <  4:1624202/151‑1 (MQ=18)
                                                                                  ccGAAATTAAATATTCTGCCCATGCGGGTAAGTATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACttt                                                                  <  8:1634356/151‑1 (MQ=255)
                                                                                   cGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTg                                                                                               <  4:1254487/121‑1 (MQ=21)
                                                                                   cGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCACAGGCAACGACCCg                                                                                                                     <  2:1555219/99‑1 (MQ=255)
                                                                                    gCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttcc                                                                 >  7:1890163/1‑148 (MQ=37)
                                                                                     cAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCttcctt                                                               >  3:1903535/1‑147 (MQ=255)
                                                                                       aTCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTGCata                                                              >  7:794/1‑148 (MQ=255)
                                                                                        cgatattatTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCaaaa                                                             <  6:2210839/145‑1 (MQ=16)
                                                                                           aaaTTTTCTGCCCTTGGGGGTATGTATGATATGCTTCTACTATGGTTCGCCTCGTGCGTCTGCTATTTTGGGTTGCCGCTGGAAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGcc                                                          <  8:510812/151‑1 (MQ=2)
                                                                                            aaTATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGccc                                                         >  2:685162/1‑151 (MQ=38)
                                                                                              tatTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTACaaaaataatc                                                       >  6:2058261/1‑145 (MQ=255)
                                                                                                ttctgagcttGCGGGGAAGTATTAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCgg                                                      <  8:2155612/141‑1 (MQ=11)
                                                                                                 tCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGa                                                     <  1:2001009/150‑1 (MQ=16)
                                                                                                  cTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTACataaaaaaagatat                                                   >  8:796844/1‑141 (MQ=255)
                                                                                                    gCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGaa                                                    <  4:401744/148‑1 (MQ=38)
                                                                                                     cccATGCGGGGATGGATGAGACCCTGCGCCCAAGGTGCGCGGCGACCGCCTGCTATAGCGGGTGGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATtct                                                <  7:1600195/151‑1 (MQ=11)
                                                                                                      ccgctcttcgtatcttGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATtctc                                               <  2:111006/136‑1 (MQ=11)
                                                                                                       cATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCa                                                                                                                                                    <  4:1428753/49‑1 (MQ=255)
                                                                                                       cATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCa                                              <  1:2171838/151‑1 (MQ=255)
                                                                                                         tGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCa                                              <  7:1533844/149‑1 (MQ=38)
                                                                                                          ggggggttggttGATAAGCTTCGACCAAGGTTCGCCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAAc                                           <  3:1569402/140‑1 (MQ=11)
                                                                                                          gCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCaaa                                            <  4:1355130/150‑1 (MQ=255)
                                                                                                          gCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAAc                                           <  3:636955/151‑1 (MQ=255)
                                                                                                            ggggAAGGATTATAAGCTTCGACAAAGGTTCGGCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGa                                         <  1:2066096/151‑1 (MQ=255)
                                                                                                             gggAAGGATGAGAAGCTTCGACCAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGaa                                        <  4:2086921/151‑1 (MQ=255)
                                                                                                             gggAAGGAGGAGAAGCGTCGACTACGGTGCGGCTCGAGCGTCTGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAATGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGaa                                        <  7:356133/151‑1 (MQ=14)
                                                                                                                aaGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATcc                                     <  6:153891/151‑1 (MQ=255)
                                                                                                                  ggATGATAAGGTTCTACTAAGGTTCGTCTCTAGCGTCATCTTGTTTTGTTTGCCGCAGTCAACGACCAGAAGGTCTAAGCTCGCAGCTCTGAGTAATCATTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCg                                    <  6:572898/150‑1 (MQ=9)
                                                                                                                         acgCTTCGAGAAGGGTCCGCCGCGAGCGCCAGCGAGAGAGCGTTGCCGCAGCCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCaaa                            <  7:1017599/149‑1 (MQ=11)
                                                                                                                           gCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCACAGGAATTCTCAAACGAATCCGCAATCaa                             >  4:750878/1‑148 (MQ=255)
                                                                                                                           gCTGCGACGACGGTTCGGCGCGCGCGTCCGCGAGAGAGCGTGGCCGCTGGCAACGTCCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACGCTCACTTTCAAACGCCCCGTAATTCTCAAACGAATCCGCAATCAAata                          <  7:217122/151‑1 (MQ=14)
                                                                                                                                gACGAAGGTTCGACTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTg                    <  1:163279/151‑1 (MQ=255)
                                                                                                                                 acaaaGGTTCGCCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGc                   <  8:1786295/148‑1 (MQ=255)
                                                                                                                                     aGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTACTTAAGCACCCGAATTCTGAAACGCATCCGCAATAAATACTTATGCCCa                >  5:1236309/1‑150 (MQ=14)
                                                                                                                                       gTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCACCGGAATTCTCAAACGACTACGGAATCAATAAATCTGGCCCAGc             >  5:1045776/1‑151 (MQ=14)
                                                                                                                                          cGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCtt                                                                   >  2:1187420/1‑95 (MQ=25)
                                                                                                                                          cGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCggg          >  3:482738/1‑151 (MQ=255)
                                                                                                                                           ggCGCGCGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCgggg         <  8:1474882/151‑1 (MQ=255)
                                                                                                                                                    gTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTTCTAAAGGCCCGCAAAACTCTTACGACTCCGCAATCAACGAGTATGCCCCAGCGGGGaacgcgg  >  3:221272/1‑144 (MQ=11)
                                                                                                                                                    gTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCACCACCCCCATCACGTTCAAACGCCGCGGATTTCTCAACCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATg  >  5:2296099/1‑149 (MQ=14)
                                                                                                                                                     |                                                                                                                                                    
CCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATG  >  DH1/2591365‑2591660

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: