Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,169,5480GT28.0% 472.1 / 79.1 246P671Q (CCG→CAG) rsxCelectron transporter RsxC
Reads supporting (aligned to +/- strand):  ref base G (26/149);  new base T (17/51);  total (46/200)
Fisher's exact test for biased strand distribution p-value = 9.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

GCGGCAACTGCCGCTTTGCGCGGATCAACCTGTTCTTCTGGTACCGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGG  >  DH1/2169398‑2169698
                                                                                                                                                      |                                                                                                                                                      
gcgGCAACTGCCGCTTTGCGCTGATCAACCTGTTCTTCTGGTACCGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCg                                                                                                                                                        <  7:203145/151‑1 (MQ=255)
gcgGCAACTGCCGCTTTGCGCGGATCAACCTGTTCTTCTGGTACCGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCg                                                                                                                                                        <  5:474361/151‑1 (MQ=255)
 gggCTACTGGCGCTTTGTGATGATCAAGTTGTTGTTCTGGTACGTTTTTATCCTGTTGGTGTTCGTGTTTGATTGCTTTGGCGCGTGGAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCgg                                                                                                                                                       <  2:455169/150‑1 (MQ=19)
  ggCAACTGCCGCTTTGCGCGGATCAACCTGTTCTTCTGGTACCGCATTAGCCTGTTGCTGTTACAGCTTAAGAGCTGTGGACAGTGCTATAGGCGCTTGGAGGGCTGCTTTGCGCGGATCTCGCGGTTATTATGGCTCCGCATTAGCCGGt                                                                                                                                                      >  1:1479634/1‑151 (MQ=18)
    cAACTGCCGCTTTGGGGGGATCATCCTGTTGTTCTGGTACCTCATTAGCCTGTTGCTGTTCGAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCTGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTg                                                                                                                                                    <  7:1252383/151‑1 (MQ=255)
      aCTGCCGCTTTGCGCGGATCAACCTGTTCTTCTGGTACCGCATTAGCCGGTTGCTGTTCCAGATTGAGTGCTTAGGCGCGTGGAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCGACCTGTTCTTCTGGCTGGGCATTCGCCGGtgggg                                                                                                                                                  >  7:1534371/1‑147 (MQ=17)
         gCCGCTTTGCGCGGATCAACCTGTTCTTCTGGTACCGCATTAGCCTGTTGCTGTTACAGCTTGCGAGATTTGGCGCGTGCAATAGCAGCTTCGACGGCTGCTTTGATCGGATCAACAAGTTCTTCTGGGTCCGCATTAGCGCGTTACTggg                                                                                                                                               >  8:1329224/1‑149 (MQ=17)
         gCCGCTTTGCGCGGAGCAACCTGTTCTTCTGGTACCGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCTCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGt                                                                                                                                                <  5:417433/150‑1 (MQ=255)
         gCCGCTGTGCGCGGATCCACCTGCTCTTCTGGTACCGCATTAGCCTGTTGCTGTGCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGtt                                                                                                                                               <  3:927883/151‑1 (MQ=255)
           cgcgTTGCGCGGATCAACCTGTTCTTCTGGTACCGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTcc                                                                                                                                             <  2:141280/151‑1 (MQ=255)
                 gcgcGGATCAACCTGTTCTTCTGGTACCGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTg                                                                                                                                       <  8:725297/151‑1 (MQ=255)
                 gcgcGGATCAACCTGATCTTCTGGTACCGCATTAGCCTGTTGCTGGTCCAGCTTGCGGGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGttt                                                                                                                                        <  5:1620453/150‑1 (MQ=255)
                 gcgcCGAGCAACCTGCTCTTCTGGTACCGCATTACCCTGTTGCTGTTCCAGCTTGCGGGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTCTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGttt                                                                                                                                        <  8:1105162/150‑1 (MQ=255)
                  cgcgGATCAACCTGGTCTTCTGGTACCGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCGTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTgc                                                                                                                                      <  7:360054/151‑1 (MQ=255)
                           aCCTGTTCTTCTGGTACCGGATTAGCCTGTTGCTGTTCGTGCTTTCGTGCGTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTgg                                                                                                                             <  5:1041984/151‑1 (MQ=255)
                            ccTGTTCTTCTGGTACCGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCAGCTAGGACGGCAGCTATTCGCGGATCCAGAGGGTATGATCGATCGGCATTAGCCTGTTGCTGTTCCAGATAGCGAGCTTTggg                                                                                                                            >  7:1319846/1‑150 (MQ=2)
                               gTTCTTCTGGTACCGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCt                                                                                                                           <  3:774844/149‑1 (MQ=255)
                               gTTCTTCTGGTACCGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCg                                                                                                                         <  8:826588/151‑1 (MQ=255)
                               gTTCTTCTGGTACCGCATTAGCCTGTTGCTGGTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTc                                                                                                                          <  3:1075693/150‑1 (MQ=255)
                               gTTCTTCTGGTACCGCATCCTCCTGGTGCTGTTCGTGCTTGGTTGCTTTGGCGCGTGCAATAGCAGCGTCGACGGGAGCTTTGCGCGGTTCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGGGCGCTTTGGCtgg                                                                                                                         <  2:607891/151‑3 (MQ=255)
                               gTTCTTCTGCGACGGCATTAGCCTGGTGCTGGTCCAGCTTGCGGGCTTTGGCGCGTGCAATAGCCGCTGCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCt                                                                                                                           <  6:988681/149‑1 (MQ=25)
                                gtcttctGGCACCGCATGAGCGTGCTGCTGTGCCTGCTTGCGTGCGTTGGCGCGTGCAATAGCCGCGTCCACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGCCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCgg                                                                                                                        <  5:1445505/150‑1 (MQ=25)
                                 tcttctGGTACCGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTCGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCgg                                                                                                                        <  8:631817/150‑1 (MQ=255)
                                  cttctGGTACCGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGACGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGc                                                                                                                      >  7:64812/1‑151 (MQ=255)
                                    tctGGTACCGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCa                                                                                                                     <  3:597252/150‑1 (MQ=255)
                                    tctGGTACCGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCa                                                                                                                     <  3:1717194/150‑1 (MQ=255)
                                    tctGGGACCGGATTAGCCTGTTGCTGGTCCCGCTTGCGGGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCa                                                                                                                     <  1:545899/150‑1 (MQ=255)
                                       ggcaCCGCTTTAGCCTGTTGCTGTTCGTGGTTTCTTGCTTTTGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAg                                                                                                                 <  6:1103197/148‑1 (MQ=25)
                                       ggTACCTCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCTACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCaa                                                                                                                    <  4:764141/148‑1 (MQ=255)
                                       ggTACCGCATTAGCCTGTTGCTGGTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCaa                                                                                                                    <  2:1493773/148‑1 (MQ=255)
                                         gtccGCTTTAGCCTGGTGCTGTTCGCGCTTTGGGGCTTTGGCGCGGGCAATAGCCGCTGCTACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCg                                                                                                               <  2:233587/149‑1 (MQ=18)
                                            ggCATTAGCCTGTTGCTGTTCGTGCTTGGGGGCTTTGGCGCGTGCAATAGGCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCt                                                                                                            <  2:1551029/150‑1 (MQ=255)
                                             ggaTTAGGCTGTTGCTGTTCGAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTGCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCtt                                                                                                           <  3:364417/149‑1 (MQ=255)
                                                 tAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCCACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGATTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGAc                                                                                                       >  3:1831883/1‑151 (MQ=255)
                                                 acgCCTGTTGCTGGTCCAGCTTGCGGGCGTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTCTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGAc                                                                                                       <  1:1692313/149‑1 (MQ=255)
                                                  aGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGAcg                                                                                                      <  4:244981/151‑1 (MQ=255)
                                                  aGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGAcg                                                                                                      <  5:1338449/151‑1 (MQ=255)
                                                   gCCTGGTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGAcgg                                                                                                     <  3:481578/151‑1 (MQ=255)
                                                   cccTGGTGCTGGTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGAcgg                                                                                                     <  1:677689/150‑1 (MQ=255)
                                                    ccTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGAcggc                                                                                                    <  4:1523803/151‑1 (MQ=255)
                                                    ccTGTTGCTGGTCCTGCTTTATTGCTTTGGCGCGTGCAATAGCCGCGTCTACTGCAGCTTTGGGCGGATCAACCTGTTCTTCTGGATCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGAcggc                                                                                                    <  4:931434/151‑1 (MQ=17)
                                                      tGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGATTTACGAGGATCACCAAGTGGTTCTGGGTCCTGATTTGCCTGTTGCTGTTCCAGTTTGCGGGGTTTGTCAAGGGAAATAGCTCCTTCAAcggc                                                                                                    >  4:362820/1‑149 (MQ=2)
                                                       gTTGCTGTTCCCGGTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCTACGGCCGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGAcggcgg                                                                                                  <  3:1368963/150‑1 (MQ=255)
                                                       gTTGCTGTTCCCGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGAcggcg                                                                                                   <  5:369339/149‑1 (MQ=255)
                                                        gtGCTGGTCGAGCTTGCGTGCGTTCGCGCGGGCAATAGCCGCGGCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCt                                                                                                <  5:1020288/150‑1 (MQ=17)
                                                         tGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAACATCTTAGGCGCCTCCTTTCCGCCAATCCACGTGGTCTTATGGGTCCGCATTAGCCTGTGGCTGTGAAATATTGAGAGCTTTGGCTCGGGCAATGGGGGcggcggcggggggtg                                                                                               >  6:1364318/1‑137 (MQ=2)
                                                             gtgcctcgttGCGGGCTTTGGCGCGGGCAATAGCCGCTGCGACTGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTgc                                                                                            <  6:1735895/142‑1 (MQ=17)
                                                             gTTCCACCTTGCGTGCTTTCGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTgc                                                                                            <  1:94014/150‑1 (MQ=255)
                                                                ccAGCTTGCGTGCTTTTGCGCGTGCAATAGCCGCTTCTACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGcgcg                                                                                         <  6:216944/150‑1 (MQ=255)
                                                                 cAGGTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGa                                                                                       <  7:244307/151‑1 (MQ=255)
                                                                  cgCTTGCGGGCGTTGGCGCGGGCAATAGCCGCTGCGACCGCCGCTGTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGAt                                                                                      <  1:1363723/150‑1 (MQ=14)
                                                                   gcttgcttgcttTGGCGGGTGCAATAGCCGCTTCGTCGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATc                                                                                     <  3:1278705/151‑1 (MQ=255)
                                                                   gCTTGCGTGCGTTGGCGCGTGCAATAGCCGCTGCGACTGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATc                                                                                     <  7:965711/151‑1 (MQ=255)
                                                                   gCTTGCGGGCTTGGGCGCGGGCAATAGCCGCTTCTACTGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATc                                                                                     <  6:619316/151‑1 (MQ=255)
                                                                    cTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGAGCGGATCAACCGGTTCTTCTCGATCCGCATTAGCCGGTTGCTGTTCCCGTTGGCGCGCTTTTTCTCGGGCAATAGCGGCTTCGCCGGCGGCTTGGCGGGGCTCg                                                                                    >  1:1232313/1‑151 (MQ=255)
                                                                     ttGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGAAGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCAGATTGCGCGCTTTGGCTCGGGCCATAGCGGCTGGGACGGCGGCTGGGCGGGGCTCgc                                                                                   >  3:247151/1‑150 (MQ=21)
                                                                        cccccgcccccccgcgCAATACCCGCGGCCACGGCGTCTTTGCGCGGATCAACCGGGGCTTCTGGATCGGCCTATGCCTGTTGCTGTCCCAGCTTGCGGGCTCTGCCCCGGGCAATAGCGGCTTCTACGGCGGCTTTGCGCGGATCGACCt                                                                                <  7:1207792/136‑1 (MQ=0)
                                                                        cGTGCTTGCGCGCGTGCAATAGCCGCTTCGACTGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGGGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGAcc                                                                                 <  2:1127736/150‑1 (MQ=255)
                                                                         gggcctgggCGCGTGCAATAGCCGCTGCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCtg                                                                               <  6:700122/144‑1 (MQ=255)
                                                                           gcgtttgcgcGTGCAATAGCCGCTGCGACGGCAGCTGTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCtgtt                                                                             <  3:1123139/145‑1 (MQ=255)
                                                                           gcgttGGCGCGTGCAATAGCCGCTTCGACTGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCtgtt                                                                             <  5:795801/148‑1 (MQ=255)
                                                                            cTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCGAACAGGTGTTCTGGCTCCGCATTAGCCAGTTGCTGTTCCGGTTTGCGCGGGTTGGCTCGGGCAATAGCGTCTTCGACGGCGGCTTTGCGCGGATCGAcgggggg                                                                            >  4:559435/1‑145 (MQ=14)
                                                                             gttGGCGCGGGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGTGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCtgttgt                                                                           <  3:678381/150‑1 (MQ=255)
                                                                                     gTGCAATAGCCGCTTCGACTGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTc                                                                   <  7:706037/151‑1 (MQ=255)
                                                                                       gCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCTGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCgg                                                                 <  5:570653/151‑1 (MQ=255)
                                                                                       gCAATAGCCGCTGCAACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGGGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCgg                                                                 <  4:163668/151‑1 (MQ=255)
                                                                                        cAATAGCCTCTTCGACTGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCTGTTGCTGTTCCCGTTTGCGGGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGGc                                                                <  3:80399/151‑1 (MQ=255)
                                                                                        cAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCGACCTGTTCTTCTGGCTCCGCATTAGCCTGTTGCTGTTCCCGGTTGCGCGCTTTGGCGCGGGCAATAGATGCTTCAACGGCAGATTTGCGGGGATAGAGCGGTTGTGCTGCGTCggg                                                                >  2:826472/1‑150 (MQ=11)
                                                                                         aaTAGCCGCTTCGACGGCAGCGGTGCGCGGAGCCACCTGCTCTTCTCGCTCCGCATTAGCCGGTCGCGGCGCCCGTCTGCGCGCTTTGGCGCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGGCa                                                               <  6:382367/151‑1 (MQ=17)
                                                                                            aGCGGCGTCGACTGGAGCGTTTGGCGGATCAACCTGGTCTTCTGGCTCCGGATTAGCCTTTTGCTGTTCCCGTTTTCGGGCCTTGGCGCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGGCAtt                                                             <  5:828093/150‑1 (MQ=14)
                                                                                             gCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCTGTTGCTGTTCCCGGTTGCGGGCTTTAGCTCGGACGAGAGGGGGTTAGGCGGCGGCTTTGCGCGGATCTACATGTGGTTCTGGCTCGGCAtt                                                             >  2:1634733/1‑149 (MQ=14)
                                                                                              ccGCGTCTACTGCCGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCGCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGGCATTAGc                                                          <  6:1201010/151‑1 (MQ=255)
                                                                                               cgcggcaactgCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCTGTTGCTGTTCCCGTTTGCGTGCTTTGGCGCGGGCAATAGCCGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGGCATTAGcc                                                         <  1:284461/141‑1 (MQ=11)
                                                                                                  gtccaccgccgcgttGCGCGGAGCAACCTGTTCTTCTGGCTCCGCATTAGCCGGGTGCTGTTCCCGGTTGCGCGCGTTGGCGCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGGCATTAGCCTGt                                                      <  8:1208474/138‑1 (MQ=14)
                                                                                                  ggcGACGGCAGCTTTGCGCGGAGCAACCTGTTCTTCTGGCTCCGCATTAGCCGGGTGCTGTTCCCGTTTGCGCGCTTTGGCGCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGGCATTAGCCTGt                                                      <  3:987775/149‑1 (MQ=255)
                                                                                                      aCGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGc                                                   <  6:841374/150‑1 (MQ=255)
                                                                                                         gccgagctGCGCGTAGCCACCTCGTCTTCGGCCTCCGCAACAGCCTGGTGCTGGGCGCGTTTGCTCGCGTTCGCGCGGGCAATAGCCGCTTCCACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGtt                                               <  7:765789/144‑1 (MQ=9)
                                                                                                         gccgCTTTGCTCTTATTAAACTGTTCTTCTGCCCCCGCATTAGCCGTTTTCTCTTCCCGTTTGTGCGCTTTCGGTCGGGCTATAGAGGCTTCTACTGCCGCTTTGCGCGGATCCACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGtt                                               <  2:1487917/148‑1 (MQ=255)
                                                                                                          cagcttttgtatgATGAACCTGTTCTTCTGGGACCGCATTAAACGGTTGGTGTGGCCGGTTGCGTGCTGTGGCGCTGTCAATAGCCGCTTCTACTGCCGCTTTGCGCGGATCCACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGTTc                                              <  1:413026/139‑1 (MQ=9)
                                                                                                           cgCTGTGCGCGGAGCAACCTGTTCTTCTGGCTCCGCATTAGCCTGTTGCTGTTCCGGTTTTCGCGCTTTGGCGCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGTTcc                                             <  3:540248/150‑1 (MQ=14)
                                                                                                               ttggggttagcAACCTGTTCTTCTGGCTCCGCATTATCCTGTTGCTGTTCCCGTTTTTTGGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCTACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGTTCCAGc                                          <  4:482355/140‑1 (MQ=14)
                                                                                                               ttGCGCTGATCAACCTGTTATTCTGGCTGCCCATTAGCCGGTTGCTGTTCCGCTCTGCGGGCGTGCGCGCGGGCAATAGCGGCTTCTACTGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGTTCCAg                                           <  3:358799/149‑1 (MQ=14)
                                                                                                                tGCTCGGAGCAATCTGTTCTTCTGGTTCCGCATTAGCCGGTTGCTGTTCCCGTTTGTTGGCTTTTGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCTACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGTTCCAGCt                                         <  2:636047/150‑1 (MQ=14)
                                                                                                                 gcgcTGAGCAACCTGATGTTCTGGCTCCGGATTAGCCTGGTGCTGTTCGCGTTTGCGGGCTTTGGCGCGGGCAATAGCGGCTGCTACTGCGGCTTTGCGCGGATCTACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGTTCCAGCTTg                                       <  8:1782356/151‑1 (MQ=14)
                                                                                                                  cgcgGATCCACGTGGTCTTCTGGCTCCTCAATAGCCTGCTGCTGTTCGCGTTTGCTCGCTTGGGCTCGGGCAATAGCGGCTTTGACGGCCGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGTTCCAGCTTGc                                      <  1:181107/151‑1 (MQ=14)
                                                                                                                       aTCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACATGTGGTTCTGGCTCGGCATTTGCCTGATGCTGTTACAGCTTGAGCGCt                                 >  7:1635600/1‑151 (MQ=14)
                                                                                                                        gctacttcttcttcTGGCTCCGCATTAGCCTGTTGCTGTTCGGGTTTGGTGGCTTTGGCGCGGGCAATAGCGGCTTCTACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCt                                 <  2:921706/142‑1 (MQ=14)
                                                                                                                         cAACCTGTTCTGCTGGCTCCGCATTAGGCTGTTTCTGATGCCGGTTGCGCGCTTTGGGGCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCttt                               <  1:1626921/151‑1 (MQ=14)
                                                                                                                          aaCCTGTTCTTCTGTCTCCGCATTAGCGTGTTGCTGTTGCCGTTTGGGGGCTTTGGCGCGGGCAATAGCGGCTTCTACGGCGGCTTTGCGCGGATCTACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTg                              <  8:387265/151‑1 (MQ=14)
                                                                                                                            ccTCTTCTTCTCCCTCCGTATCAGTCGGTGTCTCTCCTCGTTTGTGTGCTTTGGCGCGGGCAATAGCGGCTTCGACGGCTGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGc                            <  6:810989/151‑1 (MQ=14)
                                                                                                                             gtGTTCTTCTGGCTCGGCCATAGCCGGTGGCTGGTCCGGTGTGCCCGCGTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCa                           <  1:770764/150‑1 (MQ=14)
                                                                                                                               gcagttctTGCTCCGCTTTACCCTTGTGCTGGTCGGGTTTGTTGGCGTGTGCGCGGGCAATAGCCGCTGCTACTGCCGCTTTGCGCGGATCTACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCa                           <  2:1443169/145‑1 (MQ=2)
                                                                                                                               gTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGATTGCGCGCTTTTGCACTCGAAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTAGTTCTGGGTCCGCATTCGCCTGTTGCTGGTCGGGCGTGCGGGCGTGGGCACg                         >  5:547288/1‑151 (MQ=11)
                                                                                                                                  attcggccTCGGCATTAGCCTGCTGCTGGTCCCGTTTGCGCGCGTTCGCGCGGGCAATAGCCGCGTCCACGGCCGCTTTGCGCGGATCCACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGc                      <  1:1097620/144‑1 (MQ=14)
                                                                                                                                   ttcggcctcgggatgaGCCTGGTGCTGGGCCGGTTTGGGGGCGTGGGCGCGGGCAATAGCGGCTGCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGc                      <  7:1061327/135‑1 (MQ=11)
                                                                                                                                       ggctcgtgaTTAGCCTTTTGCTGTTCGCGTTTGCGGGCTTGCGCGCGGGCAATAGCGGCTTCTACGGCGGCTTTGCGCGGATCTACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCaa                    <  4:1240786/140‑1 (MQ=14)
                                                                                                                                        gCTCCGCTTGACCCTGTTGCTGGGCGCGTTTGCGCGCGTGCGCGCGGGCAATAGCCGCTTCCACGGCGGCTTTGCGCGGATCCACCTGTTCTTCTGGCTCGGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGc                <  3:1663217/151‑1 (MQ=14)
                                                                                                                                              cATTTGGCGGTTGCTGTTCTAGTTTTCGAGCTTTTGTGTTTGCTATATCGTCTTTGTGGGGTGCTTTGCGCGGATCTACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCtt           <  6:1460383/150‑1 (MQ=9)
                                                                                                                                                 tAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCTGCTTCGACGGCAGGTTTCCGCGGATCCCCCTGTTGTTCTGGGTCGGGATTATCCTGTTGCTCTTCCAGCTTGCGTGATTAGTCACGGTGAATAGCTTCTTCTAc       >  8:1581491/1‑151 (MQ=255)
                                                                                                                                                  cgCCTGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCCGCTTCGACTGCCGCTTTGCGCGGATCTACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGAc       <  5:970033/149‑1 (MQ=14)
                                                                                                                                                      |                                                                                                                                                      
GCGGCAACTGCCGCTTTGCGCGGATCAACCTGTTCTTCTGGTACCGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGGCATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGG  >  DH1/2169398‑2169698

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