Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,983,0080AG52.2% 61.8 / inf 522Q285Q (CAA→CAG) sgcXendoglucanase
Reads supporting (aligned to +/- strand):  ref base A (126/122);  new base G (134/137);  total (261/261)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 523 total aligned reads displayed.

ATCACTCCATCGTGTGACACCCCTGATTTACACGATTATTCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCTATACCCATTCTCCGGCGGAGGTTG  >  DH1/3982859‑3983156
                                                                                                                                                     |                                                                                                                                                    
 tCACTCCATCGTGTGACACCCCTGATTTACACGATTATTCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTTCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAAc                                                                                                                                                    >  2:148098/1‑151 (MQ=255)
  cACTCCATCGTGTGACACCCCTGATTTACACGATTATTCCTAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGAc                                                                                                                                                    <  7:415269/150‑1 (MQ=255)
   aCTCCATCGTGTGACACCCCTGATTTACACGATTATTCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACgg                                                                                                                                                  <  6:1402358/151‑1 (MQ=255)
        aTCGTGTGACACGCCTGATTTACAGGATTATTCCGAGGTCAGTATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCtt                                                                                                                                             <  7:1034935/151‑1 (MQ=255)
        aTCGTGTGACACCCCTGATTTACACGATTATTCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCtt                                                                                                                                             <  3:1223685/151‑1 (MQ=255)
        aTCGTGTGACACCCCTGATTTACACGATTATTCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCtt                                                                                                                                             >  5:652513/1‑151 (MQ=255)
        aTCGTGTGACACCCCTGATTTACACGATTATTCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCtt                                                                                                                                             >  1:323950/1‑151 (MQ=255)
         tCGTGTGACACCCCTGATTTACACGATTATTCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACAGATGTTGGAACAGACGGCTCTTg                                                                                                                                            >  1:1951184/1‑151 (MQ=255)
          cGTGTGACACCCCTGATTTACACGATTATTCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTg                                                                                                                                            >  4:1940927/1‑150 (MQ=255)
           gtgtGACACCCCTGATTTACACGATTATTCCGAGGTGAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGaa                                                                                                                                          <  6:1825573/151‑1 (MQ=255)
           gtgtGACACCCCTGATTTACACGATTATTCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGaa                                                                                                                                          <  6:1621224/151‑1 (MQ=255)
           gtgtGACACCCCTGATTTACACGATTATTCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGa                                                                                                                                           >  8:116498/1‑150 (MQ=255)
             gtgACACCCCTGATTTACACGATTATTCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAaca                                                                                                                                        <  2:1623284/151‑1 (MQ=255)
               gACACCCCTGATTTACACGATTATTCCGAGCTCAGGATTAATCAGGGCGTTGTGATCACATGCCTGAACTACCATGGTCGGGCAACGCTCTCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAaaac                                                                                                                                       >  5:745683/1‑150 (MQ=255)
                acacCCCTGATTTACACGATTATTCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGCATTAATCACGCCTCCTCGCCTGATACGGATCTTGGAACAGACGGCTCATGAAcaca                                                                                                                                      >  6:1971683/1‑150 (MQ=255)
                  acCCCTGATTTACACGATTATTCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACacaac                                                                                                                                    >  3:1906294/1‑150 (MQ=255)
                  acCCCTGATTTACACGATTATTCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACacaac                                                                                                                                    >  3:1383716/1‑150 (MQ=255)
                        gATTTACACGATTATTCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCgg                                                                                                                             >  4:667173/1‑151 (MQ=255)
                        gATTTACACGATTATTCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCgg                                                                                                                             <  2:1312479/151‑1 (MQ=255)
                             agacGATTATTCCGAGGTCAGGATTAATGAGGGCGTTGGGATCACCTGCCTGAACTAGGATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCGTCGCCTGATAGGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAg                                                                                                                        <  4:329168/149‑1 (MQ=255)
                             acacGATTATTCCGAGGTCAGGATTAAGCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAg                                                                                                                        <  7:708980/151‑1 (MQ=255)
                               acGATTATTCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCg                                                                                                                      <  1:244488/151‑1 (MQ=255)
                                    tattCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCgaga                                                                                                                   >  4:1044311/1‑149 (MQ=255)
                                      ttCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTg                                                                                                               >  3:1784212/1‑151 (MQ=255)
                                       tCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTgg                                                                                                              >  5:120779/1‑151 (MQ=255)
                                        gcGAGGTCTGTATTAATGAGGGCGTTGGGATCTCGTGCCTGAACTACCATGGTGGGGGAAGGCTGGCGGGATTATTCACGGCGCGTCGCCTGATACGGATGTTGGAACAGTCGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTgg                                                                                                              <  2:739327/149‑1 (MQ=255)
                                         cGAGGTCAGGATTAATCAGGGCGTTCGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCt                                                                                                            >  1:1281052/1‑151 (MQ=255)
                                                 ggATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGt                                                                                                    <  3:1097033/151‑1 (MQ=255)
                                                 ggATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGAGAGAAGTGGCTCCCGGCGt                                                                                                    >  6:897943/1‑151 (MQ=255)
                                                   aTTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTg                                                                                                   >  1:797402/1‑150 (MQ=255)
                                                         cAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTACTCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACa                                                                                             <  8:159168/150‑1 (MQ=255)
                                                         cAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACa                                                                                             <  7:636097/150‑1 (MQ=255)
                                                              cGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAAcc                                                                                       >  4:909564/1‑151 (MQ=255)
                                                              cGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAAcc                                                                                       >  2:1089496/1‑151 (MQ=255)
                                                               gTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCg                                                                                      >  2:645188/1‑151 (MQ=255)
                                                                 tGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGGGTGATAACAGAATCCGGt                                                                                    >  2:359803/1‑151 (MQ=255)
                                                                   ggATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTta                                                                                  <  8:456745/151‑1 (MQ=255)
                                                                    gATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTtat                                                                                 <  2:223705/151‑1 (MQ=255)
                                                                      ccACCTGCCTGAACTACCATGGGCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTatat                                                                               <  6:662234/150‑1 (MQ=255)
                                                                       cACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACACAAACCGGTTATATc                                                                              >  1:1603274/1‑151 (MQ=255)
                                                                       cACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATc                                                                              >  2:445721/1‑151 (MQ=255)
                                                                       cACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATc                                                                              <  2:1197589/151‑1 (MQ=255)
                                                                       aaCCTGCCTGAACTACCATGGTCGGGGAACGCTGTCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATc                                                                              <  7:1301035/150‑1 (MQ=255)
                                                                        aCCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATcc                                                                             <  1:1900904/151‑1 (MQ=255)
                                                                        aCCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATcc                                                                             >  8:1329405/1‑151 (MQ=255)
                                                                         cctgcctgAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCa                                                                            <  3:752424/151‑1 (MQ=255)
                                                                         cctgcctgAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCa                                                                            <  3:361958/151‑1 (MQ=255)
                                                                            gcctgAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGGACAAACGGCTCTTGAACACAACATTCCGGTGCCGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGt                                                                         <  1:54795/151‑1 (MQ=255)
                                                                             cctgAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCACCCGGCGTGATAACAGAAACCGGTTATATCCAGGtt                                                                        <  4:1895948/151‑1 (MQ=255)
                                                                               ggAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGa                                                                      <  8:398352/150‑1 (MQ=255)
                                                                                gAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTg                                                                       <  2:976198/149‑1 (MQ=255)
                                                                                gAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTg                                                                       <  3:610408/149‑1 (MQ=255)
                                                                                  aCTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCa                                                                   <  1:1673070/151‑1 (MQ=255)
                                                                                  aCTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCa                                                                   <  5:1880087/151‑1 (MQ=255)
                                                                                   cTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCac                                                                  >  1:1070361/1‑150 (MQ=255)
                                                                                    tACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAgg                                                                 >  3:540583/1‑151 (MQ=255)
                                                                                      ccATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGAt                                                               <  5:312053/151‑1 (MQ=255)
                                                                                       cATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATg                                                              <  1:1789124/151‑1 (MQ=255)
                                                                                        aTGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACACACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATgg                                                             >  6:797170/1‑151 (MQ=255)
                                                                                          ggTCGGGGAACGCTGGCCGGAGTAATCACGCCTCCTCGGCTGATACGGAAGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTGATATCCAGGTTGAGCAGGATGGt                                                            <  8:1996807/150‑1 (MQ=255)
                                                                                             cGGGTAACGCTGTCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCATGTTGAGCAGGATGGTAttn                                                        <  4:1223469/151‑2 (MQ=255)
                                                                                                ggAACGCTGGGAGGATTAATAAGGCCTCCTCTCCTGATTCGGATGTTGGAACAAACTGCTCTTGTTCACAAGATTCCGGTGCAGCGAGAAGTGGCTCCTGGCGTGATAACTGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCt                                                     <  6:124635/151‑1 (MQ=255)
                                                                                                 gAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTg                                                    <  1:1388812/151‑1 (MQ=255)
                                                                                                 gAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTg                                                    <  8:1881561/151‑1 (MQ=255)
                                                                                                  aaCGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATGCAGGTTGAGCAGGATGGTATTCCCTgc                                                   <  7:982546/151‑1 (MQ=255)
                                                                                                   aCGCTGGCAGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAACCCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTgcg                                                  >  2:1586287/1‑151 (MQ=255)
                                                                                                     gCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGcc                                                >  7:1624472/1‑151 (MQ=255)
                                                                                                     gCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGcc                                                <  1:737562/151‑1 (MQ=255)
                                                                                                        ggCCGGATTAATCACGCCTCCTCTCCTGATACCGATGTTGGACCAGACGGCTCTTGAACACAACATTCCTGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAg                                              <  6:365829/150‑1 (MQ=255)
                                                                                                        ggCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGt                                             <  4:990532/151‑1 (MQ=255)
                                                                                                            ggATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGtctct                                         >  7:1880289/1‑151 (MQ=255)
                                                                                                            ggATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGtctct                                         >  4:536582/1‑151 (MQ=255)
                                                                                                              aTTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTctctct                                       >  3:200838/1‑151 (MQ=255)
                                                                                                                tAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTa                                      <  1:630444/150‑1 (MQ=255)
                                                                                                                  aTCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTc                                   >  4:1617640/1‑151 (MQ=255)
                                                                                                                   tCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTcc                                  >  7:1968351/1‑151 (MQ=255)
                                                                                                                      cGCCTCCTCGCCTGATACGGAGGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTg                               <  4:1271006/151‑1 (MQ=255)
                                                                                                                       gCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGt                              <  3:1592111/151‑1 (MQ=255)
                                                                                                                        cctcctCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTc                             >  5:1994810/1‑151 (MQ=255)
                                                                                                                         ctcctcGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCg                            >  5:843377/1‑151 (MQ=255)
                                                                                                                         ctcctcGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCg                            >  7:826400/1‑151 (MQ=255)
                                                                                                                             tcGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCtat                        >  4:1434845/1‑151 (MQ=255)
                                                                                                                              cGCCTGATACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCtat                        >  5:1770748/1‑150 (MQ=255)
                                                                                                                               gCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCTATAc                      <  6:1735881/151‑1 (MQ=255)
                                                                                                                               gCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCTATAc                      >  7:1785599/1‑151 (MQ=255)
                                                                                                                                 cTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCTATAccc                    >  5:1105154/1‑151 (MQ=255)
                                                                                                                                    attcGGTTGTTGGTACAGACGTCTCTTGAAGAGAACATTGCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACGGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCTATACCCAt                  <  6:453345/147‑1 (MQ=255)
                                                                                                                                    aTACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCTATACCCAtt                 <  7:1100841/151‑1 (MQ=255)
                                                                                                                                    aTACGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCTATACCCAtt                 >  4:1791508/1‑151 (MQ=255)
                                                                                                                                      aCGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCTATACCCATtct               >  1:472132/1‑151 (MQ=255)
                                                                                                                                      aCGGATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCTATACCCATtc                <  5:1621960/150‑1 (MQ=255)
                                                                                                                                       cGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCTATACCCATtctc              <  3:875067/151‑1 (MQ=255)
                                                                                                                                         gATGTTGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCTATACCCATTCTCcg            <  2:1949751/151‑1 (MQ=255)
                                                                                                                                         gATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCTATACCCATTCTCcg            >  3:144244/1‑151 (MQ=255)
                                                                                                                                            gTTGGAACAAACGGCTCTTGAACACAACATTCGGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAAGCGGTTATATCCTGGTTGAGCAGGATCGTATTCCCTGCGCCTGTCTCTCTCTTACTTGTCGCTATACCCATTCTCcggcg         <  6:1140662/151‑1 (MQ=255)
                                                                                                                                              tGGAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCTATACCCATTCTCcggcgg        <  2:1876100/150‑1 (MQ=255)
                                                                                                                                               ggAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCTATACCCATTCTCCGGCGGa       <  7:1210165/150‑1 (MQ=255)
                                                                                                                                               ggAACAGACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCTATACCCATTCTCCGGCGGAg      <  5:1524480/151‑1 (MQ=255)
                                                                                                                                               ggAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCTATACCCATTCTCcggcgg        >  8:636339/1‑149 (MQ=255)
                                                                                                                                                 aaCAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCTATACCCATTCTCCGGCGGAGGt    >  1:96435/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATCACTCCATCGTGTGACACCCCTGATTTACACGATTATTCCGAGGTCAGGATTAATCAGGGCGTTGGGATCACCTGCCTGAACTACCATGGTCGGGGAACGCTGGCCGGATTAATCACGCCTCCTCGCCTGATACGGATGTTGGAACAAACGGCTCTTGAACACAACATTCCGGTGCAGCGAGAAGTGGCTCCCGGCGTGATAACAGAAACCGGTTATATCCAGGTTGAGCAGGATGGTATTCCCTGCGCCAGTCTCTCTATTCCTTGTCGCTATACCCATTCTCCGGCGGAGGTTG  >  DH1/3982859‑3983156

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑