Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,448,5260AC25.5% 775.2 / 85.2 418intergenic (‑86/+70)malZ/proYalpha‑glycosidase/proline‑specific permease ProY
Reads supporting (aligned to +/- strand):  ref base A (234/76);  new base C (4/102);  total (239/179)
Fisher's exact test for biased strand distribution p-value = 2.93e-42
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 435 total aligned reads displayed.

AGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGAAACC  >  DH1/3448376‑3448676
                                                                                                                                                      |                                                                                                                                                      
cgCAGTTGATCTTTGCTTTGTTTAACAACTGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggaggagga                                                                                                                                                        <  2:1410161/150‑1 (MQ=255)
aGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggaggagga                                                                                                                                                        <  6:1554450/151‑1 (MQ=255)
 ggaGTTGTTGTTTGGTTTGTTTAAGAAATGGGGGCACGGGGTGGTGCGATGGATTTAAGATGATATTTCGCCTGTGATGAAAATCTGATGAGGTTGGGATAGCCTAATGCAGGATGAGTCATGATGTGAAGGATATTTGGGGGAGGAGGAc                                                                                                                                                       <  6:259892/149‑1 (MQ=255)
   aGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACgg                                                                                                                                                     >  7:708796/1‑151 (MQ=255)
       gATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATGATAATTCCCCTGTGATGAAAATCTGATCAGCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGc                                                                                                                                                 >  6:581108/1‑151 (MQ=255)
           tttgctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGAt                                                                                                                                             >  2:1441561/1‑151 (MQ=255)
            ttgctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCg                                                                                                                                                <  4:72122/147‑1 (MQ=255)
              gctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATaa                                                                                                                                           >  6:1215689/1‑150 (MQ=255)
               gttttatatacacctGGGGGCACCGGCTGGTGCCATCCATTTAAAAACATACATCCCCTGTGATCAAAAGCTTATCCCCTCCCCTTACCCATATGCAGGATCACTCATCATCTGACGCATAGTTGGGGGCGGCGGCCGGGGGGGGATAACt                                                                                                                                         <  2:263656/137‑1 (MQ=255)
               ctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACt                                                                                                                                         >  3:1031649/1‑151 (MQ=255)
                   gTTTAACAAATGGGGGAACAGGCACGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCCCCTTTCCATAGCCTAATGCAGGATCACTCATCATCTGAAGCATATTTGGGGGAGGAGGACGGGGGCGCATAACTTTCg                                                                                                                                     >  4:1735191/1‑151 (MQ=255)
                    tttAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCgg                                                                                                                                    >  5:1525610/1‑151 (MQ=255)
                    tttAACAAATGGGGGCACCGGCAGGTGCAATTCATTTAAAATCATAAATCCGATTTGATGAAAATCTGATCACCTTGCAATAGCCTAATGCATGATCACTAATCATCTGAAGGATATTTGGGGGAGGAGGACGCGGGCGGATAACTTTCg                                                                                                                                     <  1:691317/150‑1 (MQ=255)
                     ttAACAAATGGGGGCCCCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATACCTTTCgg                                                                                                                                    >  5:1738934/1‑150 (MQ=255)
                     ttAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCgg                                                                                                                                    >  7:1293839/1‑150 (MQ=255)
                      tAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGa                                                                                                                                   >  4:761308/1‑150 (MQ=255)
                        acacctGGGGGCAGCGGCAGGTGGCATTCATTTAACATCATAATTCCCGTGTTATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGATGGTTTTGGGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATa                                                                                                                                <  7:366330/146‑1 (MQ=255)
                         caatggggggCACCGGCAGGTGCCATGAATTTAACATCATAATTCCCCTGTTATTAAAATCTGATCACCTTGCCATAGCCTAATGCCGGATCACTCATCATCTGAATTATATTTGGGGGAGGCGGACGGGGGCGGATAACTTTCGGAATAg                                                                                                                               <  3:968165/146‑1 (MQ=255)
                         cAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGCTGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAAt                                                                                                                                 >  3:1103170/1‑149 (MQ=255)
                         cAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAAt                                                                                                                                 >  1:1604759/1‑149 (MQ=255)
                          aaaTGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGc                                                                                                                              >  2:1532234/1‑151 (MQ=255)
                          aaaTGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGc                                                                                                                              >  8:1807314/1‑151 (MQ=255)
                           aaTGGGGGCACCGGCAGGTGCCATGCATTTAAAATCATAATTCCCCTGATATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGTGGGAGGAGGACGGGGGAGGATAACTTTCGTAATAGcc                                                                                                                             >  2:1696588/1‑151 (MQ=255)
                            aTGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCg                                                                                                                            >  7:1379845/1‑151 (MQ=255)
                              ggggcCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGTATATTTGGGGGCGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGa                                                                                                                          <  2:772092/151‑1 (MQ=255)
                                 ggCACCGGGAGGTGCCATTCATTTAACATAATAATTCCCATGTTATGAAAATCTGATCACCTTGCCATAGCCTAATGGAGGATCACTCATCATCTGTAGTATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCa                                                                                                                       <  7:151178/151‑1 (MQ=255)
                                 ggCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATc                                                                                                                        >  2:816832/1‑150 (MQ=255)
                                  gCACCGGCAGTTTCCATTCATTTAAAATCATATTTCCCCTTTTATTAAAATCTTATCACCTTTCCATTTCATAATGCATGTTCACTCATATTCTTAAGTATTTTTGGTGGTGGAGGCGGGGGGCGGATAACTTTCGGAATAGCCGGATCAg                                                                                                                      <  4:1408892/151‑1 (MQ=255)
                                  gCACCGGCAGGTGCCATTCATTAAACATCATAATTCCCCTGTGATCAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGTATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAg                                                                                                                      <  7:1482630/151‑1 (MQ=255)
                                  gCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATc                                                                                                                        >  2:1145693/1‑149 (MQ=255)
                                   ccccGGCCGTTGCAATGCATTTAACATCAAAATTCCCCTGTGATTAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAgg                                                                                                                     <  8:247371/149‑1 (MQ=255)
                                    acgggCAGGTGCAATTAATTTAACATCATAATTCCCCTGTGATCAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTGGGGGGAGCAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGa                                                                                                                    <  7:118413/148‑1 (MQ=255)
                                    aCCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAgg                                                                                                                     >  1:1308461/1‑150 (MQ=255)
                                     ccGGCAGGTGCCATTCATTTAACATCATAATTCCCCTGAGATCAAACTCTGAACACCTTGCCATACCCTCATGCAGGATCCCTCCTCATCTGACGGATATTGGGGGGCGGCGGCCGGCCGCGGATAACTTTCGGAATAGCCGGATCAGGAt                                                                                                                   <  5:529317/151‑1 (MQ=255)
                                        gCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtc                                                                                                                >  4:1292620/1‑151 (MQ=255)
                                        gCAGGGGCCATGCATTAAACATCATCCTTCCCCTGTCATGAAACTCTGATCCCCTTGCCATAGCCTACTGCCGGATCACCCATCATCTGACGGATATTTGGGGGCGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtc                                                                                                                <  4:1633021/151‑1 (MQ=255)
                                          agtctccATGCATTAAACATAATAATTCCCCTGTGATGAAAATCTGATCACCTTCCCATAGCCTAATGCAGGACCACTCATCATCTGAAGGATATTTGGGGGAGGCGCACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcct                                                                                                              <  5:1316148/146‑1 (MQ=255)
                                           ngTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcctc                                                                                                             >  4:1550298/2‑151 (MQ=255)
                                              gCCATGCTTTTAACATCATAATGCCCCTGTGATCAAAATCTGATCCCCTGGCCATAGCCTAATGCAGGACCACTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTg                                                                                                          <  1:1032994/151‑1 (MQ=255)
                                               ccATGCATATAACCTCATAATTCCCCTGTCAATAAAATCTCTCCACCTTGCCCTACCCTACTGCACGCTCCCTCATCATCTGACGGATATTTGGGGGAGGCGGCCGCGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTg                                                                                                          <  2:965419/150‑1 (MQ=255)
                                                   gCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTACTGAACCg                                                                                                     >  4:957043/1‑151 (MQ=255)
                                                      tttAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCt                                                                                                  >  5:196350/1‑151 (MQ=255)
                                                        tAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCAGAACCGGCtt                                                                                                 >  4:1618695/1‑150 (MQ=255)
                                                        aaacaacAAAATTCCCGTGTGATGAAAATCTGATCACCTTGCCATAGCGTAATGCAGCATCCCTCATCATCTGAAGGATTTTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTa                                                                                                <  8:781049/150‑1 (MQ=255)
                                                         aacaacATAATCCCCCTGAGATCAAAATCTGATCACCTTCCCATAGCCGAAGCCACGACCACGCATCAGCTGAAGGATATTGGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAc                                                                                               <  7:1662766/151‑1 (MQ=255)
                                                         aaCATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGCACCGGCTTAc                                                                                               >  2:54110/1‑151 (MQ=255)
                                                         aaCATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAc                                                                                               >  3:1193896/1‑151 (MQ=255)
                                                         aaCATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGACGGATCAGGATTTACTCCTGAACCGCCTTa                                                                                                >  5:1159969/1‑150 (MQ=255)
                                                            atcatAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAc                                                                                               >  6:798274/1‑148 (MQ=255)
                                                            aacaaaaTCCCCCTGTCAAGAAAATCGCCTCACCGTCCCATCCCCTAACCCCGGAGCACTCCTCATCTGAAGGATATTGGGGGGCGGCGCCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGa                                                                                            <  3:1539837/146‑1 (MQ=255)
                                                                 aaTGCCCCTGTTATCAAAATCTGATCACCCTCCAATAGCCTAATCCAGGCTCACTCATCATCTGAAGGATTTGTGGGGGCGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGc                                                                                        <  6:946301/150‑1 (MQ=255)
                                                                  aTTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCa                                                                                       >  6:1705176/1‑150 (MQ=255)
                                                                    tCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTAGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAgg                                                                                    >  8:1753039/1‑151 (MQ=255)
                                                                       ccTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGAtt                                                                                 >  6:133846/1‑151 (MQ=255)
                                                                          gtgaataaaaTCTCATCACCTTCCCATACCCTAATGCATCACCACTCATCATCTGAATCATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTAc                                                                               <  8:1244960/144‑1 (MQ=255)
                                                                            tatgaAAATCTGATCACCTTCCCATCGCCTACTGCAGGATCACTCATCATCTGAAGGTTATTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTg                                                                             <  3:212780/149‑1 (MQ=255)
                                                                              tgaAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGtt                                                                          >  6:648488/1‑151 (MQ=255)
                                                                               aaaaaTCTGATCACCTTCCCATACCCTAATGCAGGACCACTCATCTGCTGACGGTTATTGGGGGGAGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGttt                                                                         <  2:438732/150‑1 (MQ=255)
                                                                                  aaTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTc                                                                       >  3:706029/1‑150 (MQ=255)
                                                                                   aTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTc                                                                       >  2:1049872/1‑149 (MQ=255)
                                                                                       gATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcc                                                                   <  4:420580/149‑1 (MQ=255)
                                                                                       gATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcc                                                                   >  4:1042588/1‑149 (MQ=255)
                                                                                          caacTTGCCATAGCCTAATGCACGACCACTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCt                                                               <  1:1142485/147‑1 (MQ=255)
                                                                                          cACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTAATGGATTTCAGCCAGATg                                                              >  2:655308/1‑151 (MQ=255)
                                                                                            ccTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGAc                                                            >  3:555796/1‑151 (MQ=255)
                                                                                            ccTTGCCATAGCCTAATCCAGCATCACGCATCATCTGAAGGATAATTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGAc                                                            <  5:1689574/151‑1 (MQ=255)
                                                                                            ccTTGACATATCCTATTGCAGTATCAATCATCATCTTAATTTTATTTGGGTGCGGTGGAGGGGGGGGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGAc                                                            <  2:881878/151‑1 (MQ=255)
                                                                                             cTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACg                                                           >  8:410439/1‑151 (MQ=255)
                                                                                               tGCCATAGCCTAATGCAGTCTCACTCATCATCTGCAGGATAGTTGGGGGAGGAGGCCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACg                                                           <  2:1294332/149‑1 (MQ=255)
                                                                                                gCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATc                                                        >  1:774797/1‑151 (MQ=255)
                                                                                                gCCATAGCCTAATGCAGGAGCACTAATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGAt                                                         <  1:219428/150‑1 (MQ=255)
                                                                                                cccctAGCCTACTCCAGGATCTCTCTTCATCTTTAGTATTTTGGGGGGAGGCGCCCGGGGGGGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATc                                                        <  8:870839/147‑1 (MQ=255)
                                                                                                   aTAGCCTAAGGCAGCATCACTCATCATCTGACGGATATTGGGGGGCGCAGCACGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTg                                                     <  5:1493492/151‑1 (MQ=255)
                                                                                                    tAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTgg                                                    >  7:4890/1‑151 (MQ=255)
                                                                                                     aGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggc                                                   >  8:569806/1‑151 (MQ=255)
                                                                                                        cTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcggc                                                >  3:377999/1‑151 (MQ=255)
                                                                                                         tAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGCTCGTGgcggcg                                               >  1:437519/1‑151 (MQ=255)
                                                                                                          aaTGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGt                                              >  1:368559/1‑151 (MQ=255)
                                                                                                          aaTGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGt                                              >  7:990264/1‑151 (MQ=255)
                                                                                                            cccacgATCACTCATCATCTGCCGGATAGTTGGGGGCGGAGCCCGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGttt                                            <  6:289364/146‑1 (MQ=255)
                                                                                                             gCACGATCACTCATCATCTGACGTATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTa                                           <  8:481837/151‑1 (MQ=255)
                                                                                                               aGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTaa                                          >  1:290883/1‑150 (MQ=255)
                                                                                                                  aTCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCCGCGTTTAAACAt                                      >  5:756127/1‑151 (MQ=255)
                                                                                                                    cACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATcc                                    >  4:590501/1‑151 (MQ=255)
                                                                                                                     aCTCATCATCTGAAGTATCTTTGGGGGCGGAGCACGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCa                                   <  5:689653/151‑1 (MQ=255)
                                                                                                                       tcatcaTCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGc                                 >  4:1339756/1‑151 (MQ=255)
                                                                                                                        catcatCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGCTTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGcc                                >  3:1576866/1‑151 (MQ=255)
                                                                                                                        catcagCTGAAGGATATTGGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGcc                                <  1:1416874/151‑1 (MQ=255)
                                                                                                                           catcTGAAGGATTTTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAAt                             <  7:1100579/151‑1 (MQ=255)
                                                                                                                           catcTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAAt                             >  1:1697575/1‑151 (MQ=255)
                                                                                                                             tcTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGCTCGTGGCGGCGTTTAAACATCCAGCCAATc                            >  3:1407072/1‑150 (MQ=255)
                                                                                                                             tcTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATc                            >  5:315474/1‑150 (MQ=255)
                                                                                                                                gAAGGATATTTGGGGGCGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATc                            <  2:512716/147‑1 (MQ=255)
                                                                                                                                gAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGCCGGCGTTTAAACATCCAGCCAATCAAc                         >  5:253430/1‑150 (MQ=255)
                                                                                                                                       atTTGGGGGAGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCTATCAACAGCacaac                 <  4:340547/151‑1 (MQ=255)
                                                                                                                                           gggggCGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATc             <  8:1440457/151‑1 (MQ=255)
                                                                                                                                            ggggCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATcc            <  8:1399255/151‑1 (MQ=255)
                                                                                                                                               gaggaggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAcg         >  7:188030/1‑151 (MQ=255)
                                                                                                                                                    ggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGaa     >  7:244511/1‑150 (MQ=255)
                                                                                                                                                      aCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGAAAcc  >  5:1437912/1‑151 (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 < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: