Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,448,5260AC28.0% 588.1 / 85.8 318intergenic (‑86/+70)malZ/proYalpha‑glycosidase/proline‑specific permease ProY
Reads supporting (aligned to +/- strand):  ref base A (180/49);  new base C (3/86);  total (183/135)
Fisher's exact test for biased strand distribution p-value = 4.97e-38
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 334 total aligned reads displayed.

AGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGAAACC  >  DH1/3448376‑3448676
                                                                                                                                                      |                                                                                                                                                      
aGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggaggagga                                                                                                                                                        <  3:1054640/151‑1 (MQ=255)
aGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggaggagga                                                                                                                                                        <  3:378628/151‑1 (MQ=255)
    gTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACggg                                                                                                                                                    <  3:809085/151‑1 (MQ=255)
         tgtttgctttgTTTAACTAATGGGGTCACCGTCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGTTCACCTTTCCATAGCCTAATGCAGTATCACTCATCATCTGAAGGATATTTGGGGGAGTAGGACGGGGGCg                                                                                                                                                <  6:173352/148‑1 (MQ=255)
                 ttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACttt                                                                                                                                       >  3:852468/1‑151 (MQ=255)
                      aaaCAAATGGGGCCCCCGGCAGGTGCCATGAATTAAACATCATAATTCCCCTTTCATCAAAATCTCATCACCTTGCCATAGCCTAATGCAGCATCACTCATCATCTGAAGGATTTTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGaa                                                                                                                                  <  2:591551/150‑1 (MQ=255)
                       aaCAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAAt                                                                                                                                 >  1:42814/1‑151 (MQ=255)
                       aaCAAATCGGGGCACCGGCAGGTGCCATCCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGCATATTTGGGGGAGGACGACGGGGGCGGATAACTTTCGGAAt                                                                                                                                 >  4:124384/1‑151 (MQ=255)
                           aaTGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGcc                                                                                                                             >  4:400057/1‑151 (MQ=255)
                            aTGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGACg                                                                                                                            >  7:779642/1‑151 (MQ=255)
                            aTGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATACCCg                                                                                                                            >  8:281933/1‑151 (MQ=255)
                              gggggCCCGGGCAGGGGCCATGCATTTAACATCATAATTCCCGTGTGATGAAATTCTGATCACCTTGCCATAGCCTATTGCAGGATCACTCATCATGTGAAGTATTTTTGGGGGCGGAGGACGGGGGCGGATAACTTTCGGAATAGCCgg                                                                                                                           <  1:1133334/150‑1 (MQ=255)
                              gggggCCCCGGCAGGTGGCATGCATTTAACATCATAATTCCCCTGTTATTAAAATCTCATCACCTAGCCATAGCCTCATGCTGGATCAGTAATCATCTGAAGGATTTTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCg                                                                                                                            <  6:808/149‑1 (MQ=255)
                              ggggcCACCGGCAGGTGCCATGCATTTAACATAATAATGCCCCTGTGATGAAAATCTCCTCACCTTGCCATCGCCTAATGCAGGATCACTCATCATCTGAAGTATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCgg                                                                                                                           <  1:892548/150‑1 (MQ=255)
                                 ggCACCGGCTGGTGCCATTAATATAACATCATAATTCCCCTGTTATTAAAATCTGCTCCCCTTGCCATAGCCTAAGGCCGCATCACTCATCATCTGCAGTATATTGGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCa                                                                                                                       <  2:683519/151‑1 (MQ=255)
                                  gCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATc                                                                                                                        >  3:90560/1‑149 (MQ=255)
                                   ccccGGCTGGTGCCATGCATTAAACACCATAATTCCCCTGTGATCAAAATCTGATCCCCTTGCCATAGCCTAATGCAGCATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAgg                                                                                                                     <  2:1140183/149‑1 (MQ=255)
                                      cGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATCAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAACGATATTTGGGGGAGGAGGACGGGGGCGGATAACGTTCCGAATAGCCGGATCAGGAtt                                                                                                                  >  2:412577/1‑151 (MQ=255)
                                        gCATGTGCCATGCATTTAACATCATAATTCCCCTGTTATAAAAATCTGATCACCTTCCCATACCATAATGCAGGATCACTCATCATCTTAATGATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtc                                                                                                                <  4:80953/151‑1 (MQ=255)
                                          aGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcct                                                                                                              >  6:231168/1‑151 (MQ=255)
                                          aGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcct                                                                                                              >  1:776956/1‑151 (MQ=255)
                                            gTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcctcc                                                                                                            >  3:1136343/1‑151 (MQ=255)
                                              gcaattcATTAAACATCAAAATTCCCCTGTCATCAAAATCTCATCACCGTGCCATCCCCTAATGCAGCCTCACTCATCATCTGAAGGTTATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTcctcct                                                                                                           <  6:302494/144‑1 (MQ=255)
                                              gCCATGCATTTAAAATAAAACTTCCCCTGTGATCAAAATCTGCTCCCCTTGCCATCCCCTAATGCAGCATCACTCATCATCTGACGCATATTGGGGGGCCGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTg                                                                                                          <  1:565750/151‑1 (MQ=255)
                                              gCCATCCATTAAACATCATAATGCCCCTGTCATGAAAATCTGATCCCCTTCCCATAGCCTCATGCAGGATCACTCATCATCTGCCGGCTCGTGGCGCGCCCCGCCCGCGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTg                                                                                                          <  7:1003862/151‑1 (MQ=255)
                                                cATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGa                                                                                                         >  5:675486/1‑150 (MQ=255)
                                                  tGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGAGTTCCTCCTGAAcc                                                                                                      >  2:300096/1‑151 (MQ=255)
                                                   tcATTTAACATAATTATTCCCTTTTGATTAACATCTGATCACCTTGACTTATCCTAATGGCTTATCACTCTTCATCTGATGTTTATTGGGGGGCGGCGGCCGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAAcc                                                                                                      <  4:897139/149‑1 (MQ=255)
                                                   taaTTTAACATCATAATTCCCGTGTCATTAAAATGTCATCACCTTGCCATCGCCTAATTCAGCATCCCTCTTCAGCTGACGGTTATTTGGGGGCGGCGGCCCGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAAcc                                                                                                      <  5:1355387/148‑1 (MQ=255)
                                                   gCATTTAACATCATAAGTCCCCTGTGATGAAAATCTGACCCCCTTCCCATAGCCTAATGCAGGAGCACTCATCATCTGAAGGATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCg                                                                                                     <  5:137014/151‑1 (MQ=255)
                                                     aTTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGc                                                                                                   >  2:900566/1‑151 (MQ=255)
                                                     aTTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGc                                                                                                   >  7:947460/1‑151 (MQ=255)
                                                     aTTTAACAATATAATTCCCCTGTTTTTAAAATCTGAGCAACTTGCCATACCCTAATGCAGGATCACTCATCATCTGATGGATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGc                                                                                                   <  8:383566/151‑1 (MQ=255)
                                                        tAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCtt                                                                                                 >  5:578256/1‑150 (MQ=255)
                                                         aaCATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTCAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAc                                                                                               >  1:375540/1‑151 (MQ=255)
                                                            atcatAATTCCCCTGTTAAGAAAATCTCCTCACCTTCCAATCCCCTACTCCAGGAGCACTCATCATCTGAAGGATATTTGGGGGCGGCCGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGa                                                                                            <  4:7010/151‑1 (MQ=255)
                                                            atcatAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAc                                                                                               >  7:1185750/1‑148 (MQ=255)
                                                            atcatAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAc                                                                                               >  7:150909/1‑148 (MQ=255)
                                                            atcatAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAc                                                                                               >  2:143105/1‑148 (MQ=255)
                                                                 aaTTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGc                                                                                        >  3:669273/1‑150 (MQ=255)
                                                                   ttCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCa                                                                                       >  5:304937/1‑149 (MQ=255)
                                                                    gccccTGTGATGAAAATCTCCTCCCCTTCCCAGAGCCTAATCCACGATCACTCATCATCTGAAGGATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAgg                                                                                    <  1:347842/150‑1 (MQ=255)
                                                                     ccccTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAgg                                                                                    >  5:780219/1‑150 (MQ=255)
                                                                          gtgATTAAACTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGATGGATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTa                                                                                <  2:810993/149‑1 (MQ=255)
                                                                           tgatgaAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTg                                                                             >  5:406946/1‑151 (MQ=255)
                                                                            gatgaAAAGCTGATCACCTTCCCATCGCCAAATGCAGCATCCCTCATCATCTGACGGATATTTGGGGGCGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTg                                                                             <  3:598895/150‑1 (MQ=255)
                                                                                aaaaTCTGATCCCCGTCCCAAAGCCTCAGGCCGGAGCCCGCCTCATCTGAAGGATCTTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGtttt                                                                        <  5:901335/151‑1 (MQ=255)
                                                                                aaaaTCTGATCACCTGGCCATAGCCTAATGCAGGATCACTCATCATCTGACGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGtttt                                                                        <  7:1084316/151‑1 (MQ=255)
                                                                                  aaTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTc                                                                       >  8:1006895/1‑150 (MQ=255)
                                                                                   aTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTc                                                                       >  3:180077/1‑149 (MQ=255)
                                                                                     cTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcc                                                                   >  7:503908/1‑151 (MQ=255)
                                                                                      ggagcACCTGCCCATCGCATACTGCCGTATCACTCAGCATCTGACGCATATTGGGGGGCGGCGGCCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcca                                                                  <  3:1100775/147‑1 (MQ=255)
                                                                                       gaacACCTTGCCATACCCCACTCCAGGATCACTCATCATCTGAAGGATTTTGGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcca                                                                  <  2:532436/147‑1 (MQ=255)
                                                                                       gATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcc                                                                   >  8:1285029/1‑149 (MQ=255)
                                                                                         tCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCt                                                               >  1:1285506/1‑151 (MQ=255)
                                                                                         tCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCt                                                               >  6:111117/1‑151 (MQ=255)
                                                                                          cACCTTGCCATCGCCTAATGCAGGATCACTCATCTTCTGAAGTATTTTTGGGGGCGGAGGAGGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTg                                                              <  6:477046/151‑1 (MQ=255)
                                                                                             cTTGCCATAGAATAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCCGATAACTTTAGGAATAGCCGGATCATGATAACCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCAAGCTGACg                                                           >  4:62441/1‑151 (MQ=255)
                                                                                              ttGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACg                                                           >  7:30159/1‑150 (MQ=255)
                                                                                               tcccATAGCCTAATGCAGGAGCACTCATCAGCTTAAGGATATTTGGGGGCGGAGGAGGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGa                                                          <  4:27407/148‑1 (MQ=255)
                                                                                               acccATACCCTAAGCCAGGACCACTCAGCATCTGAAGGAGCTTGGGGGCCGGCGGCCGGGGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGAt                                                         <  7:382877/149‑1 (MQ=255)
                                                                                                 ccATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCATCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGATGACGATCg                                                       >  5:245553/1‑151 (MQ=255)
                                                                                                  cATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGt                                                      >  3:170311/1‑151 (MQ=255)
                                                                                                  cATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGt                                                      >  8:1132295/1‑151 (MQ=255)
                                                                                                   atccccTAATGCACGACCACTCATCATCTCAAGGATATTTGGGGGCGGAGGACGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTg                                                     <  6:773139/147‑1 (MQ=255)
                                                                                                   aTAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTg                                                     >  3:977043/1‑151 (MQ=255)
                                                                                                     accctactcccccaTCACTCCTCCTCTGACGGTTCTTTGGGGGCGGCGGCCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggc                                                   <  4:170226/138‑1 (MQ=255)
                                                                                                      gCCTAAGGCAGGATCACTCATCCTCAGAAGGATTTTTGGGGGCGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcg                                                  <  3:562168/151‑1 (MQ=255)
                                                                                                      cccTCATCCAGGAACACCCATCATCTGCAGGATTTGTGGGGGCGCCGGCCGCGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcg                                                  <  5:1128464/150‑1 (MQ=255)
                                                                                                      cccTAATGCCGGAGCCCTCATCATCTGAAGGATATTGGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcg                                                  <  3:1201602/150‑1 (MQ=255)
                                                                                                         tAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGgcggcg                                               >  7:820731/1‑151 (MQ=255)
                                                                                                         tAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGgcggcg                                               >  5:957655/1‑151 (MQ=255)
                                                                                                          actgcagccgccctcctcAGCGGAAGGAGATTGGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcggc                                                <  1:1209541/137‑1 (MQ=255)
                                                                                                            tGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGttt                                            >  6:353460/1‑151 (MQ=255)
                                                                                                             ccAGGACCACGCAACATCTGAAGGAGATGGGGGCGCCGCGCCCCGGGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTa                                           <  3:417575/150‑1 (MQ=255)
                                                                                                               acgATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTaaa                                         <  3:62074/149‑1 (MQ=255)
                                                                                                                ggATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAAc                                        >  1:1352915/1‑151 (MQ=255)
                                                                                                                  agcACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACAt                                      <  5:1281005/149‑1 (MQ=255)
                                                                                                                     ccgccgcatcTGCAGGATATTGGGGGGCGGAGGCCGCGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCa                                   <  5:548717/145‑1 (MQ=255)
                                                                                                                     aCTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATcc                                    >  8:144356/1‑150 (MQ=255)
                                                                                                                      cTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATcc                                    >  8:1239577/1‑149 (MQ=255)
                                                                                                                       acatcatCTGACGGCTCTTTGGCGGACGCCGACGGCCCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGc                                 <  2:7397/150‑1 (MQ=255)
                                                                                                                        catcttCTGAAGGATATTTGGGGGAGGAGGAGGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGc                                 <  4:982442/150‑1 (MQ=255)
                                                                                                                        cagcatcTGACGGATATTGGGGGGCGGAGGACGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAg                                  <  5:206797/146‑1 (MQ=255)
                                                                                                                        cagcatcTCAAGGATATGTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGcc                                <  5:1170152/148‑1 (MQ=255)
                                                                                                                         atcatcTTAATTATATTTGGGGGTGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCa                               <  6:568796/151‑1 (MQ=255)
                                                                                                                             tcttaagtttatTTGGGGGAGGAGGCGGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCa                               <  8:565679/138‑1 (MQ=255)
                                                                                                                              cTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCa                           >  1:1089285/1‑150 (MQ=255)
                                                                                                                               tGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGACAATCAAc                         >  8:736870/1‑151 (MQ=255)
                                                                                                                                    gATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGcac                    >  6:1264820/1‑151 (MQ=255)
                                                                                                                                       atTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCacaac                 >  1:1059331/1‑151 (MQ=255)
                                                                                                                                        tttGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCacaac                 >  7:1361137/1‑150 (MQ=255)
                                                                                                                                        attGGGGGAGGAGGACGGTGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTCGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCCATCAACAGCacaac                 >  4:1042393/2‑150 (MQ=255)
                                                                                                                                         ttGGGGGCGGAGGCCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCacaaca                <  4:57344/150‑1 (MQ=255)
                                                                                                                                         gtGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCacaaca                <  8:673759/149‑1 (MQ=255)
                                                                                                                                          tGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCAcaacaa               <  2:366688/150‑1 (MQ=255)
                                                                                                                                            ggggAGGAGGCCGGGGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATcc            <  2:1285831/151‑1 (MQ=255)
                                                                                                                                            ggcgaggaggaCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATcc            <  7:1106060/148‑1 (MQ=255)
                                                                                                                                                aggaggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAcgc        >  6:1165402/1‑151 (MQ=255)
                                                                                                                                                      aCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTNCAGANCGNNNNNNNTNNNNNTTNTCANCCNNCTNNNNATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGAAAcc  >  8:894080/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: