Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,281,5430AC38.8% 313.5 / 31.9 152V288G (GTG→GGG) cusAcation efflux system protein CusA
Reads supporting (aligned to +/- strand):  ref base A (79/14);  new base C (0/59);  total (79/73)
Fisher's exact test for biased strand distribution p-value = 4.52e-29
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 156 total aligned reads displayed.

AGGTTGTCGATGGCGCGGTCAATGAGCTGGCTGCGATCGTATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCTTTTAAAACGATGTGATTAAAGTCGTCGAGCGTT  >  DH1/3281393‑3281692
                                                                                                                                                       |                                                                                                                                                     
aGGTTGTCGATGGCGCGGTCAATGAGCTGGCTGCGATCGTATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATcacca                                                                                                                                                       <  2:972232/151‑1 (MQ=255)
  gTTGTCGATGGCGCGGTCAATGAGCTGGCTGCGATCGTATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCaccacc                                                                                                                                                     <  2:1282428/151‑1 (MQ=255)
    tGTCGATGGCGCGGTCAATGAGCTGGCTGCGATCGTATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCAccc                                                                                                                                                    <  2:556172/150‑1 (MQ=255)
     gTCGATGGCGCGGTCAATGAGCTGGCTGCGATCGTATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCAcccc                                                                                                                                                   <  1:1266350/150‑1 (MQ=255)
       cGATGGCGCGGTCAATGAGCTGGCTGCGATCGTATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACcccgcc                                                                                                                                                >  2:795098/1‑151 (MQ=255)
        gATGGCGCGGTCAATGAGCTGGCTGCGATCGTATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACcccgccc                                                                                                                                               >  1:275445/1‑151 (MQ=255)
             cgcgGTCAATGAGCTGGCTGCGATCGTATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCAc                                                                                                                                          >  2:1122820/1‑151 (MQ=255)
               cgGTCAATGAGCTGGCTGCGATCGTATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCActt                                                                                                                                        <  1:229026/151‑1 (MQ=255)
                    aaTGAGCTGGCTGCGATCGTATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCActtcgcct                                                                                                                                   >  1:1210686/1‑151 (MQ=255)
                     aTGAGCTGGCTGCGATCGTATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCActtcgcctt                                                                                                                                  >  1:661553/1‑151 (MQ=255)
                      tGAGCTGGCTGCGATCGTATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCActtcgccttc                                                                                                                                 >  2:1193126/1‑151 (MQ=255)
                       gAGCTGGCTGCGATCGTATTTTGTAACTATCTCCACGCCTTCCGGCAGCCTCCTTTTCAGCGTTTCCAGTTTGTCCTCCACGGCGGCGATCCCTTCTCGGGCGTTTTTGCGGGATCGCAGGATCACCCCCCCGCCCGCCActtcgccttc                                                                                                                                 <  2:275444/150‑1 (MQ=255)
                         gctggctgCGATCGTATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACttcgccttcgcc                                                                                                                              >  1:1818582/1‑151 (MQ=255)
                             gctgCGATCGTATGTTGTACCTATCTCCACGCCTTCGGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGTGGGCGATCACTTCTCGGGCGTTTTGGCCGGATCGCAGGATCACCCCCCCGCCCGCCACTTCGCCTTCGCCGttt                                                                                                                          <  2:661554/151‑1 (MQ=255)
                               tgCGATCGTATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGttt                                                                                                                          >  1:749039/1‑149 (MQ=255)
                                gCGATCGTATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGCTTCCAGTTTGTCCTTAACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCCCCGTTTAGt                                                                                                                       >  1:2013815/1‑151 (MQ=255)
                                     cGTATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGc                                                                                                                  >  2:2306863/1‑151 (MQ=255)
                                      gTATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGc                                                                                                                  >  1:2056837/1‑150 (MQ=255)
                                      gTATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGc                                                                                                                  >  1:308709/1‑150 (MQ=255)
                                       tATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCa                                                                                                                 >  1:500371/1‑150 (MQ=255)
                                       tATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCa                                                                                                                 >  2:1047668/1‑150 (MQ=255)
                                        aTGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAAt                                                                                                               >  2:1672526/1‑151 (MQ=255)
                                             gtAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGc                                                                                                                  >  1:2210940/1‑143 (MQ=255)
                                              tAACTATCTCCACGCCTTCCGGCAGACTATTTTTCAGCGTTTCCAGTTTGTCCTTCCCGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGcccc                                                                                                          <  2:2025129/150‑1 (MQ=255)
                                                  tATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGgcgc                                                                                                     >  2:584599/1‑151 (MQ=255)
                                                       ccACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCAtctc                                                                                                >  1:608841/1‑151 (MQ=255)
                                                        cccGCCTTCCGGCAGACTACTTTTCAGCTTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCTTTTTTCCCGGATCGCAGGATCCCCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTcc                                                                                               <  2:500371/149‑1 (MQ=255)
                                                         aCGCCTTCGGGCAGACTATTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCCCCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCg                                                                                              <  1:2306854/151‑1 (MQ=255)
                                                             cTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGcc                                                                                          >  1:898447/1‑151 (MQ=255)
                                                                ccGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAAt                                                                                       >  1:2178686/1‑151 (MQ=255)
                                                                   gCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTg                                                                                    >  1:1214074/1‑151 (MQ=255)
                                                                    cAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTgg                                                                                   >  1:801308/1‑151 (MQ=255)
                                                                    cAGACTACTTTTCAGCGTTTCAAGTTTGTCCTTCACGGCGCCGACCACTCCTCGGGCGTTTTTGCCGGACCGAAGGATCACCCCCCCGCCCGCCACTTCGCCTTCGCCGGTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTgg                                                                                   <  2:898448/151‑1 (MQ=255)
                                                                      gACTACTTTTCAGCTTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCTTTTTTGCCGGATCGCAGGATCCCCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCa                                                                                         <  2:190636/143‑1 (MQ=255)
                                                                      gACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCa                                                                                         >  1:190637/1‑143 (MQ=255)
                                                                       actcctTTTCAGCTTTTCCAGTTTGTCCTTCCCGGCGGCAATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCCCCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGAcc                                                                                <  2:1355792/151‑1 (MQ=255)
                                                                        tttctTTCCAGCGTTTCCAGTTTGTCCTTCGCGGCGGCGATCACTTCTCGGGCGTTTTGGCGGGATCGCGGGATCACCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCt                                                                                     <  2:551669/142‑1 (MQ=255)
                                                                        ctactTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGATTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCt                                                                                     >  1:551668/1‑145 (MQ=255)
                                                                         tactTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCtt                                                                              >  2:1962574/1‑151 (MQ=255)
                                                                          tctTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCCCTTCTCGGGCGTTTTTGCCGGATCGCAGGATCCCCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCtt                                                                              <  2:801309/149‑1 (MQ=255)
                                                                          actTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTc                                                                             >  1:1369082/1‑151 (MQ=255)
                                                                           ctTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCg                                                                            >  1:293534/1‑151 (MQ=255)
                                                                             tttCAGCTTTTCCAGTTTGCCCTTCCCGGCGGCGATCACTTCTGGGGCGTTTTTGCCGGATCGCAGGATCCCCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTc                                                                             <  2:2178694/148‑1 (MQ=255)
                                                                                cAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCTGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAAc                                                                        >  1:2622056/1‑150 (MQ=255)
                                                                                 aGCGTTTCGAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGAGCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTTGTTCGGCAATGCCCCCCCGCATCTCCGGGCCAATCTGGACCTTCGCAAcc                                                                       >  1:2136029/1‑149 (MQ=255)
                                                                                     tttCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATcgcg                                                                  >  1:2345904/1‑151 (MQ=255)
                                                                                         cAGTTTGTCCTTCACGGCGGCGATCCCTTCTCGGGCTTTTTTGCCGGATCGCAGGATCCCCCCCCCCCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGa                                                              <  1:1962568/151‑1 (MQ=255)
                                                                                          aGTTTTTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGACCGCAGAATCCCCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGAt                                                             <  2:293533/151‑1 (MQ=255)
                                                                                           gTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGAt                                                             >  2:2247261/1‑150 (MQ=255)
                                                                                           gTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGAt                                                             >  2:157632/1‑150 (MQ=255)
                                                                                            tttGTCCTTCACGGCGGCGATCCCTTCGCGGGCGTTTTGGCCGGATCGCAGGACCACCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCg                                                                            <  1:958790/134‑1 (MQ=255)
                                                                                            tttGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATa                                                            >  2:2630152/1‑150 (MQ=255)
                                                                                              ggTCCTTCACGGCGGCGATCACTTCTCGGGCTTTTTTGCCGGATCGCAGGATCCCCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATaaa                                                          <  2:1166332/149‑1 (MQ=255)
                                                                                               gTCCTTCCCGGCGGCGCTCCCTTCTCGGGGGTTTTTGCCGGACCGCAGGATCCCCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGAt                                                             <  1:581560/146‑1 (MQ=255)
                                                                                                 ccTTCCCGGGGGCGATCCCTTCCCGGGCTTTTTTGCCGGATCGCAGGATCACCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACgg                                                       <  2:608842/150‑1 (MQ=255)
                                                                                                 ccTTCACGGCGGCGATCCCTTCTCGGGCTTTTTTGCCGGATCGCAGGATCACCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACggg                                                      <  2:1818587/151‑1 (MQ=255)
                                                                                                   ttCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGc                                                     >  2:1978674/1‑150 (MQ=255)
                                                                                                    ccACGGCGGCGATCTCTTCCCGGGTGTTTTGGCCGGAGCGCATGACCACCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCAc                                                   <  2:219310/150‑1 (MQ=255)
                                                                                                    acACGGCGGCGCTCACTTCTCGGGCGTTTTGGCGGGATCGCAGGACCCCCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCa                                                    <  1:621434/149‑1 (MQ=255)
                                                                                                     cagggcggcGATCACTTCTCGGGCGTTTTTCCCGGATCGCAGGACCCCCCCCCCGCCCGCCACTTCGCCTTCACCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACg                                                  <  1:2630136/148‑1 (MQ=255)
                                                                                                     cACGGCGGCGTTCACTTCTGGGGCTTTTTTGCCTGATCGCAGGATCACCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATACCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACg                                                  <  2:878502/151‑1 (MQ=255)
                                                                                                     cACGGCGGCGATCCCTTCTCGGGCTTTTTTGCCGGATCGCAGGTCACCCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACg                                                  <  1:2277694/151‑1 (MQ=255)
                                                                                                     cACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCCCCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCAc                                                   <  2:2345914/150‑1 (MQ=255)
                                                                                                     cACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACg                                                  >  2:5173/1‑151 (MQ=255)
                                                                                                     cACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACg                                                  >  1:736919/1‑151 (MQ=255)
                                                                                                       cggcggCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGcc                                                <  1:2510461/151‑1 (MQ=255)
                                                                                                       ccatccgatccattcctcGGCAGTTTTTGCCGGTTCGCAGGCCCCCCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGcc                                                <  1:776705/136‑1 (MQ=255)
                                                                                                        ggcggcGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGcc                                                >  1:1243944/1‑150 (MQ=255)
                                                                                                           ggcgATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCAttt                                            >  2:802254/1‑151 (MQ=255)
                                                                                                            gcaaTCCCTTCTCGGGCGTTTTTGCCGGATCGCGGGATCCCCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCAtttt                                           <  2:668112/148‑1 (MQ=255)
                                                                                                             cgATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGcc                                                >  2:776706/1‑145 (MQ=255)
                                                                                                               aTCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCAc                                        >  2:243769/1‑151 (MQ=255)
                                                                                                                  ccTTCTCGGGCTTTTTTGCCGGATCCCGGAACCCCCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTgg                                     <  1:1115033/150‑1 (MQ=255)
                                                                                                                  aCTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTgg                                     >  2:342733/1‑151 (MQ=255)
                                                                                                                  aCTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTgg                                     >  2:1115034/1‑151 (MQ=255)
                                                                                                                   cTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACt                                       <  1:1030271/148‑1 (MQ=255)
                                                                                                                   cTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACt                                       >  2:1030272/1‑148 (MQ=255)
                                                                                                                    ttCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCt                                   >  2:2140135/1‑151 (MQ=255)
                                                                                                                       tcGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCtttt                                >  1:1426243/1‑151 (MQ=255)
                                                                                                                       attctcGTTTTCGCGGGACCGAAGGATCCCCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCtttt                                <  2:1550298/146‑1 (MQ=255)
                                                                                                                        cGGGCTTTTTTGCCGGATCGAAGGACCACCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCttt                                 <  1:157632/149‑1 (MQ=255)
                                                                                                                        cGGGCGTTTTTGCCGGATCGCAGGATCCCCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCttt                                 <  2:1214076/149‑1 (MQ=255)
                                                                                                                        cGGGCGTTTTTGCCGGATCGCAGGACCCCCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCtttt                                <  2:736920/150‑1 (MQ=255)
                                                                                                                        cGGGCGTTTTTGCCGGATCGCAGGACCACCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCtt                                  <  1:1029961/148‑1 (MQ=255)
                                                                                                                         gggCTTTTTTGCCGGATCGCAGGATCCCCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCTTTTaa                              <  2:1426248/151‑1 (MQ=255)
                                                                                                                             gTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCTTTTAAAACg                          >  2:1992054/1‑151 (MQ=255)
                                                                                                                              tctttGCCGGATCGTAGGATCACCACACCGCCCTCTACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCATTCTGGACCTTCGCAACATCGCGCAGATAAACGGCCACGCCATTTTCACTGGCTTTTAAAACg                          >  1:2485538/3‑150 (MQ=255)
                                                                                                                               ttttGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCTTTTAAAACGAt                        >  1:1648496/1‑151 (MQ=255)
                                                                                                                                tttCCCGGTCCGCAGGACCCCCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCTTTTAAAACGAt                        <  1:675446/150‑1 (MQ=255)
                                                                                                                                 ttGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCTTTTAAAACGAt                        >  1:1112700/1‑149 (MQ=255)
                                                                                                                                  tGCCGAATCGCAGGACCCCCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCTTTTAAAACGAtgt                      <  1:1472447/150‑1 (MQ=255)
                                                                                                                                    ccGGACCGCAGAGTCCCCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCTTTTAAAACGATGTGAt                   <  2:2136037/151‑1 (MQ=255)
                                                                                                                                        cgaccAGGTTCCCCACCCGCCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCGACTCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCTTTTAAAACGATGTGATTaaa               <  2:2485550/147‑1 (MQ=255)
                                                                                                                                        aTCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCTTTTAAAACGATGTGATTaa                >  1:876683/1‑150 (MQ=255)
                                                                                                                                          cGCAGGACCACCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCTTTTAAAACGATGTGATTAAAgt             <  1:2157678/151‑1 (MQ=255)
                                                                                                                                           gCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCTTTTAAAACGATGTGATTAAAgtc            >  1:1673295/1‑151 (MQ=255)
                                                                                                                                            gagaaTCACCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGTCATTTTCACTGGCTTTTAAAACGATGTGATTAAAgtcg           <  1:802253/147‑1 (MQ=255)
                                                                                                                                             aGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCTTTTAAAACGATGTGATTAAAgtcgt          >  2:2633788/1‑151 (MQ=255)
                                                                                                                                              ggATCACCCCCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCTTTTAAAACGATGTGATTAAAgtcgtc         <  1:5173/151‑1 (MQ=255)
                                                                                                                                                      caccCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCTTTTAAAACGATGTGATTAAAGTCGTCGAGCGtt  >  2:1415804/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
AGGTTGTCGATGGCGCGGTCAATGAGCTGGCTGCGATCGTATGTTGTAACTATCTCCACGCCTTCCGGCAGACTACTTTTCAGCGTTTCCAGTTTGTCCTTCACGGCGGCGATCACTTCTCGGGCGTTTTTGCCGGATCGCAGGATCACCACCCCGCCCGCCACTTCGCCTTCGCCGTTTAGTTCGGCAATGCCCCGGCGCATCTCCGGGCCAATCTGGACCTTCGCAACATCGCGCAGATAAACGGGCACGCCATTTTCACTGGCTTTTAAAACGATGTGATTAAAGTCGTCGAGCGTT  >  DH1/3281393‑3281692

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: