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% 1966.1 / 19.8 436coding (1002/1434 nt)bglA6‑phospho‑beta‑glucosidase
Reads supporting (aligned to +/- strand):  ref base T (4/1);  new base . (183/244);  total (190/246)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

Only 100 of 452 total aligned reads displayed.

CACTTTGTCGTAAGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTAAATCCGCGACGCTCCCAC  >  DH1/838944‑839244
                                                                                                                                                      |                                                                                                                                                      
cACTTTGTCGTAAGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATTCGGGTTTGGTACGCTGCCt                                                                                                                                                        >  5:968825/1‑151 (MQ=255)
  cTTTGTCGTAAGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTCACATACGGGTTTGGTACGCTGCCtc                                                                                                                                                       >  7:1238841/1‑149 (MQ=255)
     tGTCGTAAGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCtcaaa                                                                                                                                                    >  5:1175111/1‑146 (MQ=255)
        cGTAAGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTTCAAAAcc                                                                                                                                                <  4:658758/151‑1 (MQ=255)
        cGTAAGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCa                                                                                                                                               <  1:1436328/151‑1 (MQ=255)
            aGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAgaga                                                                                                                                           <  3:409539/151‑1 (MQ=255)
            aGCGCCAAAACCGATTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAgaga                                                                                                                                           <  2:1000846/151‑1 (MQ=255)
               gcCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCag                                                                                                                                              >  5:1564813/1‑145 (MQ=255)
               gcCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGatcg                                                                                                                                        <  2:683265/151‑1 (MQ=255)
                    aaCCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGatcgcatcg                                                                                                                                   <  1:1705778/151‑1 (MQ=255)
                       cGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCg                                                                                                                                >  2:38102/1‑151 (MQ=255)
                        gTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCgg                                                                                                                               <  6:631662/151‑1 (MQ=255)
                        gTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGATGCCTCAAAACCAGAGATCGCATCGCCgg                                                                                                                               >  2:433307/1‑151 (MQ=255)
                         ttttCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGt                                                                                                                              >  5:1304983/1‑151 (MQ=255)
                            tCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTTCAAAACCAGAGATCGCATCGCCGGTgc                                                                                                                            >  7:1355114/1‑151 (MQ=255)
                             cGACAATAAACAGCGGACTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCAAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGACAACACATCGcctgacccc                                                                                                                          >  8:1282419/1‑144 (MQ=255)
                                cAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCt                                                                                                                       >  5:539421/1‑151 (MQ=255)
                                cAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGCATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCt                                                                                                                       <  5:1729729/151‑1 (MQ=255)
                                  aTAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTc                                                                                                                     >  8:671634/1‑151 (MQ=255)
                                    aaaCAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCg                                                                                                                    <  7:764049/150‑1 (MQ=255)
                                      aCAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCACAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGcc                                                                                                                 >  4:502305/1‑151 (MQ=255)
                                      aCAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATAAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCAGCGCCGGTGCCGCCTTCGGcc                                                                                                                 >  3:663194/1‑151 (MQ=255)
                                       cAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCt                                                                                                                <  4:843217/151‑1 (MQ=255)
                                         gCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTc                                                                                                              <  4:1041014/151‑1 (MQ=255)
                                          cGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCa                                                                                                             <  6:275779/151‑1 (MQ=255)
                                           ggCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCAc                                                                                                            >  3:25385/1‑151 (MQ=255)
                                             ccTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTg                                                                                                          >  5:1013559/1‑151 (MQ=255)
                                               tctGATAACGCTAATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGc                                                                                                         >  4:1577262/1‑150 (MQ=255)
                                                  gATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCa                                                                                                        <  5:758206/148‑1 (MQ=255)
                                                  gATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTg                                                                                                     >  3:1094463/1‑151 (MQ=255)
                                                  gATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTGTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCAtt                                                                                                      <  5:650660/150‑1 (MQ=255)
                                                   aTAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGTATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTgg                                                                                                    <  6:1343351/151‑1 (MQ=255)
                                                     aaCGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTc                                                                                                  <  6:453073/151‑1 (MQ=255)
                                                      ncGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTc                                                                                                  >  5:1671719/2‑150 (MQ=255)
                                                       cGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCAt                                                                                                >  3:527011/1‑151 (MQ=255)
                                                        gCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATg                                                                                               <  8:415360/151‑1 (MQ=255)
                                                        gCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATg                                                                                               <  2:1387885/151‑1 (MQ=255)
                                                          tCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGt                                                                                              >  3:1415228/1‑150 (MQ=255)
                                                           aaTACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCGGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTaa                                                                                            <  1:1060198/150‑1 (MQ=255)
                                                            aTACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAAt                                                                                           >  7:602061/1‑151 (MQ=255)
                                                             tACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTCGGTACGCTGCCTTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAAt                                                                                           >  2:1610593/1‑151 (MQ=255)
                                                              aCAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGCTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAg                                                                                         <  6:1530350/151‑1 (MQ=255)
                                                                aGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCt                                                                                       <  1:678186/151‑1 (MQ=255)
                                                                 gTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTg                                                                                      >  8:1369868/1‑151 (MQ=255)
                                                                 gTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTg                                                                                      <  2:845948/151‑1 (MQ=255)
                                                                  ttCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTg                                                                                      >  1:1173825/1‑150 (MQ=255)
                                                                    cGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGaa                                                                                    <  7:1139123/150‑1 (MQ=255)
                                                                    cGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCCCCGGTGCCTCCATCGGCCTTCACTGCATTGGTCATGTAATAGCTGaac                                                                                   >  8:314951/1‑150 (MQ=255)
                                                                      cAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAAcc                                                                                 <  4:1470896/151‑1 (MQ=255)
                                                                      aaaaGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAAcc                                                                                 <  1:1016319/150‑1 (MQ=255)
                                                                       aaaGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAAcc                                                                                 >  1:787819/1‑150 (MQ=255)
                                                                         aGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCa                                                                                >  7:1746557/1‑149 (MQ=255)
                                                                          gTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGa                                                                             <  2:1168599/151‑1 (MQ=255)
                                                                            gCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGAt                                                                            <  5:432394/150‑1 (MQ=255)
                                                                            gCATAGCGCAGACCTACTGGATCAATCTCCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGAt                                                                            >  5:1607313/1‑150 (MQ=255)
                                                                              aTAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAAt                                                                         <  5:837014/151‑1 (MQ=255)
                                                                                   gcAGACCTACTGGATCAATCTGCCAGCCCGAGTCAGATGCTTTAACATACTGGTTTGGTACGCTTCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCa                                                                     <  2:660826/150‑1 (MQ=255)
                                                                                   gcAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAg                                                                    <  2:555082/151‑1 (MQ=255)
                                                                                    cAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAgg                                                                   >  1:467310/1‑151 (MQ=255)
                                                                                     aGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGt                                                                  <  3:1025310/151‑1 (MQ=255)
                                                                                      gACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTg                                                                 >  1:1621664/1‑151 (MQ=255)
                                                                                        ccTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGcc                                                               >  6:482427/1‑151 (MQ=255)
                                                                                          aaCTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCtt                                                             <  6:787451/150‑1 (MQ=255)
                                                                                           aCTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCt                                                              <  8:8731/149‑1 (MQ=255)
                                                                                            cTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATAGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTc                                                            >  4:326207/1‑150 (MQ=255)
                                                                                              ggATCAATCTGCCAGCCCCAGTCAGTTGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCa                                                           <  6:302249/149‑1 (MQ=255)
                                                                                                  cAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTGCGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACgcag                                                      <  2:1696714/150‑1 (MQ=255)
                                                                                                   aaTCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACgcagc                                                     >  5:346008/1‑150 (MQ=255)
                                                                                                   aaTCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACgcagc                                                     >  2:873307/1‑150 (MQ=255)
                                                                                                       tGCCCGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATc                                                >  1:1498594/1‑151 (MQ=255)
                                                                                                        gCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATcc                                               <  2:1040514/151‑1 (MQ=255)
                                                                                                          cAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATGGCATCGCCGGTGGCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATGCAg                                             <  2:1547397/151‑1 (MQ=255)
                                                                                                           aGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGGCTTCACTGCATTGGTCATGTAATAGCTGAAACAAAGATAATCGCAGGTGCCTTCACGCAGCACATcc                                               >  6:1104817/1‑148 (MQ=255)
                                                                                                            gCCCCAGTCAGGTGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGAt                                           <  3:1305502/151‑1 (MQ=255)
                                                                                                            gCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGAt                                           <  1:1473559/151‑1 (MQ=255)
                                                                                                             ccccAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATc                                          <  6:1248059/151‑1 (MQ=255)
                                                                                                             ccccAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATc                                          <  8:1454494/151‑1 (MQ=255)
                                                                                                              cccAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCg                                         <  3:924665/151‑1 (MQ=255)
                                                                                                                cAGTCTGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGcc                                       <  7:399255/151‑1 (MQ=255)
                                                                                                                cAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGCTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCACATCGcc                                       >  8:1567586/1‑151 (MQ=255)
                                                                                                                    cataTGCTTTAACATACGGGTTTGGTACGCTGACTAAAAAACAGAGATCGCATCGCCGGTGCCGACTTCGGCCTTCAATGCATTCGTCATGTAATACCTGAAACCAAGATAATCGAAAGTGACTTCACGCAGCACATCAAGATCGcagac                                    >  8:1096097/4‑146 (MQ=255)
                                                                                                                    cAGATGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCt                                   >  1:855835/1‑151 (MQ=255)
                                                                                                                       aTGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTcc                                <  6:1675337/151‑1 (MQ=255)
                                                                                                                       aTGCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTcc                                <  5:1452506/151‑1 (MQ=255)
                                                                                                                         gCTTTAACATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCa                               <  4:1752052/150‑1 (MQ=255)
                                                                                                                            ttAACATACGGGTTTGGTACTCTTCGTCATACCCAGAGAGCGCATCTCCGGTGCCGCCTTCGGCCTTCACTGGATTGGTCATGTAATAGCTGAAACCAAGATAATCTCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCAtt                             <  1:420700/149‑1 (MQ=255)
                                                                                                                              aaCATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTg                          >  8:58419/1‑150 (MQ=255)
                                                                                                                              aaCATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTg                          <  8:1225200/150‑1 (MQ=255)
                                                                                                                              aaCATACGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTg                          <  7:1414567/150‑1 (MQ=255)
                                                                                                                                 aTACGGGTGTGGTACGCTGCCTCAAAACCAGAGATCGCCTCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGt                      <  1:576867/151‑1 (MQ=255)
                                                                                                                                    cGGGTTTGGTACGCTGCCTCAAATCCAGAGATCGCATCGGCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGACGTTaa                   <  8:1031793/151‑1 (MQ=255)
                                                                                                                                    cGGGTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTaa                   <  7:100424/151‑1 (MQ=255)
                                                                                                                                     gggTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTaaa                  <  4:1366647/151‑1 (MQ=255)
                                                                                                                                     gggTGTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTaaa                  <  4:1299669/151‑1 (MQ=255)
                                                                                                                                       gTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTAAATc                <  7:950331/151‑1 (MQ=255)
                                                                                                                                       gTTTGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTAAATc                <  7:338585/151‑1 (MQ=255)
                                                                                                                                          tGGTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTAAATCcgc             >  7:660240/1‑151 (MQ=255)
                                                                                                                                           ggTACGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTAAATCcgc             <  8:1829445/150‑1 (MQ=255)
                                                                                                                                              aCGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTAAATCCGCGAc          <  7:1601236/150‑1 (MQ=255)
                                                                                                                                              aCGCTGCCTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTAAATCCGCGAc          <  2:1781835/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CACTTTGTCGTAAGCGCCAAAACCGTTTTCGACAATAAACAGCGGCCTCTGATAACGCTCATACAGTTCGCAAAGTGCATAGCGCAGACCTACTGGATCAATCTGCCAGCCCCAGTCAGATGCTTTAACATACGGGTTTGGTACGCTGCCTTCAAAACCAGAGATCGCATCGCCGGTGCCGCCTTCGGCCTTCACTGCATTGGTCATGTAATAGCTGAAACCAAGATAATCGCAGGTGCCTTCACGCAGCACATCCAGATCGCCGTCTTCCATTTTGATGTTAAATCCGCGACGCTCCCAC  >  DH1/838944‑839244

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: