Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,448,5260AC28.6% 578.4 / 84.8 345intergenic (‑86/+70)malZ/proYalpha‑glycosidase/proline‑specific permease ProY
Reads supporting (aligned to +/- strand):  ref base A (186/58);  new base C (4/94);  total (190/155)
Fisher's exact test for biased strand distribution p-value = 6.10e-38
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 360 total aligned reads displayed.

AGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGAAAC  >  DH1/3448376‑3448675
                                                                                                                                                      |                                                                                                                                                     
 gCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGAc                                                                                                                                                      <  6:1098059/151‑1 (MQ=255)
   aGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACgg                                                                                                                                                    <  5:607711/151‑1 (MQ=255)
    gTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACggg                                                                                                                                                   >  3:135267/1‑151 (MQ=255)
    gTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACggg                                                                                                                                                   <  4:966143/151‑1 (MQ=255)
    gTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACggg                                                                                                                                                   >  1:903068/1‑151 (MQ=255)
      tGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACgggg                                                                                                                                                  <  7:160954/150‑1 (MQ=255)
      tGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCGAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACggggg                                                                                                                                                 >  6:1406570/1‑151 (MQ=255)
        aTCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGc                                                                                                                                                <  7:1210635/150‑1 (MQ=255)
         tCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCgg                                                                                                                                              >  4:263723/1‑151 (MQ=255)
            ttgcttttTTTAAAAATTGGGGGCACGGGCTGGTGCCATTCATTTAAAATAATAATACCCCTTTTATTAAAATCTGATCACCTTGCCATAGCCTAATGCAGCATCACTCATCATCTGAAGTATATTTGGGGGCGGAGGCCGGGGGCGGAt                                                                                                                                            <  8:1406704/150‑1 (MQ=255)
              gctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAAc                                                                                                                                         >  4:798749/1‑151 (MQ=255)
              gctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAAc                                                                                                                                         >  7:559054/1‑151 (MQ=255)
                tttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGAAAACtt                                                                                                                                       >  3:432653/1‑151 (MQ=255)
                 ttgTTAAACACATGGCGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTCATTAAAATCTCCTCACCTTGCCATAGCCTAATGCACGCTCACTCATCATCTGAAGGATATTTGGGGGCGGCGGCCGGCGGCGGATAACtt                                                                                                                                       <  4:1241496/150‑1 (MQ=255)
                   ttttAACAATTGGGTGCACCGGCAGGTGCAATTATTTTAACATCTTAATCCCCTTTTTATCAAAATATGTTCATGTTCCCTTAGTCTAATGCATCCTCACTCATGATCTTAAGTATTTTTGGGGGAGTTTGACGGGGGCGGATAACTTTc                                                                                                                                     <  1:858097/149‑1 (MQ=255)
                      tAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGa                                                                                                                                  >  1:1503218/1‑150 (MQ=255)
                         cAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAAt                                                                                                                                >  2:1357593/1‑149 (MQ=255)
                         cAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATac                                                                                                                              >  6:54733/1‑150 (MQ=255)
                          aaaTGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGc                                                                                                                             >  6:1022494/1‑151 (MQ=255)
                               gcggCCCCGGCCGGGGCCATGCATTAAACATCATAATTCCCCGGTGATCAAACTCAGACCCCCTGCCCATCCCCTAATGCAGGATCACTCAACATCTCAAGGATATGTGGGGGCGGAGGCCGGGCGCGGATAACTTTCGGAATAGCCGGAt                                                                                                                        <  5:697387/149‑1 (MQ=255)
                                 ggCAGCGGCCCGTGCCATTAATTTAACATCATAATCCCCCTGTGATGAAAATCTTCGCCCCTGGCCATCGCCTAATGGCGGATCACTCATGATCTGATGGATCTTGGGGGGCGGAGGCGGGGGGCGGATAACTTTCGGAATAGCCGGATCa                                                                                                                      <  3:1203330/151‑1 (MQ=255)
                                  ccACCGCCAGGTGCCACCCATTAAACATCAAAATTCCCCTGTCATCAAAATCTCATCACCTTCCCATACCCTAATCCAGCACCACTCATCATCTGAAGTATATTTGGGGGAGGCGGCCGGGGCCGGATAACTTTCGGAATAGCCGGATCAg                                                                                                                     <  6:1224291/150‑1 (MQ=255)
                                   cACCGGCAGGTGCCATTAATTAAAAAAAATAAAACCCCTTTTATCAAAATCTCATCAACTTCCCATCCCCTAATGCAGGATCACTCTTCATCTCATGTATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAgg                                                                                                                    <  7:918886/151‑1 (MQ=255)
                                    acgggCAGGTGCCATGAATTTAAAATCATAATTCCCCTGTTATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGTATCACTCATCATCTGAAGTATTTTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGa                                                                                                                   <  7:1482449/148‑1 (MQ=255)
                                       ggCAGGCGCCATGCATTAAACATCAAAACCCCCCTGTCATCAAAATCTGACCCCCTTCCCATCGCCTAATGCAGGATCACTCATCATCTCCAGCATATTTGGGGGCGCAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGAttt                                                                                                                <  2:694242/151‑1 (MQ=255)
                                          cggggCCATTAATTTAAATTAATAATTCCCCTGTTATGAAAATCTGAGCACCTTGCCATCGCCCTATGCAGGATCTGTCATCATCTGACGGATTTTTGGGGGCGGAGGAGGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcct                                                                                                             <  1:968064/147‑1 (MQ=255)
                                          aGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcct                                                                                                             >  8:941677/1‑151 (MQ=255)
                                          aGGTGCCATGCATTTAAAATCATAATTCCCCTGTGATCAAAATCTGATCACCTTGCCATCCCCTAATGCTCGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcct                                                                                                             <  8:493041/151‑1 (MQ=255)
                                           ggTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcctc                                                                                                            >  8:1522077/1‑151 (MQ=255)
                                            gTGCCATGCTTTTAACATCATAAGTCCCCTGTGATGAACATCTTCTCACCTTGCCATAGCCTACTGCAGGCGCACTCATCATCTGCAGGATTTTTGGGGTAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcctcc                                                                                                           <  4:887938/151‑1 (MQ=255)
                                            gTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcctcc                                                                                                           >  3:1398406/1‑151 (MQ=255)
                                              gCCATGCATTTAACATCATAATTCCCCTGTCATGAAAATCTGATCACCTTGCCATAGCATAATGCAGGATCACTCATCAACTGAAGGATATTTGCGGCAGCAGCACGGGGGCGCATAACTTTCGGAATAGCCGGATCAGGATTtcctaatc                                                                                                         >  7:314851/1‑147 (MQ=255)
                                                    cATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCgg                                                                                                   >  8:1235842/1‑151 (MQ=255)
                                                       ttAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCtt                                                                                                >  3:675970/1‑151 (MQ=255)
                                                        tAATATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCtt                                                                                                >  2:1461185/1‑150 (MQ=255)
                                                        tAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCtt                                                                                                >  5:1030491/1‑150 (MQ=255)
                                                           aatcatAATTCCCCTGTGATGAAACTCTGATCACCCTGCCATACCCTAATCCAGGACCACTCCTCCTCTGAAGGATATTTGGGGGAGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAg                                                                                            <  7:1050447/150‑1 (MQ=255)
                                                             tcatAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGa                                                                                           >  3:1037085/1‑150 (MQ=255)
                                                              aaaaaatCCCCTGTGATGAAAATCTCAACCCCGTGCCATACCCAAATGCCGCACCACGCATCACCTGACGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAAc                                                                                         <  3:940016/145‑1 (MQ=255)
                                                               atcctcccccTGTGATGAAACTCTCATCACCTTCCCATACCCTAATGCAGCATCACTCATCATCTGACGGATATTGGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACg                                                                                        <  8:414023/145‑1 (MQ=255)
                                                               atAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGGTCAGGATTTCCTCCTGAACCGGCTTACAGAACg                                                                                        >  1:1387598/1‑151 (MQ=255)
                                                               atAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACg                                                                                        >  4:858398/1‑151 (MQ=255)
                                                                 aaTTCCCCTGTTCTTAAAATCTGCTCCCCTTCCCATAGCCTAATGCAGGATCACTCATCATCTGACGGAGAGTTGGGGGCCGAGGCCGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGc                                                                                       <  2:1064072/150‑1 (MQ=255)
                                                                  aTTCCCATGTGATGAAAATCTGAGCCCCTTCCCATAGCCTAATCCAGGATCCCTCTTCATCTGAAGGATATTGGGGGGAGGAGGACGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCaa                                                                                     <  2:395510/151‑1 (MQ=255)
                                                                   ttCCCCTGTGATGAAAATCTGATCACATTGCAATAGCCTAATGCAGGATCACTCATCATCTGAAGTATATTTGGGGGAGGAGGACGGGGGCGGATAAATTACGAATTAGCCGGATCAGCATTTCATCATGAACCGGCTTAAAGAAGGCa                                                                                      >  8:820946/1‑149 (MQ=255)
                                                                   atCCCCTGTGATCAAAAGCTGATCCCCTTGCCATCCCCTATTCCCGGATCACTCATCATCTGACGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAg                                                                                    <  8:670791/150‑1 (MQ=255)
                                                                    tCCCCTGTCATGAAAACCTGATCACCGTCCCATCGCCTAATGCAGGATCACTCATCATCTGAAGGATCTGTGGGGGAGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAg                                                                                    <  5:326879/150‑1 (MQ=255)
                                                                     cccgTGTGATCAATATCTGATCACGTTGCCATAGCCTAATGCAGGATCACTCATCTTGTGACGGATTTTTGGGGGCGGTGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAgg                                                                                   <  2:1177313/150‑1 (MQ=255)
                                                                       ccTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGAtt                                                                                >  1:1047772/1‑151 (MQ=255)
                                                                       catgttATCAAAATCTCATCACCTAGGCACAGCCTAATGCAGGATCACTCATCATCTGAATTATATTTGGGGGAGGTGGTCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGAtt                                                                                <  5:85602/149‑1 (MQ=255)
                                                                         tgtgATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTAc                                                                              >  6:716929/1‑151 (MQ=255)
                                                                          gttatgaAAATCTGAGCACCTTGCAATAGCATAATGCAGTAGCAGTCATGATCTTAAGTATTTTTGGGGGCGGAGGAGGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGGTTACAGAACGCAAGGATTa                                                                               <  4:1120814/146‑1 (MQ=255)
                                                                           ggatcaAATTCTGCTCACCTTGCCCTCGCCTAATCCAGGATCACTCATCATCTTAACGATTTTTGGGGGCGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTg                                                                            <  6:966153/150‑1 (MQ=255)
                                                                              ataaaaTATGATCACCTTGCCATATCCTAATGCAGCATCACTCATCATCTTAAGTATATGTGGGGGAGGAGCCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGtt                                                                         <  3:333188/149‑1 (MQ=255)
                                                                                aaaaTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGtttt                                                                       >  1:605971/1‑151 (MQ=255)
                                                                                 aaaTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGCTTTTc                                                                      >  8:1115872/1‑151 (MQ=255)
                                                                                  aagcTGATCACCTTCCCATACCCTAATGCAGGAGCACTCATCATCTCAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTca                                                                     <  5:127799/148‑1 (MQ=255)
                                                                                  aaTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTc                                                                      >  5:430312/1‑150 (MQ=255)
                                                                                   aTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTc                                                                      >  1:476650/1‑149 (MQ=255)
                                                                                       tctCACCTTGCCATAGCATATTGCAGGATCACTCATCTTCTGATGTTTTTTTGGGGGTGGCGGACGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcca                                                                 <  6:580046/148‑1 (MQ=255)
                                                                                       gagcagcTTGCCATAGCCTATGGGCGGATCAGACATCTTCTGAAGCTTATTTGGGGGAGGAGGCGGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagccag                                                                <  7:262784/145‑1 (MQ=255)
                                                                                       gATCACCTTGCCATCGCCTTCTGCAGGATCATTCATCATCTTATTGATTTTTGGGGGAGGAGGCGGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcca                                                                 <  2:519952/150‑1 (MQ=255)
                                                                                         tCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCt                                                              >  5:1415048/1‑151 (MQ=255)
                                                                                            ccTTCCCATACCCTAATCCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGAc                                                           <  4:576129/151‑1 (MQ=255)
                                                                                               tGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCTGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGAt                                                        >  3:537728/1‑151 (MQ=255)
                                                                                               tGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGAt                                                        >  7:943997/1‑151 (MQ=255)
                                                                                               tGCCATACCCTAATGCAGCATAAGTCATCATCTGAAGGATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACg                                                          <  3:559549/149‑1 (MQ=255)
                                                                                               agCCTTATCCTAATGGAGGATCCCTCAGCATCTGATGGATATTTGGGGGAGCAGGCCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGa                                                         <  6:915355/149‑1 (MQ=255)
                                                                                                 ccATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCg                                                      >  7:1350608/1‑151 (MQ=255)
                                                                                                  cATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGt                                                     >  6:363311/1‑151 (MQ=255)
                                                                                                   aTAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTg                                                    >  3:1538039/1‑151 (MQ=255)
                                                                                                   aTACCCTAATGCAGGATCACTCATCATCTGAAGGATATTGGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTg                                                    <  7:654893/151‑1 (MQ=255)
                                                                                                    aaGCCTAACGCACCAGCACTCATCAGCTGACGGATAGGGGGGGCCGGCGGCCGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTgg                                                   <  7:1130024/150‑1 (MQ=255)
                                                                                                       ccTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcgg                                                >  1:286433/1‑151 (MQ=255)
                                                                                                         tAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGACAGCTGACGATCGTGgcggcg                                              >  7:50085/1‑151 (MQ=255)
                                                                                                          aaTGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGt                                             >  1:1207840/1‑151 (MQ=255)
                                                                                                             gCAGGATCACGCATCAGCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTa                                          <  3:611565/151‑1 (MQ=255)
                                                                                                               agcacccctccacatcTGAAGGATATTTGGGGCCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTaaa                                        <  6:619715/139‑1 (MQ=255)
                                                                                                                    cACTGATCTTCTGATGCATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATcc                                   <  2:634363/151‑1 (MQ=255)
                                                                                                                     aCTCATCATCTGACGGATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCa                                  <  5:788730/151‑1 (MQ=255)
                                                                                                                     aCTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATcc                                   >  1:452830/1‑150 (MQ=255)
                                                                                                                         atcatcTGCAGGATATTTGGGGGCGGAGGCCGCGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCa                              <  2:114836/151‑1 (MQ=255)
                                                                                                                              gtgtcggagtgttGGGCGCGGAGGAGGGGGGCGGAGAACTTTCGGAATAGCCGGATCAGGATTTCCTCGTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGATGACGATCGTGGGGGCGTTTAAACATCCAGCCATTCaa                         <  4:324150/140‑1 (MQ=255)
                                                                                                                              cTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGACCCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCa                          >  6:1485751/1‑150 (MQ=255)
                                                                                                                              cTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCa                          >  8:604515/1‑150 (MQ=255)
                                                                                                                                gAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACCTCCAGCCAATCAAc                        >  1:1325525/1‑150 (MQ=255)
                                                                                                                                gAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAAc                        >  4:156179/1‑150 (MQ=255)
                                                                                                                                gAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAAc                        >  5:1386222/1‑150 (MQ=255)
                                                                                                                                 aaGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGACCGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAAc                        >  1:761183/1‑149 (MQ=255)
                                                                                                                                      tatTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGCACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCCGCCAATCAACAGcaca                  >  2:798814/1‑150 (MQ=255)
                                                                                                                                      tatTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGcaca                  >  8:709581/1‑150 (MQ=255)
                                                                                                                                        tttGGGGGCGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCacaac                <  4:1411106/150‑1 (MQ=255)
                                                                                                                                         ttcggggAGGAGGACGGGGGCCGATAACTTTCGGAATAGCCGGATCAGCATTTCCTCCTGAACCGGCATACAGAACGCAAGGATTAGAGGTTTTCAGCCAGCAAACGATCTAGGCGGCGTTTAAACATCCAGCCAATCAACAGCacaaaa               >  4:1042753/4‑148 (MQ=255)
                                                                                                                                              ggccgaggcCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAc         <  2:252002/147‑1 (MQ=255)
                                                                                                                                               gaggaggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAcg        >  6:1122328/1‑151 (MQ=255)
                                                                                                                                                  gccgaCGGCCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGa     <  4:437789/148‑1 (MQ=255)
                                                                                                                                                   aggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCTAGGATTACTGGTTTTCAGCCAGCTAACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCCCCCGa     >  8:1074927/1‑150 (MQ=255)
                                                                                                                                                   aggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGa     >  1:1135506/1‑150 (MQ=255)
                                                                                                                                                   aggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGa     >  4:523125/1‑150 (MQ=255)
                                                                                                                                                     gaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGAAAc  >  8:237282/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGAAAC  >  DH1/3448376‑3448675

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