Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 2,687,814 G→A S21N (AGC→AAC)  purB → adenylosuccinate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,687,8140GA96.6% 1943.7 / 31.9 557S21N (AGC→AAC) purBadenylosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base G (12/7);  new base A (277/260);  total (290/267)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 575 total aligned reads displayed.

GGACGTCTGCAACTGATGTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAGCGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCTTTTGCTGCCGACGCAATCGGTTACCTTGATGCAAT  >  DH1/2687666‑2687961
                                                                                                                                                    |                                                                                                                                                   
ggACGTCTGCAACTGATGTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAacg                                                                                                                                                   <  4:64435/151‑1 (MQ=255)
 gACGTCTGCAACTGATGTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAacgc                                                                                                                                                  >  4:943217/1‑151 (MQ=255)
  aCGTCTGCAACTGATGTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAacgcg                                                                                                                                                 <  1:2008608/151‑1 (MQ=255)
  aCGTCTGCAACTGATGTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAacgcg                                                                                                                                                 >  8:593157/1‑151 (MQ=255)
   cGTCTGCAACTGATGTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAacgcgc                                                                                                                                                >  3:1137161/1‑151 (MQ=255)
    gTCTGCAACTGATGTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCt                                                                                                                                               <  3:1510979/151‑1 (MQ=255)
     tCTGCAACTGATGTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTg                                                                                                                                              >  5:487260/1‑151 (MQ=255)
      cTGCAACTGATGTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTgc                                                                                                                                             >  8:422249/1‑151 (MQ=255)
        gCAACTGATGTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTgcg                                                                                                                                            <  2:1705831/150‑1 (MQ=255)
         cAACTGATGTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGcgcg                                                                                                                                          >  6:401009/1‑151 (MQ=255)
           aCTGATGTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCggg                                                                                                                                        <  2:107034/151‑1 (MQ=255)
            cTGATGTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCggg                                                                                                                                        >  1:106308/1‑150 (MQ=255)
            cTGATGTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCggg                                                                                                                                        >  5:1633431/1‑150 (MQ=255)
            cTGATGTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCggg                                                                                                                                        >  4:1274138/1‑150 (MQ=255)
             tGATGTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGAt                                                                                                                                      >  7:912863/1‑151 (MQ=255)
              gATGTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGa                                                                                                                                       <  8:225026/149‑1 (MQ=255)
              gATGTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGAt                                                                                                                                      <  6:1839392/150‑1 (MQ=255)
                 gTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTc                                                                                                                                  <  7:203854/151‑1 (MQ=255)
                 gTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTc                                                                                                                                  >  3:688775/1‑151 (MQ=255)
                  ttttCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTc                                                                                                                                  >  4:1571159/1‑150 (MQ=255)
                    ttCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAg                                                                                                                                <  5:1952582/150‑1 (MQ=255)
                     tctcGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGc                                                                                                                               <  7:796733/150‑1 (MQ=255)
                        cGTNNTCGCNTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAg                                                                                                                                <  7:1650832/146‑1 (MQ=255)
                        cGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGaa                                                                                                                            <  6:1074366/150‑1 (MQ=255)
                        cGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAat                                                                                                                           <  4:264158/151‑1 (MQ=255)
                        cGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGCTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGaa                                                                                                                            <  1:1623694/150‑1 (MQ=255)
                         gTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAat                                                                                                                           >  8:1918444/1‑150 (MQ=255)
                           aaTCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAGCGCGCTGCGCGGGATTTTCAGCGAATATg                                                                                                                        >  7:1284117/1‑151 (MQ=255)
                           aaTCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATg                                                                                                                        <  5:38940/151‑1 (MQ=255)
                           aaTCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGAGCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAatat                                                                                                                         <  4:1464066/150‑1 (MQ=255)
                            aTCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATgg                                                                                                                       >  6:745293/1‑151 (MQ=255)
                              cGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGtt                                                                                                                     <  5:832907/151‑1 (MQ=255)
                              cGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGt                                                                                                                      <  1:1803848/150‑1 (MQ=255)
                              cGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGt                                                                                                                      <  3:37969/150‑1 (MQ=255)
                              cGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGt                                                                                                                      <  5:585862/150‑1 (MQ=255)
                                 cTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGc                                                                                                                  >  8:1219470/1‑151 (MQ=255)
                                 cTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGc                                                                                                                  >  5:1526320/1‑151 (MQ=255)
                                   gACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCGCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTg                                                                                                                <  1:255813/151‑1 (MQ=255)
                                     ccACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGaa                                                                                                              <  3:482294/151‑1 (MQ=255)
                                         tCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCCACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTc                                                                                                          >  1:1914342/1‑151 (MQ=255)
                                          cAGGCAAAAGAAATTCTTGCTCATTTAACCCCTGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTc                                                                                                          <  2:1309724/150‑1 (MQ=255)
                                          cAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAGCGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTcc                                                                                                         >  8:686984/1‑151 (MQ=255)
                                          cAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTcc                                                                                                         <  3:1596321/151‑1 (MQ=255)
                                          cAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTcc                                                                                                         <  7:637583/151‑1 (MQ=255)
                                             gCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCgtg                                                                                                      <  4:1757835/151‑1 (MQ=255)
                                               aaaaCAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCgtgt                                                                                                     >  3:1705595/1‑150 (MQ=255)
                                                aaaCAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCTTGTAc                                                                                                   >  1:401915/1‑151 (MQ=255)
                                                 aaCAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACa                                                                                                  <  1:602341/151‑1 (MQ=255)
                                                  aCAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACa                                                                                                  >  2:221559/1‑150 (MQ=255)
                                                  aCAAATTCTTGCTAATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCACCTGTCGATGGACGCTACCGCGATAAAGTCAGCGCGCTGCGCGGGATTTTAAGCCAATATGGTTTGATTAAATTCCGTGTACa                                                                                                  >  6:1963694/1‑150 (MQ=255)
                                                   cAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACa                                                                                                  >  4:825504/1‑149 (MQ=255)
                                                         cTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGa                                                                                            >  6:514311/1‑149 (MQ=255)
                                                             gtCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACg                                                                                      <  3:601540/150‑1 (MQ=255)
                                                              tCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGt                                                                                     >  5:1231860/1‑151 (MQ=255)
                                                                aTTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTg                                                                                   >  7:1190915/1‑151 (MQ=255)
                                                                aTTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTg                                                                                   >  1:1418068/1‑151 (MQ=255)
                                                                    aaCCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGATTCCCCTGTCGATGGACGCAACGGCGATAAAGTCAACGCGCTGCGCGGGATTATCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTg                                                                               >  4:517302/1‑151 (MQ=255)
                                                                     aCCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGc                                                                              <  2:1859053/151‑1 (MQ=255)
                                                                     aCCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGc                                                                              >  3:851697/1‑151 (MQ=255)
                                                                      ccccGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCAGTGTACAAGTTGAAGTACGTTGGCTGc                                                                              >  7:1049963/1‑150 (MQ=255)
                                                                         cGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTTCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCaaaa                                                                          <  7:1513538/151‑1 (MQ=255)
                                                                            aGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACt                                                                       >  7:317307/1‑151 (MQ=255)
                                                                                gtgATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGcc                                                                   <  7:1371462/151‑1 (MQ=255)
                                                                                  gATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTGAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCcgc                                                                 >  2:91723/1‑151 (MQ=255)
                                                                                      tATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGcgcac                                                             >  6:1574183/1‑151 (MQ=255)
                                                                                      aaTGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGcgcac                                                             <  8:1493771/150‑1 (MQ=255)
                                                                                       aTGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCcgc                                                                 <  8:1116756/146‑1 (MQ=255)
                                                                                       aTGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGcgca                                                              <  3:1946244/149‑1 (MQ=255)
                                                                                         ggAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGcgcacgc                                                           >  3:1334258/1‑150 (MQ=255)
                                                                                          gAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGcgcacgc                                                           >  5:1052862/1‑149 (MQ=255)
                                                                                             ttATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCg                                                       >  2:359894/1‑150 (MQ=255)
                                                                                              tATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGAt                                                     >  4:205760/1‑151 (MQ=255)
                                                                                               aTCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATc                                                    >  6:1316869/1‑151 (MQ=255)
                                                                                               aTCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATc                                                    <  1:520269/151‑1 (MQ=255)
                                                                                                tCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTTAAATTACGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATc                                                    >  8:1855666/1‑150 (MQ=255)
                                                                                                 ccTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCa                                                   <  1:903699/150‑1 (MQ=255)
                                                                                                    cACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAgg                                                >  3:957996/1‑150 (MQ=255)
                                                                                                    cACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGa                                               <  8:136504/151‑1 (MQ=255)
                                                                                                     aCTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGaa                                              <  7:1906479/151‑1 (MQ=255)
                                                                                                     aCTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGaa                                              <  6:622789/151‑1 (MQ=255)
                                                                                                        gACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGtt                                           <  7:1779054/151‑1 (MQ=255)
                                                                                                          ccgccgTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTcc                                         <  5:1362094/151‑1 (MQ=255)
                                                                                                               gTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCtt                                    <  4:1193849/151‑1 (MQ=255)
                                                                                                                 ttCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCtttt                                  >  3:256999/1‑151 (MQ=255)
                                                                                                                  tCCCCTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCTTTtg                                 >  5:2034527/1‑151 (MQ=255)
                                                                                                                   ccccTGTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCTTTtg                                 <  6:1687989/150‑1 (MQ=255)
                                                                                                                        gTCGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCTTTtgctgc                             <  2:440392/149‑1 (MQ=255)
                                                                                                                          cGATGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCTTTTGCTGCCGAc                         <  3:1015570/151‑1 (MQ=255)
                                                                                                                            aTGGACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCTTTTGCTGCCGACGc                       <  5:142810/151‑1 (MQ=255)
                                                                                                                               gACGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCTTTTGCTGCCGACGCAAt                    <  4:1136257/151‑1 (MQ=255)
                                                                                                                                aCGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCTTTTGCTGCCGACGCAATc                   <  2:719575/151‑1 (MQ=255)
                                                                                                                                 cGCTACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCTTTTGCTGCCGACGCAATCg                  >  1:1230482/1‑151 (MQ=255)
                                                                                                                                    tACGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCTTTTGCTGCCGACGCAATCGGtt               >  5:1888496/1‑151 (MQ=255)
                                                                                                                                     aCGGCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCTTTTGCTGCCGACGCAATCGGTTa              <  8:1706988/151‑1 (MQ=255)
                                                                                                                                       ggCGATAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCTTTTGCTGCCGACGCAATCGGTTAcc            >  6:122915/1‑151 (MQ=255)
                                                                                                                                           aTAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCTTTTGCTGCCGACGCAATCGGTTACCTTg         >  2:1561114/1‑150 (MQ=255)
                                                                                                                                           aTAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCTTTTGCTGCCGACGCAATCGGTTACCTTg         <  1:772025/150‑1 (MQ=255)
                                                                                                                                            tAAAGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCTTTTGCTGCCGACGCAATCGGTTACCTTGAt       >  2:1009762/1‑151 (MQ=255)
                                                                                                                                             aaaGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCTTTTGCTGCCGACGCAATCGGTTACCTTGATg      >  8:798075/1‑151 (MQ=255)
                                                                                                                                               aGTCAACGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCTTTTGCTGCCGACGCAATCGGTTACCTTGATGCa    <  7:463736/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGACGTCTGCAACTGATGTTTTCTCGTAATCGCCTGACCACTCAGGCAAAACAAATTCTTGCTCATTTAACCCCGGAGTTGTGATCTATGGAATTATCCTCACTGACCGCCGTTTCCCCTGTCGATGGACGCTACGGCGATAAAGTCAGCGCGCTGCGCGGGATTTTCAGCGAATATGGTTTGCTGAAATTCCGTGTACAAGTTGAAGTACGTTGGCTGCAAAAACTGGCCGCGCACGCAGCGATCAAGGAAGTTCCTGCTTTTGCTGCCGACGCAATCGGTTACCTTGATGCAAT  >  DH1/2687666‑2687961

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