Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,169,6730AT27.3% 54.5 / 19.9 87R629R (CGT→CGA) rsxCelectron transporter RsxC
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base T (4/14);  total (39/48)
Fisher's exact test for biased strand distribution p-value = 5.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGC‑GGATCGACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTC‑TTCT‑GGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACC  >  DH1/2169523‑2169811
                                                                                                                                                           |                                                                                                                                            
cccaccgccgcgCTCCGCACTCGCCGCTTGCTGTTCCCGTTTGCTTTCTTTCGCTCGGGCAATAGCCGCTTCGTCTGCGGCTTTGCTC‑GTTTCGAGCTGTTG‑TC‑GCTTCGGCTCT‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCa                                                                                                                                              <  3:809067/140‑1 (MQ=9)
 agcttttctGGCTCGTGATCACCCTGCTGCTGTGCGCCTTTGTGGGCTTTCGCGCGGGCAATAGCCGCGTCTACTGCGGCTTTGCGC‑GGATCTACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCa                                                                                                                                              <  6:922887/145‑1 (MQ=2)
  ggtcttctGGGTGCGCATTAGCCGGTTGCTGTTCCCGCTTGCGCGCTTCGGCGCGGTGAATAGAGGGCTCGACTGCGGGGTTGCGC‑GGATCCACCTGTTGTTCTGGGTCCG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACg                                                                                                                                            <  1:316919/149‑1 (MQ=11)
    tcttctGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGAAGGCGGCTTTGCGC‑GGATCGACATCTTCATCGGGATCGG‑‑AA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACggg                                                                                                                                          >  7:776021/1‑151 (MQ=14)
       tctGGCTCCTCATTAGCCGGTTGCTGGTCGCGTTTGCGCGCTTTGGCGCGGGCAATAGCGGCTTCGACTGCGGCTTTGCGC‑GGATCGACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGc                                                                                                                                         <  4:528620/149‑1 (MQ=255)
       tctGGCTCCGCATTAGCCTGGTGCTGGTCCGGTTTGCGCGCGTTGGCTCGGGCAATAGCGGCGTCGACGGCGGCTTTGCGC‑GGATCGACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTACCAGCTTGCGTGCTTTG‑GCACGGGc                                                                                                                                         <  3:866242/149‑1 (MQ=14)
         tGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCTGCTTCGACGGCAGATTTACGA‑GGATCGCCAACTAGATCTGGATCGC‑‑CA‑GTAGCATGTTAATGTTCAATCTTCCGGGGTATG‑TCTCGGGcgaa                                                                                                                                      >  2:855845/1‑147 (MQ=255)
          ggCCCCGCATTAGCCGGTTGCTGCTCCCGCTTGCGCGATTTCGCGCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGC‑GGATCTACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATa                                                                                                                                     <  7:220016/150‑1 (MQ=14)
           gggccgccagaGCCGGGGGCTGCTCCCGCGAGCTCGCGTGCGCGCGGGCAATAGCCGCTTCTACTGCCGCTTTGCGC‑GGATCCACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGc                                                                                                                                   <  7:897040/141‑1 (MQ=9)
            cTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGC‑GGATCGACCTGTTGTTCTGGGTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGGGGTTTC‑GCACggggctgagcg                                                                                                                                  >  7:953102/1‑143 (MQ=14)
              ccGCATGAGCCTGTTGCTGTTCGCCTTTTCGCGCCTTCGCTCGGGCAATAGCGGCGTCTACGGCGGCTTTGCGC‑GGATCGACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCg                                                                                                                                  <  5:107050/149‑1 (MQ=14)
               gtctttAGCCTGTTGCTGTTCGGGTTTTCTGGCTTTTGCGCGGGCAATAGCCGCTTCTACTGCGGCTTTGCGC‑GGATCTACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCt                                                                                                                               <  2:560467/147‑1 (MQ=14)
               cGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTAGATCGGAAGATAGACGGTTAGACGAAGACTTTGATC‑GCATCGAAATGTTTGTCTGGTTCCG‑‑CT‑TTAAAGTATTAATATTCCAGCTTGCGGGGGTTG‑GCGCGGGcgggggggggg                                                                                                                               >  2:525730/1‑141 (MQ=255)
                  aTTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGC‑GGATCGACCTGTTGTTCTGACTACG‑‑GA‑TTAGCCTCTTGCTGTTCAAGCTTGAGTTCTATG‑GGCAGGGCAATCGCggatgcg                                                                                                                            >  4:156670/1‑146 (MQ=14)
                     agccttttGCTGTTCGCGTTTGCGCGCGTTGGCTCGGGCAATAGCCGCTTCTACGGCCGCTTTGCGC‑GGATCTACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGAcg                                                                                                                         <  5:399803/145‑1 (MQ=14)
                     aGCCTGGTGCTGGTCCCGTTTGGTCGCTTTGGCGCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGC‑GGATCCACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGAcg                                                                                                                         <  7:1017105/151‑1 (MQ=14)
                     aGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCTGCTTCGACGGCAGATTTACGA‑GGATCAGCAACTTGTTCTAGATCGG‑‑AA‑TAACCAAGTTGATGTTAAAACTTTCGTGCTGTG‑GCTAGCTCTATATCGGCTTCGaat                                                                                                                         >  4:302233/1‑149 (MQ=255)
                      cccttctgctgggcgcGTTTGCGGGCGTTCGCGCGGGCAATAGCCGCTTCCACGGCCGCTTTGCGC‑GGATCTACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGAcgg                                                                                                                        <  3:1004418/136‑1 (MQ=14)
                      cccGGTCGCTGTTCCATTTTTCGCGCTTTGGCTCGGGCAATAGCCGCTTCTACGGCGTGACTGCCG‑TGAACCCGCTGTTTCTCCCGATCTG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGAcgg                                                                                                                        <  3:447294/150‑1 (MQ=9)
                        cggttgcggcggccGTGTGCGGGCTTCGGCGGGGGCATTAGCGGCTGCGACGGCCGCTTTGCGC‑GGATCGACCTGTTCTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGAcggcg                                                                                                                      <  5:320800/139‑1 (MQ=11)
                           ttGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGAGC‑TCGTAGGCGGCGTCTTCT‑GCAGGC‑‑CAGTGACGCGGTAGCTGTTCCGGCTTGCGGGCTTTG‑GCACGGTCAATAGAGGCTTGGCcggggggg                                                                                                                   >  1:804624/1‑149 (MQ=2)
                             gCTGTTCCGGTTTTCGGGCGTTTGCGCGGGCAATAGCCGCGTCTACGGCCGCTTTGCGC‑GGATCTACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCttt                                                                                                                 <  8:150467/151‑1 (MQ=14)
                             gCTGTGCGCCTTTGCGGGCGTGCGCGCGGGCAATAGCCGCTTCGACGGCCGCTTTGCGC‑GGATCCACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCttt                                                                                                                 <  5:74546/151‑1 (MQ=14)
                               tGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGATTTGCGA‑GGATCGACCTGTTGTTCTGGGTCGG‑‑CA‑TTAGCCTCTTGATGTTCAAGATTGCGTGCTTAG‑GCACGGGCAATAGCTGCTTCAACGGCGGAtttc                                                                                                                >  6:488783/1‑149 (MQ=14)
                                 gtCCCGCTTGCGCGCTTGGGCGCGGGCAATAGCCGCGTCGACTGCCGCTTTGCGC‑GGATCAACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTgcgc                                                                                                             <  7:879355/150‑1 (MQ=14)
                                     cgattGCGCGCTTTGGCTCGTGCAATAGCTGCTTCGACTGCGGCTTTGCGC‑GGATCTACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGAt                                                                                                         <  4:345321/148‑1 (MQ=255)
                                     cGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGC‑GGATCGACCTGTTGTTCTGGTTCGG‑‑CA‑TTAGCGTGTTGATGTTCCCGCTTGGGAGCTGTG‑GGCCGTGCGGTGGCGGGGTGGGCGGGGGCTTTGCGGGGAt                                                                                                         >  7:115520/1‑151 (MQ=11)
                                      gTTTGCGCGCTTTGGCTCGGGCAATAGCTACATCGCACGC‑GCAGTGGGCGGGAGCGAAGTCTGTAGGTGGATCCG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGAGAGATTTG‑GCGCGGGCGGTGGCGGGGTcggcggcggggggggggggggg                                                                                                        >  5:222662/1‑138 (MQ=2)
                                      gTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGATTTGCGC‑GGATCGCCCTGGTGTTCTGGGTCGG‑‑CA‑TTAGCCTGTTGGTGTTCGAGCTGGCGTGCTATG‑GCCCGGCCCATAGCTGCTTCGACGGCGGGTTTCCGGGGATc                                                                                                        >  7:837755/1‑151 (MQ=14)
                                       gttGCGCGCGTGGGCGCGGGCAATAGCGGCTGCTACGGCGGCTTGGCGC‑GGAGCGACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCg                                                                                                       <  3:204066/150‑1 (MQ=14)
                                       attgcgaggtgttccTCGGTCATTACCTGCGTCGACGTCCGCTTTGCGC‑TTATCGACGTGCTTGTCTGGTACCG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATc                                                                                                        <  6:420804/136‑1 (MQ=2)
                                           cgcgcTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGGTTTGCGC‑GTATCGACCTGTGCGTGGGGATCTG‑‑CA‑TTAGCGGCTTGCTGTTCCAGCTTGCGGGCTTTG‑GCTCGGGCAATAGCAGCTTCGACGGCGGGTTTGCggggcgggcgcg                                                                                                   >  1:69436/1‑143 (MQ=2)
                                           cgagggtgcgCTCGGGCTCTAGCGGCTTCGTCGGCTGCTGCGCTC‑TGATCCACCTGTTGTTCTGGCTCCG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCt                                                                                                   <  3:462891/142‑1 (MQ=11)
                                            tggcttgtgcgcGGGCAATAGCGGCTTCTACTGCCGCTTTGCGC‑GGATCCACCTGTTGTTCTGGTTCGG‑‑CA‑TTAGCCTGTTGCTGTTCGAGCTTGCGTGCTTTG‑GCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCt                                                                                                   <  8:1014758/139‑1 (MQ=11)
                                            gggcctgcgcgcGGGCAATAGCCGCGTCCACGGCCGCTTTGCGC‑GGATCGACCTGTTGTTCTGGCACGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCGTTG‑GCCCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTg                                                                                                  <  5:634399/140‑1 (MQ=11)
                                             tgcTTTTGCGCGGGCAATAGCCGCTGCTACTGCTGCTTTGCGC‑GTATCTCCCTGTTCTTCTGGGTCGG‑‑CA‑TTAGCCTGTTGCTGTTCGAGGTTGCGTGCTTTT‑GCTCGGGCAATAGCGGCTTCGACTGCGGCTTTGCGCGGATCGACCTg                                                                                                  <  4:390159/149‑1 (MQ=11)
                                               ccgtGGCTCGGGCAATAGCGGCTTCTAGTGCGGCTTTGCGC‑GGAGCGACCTGGTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCGCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGttc                                                                                               <  5:958964/148‑1 (MQ=14)
                                               cTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGGTTTGCGC‑GGATCGCCATGTTTTTCTGGGTCGG‑‑CA‑TTATCATGTAGCTGTTCCAGCTTTCTGGATTTG‑GCTCGCGCCATAGCTGCTTCGACGGCGGGTTTGCGAGGATcaaaaaattt                                                                                               >  3:326497/1‑142 (MQ=2)
                                                 ttcgcgagggCAATAGCCGCGTCCACTGCCGGTTTGCGC‑TGATCTATCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGGTGCTGTGCCTGCTTGCGTGCTTTG‑GCTCGGGCAATAGCGGCTTCTACGGCGGCTTTGCGCGGATCGACCTGttc                                                                                               <  2:625367/142‑1 (MQ=14)
                                                  tcgcgcGGGCAATAGGGGCTGCGACGTCGGCTTTGCTC‑TGATCTACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCTGCTTGCTTGCTTTG‑GCCCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGttc‑tt                                                                                            <  3:923160/145‑1 (MQ=14)
                                                  tccctcgcccaagcgCGGCCTCGCCTGCCGCTTTGCGC‑GGATCTACCTGTTCTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCTCGTGCAATAGCTGCTTCGACGGCGGCTTTGCGCGGATCGACCTGttc‑tt                                                                                            <  4:168751/136‑1 (MQ=11)
                                                  tGGCGCGGGCACGCGCGGCTTGGACGGCAGCCGTGGGC‑TGATCCCAGCGTTCTTCTGGCTCCG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGttc‑tt                                                                                            <  5:608636/150‑1 (MQ=11)
                                                    gcgcGGGCAATAGCCGCGGCAACTGCGGCTTTGCGC‑GGATCTACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGGTGCTGTTCCAGCTTGCGTGCTTTG‑GCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTtc‑ttct                                                                                          <  2:1005225/147‑1 (MQ=11)
                                                    gCTCGGGCAATAGCGGCTTCGAGGGCGGCTGAGCGC‑GGATCGCCCGGTTGCTCTGGCTCGG‑‑GA‑TTAGCCGGTTGCTGTTCCAGCTTGCGGGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTC‑TTCT‑g                                                                                        <  3:156618/151‑1 (MQ=14)
                                                     cgcGGGCAATAGCCGCGTCCACGGCCGCTTTGCGC‑GGATCCACCTGTTGTTCTGGCGCGG‑‑CA‑TTCGCCGGTGGCCGGTCCAGCTTTCGGGCGTTG‑GCTCGGGCAATAGCGGCTTCCACGGCGGCTTTGCGCGGATCGACCTGTTC‑TTCT‑gg                                                                                       <  5:15173/149‑1 (MQ=11)
                                                      gcGGGCAATAGCCGCGTCAACGGCCGCTTTGCGC‑GGATCAACCTGTTCTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTGAAGCTTGCTTGCTTTT‑GCTCGGGCAATAGCGGCTTCTACTGCGGCTTTGCGCGGATCGACCTGTTC‑TTCT‑GGt                                                                                      <  8:378585/150‑1 (MQ=11)
                                                       cGGGCAATAGCCGCTTCCACGGCCGCTTTGCGC‑GGATCGACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTC‑TTCT‑gg                                                                                       <  7:516703/149‑1 (MQ=255)
                                                         ggCAATAGCTGCTTCAACTGCCGCTTTGCGC‑GGAGCTACCTCTTGTTCTGGATCCG‑‑CA‑TTAGCCTGTTGCGGTTCCAGCTTGCGTGCTTTC‑GCTCGGGCAATAGCGGCTTCTACTGCCGCTTTGCGCGGATCGACCTGTTC‑TTCT‑GGTTcc                                                                                   <  8:733571/151‑1 (MQ=11)
                                                          gcattttgcgCTTCTCCGGCTTCTTTTCGC‑TGATCGACGTGTTGTTGTGTTTGTG‑‑CA‑TTTGCCTGTTGCTGGTCGTGTTTTTGTGCTTTT‑GCTCGGGCAATAGCGGCTTGTAGTGCGGCTTTGCGCGGATCGACCTGTTC‑TTCT‑GGTTCCg                                                                                  <  6:679653/142‑1 (MQ=2)
                                                           ccatccccgCTTCCACCGCCGCTTTGCGC‑TTATCTAAATTGTGTTGTGGACCCG‑‑CA‑TTAGCCTGTTGCGGTTCCAGCTTGCTTGCTTTC‑TCTAGGGCAATAGCCGCTGCTACTGCGGCTTTGCGCGGATCGACCTGTTC‑TTCT‑GGTTCCGc                                                                                 <  4:142292/143‑1 (MQ=2)
                                                                gcggcgTCCACGGCGTGTTTGCGC‑GGATACACCTGTTTTTCTGTCCACG‑‑AA‑TCTGGTGGTTGCTGGGCGTGGTTGTTGGCTTGT‑GCGCGGGCAATAGCGGCGTCTACTGCGGCTTTGCGCGGATCTACCTGTTC‑TTCT‑GGTTCCGCATTa                                                                             <  8:672367/150‑1 (MQ=255)
                                                                            gcggcTTTGCGC‑GGATCGACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTC‑TTCT‑GGTTCCGCATTAGCCTGTTGCTGtt                                                                >  4:535888/1‑151 (MQ=40)
                                                                                cTTTGCGC‑GGATCGACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCTATAGCGGCTTCGACGGCAGATTTGCGTGGATCACCCACTTC‑TTCT‑GCTTCCTTAGTAGCCTGAAGGTGTTcc                                                              >  7:394965/1‑149 (MQ=14)
                                                                                 tttGCGC‑GGATCGACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCCCGGGCGATAGCTGCTTCGACGGCGGGTTTGCGGGGATCCCCAACGTT‑TTCT‑GGTTCCTTAGTATCCTGATGATGTTCCAgg                                                           >  3:807043/1‑150 (MQ=14)
                                                                                  ttGCGC‑GGATCGACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGAGCTTTG‑GCACGGGCAATAGCGGCTTCCACGGCGGGTTTCCGCGGATCGCCAACTGT‑TTCT‑GGGTCCGGAGTAGCCTGATGTTGTGCCacg                                                           >  1:220210/1‑148 (MQ=11)
                                                                                         ggatctccgtgttgaGCTGCCTCCG‑‑GA‑CTAGCCGGCTGCTGATCGCGCTTTGTTGCGTTG‑GCTCGGGCAATAGCCGCGGCTACGGCCGCTTTGCGCGGATCGACCTGTTC‑TTCT‑GGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTgcgcg                                                    <  7:952816/142‑1 (MQ=2)
                                                                                          gATCGACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGCCCTGTTC‑TTCT‑GGGTCCGGATTATGCTGTTGCTGTTCCAGCTTgcgagg                                                   >  7:829839/1‑148 (MQ=14)
                                                                                           aTCGACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCTGCTTCGACGGCGGGTTTGCGAGGATCGCCATCGTG‑TGCT‑GCGTCCGGAGTAGCCTGTTGCTGTTCCAGCTTgcgccgt                                                  >  3:300733/1‑148 (MQ=11)
                                                                                           aTCGACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGGGCTTTG‑GCGCGGGCAATAGCTGCTTCAACGGCAGCTGTGAACCGAACCAGGAGGAC‑ATCG‑GGTAACTCAGTAGCATGTTGCTGTTCCCGCTTGCGCGCt                                                  >  3:607454/1‑151 (MQ=2)
                                                                                           aTCGACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGAGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCAGGTTTACGAGGATCGACCTCTTG‑TTCT‑GCTTCCTAATTATCCTGATGATCTTCCAGGTTccgacat                                                  >  6:224074/1‑147 (MQ=11)
                                                                                              gACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACATGTTG‑TTCT‑GGTTCGGCATTATCCTGTTGCTGTTCCAGCTTGCGCCCtttt                                               >  7:5678/1‑150 (MQ=255)
                                                                                                 ccggtcgtCTGGCTCGT‑‑GA‑TTAGCCGGGTGCTGGTCGCGCTTGCGGGCTCTGCGCA‑GGGCAATCCCGGCATAGTCGGCCGCTTCGGTCGTATCGACCTGTTG‑TTCT‑GGCTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGcg                                            <  5:641422/144‑2 (MQ=2)
                                                                                                 cTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCCCGGGCAATAGCTGCTTCGACGGCGGGTTTGCGGGGATCGCCCACTTG‑TTCT‑GGTTCCTGATTATCCTGTTGATGTTCCAGCttcccacatttgtgc                                            >  5:1022531/1‑138 (MQ=11)
                                                                                                 cTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGAGCTTTG‑GCGCGGGAAATAGCAACATCGACGTCTTCGGTAAGCGTGTCGACCGCTTG‑TTAT‑GTTTCCGGGTTAGCCGGTTGATTTTAAAACTTGCGCGCGTGGGcg                                            >  4:904487/1‑150 (MQ=2)
                                                                                                   ggtgttCTGGCTCGG‑‑CA‑TTAGCCTTGTGCTGGTCGTGGTTGGGGGCCTTC‑GCGCGGGCAATAGCGGCTGCGACGGCGGCTTTGCGCGGATCGACCTGTTC‑TTCT‑GGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACg                                          <  7:726215/149‑1 (MQ=14)
                                                                                                    ccccccccccTCCG‑‑CC‑TTCGCCGGGTGCTGTTCCAGCTTGCGTGCTTTG‑GCCCGGGCAATAGCCGCTTCTACGGCGGGATTGCAGGTATC‑ACCGGTGCGCTCC‑GGATCGGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACgg                                         <  1:457718/142‑1 (MQ=2)
                                                                                                      gttCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCCCGGGCGATAGCGGCTTCCACGGCGGGTTTGCGGGGATCGGCCTCGTG‑TTCT‑GGTTCGGGAGTTGCCTGTTGGTGTTCCTGTTTgcgtccggggtgcggtgt                                       >  1:772249/1‑136 (MQ=2)
                                                                                                      gttCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGGTTTGCGAGGATCGGCCTGTTG‑TTCT‑GGTTCCTGATTTTCCTGTTGGTGTTCCTGATTGCGGGGTTTgggcggggg                                       >  7:696489/1‑144 (MQ=11)
                                                                                                         cTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCTGCTTCGACGGCGGGTTTGCGGGGATCGGCCAGGTG‑TTCT‑GGTTCGGGATTAGCCTGTTGGTGTTCCAGCTTGCGCGCGTTGGCCCGTGCAAt                                    >  7:821677/1‑151 (MQ=11)
                                                                                                         cTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGTCCTGTTG‑TTCT‑GCTTCGGCATTAGCCTGTTGCTGTTCCAGCTTCCGCGGTTTGGCCCGGGCact                                    >  7:367385/1‑149 (MQ=11)
                                                                                                          tGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTG‑TTCT‑GGTTCCGGATTAGCCTGTTGGTGTTCCAGCTTGCGCGGTTTGGCCGGGGCAAtt                                   >  3:416423/1‑150 (MQ=11)
                                                                                                             cTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTC‑TTCT‑GGTTCCGGATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCCCGGGAAATAGct                                >  7:336341/1‑150 (MQ=255)
                                                                                                                    cA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGGTTTGCGCGGATCGCCCTGGTG‑TTCT‑GGTTCGGGATTAGCGTGTTGGGGTTCGAGGTAGCTCGCTTTGGCCCGGGGAATAGCGGCTTc                           >  3:313046/1‑151 (MQ=11)
                                                                                                                     a‑TTAGCCTGTTGCTGTTGCATCTTGGGTGCTTTG‑TCTCTTTCAATAGACGATGCGATTGCGTGGTTGCGCGGATCTACCTGTTC‑TTCT‑GGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCGGCTTCg                          <  8:42053/151‑1 (MQ=14)
                                                                                                                     a‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTC‑TTCT‑GGCTCCGCATTAGCCTGTTGCTGTTCGGCAGTGCGCGCGTGGGCACTGGcgtgaggggggtgc                          >  7:996830/1‑138 (MQ=14)
                                                                                                                        tAGCCTGTTTCTGTTCCAGCGTGCGAGCGTTC‑GCACGGCCTCTACCGCCATCGACGGCGTCCTTGCGCGGATCGAGATCTTC‑TTCT‑GGCTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCGGCTTCGa                         <  8:775088/150‑1 (MQ=14)
                                                                                                                        tAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTC‑TTCT‑GGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCGGCTTCgcc                        >  5:1032324/1‑149 (MQ=255)
                                                                                                                         aGCCTGTTGCTTTTCGAGCGTGCGGGCTGTG‑GCTCGGGCAATAGCGGCGGCGTGGGCGGCTTTGCGCGGTGCGAGCTGGCG‑ATCT‑GGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCGGCTTCGAcg                       <  3:932558/151‑1 (MQ=11)
                                                                                                                         aGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCTGCTTCCACGGCCGGTTTGCGAGGATCGACATGTTG‑TTCT‑GGTTCCTCATTATCCTGATGATGTTCCAGATTCCGAGTTTTTTCACGTTGAATAGCTGCTTctgct                       >  4:795799/1‑147 (MQ=2)
                                                                                                                         aGCCGGTTGCGGTGCCCGCTTGCCGGCTTTG‑CCGCGGGCAATCGCGGCGTCGCGGGGGGCTTTCCGCGGTTCTACCTGTTC‑TTCT‑GGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGC‑TTGGCACGGGCAATAGCGGCTTCGAcg                       <  7:348149/150‑1 (MQ=11)
                                                                                                                                ggCTGGTCCAGCTTGCGTGCTTTG‑GCGCGGGCAATAGCGGCGTCGACGGCTGCTTCGCGCGGAGCGCCCTGGTG‑TGCT‑GGGGGCGCCGTAGCGTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCtt                <  1:326636/150‑1 (MQ=14)
                                                                                                                                agctgtgcgctgctGCGCGCCTTG‑GCGCGGGCAATAGCGGCTTCTACGGCGGCTTTGCGCGGATCTACCTGTTC‑TTCT‑GGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCt                 <  3:295496/137‑1 (MQ=11)
                                                                                                                                 gCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGTAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACATGTTC‑TTCT‑GGTTCCGCATTAGCCTGTTGCTGTTCCAGATTGCGCGCTATGGCCCGGGGAATCCCGGCATCAACGGCGGGttt               >  5:712990/1‑151 (MQ=14)
                                                                                                                                 cccgTTCCAGCTTTTTTTTTTTG‑CCACGGGCAATAGCGTCTTCTACCGAGGCTTTGAGAGTATTTACCTCCTC‑TACTCTCTTCCG‑ATCAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCttt               <  2:694290/148‑1 (MQ=2)
                                                                                                                                    gagggccgacgCGTGCTCTG‑GCGCGCTCAATAGCCGCGGCCACTGCCGCTTTGCGCGGATCCACCTGTTC‑TTCT‑GGTACCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCCCGGGCAATAGCGGCTTCGACGGCGGCTTTgcg            <  1:829590/141‑1 (MQ=11)
                                                                                                                                        cAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTG‑TTCT‑GGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGTCCCGGGCAATAGCTGCTTCGACGGCTGGTTTTCGAggt        >  7:18220/1‑150 (MQ=11)
                                                                                                                                         aGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTCTTG‑TTCT‑GCTTCCGGATTAGACTGTTGCTGTGCCAGCTTGCGCCCTTTCCCTCGGGCGATAGCGGCTTCGACGGCGGTTTTtcgttgat       >  3:913983/1‑146 (MQ=11)
                                                                                                                                            ttGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTC‑TTCT‑GGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCTGCTTCGACGGCTGGTTTgcggcgctcg     >  3:517905/1‑143 (MQ=14)
                                                                                                                                               cGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTC‑TTCT‑GGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGATTTTGCACGTGAAATAGCTGCTTCTCCTGCAGTTTTTCGAGGATctccc  >  8:502999/1‑146 (MQ=11)
                                                                                                                                                           |                                                                                                                                            
CTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGC‑GGATCGACCTGTTGTTCTGGCTCGG‑‑CA‑TTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTG‑GCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTC‑TTCT‑GGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACC  >  DH1/2169523‑2169811

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