Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,448,5260AC33.0% 451.2 / 97.6 283intergenic (‑86/+70)malZ/proYalpha‑glycosidase/proline‑specific permease ProY
Reads supporting (aligned to +/- strand):  ref base A (141/48);  new base C (4/89);  total (145/138)
Fisher's exact test for biased strand distribution p-value = 3.02e-32
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 301 total aligned reads displayed.

AGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGAAACC  >  DH1/3448376‑3448676
                                                                                                                                                      |                                                                                                                                                      
aGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggaggagga                                                                                                                                                        <  7:929643/151‑1 (MQ=255)
 gCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGAc                                                                                                                                                       <  2:372171/151‑1 (MQ=255)
      tGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACgggg                                                                                                                                                   <  3:1010216/150‑1 (MQ=255)
       gATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACggggg                                                                                                                                                  <  4:972116/150‑1 (MQ=255)
       gATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGc                                                                                                                                                 >  1:274235/1‑151 (MQ=255)
       gATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGc                                                                                                                                                 >  2:1012871/1‑151 (MQ=255)
       gATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGCGGGc                                                                                                                                                 >  7:405441/1‑151 (MQ=255)
           tttgcttttTTTAACAAATGGGGGCCCCGGCAGGGGCAATTCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATCCAGGATCACTCATCATCTCAAGTATATTTGGGGGCGGAGGACGGGGGCGGa                                                                                                                                              <  6:1118102/150‑1 (MQ=255)
            ttgctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCCCCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGCGGGCGGAt                                                                                                                                             >  3:798078/1‑150 (MQ=255)
             tgctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATa                                                                                                                                            >  3:216055/1‑150 (MQ=255)
             tgctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATa                                                                                                                                            >  1:311593/1‑150 (MQ=255)
                tttcTTTAAAAACTGGGGCCCCCGGCAGGTGCAATAAATTTAAAATAATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGCGGCGGACGGGGGCGGATAACt                                                                                                                                         <  2:561010/150‑1 (MQ=255)
                 ttgTTTAACAAATGGGGGCACCGGCAGGTGCCATTCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGTTTATTTGGGGGAGGAGGACGGGGGCGGATAAc                                                                                                                                          <  1:1064011/148‑1 (MQ=255)
                   gTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCg                                                                                                                                     >  7:772734/1‑151 (MQ=255)
                     ttAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCgg                                                                                                                                    >  6:947873/1‑150 (MQ=255)
                       aaCAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAAt                                                                                                                                 >  6:483341/1‑151 (MQ=255)
                       aaCAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAAt                                                                                                                                 >  3:447707/1‑151 (MQ=255)
                        aaaactGGGGGCCCCGGCAGGTGCCATGAATAAAACATCATAATTCCCCTGTTATTAAAATCTCATCACCTTGCCATACCCTAATGCAGCATCACTCATCATCTTAAGGATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATa                                                                                                                                <  2:1122335/146‑1 (MQ=255)
                        aCAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAAt                                                                                                                                 >  8:1067453/1‑150 (MQ=255)
                          aaaTGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGc                                                                                                                              >  1:1135219/1‑151 (MQ=255)
                           aaTGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGCGGAGGACGGGGGCGGATAACTTTCGGAATAGcc                                                                                                                             >  3:1107569/1‑151 (MQ=255)
                              gtgggCACCGGCATGTGAAATGATTTTAACATCATAATACGCCTTTGATTAAAATCTGACCACCATGCCTTATACTAATGCAGGATCACTCTTCATGTGTCGGTTATTTGGGGGCGGAGGCCGGGCGGGGATAACTTTCGGAATAGCCgg                                                                                                                           <  5:150582/148‑1 (MQ=255)
                              ggggcCCCCGGCAGGTTCAATTAATTTAAAATAATAATTCCCCTTTTATCAAAATCTGATCACGTTGCCATACCCTAATGCAGGATCACTCATCTTCTTATGTATTTTTGGGGCAGGCGGCGGGGGGCGGATAACTTTCGGAATAGCCGGa                                                                                                                          <  2:333861/151‑1 (MQ=255)
                               ggggCACGGGCAGGTGACATGAATTAAAAATATTAATTCCCATGTGATCAAAATGTTATCACGTTGCCATACCATAATGCAGTATAACTCATCATATTATTTTTTTTTGGGGTAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGAt                                                                                                                         <  5:538106/151‑1 (MQ=255)
                                  gCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATc                                                                                                                        >  6:208747/1‑149 (MQ=255)
                                   cACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAgg                                                                                                                     >  1:902738/1‑151 (MQ=255)
                                    cccGGCAGGTGCAATGCATTTAACATCATAATCCCCCTGTCATGAAAAGCTCAACCCCTGGCCACCGCCTAATGCAGCAGCACTCATCATCTCAAGCATATTGGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGa                                                                                                                    <  8:1014901/150‑1 (MQ=255)
                                    aCCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAgg                                                                                                                     >  8:529744/1‑150 (MQ=255)
                                     ccGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGAt                                                                                                                   >  5:1199834/1‑151 (MQ=255)
                                         gaGGTGAGATTCATTTAACATCATAATTCCGCTGTGATCAAATTCTGTTCAGCTCGCCATAGGTAAATTCAGGTTTAGTCATATTCTCAAGTTTATTTGGGGGCGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtc                                                                                                                <  7:667239/149‑1 (MQ=255)
                                          aGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcct                                                                                                              >  5:84922/1‑151 (MQ=255)
                                           ggTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGACCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcctc                                                                                                             >  8:1120819/1‑151 (MQ=255)
                                            gggCCATTCATTTAACATCATAATTCCCCTGTGATGAAAATCTGACCCCCTGGCCATCGCCTAATGCAGGATCACTCATCATCTGAAGTATATTTGGGCGCCGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcctcc                                                                                                            <  3:241122/149‑1 (MQ=255)
                                              gCCATGCATTTAACATCATAATTCCCCTGTGATGAAACTCTGATCACCTTCCCACAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTGGGGGGCGGCGCCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTg                                                                                                          <  5:75268/151‑1 (MQ=255)
                                                cATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATCCAGGATCACTCATCATCTGAAGGATAGTTCGGGGAGTAGGCCCGGGGCGGATACCTTTCCGAATAGACCGATCAGGATTTCCTCCGGa                                                                                                         >  2:656047/1‑150 (MQ=255)
                                                   tcATATAACATCAAAAATCCCCTGCGATCAAAATCTGATCACCTTCCCATCCCCTACTGCACGAACACTCATCATCTTACGGATTTTTGGGGGCGGAGGCGCGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAAcc                                                                                                      <  4:566067/149‑1 (MQ=255)
                                                   gCATTAAACATCATAATTCCCCTGTTATGAAAAGCTGATCACCTTCCCATAGCCTAATGCAGGATCACTCATCCTCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAAcc                                                                                                      <  5:127130/150‑1 (MQ=255)
                                                   ccATTTAACATCATAATCCCCCTGTCATCAAAATCTGATCACCTTGCCATAGCCTACTGCAGGATCACTCATCCTCTCAAGCATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCg                                                                                                     <  4:502818/150‑1 (MQ=255)
                                                     aTTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGc                                                                                                   >  7:615022/1‑151 (MQ=255)
                                                     aTTTAACACCATAATTCCGCTGTCATGAAAATCTGTTCACCTTGTCATAGCCTCTTCCAGCATCACTATTCTTGTGACGGATATTGGGGGGCGCCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGc                                                                                                   <  6:973901/151‑1 (MQ=255)
                                                      tttAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGcg                                                                                                  >  4:674666/1‑150 (MQ=255)
                                                      tttAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCt                                                                                                  >  4:1142725/1‑151 (MQ=255)
                                                      tttAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCt                                                                                                  >  2:608650/1‑151 (MQ=255)
                                                       taaaCCTCATAATTCCCGTGTTATGAAAATCTGATCACCTTCCCATACCCTAATCCAGGATCCCTCATCCTCTGAAGGATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCt                                                                                                  <  6:331975/148‑1 (MQ=255)
                                                           catcatAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAc                                                                                               >  1:50155/1‑149 (MQ=255)
                                                            atcaaAATTCCCCTGGGAGGAAAATCTGAGCACCTTGCCATACCCTACTGCCGGATCACTCCTCATCTCACGGATATTTGGGGGCGGACGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGa                                                                                            <  3:1139534/151‑1 (MQ=255)
                                                               aaaaTTCCCATGTGCACAAATTCTCAACCCCTGCCCAAAGCCCAATCCAGGATCACTCCTCATCTCTAGGATTTTTGGGGGCCCCGGCGGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACg                                                                                         <  7:154478/149‑1 (MQ=255)
                                                                tAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGc                                                                                        >  1:437491/1‑151 (MQ=255)
                                                                 aaTTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGc                                                                                        >  1:1261161/1‑150 (MQ=255)
                                                                   ttCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCa                                                                                       >  1:1102614/1‑149 (MQ=255)
                                                                    tCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTACTCCTGAACAGGCTTACAGAACGCAAgg                                                                                    >  6:879003/1‑151 (MQ=255)
                                                                    gccccTGTGAGTAAAATCTGATCACCCTGCCATCGCCTAATGCCGGATCACTCATCATCTCAAGGATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAgg                                                                                    <  5:349541/150‑1 (MQ=255)
                                                                    cccccTGTCATGAAAATCTGAGCACCTGCCCATACCCTAACGCAGCACCACTCCTCATCTGACGCATATTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAgg                                                                                    <  5:985544/150‑1 (MQ=255)
                                                                     ccccTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGACCCGGCTTACAGAACGCAAgg                                                                                    >  5:1007894/1‑150 (MQ=255)
                                                                       ccTTTGATGACCATCTGATCACCTTGCCATAGCCTAATGCAGGCTCACTCATCATCTGCAGGATATTTGGGGGAGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGAtt                                                                                 <  8:725732/151‑1 (MQ=255)
                                                                           tgatgaAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCGTCCTGAACCGGGTTACAGATCGGAGGGATTACTg                                                                             >  4:406780/1‑151 (MQ=255)
                                                                            gatgaAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGACCCGGCTTACAGAACGCAAGGATTACTgg                                                                            >  5:686819/1‑151 (MQ=255)
                                                                            gacgaAAATATCCTCACCTTGCCATACCCTAATGCTTGATCAGTCATCATGTGTAGTTTATGTGTGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTg                                                                             <  6:1137069/147‑1 (MQ=255)
                                                                              tgaAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGtt                                                                          >  4:290055/1‑151 (MQ=255)
                                                                               caaaccgcgcgcACCTGGCCATCGCCGAATGCACGATCCCTCCTCCTCGGTCGGATAGTGGGCGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGttt                                                                         <  4:798621/140‑1 (MQ=255)
                                                                                aaaagctcgccACCTTGCCATCCCCTAATGCAGGACCCCTAATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGttt                                                                         <  7:849141/140‑1 (MQ=255)
                                                                                aaaaTCTGCTCACCTTGCCATAGCCTAATGCACCATCACTCATCATCTCAAGGATTTTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGtttt                                                                        <  8:117905/151‑1 (MQ=255)
                                                                                 aaaTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTc                                                                       >  5:1289972/1‑151 (MQ=255)
                                                                                 aaaTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCCTATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTc                                                                       >  8:108532/1‑151 (MQ=255)
                                                                                  aaTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGCGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGACTTACAGAACGGAAGGATACTGGTTTTc                                                                       >  5:643304/1‑149 (MQ=255)
                                                                                         ccaccgacccATCGCCTACTGCAGGAGCACTCATCAGCTGAAGGTTATTTGGGGGCGGCGGCCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCt                                                               <  5:1041465/143‑1 (MQ=255)
                                                                                         ccACCTTCCCAAAGCCGACTGCCCCATCACCCATCCGCTGCAGTATCTTTGGGCGCCGACGACGGGGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCt                                                               <  3:569469/150‑1 (MQ=255)
                                                                                           aCCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGCATATTTGGCGGAGGAGGACGGGCGCGGATAACTTTCGGAAAAGCCGGATCACGATTTCCTACTCAACCGGCTAACAGAACGCAAGCAATACAGGTTTTcagccaaatg                                                              >  2:12489/1‑146 (MQ=255)
                                                                                            ccTTGCCAGAGCCTAATCCCGGAGCCCTCCTCATCTGAATGATATGCGGGGGCGGCGGCGGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGAc                                                            <  5:239449/151‑1 (MQ=255)
                                                                                                gCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATc                                                        >  1:472065/1‑151 (MQ=255)
                                                                                                cccATAGCCTAATGCAGGAGCACTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATc                                                        <  5:785845/150‑1 (MQ=255)
                                                                                                 ccatccccTAATCCAGGAGCACTCATCATCTGAAGTATTTTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCg                                                       <  6:681450/145‑1 (MQ=255)
                                                                                                 ccAGAGCCTCCTGCAGCACCACTCATCATCTGACGTACATTGGGCGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCg                                                       <  2:264369/151‑1 (MQ=255)
                                                                                                   aaccccTAATGCCGGACCACTCCTCATCTGAAGCATATTTGGGGGAGGCGGACGGCGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTg                                                     <  8:325779/147‑1 (MQ=255)
                                                                                                   aTAGCCTAATGCAGGCACACTCATCATCTTATGGATTTTGGGGGGCGGCGCACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTg                                                     <  4:1174103/151‑1 (MQ=255)
                                                                                                     aGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGCAGGACGACCGGGCCGGAAAACTTTCGGAATAGCCCGATCAGGATTTCCTCCTCAACCGGCTTACAGAACGCAAGGATTACAGGTTTTCCGCCCGCTGCCGCTCGTggc                                                   >  3:867031/1‑151 (MQ=255)
                                                                                                      cccTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGCGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcg                                                  <  4:723015/150‑1 (MQ=255)
                                                                                                       cataatggcgtaTCATTCACCATGTTAAGTATTTTTGGGGCCGGAGTCCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAAGCGGCTTAGAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcgg                                                 <  4:235058/140‑1 (MQ=255)
                                                                                                        cTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcggc                                                >  3:146755/1‑151 (MQ=255)
                                                                                                          aaTGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCGCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGCTCGTGGCGGCGt                                              >  2:917848/1‑151 (MQ=255)
                                                                                                             gCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGttt                                            >  8:244814/1‑150 (MQ=255)
                                                                                                             gCAGGAGCACTCATCATCTCAAGGATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGttt                                            <  7:1259763/150‑1 (MQ=255)
                                                                                                               aGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTaaa                                         <  8:411029/151‑1 (MQ=255)
                                                                                                                ggATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGACCGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAAc                                        >  4:196459/1‑151 (MQ=255)
                                                                                                                    cagtcctcatcTGAAGGATATTTGGGGGAGGCGGCCGGGGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACAt                                      <  2:1187912/143‑1 (MQ=255)
                                                                                                                     ccTCATCATCTGAAGGATATGGGGGCGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCa                                   <  5:1103423/150‑1 (MQ=255)
                                                                                                                     ccTCATCAGCTGAAGGATCTGGGGGGGAGGAGCCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCa                                   <  3:760246/150‑1 (MQ=255)
                                                                                                                     aCTCATCCTCTGAAGGATATTTGGGGGCGGCGGCCGCGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCa                                   <  1:180051/151‑1 (MQ=255)
                                                                                                                      cTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATcc                                    >  7:1036058/1‑149 (MQ=255)
                                                                                                                        cctcatcTGAAGGATATTGGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGcc                                <  1:329241/149‑1 (MQ=255)
                                                                                                                            atcTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATc                            >  2:523881/1‑151 (MQ=255)
                                                                                                                            atcTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATc                            >  4:884894/1‑151 (MQ=255)
                                                                                                                             tcTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATc                            >  8:178009/1‑150 (MQ=255)
                                                                                                                                  agtatctttGCGGGCCGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCACCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGc                      <  5:659722/145‑1 (MQ=255)
                                                                                                                                      tatTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGcaca                   >  1:1190387/1‑150 (MQ=255)
                                                                                                                                         ttGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCacaaca                >  5:42174/1‑150 (MQ=255)
                                                                                                                                              ggcggcggaCGGGGGGGGATAACTTTCGGAATAGCCGGATCAGTATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAc          <  4:137034/145‑1 (MQ=255)
                                                                                                                                              ggaggcggcCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAc          <  1:1007364/151‑1 (MQ=255)
                                                                                                                                              ggaggaggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAc          <  7:1179860/151‑1 (MQ=255)
                                                                                                                                                 ggaggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCAAcgcc       >  7:1274675/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
AGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGAAACC  >  DH1/3448376‑3448676

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