Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 839,095 Δ1 bp coding (1001/1434 nt) bglA ← 6‑phospho‑beta‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH1839,0940T.98.8% 2217.2 / 18.5 488coding (1002/1434 nt)bglA6‑phospho‑beta‑glucosidase
Reads supporting (aligned to +/- strand):  ref base T (4/2);  new base . (227/253);  total (232/256)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

Only 100 of 504 total aligned reads displayed.

CACTTTGTCGTAAGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTAAATCCGCGACGCTC  >  DH1/838944‑839240
                                                                                                                                                      |                                                                                                                                                  
cACTTTGTCGTAAGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCt                                                                                                                                                    >  2:111112/1‑151 (MQ=255)
 acctgggcGTAAGCGCCAAAACCGGCTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCACAGGGCATCGCGCAGACCTACTGGATCAATCTGCCAGCCCCCGTCAGACGCTTTAACATACGGGGGTGGTACGCCGCCtc                                                                                                                                                   <  2:725773/144‑2 (MQ=255)
 aCTTTGTCGTAAGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCtc                                                                                                                                                   <  4:1369073/151‑2 (MQ=255)
   tttGTCGTAAGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCtca                                                                                                                                                  >  6:1053886/1‑148 (MQ=255)
   tttGTCGTAAGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCtca                                                                                                                                                  >  2:1411167/1‑148 (MQ=255)
    ttGTCGTAAGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCAAAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCtca                                                                                                                                                  <  3:760014/149‑3 (MQ=255)
     tGTCGTAAGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCCGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCtcaaa                                                                                                                                                >  3:1877383/1‑146 (MQ=255)
      gTCGTAAGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCaaaa                                                                                                                                              <  1:738362/150‑1 (MQ=255)
        cGTAAGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCa                                                                                                                                           <  7:1249775/151‑1 (MQ=255)
        cGTAAGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCa                                                                                                                                           <  7:355616/151‑1 (MQ=255)
        cGTAAGCGCCAAAACCGTTATCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCa                                                                                                                                           <  1:1674850/151‑1 (MQ=255)
         gTAAGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCag                                                                                                                                          <  8:1687453/151‑1 (MQ=255)
           aaccgcCAAAACCGTTTTCGACAATAAACAGCGCCCACTCAAAACGCTCATACAGTTCCCAAACTTCATAACGAAGACCTACTCGATCAATATCCCACCCCCACACACATGCTATAACATAAGCCTAAGGAACCCTCCCTCAAAACCaga                                                                                                                                         >  6:381314/4‑150 (MQ=255)
            aGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCag                                                                                                                                          >  2:572327/1‑148 (MQ=255)
             gcgcCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCag                                                                                                                                          >  6:698846/1‑147 (MQ=255)
                 cAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGatcgca                                                                                                                                  <  4:841880/151‑1 (MQ=255)
                 cAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGatcgca                                                                                                                                  <  7:1204983/151‑1 (MQ=255)
                    aaCCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGatcacatcg                                                                                                                               >  3:1205909/1‑151 (MQ=255)
                      ccGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGcc                                                                                                                             >  5:855746/1‑151 (MQ=255)
                       cGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCg                                                                                                                            <  3:130537/151‑1 (MQ=255)
                        gTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCgg                                                                                                                           >  8:55920/1‑151 (MQ=255)
                          tttCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTg                                                                                                                         >  7:1670742/1‑151 (MQ=255)
                           ttCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTg                                                                                                                         <  4:249342/150‑1 (MQ=255)
                             cGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTgccg                                                                                                                      <  8:1420240/151‑1 (MQ=255)
                             cGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTgccg                                                                                                                      <  5:1764202/151‑1 (MQ=255)
                                    aaaCAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCg                                                                                                                <  7:1357728/150‑1 (MQ=255)
                                     aaCAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGc                                                                                                              <  5:956136/151‑1 (MQ=255)
                                        aGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACCCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCtt                                                                                                           >  6:1932477/1‑151 (MQ=255)
                                              ctctGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGc                                                                                                     >  6:1987421/1‑151 (MQ=255)
                                                  gATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCa                                                                                                    <  6:1348527/148‑1 (MQ=255)
                                                   aTAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTgg                                                                                                >  6:625315/1‑151 (MQ=255)
                                                     aaCGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTc                                                                                              <  6:1228022/151‑1 (MQ=255)
                                                     aaCGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTc                                                                                              <  6:871843/151‑1 (MQ=255)
                                                     aaCGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTCGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGt                                                                                               <  3:1748662/150‑1 (MQ=255)
                                                        gCTCATACAGTTCGCAAAGTGCATCGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATg                                                                                           <  4:1843949/151‑1 (MQ=255)
                                                        gCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATg                                                                                           <  1:272697/151‑1 (MQ=255)
                                                        gCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATg                                                                                           <  2:128887/151‑1 (MQ=255)
                                                         cTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGt                                                                                          >  8:1849906/1‑151 (MQ=255)
                                                           cATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTa                                                                                         >  2:258124/1‑150 (MQ=255)
                                                            aTACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAAt                                                                                       <  7:1810057/151‑1 (MQ=255)
                                                            aTACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTGTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAAt                                                                                       >  3:1001610/1‑151 (MQ=255)
                                                              aCAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAAt                                                                                       >  6:7951/1‑150 (MQ=255)
                                                                aGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCt                                                                                   <  2:610963/151‑1 (MQ=255)
                                                                aGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCt                                                                                   <  4:344289/151‑1 (MQ=255)
                                                                 gTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTg                                                                                  <  6:1631068/151‑1 (MQ=255)
                                                                   tCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGa                                                                                 <  1:561468/150‑1 (MQ=255)
                                                                   tCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGa                                                                                 <  1:940648/150‑1 (MQ=255)
                                                                    cGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGACTTCACTGCATTGGTCATGTAATAGCTGaa                                                                                >  2:1720178/1‑150 (MQ=255)
                                                                     gCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTAAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAAc                                                                              <  8:748253/151‑1 (MQ=255)
                                                                      cAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGCTCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGACAcc                                                                             >  3:1808803/1‑151 (MQ=255)
                                                                       aacctGCATAGCGCAGACCAACTGGACCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTGTGGTACGCGGCCTCAAACCCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCa                                                                            <  4:1683164/147‑1 (MQ=255)
                                                                       aaaGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAAcc                                                                             >  1:1432567/1‑150 (MQ=255)
                                                                        aaGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCa                                                                            >  3:1401010/1‑150 (MQ=255)
                                                                         aGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAg                                                                          <  7:1585939/151‑1 (MQ=255)
                                                                          gTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGa                                                                         <  7:1887771/151‑1 (MQ=255)
                                                                          gTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGa                                                                         <  5:489473/151‑1 (MQ=255)
                                                                            gCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCTGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCTTCGCCGGTGCCGCCTTGGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGAt                                                                        <  1:1754046/150‑1 (MQ=255)
                                                                                aGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATc                                                                    <  6:259869/150‑1 (MQ=255)
                                                                                 ccgcAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATCCGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGc                                                                  <  6:499368/150‑1 (MQ=255)
                                                                                   gcAGACCTACTGGCTCGATCTGCCCGCCCCAGTCAGATGCTTTAACATACGGGTGTGGTCCGCTGCCTTCAAAACCAGAGATCGCATCGCCGGTGCGGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCa                                                                 <  2:428794/151‑1 (MQ=255)
                                                                                   gcAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGGCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAg                                                                <  3:1618007/151‑1 (MQ=255)
                                                                                    cAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAgg                                                               >  8:1030954/1‑151 (MQ=255)
                                                                                     aGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGt                                                              >  3:345373/1‑151 (MQ=255)
                                                                                         cTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAACACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCt                                                          >  8:1212877/1‑151 (MQ=255)
                                                                                          tACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCtt                                                         >  3:220059/1‑151 (MQ=255)
                                                                                           aCTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGCGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCCGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGAGCCTTc                                                        >  4:1536668/1‑151 (MQ=255)
                                                                                           aCTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAAAATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTc                                                        >  1:1010209/1‑151 (MQ=255)
                                                                                              ggATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATATGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACg                                                     >  5:1664395/1‑151 (MQ=255)
                                                                                              ggATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCAc                                                      <  4:461877/150‑1 (MQ=255)
                                                                                              ggATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACg                                                     >  4:1116631/1‑151 (MQ=255)
                                                                                               gATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACgc                                                    >  4:715227/1‑151 (MQ=255)
                                                                                               gATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACgc                                                    >  1:1958809/1‑151 (MQ=255)
                                                                                                aTCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACgc                                                    >  7:1470702/1‑150 (MQ=255)
                                                                                                 tCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACgca                                                   <  3:388432/150‑1 (MQ=255)
                                                                                                 tCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGCTGCCTTCACgc                                                    >  2:1515009/1‑149 (MQ=255)
                                                                                                  cAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACgcagc                                                 >  4:948653/1‑151 (MQ=255)
                                                                                                    aTCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGcac                                               >  5:776812/1‑151 (MQ=255)
                                                                                                        gCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATcc                                           >  1:598951/1‑151 (MQ=255)
                                                                                                        gCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATc                                            <  8:429112/150‑1 (MQ=255)
                                                                                                         ccAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCa                                          <  8:1343769/151‑1 (MQ=255)
                                                                                                         ccAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCa                                          <  2:1935292/151‑1 (MQ=255)
                                                                                                           aGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATcc                                           >  1:405028/1‑148 (MQ=255)
                                                                                                             cccctctgagttgatttAACATACGCGTTTGTTACGCTGCCTCAAAACCAGAGATCGCAGCGCGGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAAGCAAGATAATCGGAGGTGGCTTCACGCAGGTCATCCAGATc                                      <  8:1506939/137‑1 (MQ=255)
                                                                                                             ccccAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTGCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATc                                      <  7:490707/151‑1 (MQ=255)
                                                                                                               ccAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGc                                    >  7:1020694/1‑151 (MQ=255)
                                                                                                                 aGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCg                                  >  6:1148974/1‑151 (MQ=255)
                                                                                                                         gCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCAt                          <  5:697919/151‑1 (MQ=255)
                                                                                                                            ttATCATACTGTTTTGGTACGCTGCCTCAAACCGAGAGATCTCATCGCCGGTGCCGCCTTCGGCCTTCAGTGCATTGGTCATGTAATAGCTGAAACCAAGATAATGGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCAttt                        <  5:381629/150‑1 (MQ=255)
                                                                                                                            ttAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCAttt                        <  7:253707/150‑1 (MQ=255)
                                                                                                                              aaCATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTg                      <  8:296976/150‑1 (MQ=255)
                                                                                                                              aaCATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGa                     <  3:913432/151‑1 (MQ=255)
                                                                                                                               aCATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGAt                    >  3:1522777/1‑151 (MQ=255)
                                                                                                                                   aCGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGtt                 <  1:483857/150‑1 (MQ=255)
                                                                                                                                     gggTTTGGTACGCTGCCTTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTa                >  5:1353704/1‑150 (MQ=255)
                                                                                                                                      ggTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAACATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTAAAt             <  8:838914/151‑1 (MQ=255)
                                                                                                                                            gTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTAAATCCGCGa       <  1:1143974/151‑1 (MQ=255)
                                                                                                                                              aCGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTAAATCCGCGACg     >  7:731671/1‑151 (MQ=255)
                                                                                                                                              aCGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTAAATCCGCGACg     <  5:1052428/151‑1 (MQ=255)
                                                                                                                                              aCGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTAAATCCGCGACg     >  1:135553/1‑151 (MQ=255)
                                                                                                                                                gCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTAAATCCGCGACGCt   >  5:1964265/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CACTTTGTCGTAAGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTAAATCCGCGACGCTC  >  DH1/838944‑839240

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: