Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,591,5130TG28.5% 216.8 / 54.4 261D86A (GAC→GCC) ECDH1_RS12815hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (96/62);  new base G (20/43);  total (128/133)
Fisher's exact test for biased strand distribution p-value = 1.04e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 261 total aligned reads displayed.

CGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATG  >  DH1/2591367‑2591660
                                                                                                                                                  |                                                                                                                                                   
cGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGCCTTGAGCGCCAGc                                                                                                                                                 >  5:231852/1‑151 (MQ=255)
        gAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGc                                                                                                                                         >  7:949927/1‑151 (MQ=255)
               gcAGCGCTGAGTAATCCATCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTGCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCg                                                                                                                                  <  3:840383/151‑1 (MQ=255)
                tagatgtGAGGAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGc                                                                                                                                    <  4:925431/142‑1 (MQ=255)
                 aGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCa                                                                                                                                <  3:275620/151‑1 (MQ=255)
                     ggttgTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAATCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCa                                                                                                                            <  8:1172425/147‑1 (MQ=255)
                     ggttcgaaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGc                                                                                                                                                 <  6:1583077/124‑1 (MQ=255)
                      tGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAGAAGCTTCGACAAAGATACGAATCGAGACTAAGAGAGATAGAGtttcatcaggaa                                                                                                                            >  8:1404613/1‑140 (MQ=11)
                      gttctAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTATCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCaa                                                                                                                           <  4:1394100/147‑1 (MQ=255)
                       gacgaaTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAAc                                                                                                                          <  1:1766673/147‑1 (MQ=255)
                            aTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGAccc                                                                                                                     <  8:948920/151‑1 (MQ=255)
                            aTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGAccc                                                                                                                     >  5:1809736/1‑151 (MQ=255)
                                ttCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGCCAAAGGTTCGGCTGGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGCCCCGa                                                                                                                   >  6:852349/1‑149 (MQ=255)
                                  cccccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAggg                                                                                                               >  3:430309/1‑151 (MQ=255)
                                   ccccACCACCATCACTTTCAAAAGTCCCTGAACTGACAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGc                                                                                                              <  5:943284/151‑1 (MQ=255)
                                     ccaccaccaTCACTTTCAAAAGTCCCTGAACTCTCAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGa                                                                                                            <  4:486596/151‑1 (MQ=255)
                                      caccaccaTCATTTTAAAAAGTCCCTTAACTCTCAATCTAATCCGCAATAAATTATTCTGCCCATGGGGGGATGTATTATAAGCTTCGACTATGGTTCGGCTCGAGCGTCAGCTATTGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGaa                                                                                                           <  8:190620/151‑1 (MQ=14)
                                      caccaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGa                                                                                                            >  6:1736282/1‑150 (MQ=255)
                                        ccaccaTCACTTTCATATGTCCCTTAAATCTCAAGCGAATCCGCAATCTAATATTCTGCCCATGCGGGGAAGGTTTAGAAGCTTCGACTAAGGTTGGCCTCGAGGGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgc                                                                                                         <  4:1248771/151‑1 (MQ=14)
                                         caccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcg                                                                                                        >  6:774006/1‑151 (MQ=255)
                                          accaTCACTTTCAAAAGTCCCTGAACTCTCAACCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgc                                                                                                       <  2:1058017/151‑1 (MQ=255)
                                          accaTCACTTTCAAAAGTCCCTGAACTCTCAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgc                                                                                                       <  6:213047/151‑1 (MQ=255)
                                           ccaTCACTTTCAAACGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcg                                                                                                      <  4:217043/151‑1 (MQ=255)
                                           ccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcg                                                                                                      >  6:338323/1‑151 (MQ=255)
                                           ccaTCACTATCAAAAGGCCCTCACCTCTCACGCCATGCCCCACGCAAATATCCGGCCCCTCCGGGGAAGTATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcg                                                                                                      <  1:1208025/151‑1 (MQ=14)
                                             aTCACTTTCAAAAGTCCTTTAAATCTCAAGCGATTCCGCAATAAAATATTCTGGCCTTGCGGTTAATTATGATAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCa                                                                                                    <  8:1502502/151‑1 (MQ=14)
                                             aTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                     >  7:1531425/1‑150 (MQ=255)
                                             aTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                     >  1:654788/1‑150 (MQ=255)
                                               cactaccccaacgcccTGAACTCACACGGCAAGCCGCACACAAGTCTGCTGCCCCTCCCGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgc                                                                                                  <  1:1111135/138‑1 (MQ=11)
                                                ccTTTCAAAAGTCCCTCACCTCTCAAGCGAATCCGCAATCAAATATTCTGCCCCTGCGGGGAAGGATGAGAAGCCGCGACTAAGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgcg                                                                                                 <  5:1223352/150‑1 (MQ=14)
                                                agggcataaaGTCCCTGAACTCGCCATCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGAGGAGAAGCTTCGACAAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgcg                                                                                                 <  1:379900/144‑1 (MQ=11)
                                                  tttCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCt                                                                                               >  5:1729487/1‑151 (MQ=255)
                                                  tttCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCt                                                                                               >  2:1260874/1‑151 (MQ=255)
                                                   ttCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTg                                                                                              >  6:292500/1‑151 (MQ=255)
                                                    tCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgccgcgcgctcg                                                                                             >  3:951459/1‑141 (MQ=255)
                                                      aaaaGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGt                                                                                           >  6:703388/1‑151 (MQ=255)
                                                      aaaaGTCCCGTATCTCTCAACGGAGTCCGCAACCAATCTTTCGGCCCATGCGGTCATGGATGACAAGGTGCGACAACGGTTGGCCGCGAGCGGCAGCGAGATAGCGTGGCCGCAGCCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGt                                                                                           <  7:179989/151‑1 (MQ=2)
                                                       aaaGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGt                                                                                           >  8:279575/1‑150 (MQ=255)
                                                       aaaGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGt                                                                                           >  6:548398/1‑150 (MQ=255)
                                                        accccacTGAACTCTAAAGCGACTCCGAAATCAATTATTCTCCCCAGGCGGTTAAGTTTTAGAAGCTTCGACTAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTaa                                                                                         <  6:1379218/145‑1 (MQ=14)
                                                        aaGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCGATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAg                                                                                            >  3:128836/1‑148 (MQ=255)
                                                         aGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCAGCGTAAt                                                                                        >  7:147436/1‑151 (MQ=255)
                                                         aGTCCCTGAACTCCCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAAt                                                                                        <  3:1121687/151‑1 (MQ=255)
                                                             ccTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCtt                                                                                    >  5:1892312/1‑151 (MQ=255)
                                                              cTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGAGATccga                                                                                    >  1:409441/1‑148 (MQ=255)
                                                              cTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGCCTAGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGa                                                                                             >  7:128814/1‑141 (MQ=255)
                                                                gAACTGGCAACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTc                                                                                   <  5:1136327/149‑1 (MQ=255)
                                                                gAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTccc                                                                                 >  5:249114/1‑151 (MQ=255)
                                                                gAACTATCAATCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTc                                                                                   <  4:361575/149‑1 (MQ=255)
                                                                  aCTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGGCAGAGAGAGCGAGATCACGCCTCGGACCACCAGTAGGCAGAAGCGCGCAGCGCTGCGTAATCCTTccggc                                                                               >  7:730567/1‑148 (MQ=9)
                                                                   ctctCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAg                                                                                            >  7:1360101/1‑137 (MQ=255)
                                                                    tctcAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGCGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGCCCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTccccccc                                                                             >  1:751107/1‑149 (MQ=255)
                                                                    tctcAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTccccacc                                                                             >  2:230429/1‑148 (MQ=255)
                                                                     ctcAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccacc                                                                            >  3:690987/1‑151 (MQ=255)
                                                                     ctcAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccacc                                                                            >  3:1667078/1‑151 (MQ=255)
                                                                       cAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTTGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccacc                                                                            <  4:690303/149‑1 (MQ=9)
                                                                       cAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccac                                                                          >  7:1186951/1‑151 (MQ=40)
                                                                       cAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccac                                                                          >  2:461113/1‑151 (MQ=40)
                                                                        aaGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccacc                                                                         >  1:1407377/1‑151 (MQ=40)
                                                                         attccttactcAATCAAATATTCTGCCCATGCGGGTAAGTATGAGAAGCTTCGACAAAGGTTCGCCTTGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccacca                                                                        <  2:1575651/141‑1 (MQ=11)
                                                                         aGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccacc                                                                         >  1:1913837/1‑150 (MQ=40)
                                                                             aaTCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTGGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCtt                                                                                    >  8:1799125/1‑135 (MQ=255)
                                                                             aaTCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCAc                                                                    >  2:1729406/1‑151 (MQ=38)
                                                                              atcctaaaTTAAATATTCTGCGCATGGGTGGAAGGATGTTATGCTTCGACTATGGTTGGTCTCGAGCGTCAGCTAGAGAGGGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCACCCACCACCATCACt                                                                   <  6:1684113/145‑1 (MQ=11)
                                                                               tCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCtt                                                                  >  4:1717102/1‑151 (MQ=25)
                                                                                   cAATCAAATATTCTGCCAATGCGGGGAAGGATGAGAAGCATCGACAAAGTTTCGGCTCGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGAGCGCAGCGCTGAGTAATCCTTAAACCACAACCATCACttact                                                               >  6:645529/1‑147 (MQ=11)
                                                                                     aTCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCAt                                                                       >  4:1777333/1‑140 (MQ=39)
                                                                                     aTCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttgctaa                                                             >  7:339772/1‑145 (MQ=255)
                                                                                            atTCTGCCCATGGGGGGAAGGTTGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGGCAGCTAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCg                                                      <  6:986464/150‑1 (MQ=255)
                                                                                            atTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTTCtaaactaaggc                                                     >  3:564461/1‑144 (MQ=255)
                                                                                             ttCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGTCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccacc                                                                         <  2:1377555/130‑1 (MQ=17)
                                                                                              tcggCCCCTGCGGGGAAGGCGGCGAAGAGACCACCAAGGGGCGGCGCGTGGTTCTGCTAGATAGCGTGGCCCCAGGAAACGTCCCTAAGGGAGAAGCGGTCAGCGGTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGa                                                    <  3:485403/147‑1 (MQ=2)
                                                                                               cTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCttgcttaaatactcaaa                                                   >  5:1381525/1‑135 (MQ=11)
                                                                                                 gCCCATGCGGGGAAGTATGAGAAGCTTCGACAAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAAt                                                  <  4:1565072/149‑1 (MQ=255)
                                                                                                 gCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATtc                                                <  3:1793235/151‑1 (MQ=255)
                                                                                                 gCCCATGCGGGGAAGGATGAGAAGCTTCGACCACGGTTCGCCTCGAGCGCCAGCTAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAAt                                                  <  4:603029/149‑1 (MQ=255)
                                                                                                     aTGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCaa                                            <  2:944938/151‑1 (MQ=255)
                                                                                                      tGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTGGAGCGTCAGCGAGAGAGCCTTGACCCAGCCAGACACCCGAAGCGCGAAGATCGCAGCGCTGAGTAATCGTTCGCCCACAAAAAAAACTTTCACAAGCCCCGGAATtctccgc                                           >  5:543555/1‑148 (MQ=9)
                                                                                                       gCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAAc                                          <  6:580709/151‑1 (MQ=255)
                                                                                                       gCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgc                                                                                                         >  4:30040/1‑88 (MQ=255)
                                                                                                           ggAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTGTCAAAAGCCCCGGAATTCTCAAACGAAt                                      >  1:80579/1‑151 (MQ=255)
                                                                                                             aaGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATcc                                    <  7:619379/151‑1 (MQ=255)
                                                                                                               ggATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccac                                                                             >  2:700073/1‑108 (MQ=39)
                                                                                                                       cgcggcGACAAAGGGGCGCCGCGAGCGCCAGCGAGAGAGCGGTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAat                          <  5:657681/146‑1 (MQ=11)
                                                                                                                          ttCGACTAAGGTTCGGCTCGAGCGTCAACGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCCGCCCCATCACTTTCAAAAGCCCCGGCATTGTCAAACGCATCCGCAATCCAATAtt                       >  4:841323/1‑151 (MQ=14)
                                                                                                                           ccggctccggTGCGGCTCGCGCGCCCGCGAGAGCGCGTGGCCCCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTc                      <  7:190430/143‑1 (MQ=14)
                                                                                                                            cGACTAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCt                     <  8:1772105/151‑1 (MQ=255)
                                                                                                                             cacaccggggcgccTCGAGCGGCAGCGAGAGAGCGTTGCCGCCGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAACCGAATCCGCAATCAAATATTCTg                    <  5:328935/138‑1 (MQ=11)
                                                                                                                                tAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGCATCCGCAATCAAATATTCTGccc                 >  1:1648373/1‑151 (MQ=255)
                                                                                                                                caaGGTGCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGccc                 <  3:1708559/150‑1 (MQ=16)
                                                                                                                                   ggctccgcgcGAGCGCCAGCGAGAGAGCGTTGCCGCAGTCAACGACTCTAATGTCGAAGCGCGCAGCGCTCAGTAATCCTTCCACCACCATCATTACTTTCAACAGGCGCTGACTTCTTAAACGAATCCGCAATCAAATATTCTGCCCAAg              <  2:321988/142‑1 (MQ=11)
                                                                                                                                   ggTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCa                >  8:657486/1‑149 (MQ=255)
                                                                                                                                   ggTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAg              <  6:1780686/151‑1 (MQ=255)
                                                                                                                                    gTTCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGACCCGGAATTCTCAAACGAATCCGCAATCAAATAtt                       >  8:1430205/1‑141 (MQ=9)
                                                                                                                                      tCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCg            <  5:1741240/150‑1 (MQ=25)
                                                                                                                                           tCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGa        >  1:1326389/1‑149 (MQ=255)
                                                                                                                                            cGAGCGCCAGCGAGATAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCTGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAgg     <  5:1501043/151‑1 (MQ=16)
                                                                                                                                            cGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCTGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAgg     <  7:1612604/151‑1 (MQ=18)
                                                                                                                                            cGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGACCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGaa       <  5:1012789/149‑1 (MQ=25)
                                                                                                                                                 gTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCACAGCGCTGAGTAATCCTTCCCGCACAACCATCACTTTCAAAAGCACCGGAATTATCAAACGATTCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATg  >  2:1127055/1‑149 (MQ=14)
                                                                                                                                                  |                                                                                                                                                   
CGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATG  >  DH1/2591367‑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: ‑