Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,424,4860TC59.4% 357.4 / inf 531intergenic (‑1/‑350)paaA/paaZphenylacetate‑CoA oxygenase subunit PaaA/phenylacetic acid degradation bifunctional protein PaaZ
Reads supporting (aligned to +/- strand):  ref base T (109/106);  new base C (155/159);  total (264/267)
Fisher's exact test for biased strand distribution p-value = 7.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01

Only 100 of 535 total aligned reads displayed.

TTCATCCTGCACTTTGGCCAACAGAATGGCTTTACGCCGCAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGAAAGTTACACAATAGTTAAGT  >  DH1/2424336‑2424636
                                                                                                                                                      |                                                                                                                                                      
ttCATCCTGCACTTTGGCCAACAGAATGGCTTTACGCCGCAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCAtt                                                                                                                                                        >  5:1075439/1‑151 (MQ=255)
   aTCCTGCACTTTGGGCAACAGAATGGCTTTACGCCGCAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAg                                                                                                                                                     <  8:1435615/151‑1 (MQ=255)
        gCACTTTGGCCAAGAGAATGGCTTTACGCCGCAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAAGATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCt                                                                                                                                                 <  6:1844736/150‑1 (MQ=255)
        gCACTTTGGCCAACAGAATGGCTTTACGCCGCAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCt                                                                                                                                                 <  2:1042054/150‑1 (MQ=255)
             ttGGCCAACAGAATGGCTTTACGCCGCAAGGTTGGTTCGCGAGTGTTCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGc                                                                                                                                            <  3:1346846/150‑1 (MQ=255)
              tGGCCAACAGAATGGCTTTACGCCGCAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCt                                                                                                                                           <  4:1941049/150‑1 (MQ=255)
              tGGCCAACAGAATGGCTTTACGCCGCAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTc                                                                                                                                          >  5:1210400/1‑151 (MQ=255)
               gggcAACAGAATGGCTTTACGCCGCAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCg                                                                                                                                         >  4:320400/4‑151 (MQ=255)
               ggCCAACAGAATGGCTTTACGCCGCAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCg                                                                                                                                         <  8:1306633/151‑1 (MQ=255)
                  cAACAGAATGGCTTTACGCCGCAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATa                                                                                                                                      <  8:1590878/151‑1 (MQ=255)
                  cAACAGAATGGCTTTACGCCGCAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATa                                                                                                                                      <  8:1692670/151‑1 (MQ=255)
                    aCAGAATGGCTTTACGCCGCAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATATCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTTGATAGc                                                                                                                                    >  7:1950192/1‑151 (MQ=255)
                     cAGAATGGCTTTACGCCGCAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGcc                                                                                                                                   >  6:700089/1‑151 (MQ=255)
                          tGGCTTTACGCCGCAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGCTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTc                                                                                                                              >  2:504683/1‑151 (MQ=255)
                           ggCTTTACGCCGCAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTc                                                                                                                              <  6:2073005/150‑1 (MQ=255)
                               ttACGCCGGAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTg                                                                                                                           <  1:1984664/149‑1 (MQ=255)
                               ttACGCCGCAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTggg                                                                                                                         >  4:2027242/1‑151 (MQ=255)
                               ttACGCCGCAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTgg                                                                                                                          <  1:1795033/150‑1 (MQ=255)
                               ttACGCCGCAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTgg                                                                                                                          <  8:562883/150‑1 (MQ=255)
                                     cgcAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATcc                                                                                                                   >  2:785099/1‑151 (MQ=255)
                                     cgcAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATcc                                                                                                                   >  6:488278/1‑151 (MQ=255)
                                     cgcAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATcc                                                                                                                   >  1:330685/1‑151 (MQ=255)
                                      gcAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCg                                                                                                                  >  4:1732454/1‑151 (MQ=255)
                                      gcAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCg                                                                                                                  <  5:2030348/151‑1 (MQ=255)
                                      gcAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCg                                                                                                                  <  1:39153/151‑1 (MQ=255)
                                        aaGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCATTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGtt                                                                                                                >  2:1513302/1‑151 (MQ=255)
                                         aGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTg                                                                                                               <  1:1007471/151‑1 (MQ=255)
                                              ggTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCa                                                                                                           <  7:663650/150‑1 (MQ=255)
                                               gTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCaaa                                                                                                         <  7:1060774/151‑1 (MQ=255)
                                               gTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCCCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCaaa                                                                                                         <  6:2202957/151‑1 (MQ=255)
                                                 gcgcGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGc                                                                                                       <  3:1680807/151‑1 (MQ=255)
                                                    cgAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGtt                                                                                                    <  5:459648/151‑1 (MQ=255)
                                                       gTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGttctt                                                                                                 >  7:1895731/1‑151 (MQ=255)
                                                        tGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGATTCCTGGGCGATCCGTTGCTCAAAGCGttcttc                                                                                                >  6:145965/1‑151 (MQ=255)
                                                          aTCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTcttctt                                                                                              >  5:1354990/1‑151 (MQ=255)
                                                             cAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTggg                                                                                           >  4:551261/1‑151 (MQ=255)
                                                              aGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGt                                                                                          >  7:1453492/1‑151 (MQ=255)
                                                                 tACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCAc                                                                                       >  3:593781/1‑151 (MQ=255)
                                                                  aCCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACg                                                                                      >  3:2060018/1‑151 (MQ=255)
                                                                  aCCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACg                                                                                      >  4:1821620/1‑151 (MQ=255)
                                                                    cTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGtc                                                                                    >  1:1626863/1‑151 (MQ=255)
                                                                      tCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGtc                                                                                    <  1:901599/149‑1 (MQ=255)
                                                                       cAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCTGATCAATGTCTTTCTGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCAAGtcat                                                                                  <  4:72972/150‑1 (MQ=255)
                                                                        aGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGtcatca                                                                                <  5:1798303/151‑1 (MQ=255)
                                                                        aGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGtcatc                                                                                 >  2:301649/1‑150 (MQ=255)
                                                                                     acaaTCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTcc                                                                   >  6:1167727/1‑151 (MQ=255)
                                                                                      caaTCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCa                                                                  <  2:1372624/151‑1 (MQ=255)
                                                                                       aaTCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTcc                                                                   >  5:938051/1‑149 (MQ=255)
                                                                                          ctcGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGttt                                                              >  2:637448/1‑151 (MQ=255)
                                                                                            cGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGtttt                                                             <  5:1690263/150‑1 (MQ=255)
                                                                                            cGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGtttt                                                             <  5:220350/150‑1 (MQ=255)
                                                                                            cGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGtttt                                                             <  7:1574375/150‑1 (MQ=255)
                                                                                              gAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGt                                                          <  3:702947/151‑1 (MQ=255)
                                                                                                  gcgcATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGaat                                                      >  6:829668/1‑151 (MQ=255)
                                                                                                    gcATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGaatta                                                    <  8:1081211/151‑1 (MQ=255)
                                                                                                      aTGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGaattaat                                                  <  6:364026/151‑1 (MQ=255)
                                                                                                      aTGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGaattaat                                                  >  1:1425919/1‑151 (MQ=255)
                                                                                                       tgctgcCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGaattaatt                                                 >  6:1371536/1‑151 (MQ=255)
                                                                                                        gctgcCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAattaatta                                                <  4:1624046/151‑1 (MQ=255)
                                                                                                        gctgcCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAattaatta                                                <  7:2142384/151‑1 (MQ=255)
                                                                                                         ctgcCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAattaatta                                                >  3:1769505/1‑150 (MQ=255)
                                                                                                             ccAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTaaaaa                                            >  7:1266189/1‑150 (MQ=255)
                                                                                                               aaTCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAacac                                         >  8:240193/1‑151 (MQ=255)
                                                                                                               aaTCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAacac                                         >  4:974310/1‑151 (MQ=255)
                                                                                                                   tGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAAt                                     >  5:1561882/1‑151 (MQ=255)
                                                                                                                    ggCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGaa                                      <  1:2221418/149‑1 (MQ=255)
                                                                                                                    ggCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGaa                                      <  5:610384/149‑1 (MQ=255)
                                                                                                                    ggCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGTGTCACGTCATCACCTTTACGATTCCAGTTTTCGTTAATTAATTAAAAACACGaa                                      <  6:1276913/149‑1 (MQ=255)
                                                                                                                    ggCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGa                                       <  2:1787554/148‑1 (MQ=255)
                                                                                                                     gCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATca                                   <  8:2117269/151‑1 (MQ=255)
                                                                                                                     gCGGATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATc                                    >  4:1289954/1‑150 (MQ=255)
                                                                                                                        gATCAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACAt                                <  7:857953/151‑1 (MQ=255)
                                                                                                                         aTCAATGTCTTTCGGTAAGCATCTGTCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTTCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTACAGTTTTAGTGAATTAATTAAAAACACGAATCACAtt                               <  8:2014492/151‑1 (MQ=255)
                                                                                                                          tCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTg                              >  5:1446860/1‑151 (MQ=255)
                                                                                                                           cAATGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTg                              <  3:1483056/150‑1 (MQ=255)
                                                                                                                            aaTGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGtt                            >  1:590818/1‑151 (MQ=255)
                                                                                                                              tGTCTTTCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGtttt                          >  3:372893/1‑151 (MQ=255)
                                                                                                                                 cTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTAtt                       >  2:450028/1‑150 (MQ=255)
                                                                                                                                 cTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATg                       >  1:1192670/1‑151 (MQ=255)
                                                                                                                                   ttCGTTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATg                       <  2:165619/149‑1 (MQ=255)
                                                                                                                                   ttCGGTAAGCATCGGGCTTCCAGTCGTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTGAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTAt                        <  4:838269/148‑1 (MQ=255)
                                                                                                                                    tCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGaa                     <  3:1207554/150‑1 (MQ=255)
                                                                                                                                    tCGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGaa                     <  6:1615883/150‑1 (MQ=255)
                                                                                                                                     cGGTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGaaa                    <  4:2142014/150‑1 (MQ=255)
                                                                                                                                       gTAAGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGAAAGtt                 <  3:47713/151‑1 (MQ=255)
                                                                                                                                        tAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGAAAGtt                 >  5:763362/1‑150 (MQ=255)
                                                                                                                                         cagcagcGGGCATCCAGTCCTGGGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGAAAGTTa                <  4:1493725/149‑1 (MQ=255)
                                                                                                                                         aaGCATCGGGCATTCAGTCCTGTGGCTTGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGAAAGTTa                <  8:1948409/150‑1 (MQ=255)
                                                                                                                                         aaGCATCGGGCATGCAGTCCTGTGGCGCGATAGCCGTTTCCTGGGCTATCCGGTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGAAAGTTac               <  2:1944053/151‑1 (MQ=255)
                                                                                                                                         aaGCATCGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGAAAGTTa                >  6:1746159/1‑150 (MQ=255)
                                                                                                                                          aGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGAAAGTTaca              <  8:781657/151‑1 (MQ=255)
                                                                                                                                               cGGGCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGAAAGTTACACAATa         <  2:9956/151‑1 (MQ=255)
                                                                                                                                                gggCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGAAAGTTACACAATAg        <  3:225435/151‑1 (MQ=255)
                                                                                                                                                gggCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGAAAGTTACACAAt          >  2:2151198/1‑149 (MQ=255)
                                                                                                                                                gggCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGAAAGTTACACAATAg        <  3:974808/151‑1 (MQ=255)
                                                                                                                                                 ggCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGAAAGTTACACAATAGt       <  7:369655/151‑1 (MQ=255)
                                                                                                                                                 ggCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGAAAGTTACACAATAGt       <  3:1893083/151‑1 (MQ=255)
                                                                                                                                                  gCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGAAAGTTACACAATAGtt      <  4:1154924/151‑1 (MQ=255)
                                                                                                                                                  gCATCCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGAAAGTTACACAATAGtt      >  1:488471/1‑151 (MQ=255)
                                                                                                                                                      tCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGAAAGTTACACAATAGTTAAg   <  2:1094851/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
TTCATCCTGCACTTTGGCCAACAGAATGGCTTTACGCCGCAAGGTTGGTGCGCGAGTGATCCAGTTACCTTCAGGCAACATACCAACAATCTCGGAGTGCGCATGCTGCCCAATCTGGCGGATCAATGTCTTTCGGTAAGCATCGGGCATTCAGTCCTGTGGCTCGATAGCCGTTTCCTGGGCTATCCGTTGCTCAAAGCGTTCTTCTTGGGTCACGTCATCACCTTTACGATTCCAGTTTTCGTGAATTAATTAAAAACACGAATCACATTGTTTTATGAAAGTTACACAATAGTTAAGT  >  DH1/2424336‑2424636

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: ‑