Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,876,1040GC55.3% 199.0 / inf 458P77R (CCC→CGC) creCtwo‑component sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (114/90);  new base C (145/107);  total (260/198)
Fisher's exact test for biased strand distribution p-value = 7.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

Only 100 of 471 total aligned reads displayed.

AACGTTAGCCAGACGTCATTCCAGCGCGAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTG  >  DH1/3875954‑3876253
                                                                                                                                                      |                                                                                                                                                     
   gTTAGCCAGACGTCATTCCAGCGCGAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGAt                                                                                                                                                    <  2:851740/151‑1 (MQ=255)
    ttAGCCAGACGTCATTCCAGCGCGAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGAtg                                                                                                                                                   >  7:1044354/1‑151 (MQ=255)
      aGCCAGACGTCATTCCAGCGCGAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTTCCCTGAGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGACACAGATATTGGCGCGAAACGGGAGAtgtt                                                                                                                                                 >  7:159909/1‑151 (MQ=255)
      aGCCAGACGTCATTCCAGCGCGAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGAtgtt                                                                                                                                                 >  8:1484346/1‑151 (MQ=255)
       gCCAGACGTCATTCCAGCGCGAATAATCCTGTCCAACTGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGAtgttg                                                                                                                                                <  6:618732/151‑1 (MQ=255)
         cAGACGTCATTCCAGCGCGAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTa                                                                                                                                              <  1:497731/151‑1 (MQ=255)
         cAGACGTCATTCCAGCGCGAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTa                                                                                                                                              <  1:143167/151‑1 (MQ=255)
            aCGTCATTCCAGCGCGAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGc                                                                                                                                            <  4:1189392/150‑1 (MQ=255)
              gTCATTCCAGCGCGAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGa                                                                                                                                         <  5:1438792/151‑1 (MQ=255)
                cATTCCAGCGCGAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGAt                                                                                                                                        <  8:1065411/150‑1 (MQ=255)
                 aTTCCAGCGCGAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGAtt                                                                                                                                       >  7:1414111/1‑150 (MQ=255)
                  ttCCAGCGCGAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTa                                                                                                                                      >  3:534872/1‑150 (MQ=255)
                    ccAGCGCGAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAAc                                                                                                                                   <  6:228434/151‑1 (MQ=255)
                    ccAGCGCGAATAATCCAGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAAc                                                                                                                                   <  2:1766945/151‑1 (MQ=255)
                     cAGCGCGAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACg                                                                                                                                  >  2:325222/1‑151 (MQ=255)
                      aGCGCGAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGc                                                                                                                                 >  1:862933/1‑151 (MQ=255)
                       gcgcGAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGATGATTAAACGcc                                                                                                                                <  7:623630/151‑1 (MQ=255)
                         gcgAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTg                                                                                                                              >  7:313970/1‑151 (MQ=255)
                          cgAATAATGCTGTCCAACGGCTTTATTTGACGTATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGTGGATGTTGTAGCTGATTAAACGCCTgc                                                                                                                             <  3:1645005/151‑1 (MQ=255)
                          cgAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTgc                                                                                                                             >  3:632467/1‑151 (MQ=255)
                          cgAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTgc                                                                                                                             >  8:389508/1‑151 (MQ=255)
                           gAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTgcg                                                                                                                            >  6:86900/1‑151 (MQ=255)
                              taatCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGcc                                                                                                                          >  1:1009300/1‑150 (MQ=255)
                              taatCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGcc                                                                                                                          >  3:771702/1‑150 (MQ=255)
                                  ccTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTg                                                                                                                     >  6:1499073/1‑151 (MQ=255)
                                    tGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGcc                                                                                                                   >  5:523404/1‑151 (MQ=255)
                                     gTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGcc                                                                                                                   <  6:732470/150‑1 (MQ=255)
                                            ggCTTTATTTGCCGAATCGAACAATACTTTGCCGTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGtt                                                                                                           <  2:412522/151‑1 (MQ=255)
                                              cTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTgg                                                                                                         >  8:270654/1‑151 (MQ=255)
                                                ttATTTGCCGAATCGAACAATACTTTGCCCTGCGGATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTggg                                                                                                        <  4:927410/150‑1 (MQ=255)
                                                ttATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGt                                                                                                       >  6:1730615/1‑151 (MQ=255)
                                                 tATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTCCCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATCCGTTGGGTc                                                                                                      >  8:1815816/1‑151 (MQ=255)
                                                 tATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGAACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTc                                                                                                      >  5:749922/1‑151 (MQ=255)
                                                   tttGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTccc                                                                                                    >  1:680859/1‑151 (MQ=255)
                                                     tGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTcccc                                                                                                   >  5:1880625/1‑150 (MQ=255)
                                                      gCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTcccc                                                                                                   >  2:222389/1‑149 (MQ=255)
                                                       ccGAATCTAACAATACTTTGCCCTGCGCCTCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGCGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTCCCAATGAGTTGGGTcccc                                                                                                   >  4:1002352/1‑148 (MQ=255)
                                                       ccGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTcccc                                                                                                   >  6:1217320/1‑148 (MQ=255)
                                                        cGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCagag                                                                                               <  2:960369/151‑1 (MQ=255)
                                                        cGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCag                                                                                                 >  4:1359388/1‑149 (MQ=255)
                                                         gAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCagaga                                                                                              <  8:152534/151‑1 (MQ=255)
                                                            tCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGc                                                                                            >  6:1366170/1‑150 (MQ=255)
                                                             cGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCa                                                                                           <  6:167655/150‑1 (MQ=255)
                                                              gAACAATACTTTGCCCTGCGCATCGGTCATATAGCCATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCaa                                                                                          >  3:592906/1‑150 (MQ=255)
                                                               aaCAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAAt                                                                                        >  3:295864/1‑151 (MQ=255)
                                                                      atttGCCCTGCGGAGCGGTCATATTGACATGATATTGATTGCGCACTTTTTTAATGCCCCCGATATTTGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGTTCCCCAGAGAGCAAATCGGGAcg                                                                                 <  1:1743074/150‑1 (MQ=255)
                                                                        ttGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCAACGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGCActcc                                                                               >  4:382757/1‑148 (MQ=255)
                                                                         tGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGAcgcg                                                                               <  2:62686/150‑1 (MQ=255)
                                                                         tGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGAcgcg                                                                               <  8:934001/150‑1 (MQ=255)
                                                                         tGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACgcgc                                                                              >  5:990566/1‑151 (MQ=255)
                                                                           cccTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGccc                                                                            >  5:1985349/1‑150 (MQ=255)
                                                                            ccTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGATGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGcc                                                                             >  5:699195/1‑149 (MQ=255)
                                                                              tGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCt                                                                         >  2:765184/1‑151 (MQ=255)
                                                                               gcgcATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTc                                                                        >  5:118676/1‑151 (MQ=255)
                                                                                 gcATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTccc                                                                      >  1:169929/1‑150 (MQ=255)
                                                                                  cATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGc                                                                     >  1:1141242/1‑151 (MQ=255)
                                                                                  cATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGc                                                                     >  7:1224799/1‑151 (MQ=255)
                                                                                   aTCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGcc                                                                    <  3:1103817/151‑1 (MQ=255)
                                                                                         cATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAAc                                                              <  5:1807458/151‑1 (MQ=255)
                                                                                           tataGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGt                                                            >  3:331688/1‑151 (MQ=255)
                                                                                            ataGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGtt                                                           >  6:1813867/1‑151 (MQ=255)
                                                                                              aGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGc                                                         >  6:65331/1‑151 (MQ=255)
                                                                                                  aTGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGgtg                                                     <  3:1494600/151‑1 (MQ=255)
                                                                                                    gATATTCATTGAGCAGTTTGTTAATGCGAGCGATATTGGCGCGAAACGGTGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCTGATGCGTTGGGTCCTCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGgtgt                                                    <  2:1323120/150‑1 (MQ=255)
                                                                                                     atatTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCg                                                  >  6:422471/1‑151 (MQ=255)
                                                                                                      tatTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGCGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTCCGGTGTCg                                                  >  4:777973/1‑150 (MQ=255)
                                                                                                       atTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGAt                                                <  5:1058792/151‑1 (MQ=255)
                                                                                                            ttGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACg                                           >  5:891981/1‑151 (MQ=255)
                                                                                                             tGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAAGGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACg                                           <  2:1674271/150‑1 (MQ=255)
                                                                                                               cgcACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCCGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTcc                                        >  8:85743/1‑151 (MQ=255)
                                                                                                                gcACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATAGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTcac                                       >  4:662850/1‑149 (MQ=255)
                                                                                                                 cACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTcccc                                      >  2:1067143/1‑151 (MQ=255)
                                                                                                                  aCTTTGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCt                                     >  3:1906047/1‑151 (MQ=255)
                                                                                                                    tttGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTcc                                   >  4:508795/1‑151 (MQ=255)
                                                                                                                     ttGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTcc                                   <  1:1621845/150‑1 (MQ=255)
                                                                                                                      tGTTAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGt                                 >  2:1593900/1‑151 (MQ=255)
                                                                                                                        ttAATGCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACCTTGCGGTGTCGATCAACGTCCCCTCCGTTg                               >  4:1491934/1‑151 (MQ=255)
                                                                                                                          aaTGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGct                             >  3:896218/1‑151 (MQ=255)
                                                                                                                             gCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGctct                           <  6:1613618/150‑1 (MQ=255)
                                                                                                                             gCCACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGctct                           <  4:1747008/150‑1 (MQ=255)
                                                                                                                              ccACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTc                         >  3:1362496/1‑151 (MQ=255)
                                                                                                                              ccACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTc                         <  5:418668/151‑1 (MQ=255)
                                                                                                                               cACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCg                        >  5:1727085/1‑151 (MQ=255)
                                                                                                                               cACCGATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCg                        >  6:1162410/1‑151 (MQ=255)
                                                                                                                                 ccGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGcc                      >  1:1196967/1‑150 (MQ=255)
                                                                                                                                   gATATTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACg                    >  1:312990/1‑151 (MQ=255)
                                                                                                                                    atatTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGc                   >  5:1017320/1‑151 (MQ=255)
                                                                                                                                      atTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGccc                 >  7:210094/1‑151 (MQ=255)
                                                                                                                                      atTGGCGCGAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGccc                 <  1:1955542/151‑1 (MQ=255)
                                                                                                                                       ttGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCg                >  7:508639/1‑151 (MQ=255)
                                                                                                                                          gcgcGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGtt             <  1:1328156/151‑1 (MQ=255)
                                                                                                                                           cgcgAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGttt            <  7:1561847/151‑1 (MQ=255)
                                                                                                                                           cgcgAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGttt            >  7:1928652/1‑151 (MQ=255)
                                                                                                                                             cgAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTaa          <  4:29974/151‑1 (MQ=255)
                                                                                                                                              gacaCGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAAc         <  4:191443/148‑1 (MQ=255)
                                                                                                                                              gAAAGGGGGGATGTTGTAGCTGATTAAACGCGTGCGCCAGTTGCCCATGGGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACG‑GCCCTCCGTTGCTCTTCGCACGCCCGGTTTAAct        <  8:140694/151‑1 (MQ=255)
                                                                                                                                              gAAACGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAAc         >  4:1643569/1‑151 (MQ=255)
                                                                                                                                               aaaCGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCATTGCTCTTCGCACGCCCGGTTTAAct        <  2:1215296/151‑1 (MQ=255)
                                                                                                                                                 aCGGGCGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAActtc      >  7:797255/1‑151 (MQ=255)
                                                                                                                                                  cGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAActtct     >  7:1646414/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AACGTTAGCCAGACGTCATTCCAGCGCGAATAATCCTGTCCAACGGCTTTATTTGCCGAATCGAACAATACTTTGCCCTGCGCATCGGTCATATAGACATGATATTCATTGCGCACTTTGTTAATGCCACCGATATTGGCGCGAAACGGGGGATGTTGTAGCTGATTAAACGCCTGCGCCAGTTGCCCATGCGTTGGGTCCCCAGAGAGCAAATCGGGACGCGCCAGCTCCGCCAGCAACGTTGCGGTGTCGATCAACGTCCCCTCCGTTGCTCTTCGCACGCCCGGTTTAACTTCTTTG  >  DH1/3875954‑3876253

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