Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,965,8930AC25.2% 232.9 / 25.7 171intergenic (+126/‑369)aqpZ/ybjEaquaporin/transporter
Reads supporting (aligned to +/- strand):  ref base A (7/97);  new base C (21/14);  total (41/130)
Fisher's exact test for biased strand distribution p-value = 2.88e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 183 total aligned reads displayed.

ACCCTGCTGGAAAAGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGT  >  DH1/2965744‑2966043
                                                                                                                                                      |                                                                                                                                                        
aCCCTGCTGGAAAAGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAg                                                                                                                                                         <  2:1486322/151‑1 (MQ=255)
  ccTGCTGGAAAAGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAggg                                                                                                                                                       <  8:246249/151‑1 (MQ=255)
   ctgctgGAAAAGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGCgggg                                                                                                                                                      >  6:1432081/1‑151 (MQ=255)
   ctgctgGAAAAGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGCgggg                                                                                                                                                      >  6:603986/1‑151 (MQ=255)
     gctgGAAAAGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCggggggggg                                                                                                                                                    >  3:1054440/1‑150 (MQ=255)
     gctgGAAAAGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCCcggggg                                                                                                                                                    >  6:1668019/1‑150 (MQ=255)
     gctgGAAAAGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAgggg                                                                                                                                                      <  1:1683470/149‑1 (MQ=255)
       tgGAAAAGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGCgggggggg                                                                                                                                                  >  5:730599/1‑151 (MQ=255)
       gggAAAAGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGggcgg                                                                                                                                                  <  5:329242/150‑1 (MQ=255)
         gAAAAGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGGGGGggcgggc                                                                                                                                                >  5:1339806/1‑149 (MQ=255)
          aaaaGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGggggggggggggc                                                                                                                                               >  3:705367/1‑148 (MQ=255)
          aaaaGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCACGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCa                                                                                                                                              <  2:1186637/151‑1 (MQ=255)
             aGCGTGATTAATTAACGCAGTAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATc                                                                                                                                           <  1:162573/151‑1 (MQ=255)
             aGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATc                                                                                                                                           <  4:778816/151‑1 (MQ=255)
              gCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGgggggggggggggcgcc                                                                                                                                           >  7:188987/1‑144 (MQ=255)
               cGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGGGGGGGGGGccccaccc                                                                                                                                         >  3:569458/1‑145 (MQ=255)
               cGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgc                                                                                                                                         <  4:1042348/151‑1 (MQ=255)
                gTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgcc                                                                                                                                        <  4:589709/151‑1 (MQ=255)
                  gATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGCCGGGGGGGGCCaaaaaccc                                                                                                                                       >  4:984780/1‑145 (MQ=255)
                  gATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccgc                                                                                                                                      <  2:1416840/151‑1 (MQ=255)
                   attaattaACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGGTGTTATATGTTGATTCATCATCTTGTTGCGCGCCGTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccgcc                                                                                                                                     <  7:548549/151‑1 (MQ=255)
                    ttaattaaCGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGggggggggggggccccccccccc                                                                                                                                     >  1:226701/1‑138 (MQ=255)
                     taataaaCGCAGAAGCCCTGGCGTAATGTCAGGCCTTTTTTTAAGATCTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccgccg                                                                                                                                    <  6:489901/150‑1 (MQ=255)
                      aattaaCGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccgccgc                                                                                                                                   <  1:1051449/150‑1 (MQ=255)
                       attaaCGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGCGGGGGCGGGCCAACGCcccccccc                                                                                                                                 >  4:1512462/1‑151 (MQ=255)
                          aaCGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCCCGGGGGGGgccccccccgccccccccc                                                                                                                              >  2:1327268/1‑150 (MQ=255)
                          aaCGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCt                                                                                                                              <  1:441316/151‑1 (MQ=255)
                           aCGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTg                                                                                                                             <  2:731959/151‑1 (MQ=255)
                            cGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGc                                                                                                                            <  1:731647/151‑1 (MQ=255)
                             gCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCANNTTGTTGCGCTCCTNTCATGAGCTNAGCAACATATNGCCACNGGCNCAAGGAGCGCNCAGGGGGCGGNCAATCGCCGCCGCCCCCTGCa                                                                                                                           <  7:1116706/151‑1 (MQ=255)
                              cAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGTCCAATCGCAGCCGCCCCCTGCa                                                                                                                           <  2:397500/150‑1 (MQ=255)
                                gAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGGTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAccc                                                                                                                        <  7:746284/151‑1 (MQ=255)
                                  aGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGCCGGGGGGGGCCAAACcccccccccccccccccccc                                                                                                                       >  1:92820/1‑143 (MQ=255)
                                    gCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGGGGGGGCGGCGCCAACGccccccccccccccccccccc                                                                                                                     >  1:1022499/1‑150 (MQ=255)
                                      gtGGCGTAATGTCAGGCCTTTTTTAAAGATTTGTTTATTTTTTGATTATTCATGTTTTGGCGCTCATTTCTTCCTCTACCAAAAATATCGCCAGTGGCGAAAGGCGCGCGCCGGGGGGGGCCATGCGCCGCCGCCCCCTCCACCCCCGGGc                                                                                                                  <  8:176308/150‑1 (MQ=255)
                                      cTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGGCGGGGGGGGCCAAAccccccccccccccccccccccccc                                                                                                                   >  8:753254/1‑148 (MQ=255)
                                       gggCGTAATTTCTGGCCTTTTTTAAGTAAGTTGTTAAATGTTGATTCATCATCTTGTTGCCCGCATTCCATGAGCTAAGCAACATATGCCCACTGCCGCATGGACCGCGCTGGGGGCGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGct                                                                                                                 <  4:571382/150‑1 (MQ=21)
                                           gTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGGGGGGGGGGGCGGCCGccccccccccccccccccccccccccccc                                                                                                             >  5:1470772/1‑142 (MQ=255)
                                             attGTCAGGCCTTTTTTTAAGATGTTTTTATATGTTGATTCATCATCTTGTTGCGCTCTTTTCATGAGATACGCAACATATTGCCAGTGGCGCATGGAGCGCGCAGGGGGGGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGc                                                                                                           <  4:457187/148‑1 (MQ=25)
                                             aaTGTGAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGc                                                                                                           <  4:699237/150‑1 (MQ=255)
                                             aaTGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCCAAcccccccccccccccccccccccccgcaggcg                                                                                                            >  2:1010275/1‑141 (MQ=255)
                                                gTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGCGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCACGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAAc                                                                                                       <  8:1299039/151‑1 (MQ=255)
                                                  cAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCGCCCCCTGCACCCCCGGGCTCTGGCGAACaaa                                                                                                    <  2:8664/151‑1 (MQ=255)
                                                    ggCTTTTTTTTATGATGTTGTTATATGTTGATTCTTCTTGTTGTTGGTCTCCTTTCTTGAGATACGAAACTTATTGCCACTGGCGCATGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACaaa                                                                                                    <  4:121258/150‑1 (MQ=18)
                                                    ggCTTTTTTTTATGATGTTGATATATGTTGATTCATCAGCTTGTTGGCCCCCTTTCATGAGATAAGCAACATTTTGCCACGGGCGCAAGGACCGCGCCGGGGGCGGCCAATGGCGGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACaaaa                                                                                                   <  3:411491/151‑1 (MQ=21)
                                                      gttttttttAAGTTGTTGTTATTTGTTGTTTCATCATCTTTTTGCGCTCTTTTCATGAGCTAAGAAACATATTGCCACTGGCGCATGGAGCGGGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATc                                                                                                 <  2:510039/149‑1 (MQ=21)
                                                        tttttttAAGATGTTTTTTTTTGTTCTTACATCATCTTGGCGCGCGCCTTTCTTGCGCAACCAAACATATTCCCACTGGCGCACGGAGGGCGCTGGGGGCGGCCAATCCCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcg                                                                                                <  8:321131/150‑1 (MQ=12)
                                                             ttATGATGTTTTTATTTTTTTATTCTAGATGTTTTTGCGGGCCTTCCCTGGGCTACGAAACATTTTGCCCCTGGCGCAGGGAGGGCCGCGGGGGGGGCAAACCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCt                                                                                           <  5:1148072/150‑1 (MQ=9)
                                                             ttAAGATTTTGTTATAGGTTGATGCAGCTTCTTGTTGCGCGCCTTTCATCAGGTACGCAAAATATTGCCACTGGGGCACGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAACTCGCCGCt                                                                                           <  5:1647324/150‑1 (MQ=18)
                                                              tAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTTATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTgcg                                                                                        >  4:802478/1‑151 (MQ=255)
                                                               aaGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCCGGGGGGGGCCAAAACACCCCGCCCCCCGCCCCCCCCGGCGCAGGGCGAaaaaaaaaacccagac                                                                                         >  6:257358/1‑141 (MQ=255)
                                                               aaGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCGTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTgcg                                                                                        <  2:460038/151‑1 (MQ=255)
                                                                aGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAATCGCCGCCCCCCCCCGCCCCCCCCGGCGCTGGCGCCCAAAAACGCCGCGgcgc                                                                                       >  1:1441067/1‑151 (MQ=255)
                                                                  aTGTTTTTATATGTTGATTCATCATCTTGTGGCGCTCCTTTCATGAGCTAACAAAAATATTGCCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGcgcg                                                                                      <  7:426779/150‑1 (MQ=21)
                                                                     tttttATTTGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCCCTGGCGCAAGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTg                                                                                   <  6:860325/147‑1 (MQ=255)
                                                                          atatGTTTATTCATCATCTTGTTGCGCTCCTTTCATGAGCTACGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCGTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTc                                                                              <  4:1021691/150‑1 (MQ=255)
                                                                          atatGTTGATTCATCATGTTGTTGCGCTCCTTTCATGAGGAAAGCAACATATTGCAACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCg                                                                             <  1:80115/151‑1 (MQ=255)
                                                                          atatGTTGATGCAGCATCTTGTTGCGCTCCTTTCATCAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCg                                                                             <  2:897846/151‑1 (MQ=255)
                                                                          aactGTTGATTGATCTTCTTGTTGCGCTCCTTTCATCAGCTACGCAACTTATTGCCTCTGGCGCAAGGTGCGCGGTGGGGGCGGCAACGCGCGCCCGCCCCGGGCAGCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTc                                                                              <  3:17504/147‑1 (MQ=18)
                                                                              gttcttacatcatGTTGTTGCGCGCCTTGCATGCGCTACCCAACTTATGCCCACTGGCGCAGGGAGCGCGCAGGGGGCGGCCAAGCCCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACtt                                                                         <  5:29191/144‑1 (MQ=14)
                                                                              gTTCATTCATCATCTTGTTGCGCTCCTTTCATGACCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGCCAAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACt                                                                          <  2:1454094/150‑1 (MQ=255)
                                                                                tcatgcttcatcTTGTGGCGCGCCTTCCATGCGCTAAGAAACATATTGCCACTGGCGCACGGCGCGCGCCGGGGGGGGCCAATCGCCGCGGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTa                                                                        <  3:1073868/143‑1 (MQ=17)
                                                                                     catcatCTTGTTGGGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGGGCAAGGAGCGGGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCtt                                                                  <  6:132640/151‑1 (MQ=255)
                                                                                     catcatCTTGTTGCGCTCCTTTCATCAGCTAAGCAACATATTGCCACTGGCGCAAGGCGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCtt                                                                  <  1:854762/151‑1 (MQ=255)
                                                                                       acatcTTGTGGCGCGCCTTCCATGAGCTAACCAAAATTTGGCCACGGCCGCACGCAGCGCGCCGGGGGCGGCCAAGCCCCGCCGCCCCGGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTAc                                                                <  7:422633/150‑1 (MQ=14)
                                                                                        catcTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACg                                                               <  8:826676/151‑1 (MQ=255)
                                                                                        catcTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACg                                                               <  5:1664483/151‑1 (MQ=255)
                                                                                          tgttgttgCGCTCTTTTCATGAGCTAACAAAAATTTTGCCCCGGGCGCATGGCGGGCGCTGGGGGGGGCCAAGCGCCGCCGCCCCGTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACgg                                                              <  5:937167/148‑1 (MQ=17)
                                                                                              gttgCGCTCCTTTAATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTTCCCTCGACTTATCCCTTACGGCTAc                                                          <  2:92530/150‑1 (MQ=255)
                                                                                                 gggcgccTTTCATGTGGTAAGCAACATATTGCCCCGGGCGCATGGCGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCGGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCgg                                                       <  8:802286/145‑1 (MQ=14)
                                                                                                 gcgggccTTTCATGAGCTAAGCAACATATTCCCACTGGCGCAAGGCGCGGGCAGGGGGCGCCCAATCGCGGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCgg                                                       <  3:1253780/145‑1 (MQ=17)
                                                                                                 gcgcTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCggg                                                      <  8:1277891/151‑1 (MQ=255)
                                                                                                     gccTTTCATGAGCTAAGCAAAATATTGCCACTGGCGCACGGCGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCgg                                                  <  1:399731/150‑1 (MQ=255)
                                                                                                      ccTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCACTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCg                                                   <  5:1538061/149‑1 (MQ=255)
                                                                                                      ccTTGCATGAGCTACGCAACATATTGCCCCTGGCGCATGGAGCGCGCAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCggg                                                 <  6:1290992/151‑1 (MQ=255)
                                                                                                       gtttCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGgc                                                <  8:1521307/150‑1 (MQ=255)
                                                                                                         ttCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGgc                                                <  8:1551911/149‑1 (MQ=255)
                                                                                                         ttCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGgc                                                <  2:220878/149‑1 (MQ=255)
                                                                                                            aTGAGCTACGCAACATATTGCCACTGGCGCAAGCAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCGTGCACCCCCGGGCGCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                           <  5:1036912/151‑1 (MQ=255)
                                                                                                               aGCTAAGCAACATATTGCCACTGGCGCAAGGCGCGCGCGGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                        <  1:1517398/151‑1 (MQ=255)
                                                                                                               aGCTAACCAACATATTGCCACTGGCGCAAGGCGCGCGCAGGGGGCGGCCAAGCCCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                        <  1:1339553/151‑1 (MQ=255)
                                                                                                                         aatatgcccaggggCCCACGGCGCGCGCGGGGGGGGGCCAAGCCCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGt                              <  7:1418811/139‑1 (MQ=11)
                                                                                                                              tGCCCCTGGCGCAAGGAGCGCGCAGGGGGCGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaaa                          <  4:504669/150‑1 (MQ=18)
                                                                                                                               gCCACTGGCGCATGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                         <  3:956369/150‑1 (MQ=18)
                                                                                                                               gCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCATTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                         <  2:1632796/150‑1 (MQ=18)
                                                                                                                               gCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCGTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                         <  3:1609368/150‑1 (MQ=18)
                                                                                                                               gCCACGGGCGCACGGCGCGCGCCGGGGGCGGCCAAGCCCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                         <  6:582654/150‑1 (MQ=11)
                                                                                                                               cccACTGCCCCACGGCGCGCGCAGGGGGCGCCCACCCCCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcg                        <  5:7822/150‑1 (MQ=11)
                                                                                                                                     ggCGCACGGAGCGCGCTGGGGGCGGCCAATCGCCGCGGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc                   <  8:1003840/150‑1 (MQ=9)
                                                                                                                                     ggCGCAAGGAGCGCGCAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCt                  <  5:849442/151‑1 (MQ=14)
                                                                                                                                     ggCGCAAGGAGCGCGCAAGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCt                  <  8:1183500/151‑1 (MQ=14)
                                                                                                                                       cgcAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCa                <  6:1093708/151‑1 (MQ=31)
                                                                                                                                          atggcgcgggcAGGGGGGGCCCAATGGCCGCGGCCCCGTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCGCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc             <  5:698197/146‑1 (MQ=11)
                                                                                                                                          aaGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCTCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc             <  2:1439962/151‑1 (MQ=14)
                                                                                                                                           cggAGCGCGCCGGGGGCGGCCAAGCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc            <  8:133526/150‑1 (MQ=9)
                                                                                                                                                cgcgcAGGGGGCGGCCAATCGCGGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATc       <  5:546468/151‑1 (MQ=14)
                                                                                                                                                 gcgcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc      >  2:752423/1‑151 (MQ=31)
                                                                                                                                                 gcgcAGGGGGCGGCCAATCGCCGCCCCCCCCAGCCCCCCCCGGCTCTGGCGAACAAAATCCCCGCTGCGCGGTGCCCTCGACTTATACCATACCGCTACCCGGTCCGGCCCGAGGGAACATCCCTGTAAAACGCGCCCACAGCCcaccacc      >  8:963040/1‑147 (MQ=11)
                                                                                                                                                    gcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCccatccat    >  5:751601/1‑151 (MQ=14)
                                                                                                                                                      aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgt  >  1:509472/1‑151 (MQ=31)
                                                                                                                                                      |                                                                                                                                                        
ACCCTGCTGGAAAAGCGTGATTAATTAACGCAGAAGGCCTGGCGTAATGTCAGGCCTTTTTTTAAGATGTTGTTATATGTTGATTCATCATCTTGTTGCGCTCCTTTCATGAGCTAAGCAACATATTGCCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGACTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGT  >  DH1/2965744‑2966043

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: