Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,591,5130TG25.4% 246.5 / 49.1 279D86A (GAC→GCC) ECDH1_RS12815hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (106/64);  new base G (21/37);  total (149/130)
Fisher's exact test for biased strand distribution p-value = 7.08e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 280 total aligned reads displayed.

CCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATG  >  DH1/2591365‑2591660
                                                                                                                                                    |                                                                                                                                                   
  cGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGCCTTGAGCGCCAGc                                                                                                                                                 >  1:52154/1‑151 (MQ=255)
   gAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTGCAAAAGTAACGCAACGCTCAACCGAATCCCGAATAACATATTCAGCCCAACCGGGGGACCTGGAGAAGCATGCAGCAAGGTGCGGGTCGCGCGTCAGCg                                                                                                                                                >  7:113991/1‑151 (MQ=2)
   gAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAAATCTCAAGCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGACTAGAGCGACAGCg                                                                                                                                                >  2:425488/1‑151 (MQ=11)
    aaGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGACTGGAGCGGCAGCg                                                                                                                                                >  4:1421202/1‑150 (MQ=255)
     tgggcgatgccctgcgcgcGGACTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCga                                                                                                                                               <  1:1687843/135‑1 (MQ=17)
        gCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGGCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTGGAGCGCCacgctcgg                                                                                                                                            >  1:1175350/1‑143 (MQ=14)
           atcgccttcgaccttaGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCg                                                                                                                                        <  5:1035718/136‑1 (MQ=255)
                 gcAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAAATCTCAAGCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAGAAGCTTCAACTAAGGATCGTCTAGATGGTCAGCAAGAGAAAGTTGACg                                                                                                                                  >  8:762600/1‑151 (MQ=14)
                  cAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGc                                                                                                                                 >  3:631708/1‑151 (MQ=255)
                  cAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGACTTGAGCGACAGCGAGAGAGCGTTGCCGc                                                                                                                                 >  2:1122822/1‑151 (MQ=255)
                   cgaacggtctAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCa                                                                                                                                <  3:1442274/142‑1 (MQ=255)
                   aGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGCCTAGAGCGCCAGCgagatcgccagcgcgc                                                                                                                                 >  7:1825920/1‑138 (MQ=14)
                      gcTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTGGAGCGGCAGCGAGAGAGCGTTGCCGCAGGc                                                                                                                             >  7:611868/1‑151 (MQ=255)
                        tGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTAGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCa                                                                                                                            >  4:698127/1‑150 (MQ=255)
                         gAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGACTGGAGCGCCAGCGAGAGAGCGTTGCCGCAggaga                                                                                                                           >  2:547673/1‑147 (MQ=11)
                         gAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAAATCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGACTAGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAAc                                                                                                                          >  6:1229367/1‑151 (MQ=255)
                            tAATCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGACTAGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGAc                                                                                                                       >  6:1461777/1‑151 (MQ=255)
                            tAATACTTCCCCCAACACCATCAATTTCAAAAGTCCCTTAAATTTCAAACGAATCCGCAATTAAATATTCTTCCCATGGGGGGTAGGATGATAAGCTTCGACTAAGGTTCGGCTCGTTCGTGATCGATATAGCGGTGGCGCAGGCAACGa                                                                                                                        <  8:102060/150‑1 (MQ=2)
                              aTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGtt                                                                                                                                      <  3:1369912/134‑1 (MQ=255)
                              aTCCTTCCCCCACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTGGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGAccc                                                                                                                     >  1:1821253/1‑151 (MQ=255)
                                   tCCCCCACCACCATCAATTACAAAATTCCGTGAACTCTCAAGCTATTCCTCATTTAATTTTTCTGCCCTTGGGGGTAATTTTTAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCTAGATTGCGTTGCCGCTGGCAACGACCCGAAg                                                                                                                 <  6:1628919/150‑1 (MQ=255)
                                    cccccACCACCATCACTTTCAAAAGTCCCTGAACTTTCAAGCGAATCCGCAATTAAATATTCTGCCCATGCGGGTAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAggg                                                                                                               <  6:1115869/151‑1 (MQ=255)
                                    cccccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGACTAGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAggg                                                                                                               >  6:1574441/1‑151 (MQ=255)
                                     ccccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGc                                                                                                              <  1:1145224/151‑1 (MQ=255)
                                     ccccACCACCATCACTTTCAAAAGTCCCTGAACTCTCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGc                                                                                                              <  3:217812/151‑1 (MQ=255)
                                     ccccACCACCATCACTTTCAAAAGTCCCTGAACTCCCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTGCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGc                                                                                                              <  7:1765212/151‑1 (MQ=255)
                                      cccACCACCATCACTTTCAAAAGTCACTGAACTCTCAAGCAAATCCGCAATCAAATATTCTGCCAATGCGGGGAAGGATGAGAAGCTTCGACTAATGTTCGCCTAGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCg                                                                                                             >  2:1476834/1‑151 (MQ=14)
                                       ccaccaccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCg                                                                                                             >  3:999858/1‑150 (MQ=255)
                                            accaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCGGAAGGGCGAAgcgc                                                                                                       >  3:938189/1‑151 (MQ=255)
                                             ccaTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcg                                                                                                      >  3:1690486/1‑151 (MQ=255)
                                              caTCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                     >  5:1522062/1‑151 (MQ=255)
                                               caccttttCAAAAGTCCCTGAACTCCCCACCGAATCCGCAATAAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCa                                                                                                    <  1:548490/146‑1 (MQ=11)
                                                tCACTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAgcgcgc                                                                                                     >  3:33954/1‑149 (MQ=255)
                                                 cacatgaaaacGTCCATGCAGTCTCCAGCGAAGCCGCAAGCAAATACTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgc                                                                                                  <  1:1770100/144‑1 (MQ=14)
                                                  aCTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgcg                                                                                                 >  2:852139/1‑151 (MQ=255)
                                                   cTTTCAAAAGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAgcgc                                                                                                >  3:1799456/1‑151 (MQ=255)
                                                       aaaatGTCCGTGAACTATCAATCTAATCGTAAATAAAATATTCTGCCCTTGGGGGTATGTATGATAAGGTTCGACTAAGGTTGGTCTCTAGCGTCTGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAg                                                                                            <  2:902212/150‑1 (MQ=9)
                                                           ccccccTCAACGCACAACCCAAGCCCCAATCAAATATTCGGCCCACGCGGGGAAGGATGAGAAGCTTCGACCACGGTGCGACTCGAGCCCCAGCGAGAGAGCGTGCCCGCAGCCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAAt                                                                                        <  7:260334/148‑1 (MQ=1)
                                                           aGTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAAt                                                                                        >  8:1865867/1‑151 (MQ=255)
                                                            gTCCCTGAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATc                                                                                       >  3:548438/1‑151 (MQ=255)
                                                             tCCCTGAACTGTCAAGCTAATCCGAAATAAAATATTCTGCCCATGGGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCTAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATc                                                                                       <  2:887407/150‑1 (MQ=14)
                                                             gcccTGAACTCTCAACCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATcc                                                                                      <  5:1133176/150‑1 (MQ=255)
                                                              cccTGACCCCTCAAGCGACTCCCCAATCAAATATTCTGCCCACGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCt                                                                                     <  6:1841266/151‑1 (MQ=255)
                                                               ccTGAACTTTCAATCGAATCCGCAATTAAATATTCTGCCCATGCGGGTAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGGGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCtt                                                                                    <  4:1235526/151‑1 (MQ=255)
                                                                  gAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTccc                                                                                 >  5:163777/1‑151 (MQ=255)
                                                                  gAACTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTccc                                                                                 >  5:755397/1‑151 (MQ=255)
                                                                   aaCTCTCAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTcccc                                                                                >  2:159075/1‑151 (MQ=255)
                                                                      tctcAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGGCTGGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccac                                                                             >  4:326968/1‑151 (MQ=255)
                                                                         cAAGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccac                                                                          >  1:760382/1‑151 (MQ=40)
                                                                          aaGCGAATCCGCAATTAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTccagaacgga                                                                          >  6:468710/1‑142 (MQ=255)
                                                                           aGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccacc                                                                         >  3:1032541/1‑150 (MQ=40)
                                                                           aGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCccaccacc                                                                         >  4:1798090/1‑150 (MQ=38)
                                                                           aGCGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGGCTAGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTACCacaccacc                                                                         >  4:620871/1‑150 (MQ=255)
                                                                             cGAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCcca                                                                              >  2:1263313/1‑143 (MQ=38)
                                                                              gAATCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCcca                                                                              >  4:22815/1‑142 (MQ=38)
                                                                               aaTCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATcgc                                                                    >  2:270117/1‑149 (MQ=38)
                                                                               aaTCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATccc                                                                    >  6:679672/1‑149 (MQ=38)
                                                                                aTCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTAGAGCCTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGCCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCt                                                                   >  8:1001386/1‑151 (MQ=255)
                                                                                 tCCGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACGAAGGTTCGGCTAGAGCGACAGCGAGAGAGCGTTGCCGCAGGCa                                                                                                                            >  4:127945/1‑93 (MQ=255)
                                                                                   cGCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttcc                                                                >  7:393717/1‑149 (MQ=37)
                                                                                    gCAATCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCttcct                                                               >  5:1086766/1‑148 (MQ=25)
                                                                                    gCAATCAAATATTCAGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCttcc                                                                >  5:1531399/1‑148 (MQ=255)
                                                                                       aTCAAATATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCttactaa                                                             >  6:836257/1‑145 (MQ=255)
                                                                                           aaTATTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCttccttaaaga                                                         >  3:788504/1‑141 (MQ=255)
                                                                                            atatTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTACataaaaaa                                                        >  2:928372/1‑146 (MQ=255)
                                                                                             tatTCTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGa                                                                                                                        <  1:1258562/85‑1 (MQ=40)
                                                                                              atTCTGCCCATGCGGGGACGTATGAGACGCTTCGACTACGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCg                                                      <  8:1761744/150‑1 (MQ=255)
                                                                                                 gtGCCCATGCGGGGAAGGATGAGAAGCTTCGACCAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAAt                                                  <  4:1045374/150‑1 (MQ=21)
                                                                                                 cTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGATGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCCCCCCCTTC‑CGTTCCAAAGCGCCGGaaat                                                 >  3:294337/1‑149 (MQ=14)
                                                                                                 cTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTCCaaaacacaccaaa                                                   >  1:1569072/1‑141 (MQ=255)
                                                                                                 cTGCCCATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTAGAGCGACAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTACtaaaaaaatcaaa                                                   >  2:195417/1‑141 (MQ=11)
                                                                                                   gcgcATGCGGGTAAGTATGAGAAGCTTCGACTAAGGTTCGGCTCTAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAAtt                                                 <  8:1370776/147‑1 (MQ=255)
                                                                                                   gCCCATGCGGGGAAGGATGAGAAGCTTCGACAACGGTTCGCCTCGAGCGCCCGCGAGAGAGCGTGGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGaa                                                   <  1:912043/148‑1 (MQ=11)
                                                                                                    cccATGCGGGGAAGGATGACAAGCGTCCCCCAAGCTTCGCCTCTATCTTTATCTAGAGATCGTTGCAGCAGGCAAATACACTAATTTCTAATCGCTCATCGCTTAGTAATCCTTCATCCAACACCATCACTTTCAAAAGCCCCGGAATtct                                               <  2:1350926/151‑1 (MQ=255)
                                                                                                     ccttgggggtatgtaTGATACGCTTCGACTAAGTTTCGCTTCGAGCGGCAGCTAGTTTGCGTTGCCGCAGGAAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATtctc                                              <  8:226181/137‑1 (MQ=14)
                                                                                                     ccATGCGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTCCaaaacaaaacaaaaactc                                              >  7:1427838/1‑137 (MQ=11)
                                                                                                      atttcGGGGAAGGATGAGAAGCTTCGACTAAGGTTCGCCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCa                                             <  5:1738566/147‑1 (MQ=255)
                                                                                                         gCGGGGAAGGATGAGAAGCTTCGACAAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAAc                                          <  5:122838/151‑1 (MQ=255)
                                                                                                          gggggAAGGATGAGAAGCTTCGACTAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACg                                         <  5:664807/150‑1 (MQ=255)
                                                                                                           ggggTAGCATGAGACGCTGCGACTACGGTTCGCCTCGAGCGCCAGCGAGAGAGGGTGGCCGCAGGAAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGa                                        <  1:1885930/151‑1 (MQ=14)
                                                                                                            gggACGGATGAGACGCTTCGACTACGGTTCGGCTCGAGCTTCATTTTTTTATCGTTGCGGCAGGCAACGACCCGAAGGGCTAAGCGAGCAGCGCTTAGTAATCCTTCGCCCATCACATTCACTTTCAAAAGCCCCGGAATTCTCAAACGaa                                       <  6:1517323/151‑1 (MQ=9)
                                                                                                                     gagaAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAAt                               <  2:1790583/150‑1 (MQ=255)
                                                                                                                      agaAGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCTAAAAAACACGAATTCTGATACTAATCCCCAAtt                              >  5:1420014/1‑149 (MQ=11)
                                                                                                                       gaAGCTTCGACAAAGGTTCGGCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCaa                            <  3:894458/151‑1 (MQ=255)
                                                                                                                        aaGCTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCGTTCAACCACCACCATCACTTTAAAAAGCCAAGGAATTCTCAAACTAATCCGCAATCaa                            >  8:456313/1‑150 (MQ=14)
                                                                                                                           cTTCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGACAACGAAATCTCAAACTaagagccaataaaaac                         >  1:1508416/1‑136 (MQ=11)
                                                                                                                            ttCGACTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTACAAAAGAAAACGAAATCTCAAAAGCATCCGCAATCAAataat                       >  8:1086651/1‑149 (MQ=14)
                                                                                                                             tCGACTAAGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCTACTAACATAATTGCGAAGCGCTCAGCGATGATTAATACTTTTCCCACGAATTTCACTTTCTAAAGGCCCGGTATTCTCAAACGAATCCGCAATCAAATAtt                       <  6:758919/150‑1 (MQ=9)
                                                                                                                             tCGACAAAGGTTCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAatat                        <  4:832281/149‑1 (MQ=255)
                                                                                                                              cGACAAAGGTTCGGCTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTc                      <  3:1601636/150‑1 (MQ=255)
                                                                                                                              actactctcttCGGATCTAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCt                     <  6:245194/142‑1 (MQ=14)
                                                                                                                                aCTAAGGTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCCCTTTCAAAAGCACAGGAATTCTCAAACGCATACGGAATTAAATAATCTgg                   >  1:1369343/1‑150 (MQ=14)
                                                                                                                                aCTAAGGTTCGACTCGAGCGCCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTg                    <  6:328030/150‑1 (MQ=25)
                                                                                                                                 ccacggTGCGCCTCGAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCTCTCATCTCTGAGTAATCCTTCCCCCACCAACATCACTTTCAAAAGCCCCTGAATTCTCAAACGAATCCGCAATCAAATATTCTGcc                  <  4:1507098/147‑1 (MQ=11)
                                                                                                                                     ggTTCGGCTCGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCaac              >  5:535307/1‑150 (MQ=255)
                                                                                                                                              cGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTTCAAATCCCAAGCAAAGCGCAAAGCAATCCACAAACATTTATTCTGcaccacaggcgactg     >  3:1295260/1‑137 (MQ=2)
                                                                                                                                              cGAGCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAgg     >  7:1619813/1‑151 (MQ=255)
                                                                                                                                               gAGCGGCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAAGTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGa    <  5:1857759/151‑1 (MQ=255)
                                                                                                                                                 gCGTCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCACTTTCAAAAGCCCCGGAATTCTCAAACGAATCCGCAATCAAATATTCTGCCCAAGCGGGGAAGGATg  <  5:276475/151‑1 (MQ=255)
                                                                                                                                                    tCAGCGAGAGAGCGTTGCCGCAGGCAACGACCCGAAGGGCGAAGCGCGCAGCGCTGAGTAATCCTTCCCCCACCACCATCGCTTTCAAAAGCCACCGAATTCTCAAACGAATCCGCAATCAAATAATCTGCCCAAGCGGGGGAGGATg  >  7:733999/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
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: ‑