Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,965,8990CG25.5% 359.5 / 69.2 204intergenic (+132/‑363)aqpZ/ybjEaquaporin/transporter
Reads supporting (aligned to +/- strand):  ref base C (23/129);  new base G (29/23);  total (52/152)
Fisher's exact test for biased strand distribution p-value = 4.70e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 243 total aligned reads displayed.

GCTGGAAAAGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTG  >  DH1/2965749‑2966049
                                                                                                                                                        |                                                                                                                                                       
gctgGAAAAGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGggc                                                                                                                                                         <  3:1942400/151‑1 (MQ=255)
  tgGAAAAGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGCCGGGggcg                                                                                                                                                        >  2:1719751/1‑151 (MQ=255)
  tgGAAAAGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGGgcgggcgg                                                                                                                                                       >  8:706669/1‑151 (MQ=255)
    gAAAAGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGCgggggggggc                                                                                                                                                     >  4:1861578/1‑149 (MQ=255)
      aaaGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGgggggggggggggg                                                                                                                                                   >  5:1735586/1‑147 (MQ=255)
      aaaGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGggggggggggggcg                                                                                                                                                   >  2:608163/1‑147 (MQ=255)
       aaGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGCGGGGGGGGGccc                                                                                                                                                    >  6:655715/1‑149 (MQ=255)
       aaGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGggggggggggggccc                                                                                                                                                  >  4:962151/1‑146 (MQ=255)
       aaGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCCCGGGGGGGcccac                                                                                                                                                  >  3:145823/1‑148 (MQ=255)
        aGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATc                                                                                                                                                <  6:1447924/151‑1 (MQ=255)
        aGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATc                                                                                                                                                <  2:1542679/151‑1 (MQ=255)
         gCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGgggggggggggcccgac                                                                                                                                                >  1:1705172/1‑144 (MQ=255)
         gCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcg                                                                                                                                               <  7:1634186/151‑1 (MQ=255)
          cGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGAAGGGGCGGCCaaaagc                                                                                                                                              >  6:1290587/1‑147 (MQ=255)
          cGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgc                                                                                                                                              <  7:1295421/151‑1 (MQ=255)
           gTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGgggggggggggggcccccc                                                                                                                                              >  7:607566/1‑142 (MQ=255)
            tGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGCGGGGGGGGGcccaagccc                                                                                                                                              >  5:1928330/1‑144 (MQ=255)
            tGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGgggggggggggggggggggg                                                                                                                                             >  7:1687564/1‑141 (MQ=255)
              attaattaACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccgcc                                                                                                                                          <  2:1359260/151‑1 (MQ=255)
               ttaattaaCGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGggggggggggggccgcccccccc                                                                                                                                          >  1:597780/1‑138 (MQ=255)
               ttaattaaCGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGgggggggggggcgctccaccccc                                                                                                                                          >  7:232144/1‑138 (MQ=255)
               ttaattaaCGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGgggggggggggcccccccccccc                                                                                                                                          >  3:1831104/1‑138 (MQ=255)
                taattaaCGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGgggggggggggcgcccgcccccgc                                                                                                                                         >  7:1039945/1‑137 (MQ=255)
                 aattaaCGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGCCGGGGGGGGcccaaccccccccc                                                                                                                                         >  3:47337/1‑139 (MQ=255)
                 aattaaCGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGgggggggggggcccgccgccccccc                                                                                                                                        >  3:650635/1‑136 (MQ=255)
                    taaCGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                                     <  6:1764792/150‑1 (MQ=255)
                     aaCGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCGCGGGGGGGCCCAGCCCCCCCGcccccc                                                                                                                                   >  4:1021553/1‑150 (MQ=255)
                     aaCGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCt                                                                                                                                   <  4:1692347/151‑1 (MQ=255)
                        ccAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCa                                                                                                                                <  3:1778260/150‑1 (MQ=255)
                          aGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCATAGCGccccccccccgcccc                                                                                                                                >  2:203075/1‑146 (MQ=255)
                          aGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGCCGGGGGGGCCCCAACCCCGccccccccccacc                                                                                                                               >  6:649201/1‑146 (MQ=255)
                             aGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATCAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAccccc                                                                                                                           <  8:610260/151‑1 (MQ=255)
                               gCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGCGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCgg                                                                                                                         <  6:1671966/151‑1 (MQ=255)
                                 cTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCCCGGGGGGGCCACAACCCCCCCGCCCCCCCCAcccccccc                                                                                                                        >  6:845059/1‑148 (MQ=255)
                                   ggCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGGCGGGGGGGGCCAAAcccccccccccccccccccccccgccc                                                                                                                      >  2:1232706/1‑146 (MQ=255)
                                    gcttattgtacggCTTTTTTTTACGATTTTTTTATCTGTTGTTTCAGAATGTTGTTGCGCTCCTTTCATTCCCAACGCAAAATATTGCCACTGGCGCAGGGCGGGCGCCGGGGGGGGCCAACCCCCCCCGCCCCCTGCACCCCCGGGctct                                                                                                                    <  5:38133/140‑1 (MQ=12)
                                     cGTATTGTCAGCCCTTTTTTTAAGATGTTGTTATATGTTCATTCATCATCTTGTTGCGCTCCTTTCATCACCTAAGCAACATATTGCCCCTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCCCCGCCGCCCCCTGCACCCCCGGGCTCTg                                                                                                                   <  6:1892803/151‑1 (MQ=25)
                                     cGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGGGGGGGGCGGCAAAAGACGCCGCCCCCCCGAcccccccgcgc                                                                                                                     >  6:1950490/1‑143 (MQ=255)
                                        attGTCTGGCCTTTTTTTAAGATGTTTTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCGGCCGCCCCCTGCACCCCCGGGCTCTGGc                                                                                                                 <  2:670193/148‑1 (MQ=255)
                                         aTGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGCGGGGGGGGGGCCAACCCCCCCCCCCCCCCCCCCCCCCgcggcggggg                                                                                                                 >  7:1263787/1‑142 (MQ=255)
                                           gTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAAc                                                                                                             <  4:530369/151‑1 (MQ=255)
                                            tCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCCCGGGGCGGCCCAACGCCGCCCCCCCCAGCccccccccgggcaggggcac                                                                                                             >  2:44994/1‑136 (MQ=37)
                                              aGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGGCGGGGGCGCCCAAACGCAccccccccccccccccccccccccctccccaaa                                                                                                            >  4:195110/1‑135 (MQ=255)
                                                 ccTTTTTTTAATATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAACCAACATATTGCCACTGGCGCAAGGCGCGCGCAGGGGGCGCCCAACCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATc                                                                                                       <  5:1203498/151‑1 (MQ=255)
                                                 ccTTTTTTTAAGTTGTTTTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGCGCTAACCAAAATATTGCCCCTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCGGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATc                                                                                                       <  8:896818/151‑1 (MQ=21)
                                                  cTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGCGGGGGGGGCCCAAACCCGCCCCCCCCCCCCCCCCCCGGCGCAGGCGCCCaaaaaa                                                                                                       >  7:1643945/1‑148 (MQ=255)
                                                    ttttttAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgc                                                                                                     <  8:1955886/150‑1 (MQ=255)
                                                    ttttttAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgc                                                                                                     <  1:375887/150‑1 (MQ=255)
                                                      tttaacgATGTTTTTATTTGTTCATTCATCATCTTGTGGCGCCCCTTGCTGGCGCTAACCAACATATTCCCCCGGGCGCACGGCGCGCGCCGGGGGCGGCCAACCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgccg                                                                                                   <  8:1614914/144‑1 (MQ=17)
                                                         tAAGTTTTTGTTTTTTGTTGATTTTTCATTTTGTTGCGGTCTTTTCATGAGCTATGCAACATATTGCCACTGGCGCATGGAGCGGGCAGGGGGCGGCCTATGGCGGCGGCCCCGTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTg                                                                                                <  6:301508/150‑1 (MQ=17)
                                                             aTGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGGGGGGGGGGGCCAGAGGCGGCCCCCCCGGCCCCCCCCGGCGCGGGGGGACCAAAAcgcccgcgggcgg                                                                                           >  6:1959941/1‑142 (MQ=255)
                                                               gttgttTTTTGTTGATCCATCATCTTGTTGCGGTCCTTTCTTGAGCTAATAAAAATTTTGCCAGTGGCGCAAGCCGCGCGCTGGGGGGGGCAAACCCCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGt                                                                                          <  6:544377/150‑1 (MQ=17)
                                                                ttgttATATGTTTATTCATCATGTTGTTGCGCTCTTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTg                                                                                         <  2:511069/150‑1 (MQ=255)
                                                                 tgttATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCCCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGc                                                                                        <  1:230969/150‑1 (MQ=255)
                                                                            taTTCATCAGCTTGTTGCGCTCCTTCCATGAGCTAACCAACATATTGCCACTGGCGCACGGAGCGCGCTGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTAt                                                                             <  7:34284/149‑1 (MQ=18)
                                                                               tcttcatgTTGTTGCGCTCGTTTCATGTGCTAAGCAACATATTGCCACTGGCGCAAGCCGCGCGCAGGGGGCGGCGAATCGCCGCGGCCCCGTGCACCCCCGGGCTCTGGGGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATccc                                                                          <  6:1173931/147‑1 (MQ=17)
                                                                               tcatcaTCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGCGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATccc                                                                          <  8:1515389/150‑1 (MQ=255)
                                                                               tcatcaTCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATccc                                                                          <  4:465942/150‑1 (MQ=255)
                                                                                cctcatcTTCTTGCGCTCCTTCCATGAGCTAAGCAACATATTGCCACTGCCGCACGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCt                                                                         <  4:1200693/148‑1 (MQ=18)
                                                                                catcatCTTGTTGCGCGCCTTTCATCAGCTACGCAACATATTGCCACTGGCCCACGGCGCGCGCAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCtt                                                                        <  3:886016/151‑1 (MQ=18)
                                                                                   catcTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTAc                                                                      <  3:1089290/150‑1 (MQ=255)
                                                                                   catcTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTAc                                                                      <  1:1004758/150‑1 (MQ=255)
                                                                                   cagcTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGTGGGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTAc                                                                      <  6:915585/147‑1 (MQ=255)
                                                                                     tgttgttgCGCGCCTTTAATGCGTTAAGAAACATATTGCCCCTGGCGCATGGCGCGGGCTGGGGGGGGCCAATCGCGGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCGCTTACg                                                                     <  6:1558686/147‑1 (MQ=17)
                                                                                        tgtggCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTa                                                                 <  1:407941/150‑1 (MQ=255)
                                                                                            gcgcgcCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCgg                                                             <  6:1021099/150‑1 (MQ=255)
                                                                                             ggcTCCTTTCATTATCAAACAAACATTTTTCCAGGGGGCCACGGCGCGGGCAGGGGGGGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGt                                                           <  8:1233015/150‑1 (MQ=11)
                                                                                             cgcTCCTTTCATCCGATAAGCAACATATTGCCACTGGCGCACGGAGCGCGCCGGGGGCGGCCACTCGCGGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGt                                                           <  3:1315711/151‑1 (MQ=17)
                                                                                              gcTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTc                                                          <  4:809660/151‑1 (MQ=255)
                                                                                                gccTTTCATGAGCTAAGCAAAATATTGCCACTGGCGCATGGAGCGCGCTGGGGGCGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCgg                                                        <  2:1492824/150‑1 (MQ=255)
                                                                                                 ccTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCCCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCggg                                                       <  5:1140652/151‑1 (MQ=255)
                                                                                                     tcatcaGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGgcgc                                                    <  6:1500016/150‑1 (MQ=255)
                                                                                                      cATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGcgcg                                                   <  1:108933/150‑1 (MQ=255)
                                                                                                      cATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                                  <  5:1900512/151‑1 (MQ=255)
                                                                                                      cATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                                  <  2:1140215/151‑1 (MQ=255)
                                                                                                       aTGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCATGCGCCGCCGCCCCGTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                                 <  1:863145/151‑1 (MQ=255)
                                                                                                         gcgcaaaGAAACATTTTCCCCCTGGCGCACGGCGCGCGCTGGGGGCGGCCAAGGGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                                <  1:920422/145‑1 (MQ=11)
                                                                                                         gAGCTAAGCAACATATTGCCACGGGCGCAAGGAGCGCGCAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGt                                               <  5:147479/151‑1 (MQ=255)
                                                                                                         gAGCTAAGAAAAATATTGCCACTGGCGCACGGAGCGCGCAGGGGGGGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGt                                               <  2:1882238/151‑1 (MQ=255)
                                                                                                               aGCAACATATTGCCACTGGCGCAAGGCGGGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATc                                         <  2:832255/151‑1 (MQ=255)
                                                                                                                    aatatTGCCAGGGGCGCATGGACCGCGCTGGGGGGGGCCAATCGCCGCCGCCCCGGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGt                                    <  8:349421/150‑1 (MQ=11)
                                                                                                                    aatatTCCCACTGGCGCAAGGCGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGt                                    <  1:1269263/150‑1 (MQ=255)
                                                                                                                       atTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaa                                  <  2:71853/150‑1 (MQ=35)
                                                                                                                          gCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                               <  3:1394181/150‑1 (MQ=18)
                                                                                                                           ccACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                             <  1:144536/151‑1 (MQ=32)
                                                                                                                                ggCGCAAGGCGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCt                        <  8:1005068/151‑1 (MQ=14)
                                                                                                                                 gcgaacggAGCGCGCTGGGGGGGGCCAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTc                       <  3:538350/145‑1 (MQ=255)
                                                                                                                                      aggcgggcgcCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc                  <  7:1814160/145‑1 (MQ=9)
                                                                                                                                       ggAGCGGGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca                 <  3:1849836/151‑1 (MQ=14)
                                                                                                                                           cgcgcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATGCCTGTAAAACGCGCCCTCAGCCCACATc             >  2:1692955/1‑151 (MQ=14)
                                                                                                                                            gcgcAGGGGGCGGCCAATCGCCGCCGCCCCCCGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc            >  8:532779/1‑151 (MQ=14)
                                                                                                                                            gcgcAGGGGGCGGCCAATCCCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc            <  5:1323846/151‑1 (MQ=14)
                                                                                                                                                   gggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgtg       >  8:1629076/1‑151 (MQ=31)
                                                                                                                                                    ggggCGGCCAATCGCCGCCGCCCCGTGGACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTgg      <  4:1408183/151‑1 (MQ=14)
                                                                                                                                                     gggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTggg     >  6:1867229/1‑151 (MQ=31)
                                                                                                                                                        cggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTg  <  5:1801313/151‑1 (MQ=31)
                                                                                                                                                        cggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTg  <  1:756864/151‑1 (MQ=31)
                                                                                                                                                        cggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTg  >  5:697738/1‑151 (MQ=31)
                                                                                                                                                        cggcCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTg  <  5:1871205/151‑1 (MQ=14)
                                                                                                                                                        cggcCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGGCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTg  <  5:1533616/151‑1 (MQ=14)
                                                                                                                                                        |                                                                                                                                                       
GCTGGAAAAGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTG  >  DH1/2965749‑2966049

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