Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,448,5260AC27.9% 541.6 / 82.8 318intergenic (‑86/+70)malZ/proYalpha‑glycosidase/proline‑specific permease ProY
Reads supporting (aligned to +/- strand):  ref base A (170/57);  new base C (2/86);  total (173/145)
Fisher's exact test for biased strand distribution p-value = 1.07e-35
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 336 total aligned reads displayed.

AGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGAAA  >  DH1/3448376‑3448674
                                                                                                                                                      |                                                                                                                                                    
aGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggaggagga                                                                                                                                                      <  3:1327203/151‑1 (MQ=255)
 gCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGAc                                                                                                                                                     <  5:202596/151‑1 (MQ=255)
  cAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACg                                                                                                                                                    <  8:265276/151‑1 (MQ=255)
   aGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACgg                                                                                                                                                   <  7:90639/151‑1 (MQ=255)
   aGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACgg                                                                                                                                                   <  1:71820/151‑1 (MQ=255)
     ttGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACgggg                                                                                                                                                 >  4:415411/1‑151 (MQ=255)
     ttGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACggg                                                                                                                                                  <  2:537068/150‑1 (MQ=255)
      tGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCT‑‑AGGATATTTGGGGGAGGAGGACGGGGGCg                                                                                                                                              >  2:1071304/1‑151 (MQ=255)
       gATCTTTGCTTTGTTTAACAAATGGGGGCACGGGCAGGTGCCATGCATTTAACATCTTAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACggggg                                                                                                                                                <  2:876831/150‑1 (MQ=255)
        aTCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTTATTAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGc                                                                                                                                               <  8:107097/150‑1 (MQ=255)
         tgtttgctttgTTTAACAATTGGGGGCACCGGAAGTTTCAATGCATTTAACATCATAATTCCCCTGTGATTAAATTGTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCg                                                                                                                                              <  8:951375/148‑1 (MQ=255)
          ctttgctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCgg                                                                                                                                             >  2:1131578/1‑150 (MQ=255)
            ttgctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGAt                                                                                                                                           >  4:1147936/1‑150 (MQ=255)
              gctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATacc                                                                                                                                        >  6:355602/1‑149 (MQ=255)
               gtttgTTTAACAAAGGGGGTCAGCTTCAGGTGCCATGCATTTAACATCATAATTCCCTTTTGATTAAAATCTTATCACGTTGGCATATCCTAATGCTGGATCACTTATCATCTTTTGGATATTTGGGGGGGGAGGCCGGGGGCGGATAAc                                                                                                                                        <  7:528534/149‑1 (MQ=255)
               ctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACt                                                                                                                                       >  8:1445981/1‑151 (MQ=255)
                     ttAACAACTGGGGGCACCGGCAGGTGCAATTAATTTAACATCATAATTCCCCTGTTATGAAAATCTGATCACCTTGCCATACCCTAATGCAGGATCACTCATCATCTGAAGGTTATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCgg                                                                                                                                  <  1:230113/150‑1 (MQ=255)
                     taaaCAAAGGGGGGCACCGGCAGGTGCCATTCATTAAACATCATACTTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGACGGATATGTGGGGGCGGCGGACGGGGGCGGATAACTTTCgg                                                                                                                                  <  2:1412025/148‑1 (MQ=255)
                       aacacctGGGGCCCCCGGCAGTTTCCATTCATTTAACATAATTATTCCCCTGTTATCAAAATCTCATCACCTAGCCATACCCTAATGCAGGATCACTCATCATCTCAATTATTTTTGGGGGCGGCGGCCGGCGCCGGATAACTTTCGGAAt                                                                                                                               <  4:1409185/145‑1 (MQ=255)
                       aaCAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAAt                                                                                                                               >  5:500820/1‑151 (MQ=255)
                         aaaaTGGGGCCCCCGGCAGTGGCCATGCATTAAACATCATAATTCCCCTGTTATCAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATTTTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAg                                                                                                                             <  6:661099/150‑1 (MQ=255)
                          aaaTGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGc                                                                                                                            >  7:836195/1‑151 (MQ=255)
                           aaTGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATAACCTTGCCATAGCCTAATGAAGGATCACTCATCATCTGAAGGATATTTGGTGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGcc                                                                                                                           >  7:312148/1‑151 (MQ=255)
                            aTGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGTGCGGATAACTTTCGGAATAGCCg                                                                                                                          >  8:1049103/1‑151 (MQ=255)
                                gggCACCGGCAGGTTCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTTATCACCTTGCCATAGCCTAATGCAGTATCACTCATCATCTGAAGGTTATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATc                                                                                                                      <  5:565014/151‑1 (MQ=255)
                                gcgCACCGGCAGGTGCCAGTCATTTAACATCATAATTCCCCTGTGATGAAAATCAGACCACCTTGCCATAGCCTAATGCAGTATCACTAATCATCTGCAGGATATTTGGGGGCGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATc                                                                                                                      <  5:38/149‑1 (MQ=255)
                                 ggCCCCGGCATGTGCCATTCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGTATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCa                                                                                                                     <  1:1199520/151‑1 (MQ=255)
                                 ggCACGGGCAGGTGCCATTCATTTAAAATCATAATTCCCCTGTGATTAAAATCTGATCACCTTGCCATAGCATAATGCATGATCAATCATCATCTGAAGGATATTTGGGGGAGGAGGAGGGGGGCGGATAACTTTCGGAATAGCCGGATCa                                                                                                                     <  4:120740/151‑1 (MQ=255)
                                    aCCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAgg                                                                                                                   >  5:983208/1‑150 (MQ=255)
                                     ccGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGGCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGAt                                                                                                                 >  8:1364398/1‑151 (MQ=255)
                                         cAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcc                                                                                                             >  8:879229/1‑151 (MQ=255)
                                             tGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTcctcct                                                                                                         >  6:547497/1‑151 (MQ=255)
                                              gCCATGAAATTAACATCATAATTCCCCTTTGATCAAAATCTGATCCCCTTCCCATACCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTcctcct                                                                                                         <  1:562760/150‑1 (MQ=255)
                                               ccATGAATTTAACATCATAATTCCCCTGTCATTAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGATGTATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGa                                                                                                       <  4:781592/151‑1 (MQ=255)
                                                cATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGa                                                                                                       >  7:438557/1‑150 (MQ=255)
                                                   gcataaaaCATCATAATTCCCCTGTGATGAAAAACTGATCCCCTTGCCATAGCCAAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCg                                                                                                   <  1:888027/145‑1 (MQ=255)
                                                   gCATTTAACAACATAATACCCCTGTGATCAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCCCTCATCATCTGAAGGATATTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAAcc                                                                                                    <  2:329847/150‑1 (MQ=255)
                                                     aTTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtc                                                                                                              >  2:779775/1‑138 (MQ=255)
                                                       ttAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCtt                                                                                               >  1:1005087/1‑151 (MQ=255)
                                                       ttAACATCAAAATTCCCCTGTCATGAAAATCTCCCCACCTTCCCATCCCCTAATGCAGGATCACGCATCATCTTCAGGATATTTGGGGGAGGCGCCGGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCt                                                                                                <  3:7048/150‑1 (MQ=255)
                                                       aaaacaaaataaaaccccTGACAAGAAAAGCTGATCACCTTCCCATCCCCTAATGCAGCAGCACTCATCACCTCACGGATATTTGGGGCCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCtt                                                                                               <  1:170442/137‑1 (MQ=255)
                                                          aCATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAc                                                                                             >  7:27713/1‑150 (MQ=255)
                                                            atcatAATTCCCCTTTTATTAAAATCTGATCACCTTCCAATAGCCTAATGCAGTATCACTCATCTTCTGAATTATTTTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGa                                                                                          <  6:983245/151‑1 (MQ=255)
                                                            atcaaAATGCCCCTGTGATCAAAATCTGAGCACCTGGCCAAAGCCTACTGCAGCACCACTCACCAGCTCAACGATATTTGGGGGCGCAGCCCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGa                                                                                          <  3:1043897/151‑1 (MQ=255)
                                                              catAATTCCCCTGTGATGAAACTCTCAGCACCTTGCCATCCCCTCAGGCAGGCTCCCTCCTCATCTGAAGGATTTTTGGGGGCGGCGCCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAAc                                                                                        <  1:1325446/151‑1 (MQ=255)
                                                                  aTTCCCCTGTGATTAAAATCTGACCCCCTTGCCATAGCCTAATGCAGGATCACTCATCATCTTAAGGATTTTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCaa                                                                                    <  6:526607/151‑1 (MQ=255)
                                                                      cccTGTCATGAAAAGCTGATCAACATCCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGAt                                                                                <  1:641135/151‑1 (MQ=255)
                                                                       ccTGTGATCAAAACCTGATCACCTCCCCATAGCCTAATGCACGATCACTCATCATCTCAAGGATTTTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGAtt                                                                               <  8:715025/151‑1 (MQ=255)
                                                                         tgtgATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGACCGCAAGGATTAc                                                                             >  5:1072671/1‑151 (MQ=255)
                                                                            gatgaaccgcTGATCACCTTCCCATAACAAAATGCCGTATCTCTCCTCCTGTGTCGTTTTGTTGGGGGCGGGGGCCCCGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACt                                                                            <  4:1428165/140‑1 (MQ=255)
                                                                            gaggaAAATCGGAACACCTTCCCATACCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGCGGCCGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTgg                                                                          <  3:1406232/148‑1 (MQ=255)
                                                                             atgaAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGt                                                                         >  2:1158343/1‑151 (MQ=255)
                                                                             atgaAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGt                                                                         >  3:2878/1‑151 (MQ=255)
                                                                              tgaAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGACCCGGCTTACAGAACGCAAGGATTACTGGtt                                                                        >  3:124095/1‑151 (MQ=255)
                                                                              tgaAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGtt                                                                        >  8:523171/1‑151 (MQ=255)
                                                                              tgaAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGtt                                                                        >  4:1063481/1‑151 (MQ=255)
                                                                                aaaaTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGtttt                                                                      >  2:653540/1‑151 (MQ=255)
                                                                                  aaTCTCATCACCTGCCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTca                                                                    <  1:1381075/151‑1 (MQ=255)
                                                                                      ccaTCACCTTCCCAGACCCTACTGCCGGATCACTCATCATCTCAAGGATTTTTGGGGGAGGCGGCCGCGGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcca                                                                <  2:796789/149‑1 (MQ=255)
                                                                                       gaacACCTAGCCATAGCCTTATGCAGGATCACTAATCATCTGAAGGATTTTTGGGGGAGGCGGAGGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcca                                                                <  4:22922/147‑1 (MQ=255)
                                                                                       gATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcc                                                                 >  2:205325/1‑149 (MQ=255)
                                                                                         gcACCTTGCCATAGCCTAATCCAGGATCACTCAACATCTGAAGCATATTTGGGGGAGGAGGACGGAGGCGGATAACTTTCGGAATAGCCGGATAAGGATTTACTCCTGAAACCGCTTACAGAACGCAAGGATAACTGGTTTTcagccagc                                                              >  7:1189165/2‑150 (MQ=255)
                                                                                         gcACCTTCCCATAGCCTAAAGCAGGATCACTCATCATCTGAATGATTTTTGGGGGCGGCGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCt                                                             <  2:1312598/150‑1 (MQ=255)
                                                                                          cACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTg                                                            >  7:197536/1‑151 (MQ=255)
                                                                                           aCCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTgc                                                           >  5:1406107/1‑150 (MQ=255)
                                                                                                gCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATc                                                      >  2:15679/1‑151 (MQ=255)
                                                                                                gCCATAGCATAATGCAGGATCACTCATCATCTTAAGGATTTTTGGGGGCGGCCGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATc                                                      <  3:395017/151‑1 (MQ=255)
                                                                                                cccATCGCCTACTGCCGCACCCCTCATCATCTGAAGGATAGTTGGGGGCGGCGGCGGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATc                                                      <  1:397868/150‑1 (MQ=255)
                                                                                                  cATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGt                                                    >  6:482308/1‑151 (MQ=255)
                                                                                                    tAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTgg                                                  >  6:1019370/1‑151 (MQ=255)
                                                                                                      gCCTACTGCAGGACCACTAATCATCTGAAGCATATTTGGGGGCGGCGGCCCGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcg                                                <  5:121985/151‑1 (MQ=255)
                                                                                                      gCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcg                                                >  8:1155696/1‑151 (MQ=255)
                                                                                                        cTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGAGGATCGTggcggc                                              >  6:1309244/1‑151 (MQ=255)
                                                                                                          aaTGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGCTCGTGGCGGCGt                                            >  6:1183293/1‑151 (MQ=255)
                                                                                                           aTGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGtt                                           >  7:620978/1‑151 (MQ=255)
                                                                                                             gCAGGATCACTCATCATCTGAATGATTTTGGGGGGAGGAGGAGGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTa                                         <  4:518006/151‑1 (MQ=255)
                                                                                                                ggATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAAc                                      >  1:285229/1‑151 (MQ=255)
                                                                                                                   tCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATc                                   >  1:1039790/1‑151 (MQ=255)
                                                                                                                     ccTCAGCATCTGAAGGATTTTTGGGGGCGGCGGACGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCa                                 <  7:1419413/150‑1 (MQ=255)
                                                                                                                        catcatCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGcc                              >  5:1267202/1‑151 (MQ=255)
                                                                                                                            atcTGAAGGATATTTGGGGGAGGAGGACCGGGGCGGATAACTTTCGGAATAGCAGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTAGTGGTTTACAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATc                          >  2:975738/1‑151 (MQ=255)
                                                                                                                              cTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCa                         >  3:15784/1‑150 (MQ=255)
                                                                                                                               tGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAAc                       >  3:970035/1‑151 (MQ=255)
                                                                                                                                  aGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGTCGTTTAAACATCCAGCCAATCAACAGc                    >  5:749453/1‑151 (MQ=255)
                                                                                                                                  aGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGc                    >  6:157194/1‑151 (MQ=255)
                                                                                                                                   ggATTTTTGGGCGTGCAGGCCCGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGca                   <  2:284232/151‑1 (MQ=255)
                                                                                                                                     atatTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGcac                  >  8:478846/1‑150 (MQ=255)
                                                                                                                                     atatTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTCAACCGGCATACAGAACGCAAGGATTACTGCTTTTCAACCAGCTGACGATCGTGGCGGCGATTAAACATCCAGCCAATCAACATcac                  >  6:1408928/1‑150 (MQ=255)
                                                                                                                                       atTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCGGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCacaac               >  6:898804/1‑151 (MQ=255)
                                                                                                                                        tttGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCacaac               <  3:303273/150‑1 (MQ=255)
                                                                                                                                        tttGGGGGAGGAGGACGGGGGCGGATAACTTTGGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGCACGCAAGGATTACTGGTTTTCAGCCACCTGACGATCGTGGCGGCGTGTAAAGCTCCAGCCAATCAACAGCacaac               >  8:283331/1‑150 (MQ=255)
                                                                                                                                          tGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCAcaacaa             <  4:1300945/150‑1 (MQ=255)
                                                                                                                                           gggggCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACTAt            <  3:814645/150‑1 (MQ=255)
                                                                                                                                           gggggAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATc           >  4:605676/1‑151 (MQ=255)
                                                                                                                                             gggAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCa         <  8:1194572/151‑1 (MQ=255)
                                                                                                                                               gaggaggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCAGAACCGGCTTACAGCACGCAAGGCTTACTGGTTTTCAGCAAGCTGACGCTCGTGGAGGAGTTTAAACATCCACCCAATCAACACCACCACCATCCAcg       >  4:842424/1‑151 (MQ=255)
                                                                                                                                                  gaggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAcgcg     >  6:221620/1‑150 (MQ=255)
                                                                                                                                                  gaggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAcgcg     >  2:430055/1‑150 (MQ=255)
                                                                                                                                                  gaggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAcgcg     >  7:1292194/1‑150 (MQ=255)
                                                                                                                                                   aggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACACCACAACAATCCACGCGa    >  3:8879/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
AGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGCGAAA  >  DH1/3448376‑3448674

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: