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→CGArsxCelectron 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

CTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTCTTCTGGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACC  >  DH1/2169523‑2169811
                                                                                                                                                           |                                                                                                                                            
cccaccgccgcgCTCCGCACTCGCCGCTTGCTGTTCCCGTTTGCTTTCTTTCGCTCGGGCAATAGCCGCTTCGTCTGCGGCTTTGCTCGTTTCGAGCTGTTGTCGCTTCGGCTCTTTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCa                                                                                                                                              <  3:809067/140‑1 (MQ=9)
 agcttttctGGCTCGTGATCACCCTGCTGCTGTGCGCCTTTGTGGGCTTTCGCGCGGGCAATAGCCGCGTCTACTGCGGCTTTGCGCGGATCTACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCa                                                                                                                                              <  6:922887/145‑1 (MQ=2)
  ggtcttctGGGTGCGCATTAGCCGGTTGCTGTTCCCGCTTGCGCGCTTCGGCGCGGTGAATAGAGGGCTCGACTGCGGGGTTGCGCGGATCCACCTGTTGTTCTGGGTCCG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACg                                                                                                                                            <  1:316919/149‑1 (MQ=11)
    tcttctGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGAAGGCGGCTTTGCGCGGATCGACATCTTCATCGGGATCGG‑‑AATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACggg                                                                                                                                          >  7:776021/1‑151 (MQ=14)
       tctGGCTCCTCATTAGCCGGTTGCTGGTCGCGTTTGCGCGCTTTGGCGCGGGCAATAGCGGCTTCGACTGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGc                                                                                                                                         <  4:528620/149‑1 (MQ=255)
       tctGGCTCCGCATTAGCCTGGTGCTGGTCCGGTTTGCGCGCGTTGGCTCGGGCAATAGCGGCGTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTACCAGCTTGCGTGCTTTGGCACGGGc                                                                                                                                         <  3:866242/149‑1 (MQ=14)
         tGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCTGCTTCGACGGCAGATTTACGAGGATCGCCAACTAGATCTGGATCGC‑‑CAGTAGCATGTTAATGTTCAATCTTCCGGGGTATGTCTCGGGcgaa                                                                                                                                      >  2:855845/1‑147 (MQ=255)
          ggCCCCGCATTAGCCGGTTGCTGCTCCCGCTTGCGCGATTTCGCGCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCTACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATa                                                                                                                                     <  7:220016/150‑1 (MQ=14)
           gggccgccagaGCCGGGGGCTGCTCCCGCGAGCTCGCGTGCGCGCGGGCAATAGCCGCTTCTACTGCCGCTTTGCGCGGATCCACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGc                                                                                                                                   <  7:897040/141‑1 (MQ=9)
            cTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGGTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGGGGTTTCGCACggggctgagcg                                                                                                                                  >  7:953102/1‑143 (MQ=14)
              ccGCATGAGCCTGTTGCTGTTCGCCTTTTCGCGCCTTCGCTCGGGCAATAGCGGCGTCTACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCg                                                                                                                                  <  5:107050/149‑1 (MQ=14)
               gtctttAGCCTGTTGCTGTTCGGGTTTTCTGGCTTTTGCGCGGGCAATAGCCGCTTCTACTGCGGCTTTGCGCGGATCTACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCt                                                                                                                               <  2:560467/147‑1 (MQ=14)
               cGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTAGATCGGAAGATAGACGGTTAGACGAAGACTTTGATCGCATCGAAATGTTTGTCTGGTTCCG‑‑CTTTAAAGTATTAATATTCCAGCTTGCGGGGGTTGGCGCGGGcgggggggggg                                                                                                                               >  2:525730/1‑141 (MQ=255)
                  aTTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGACTACG‑‑GATTAGCCTCTTGCTGTTCAAGCTTGAGTTCTATGGGCAGGGCAATCGCggatgcg                                                                                                                            >  4:156670/1‑146 (MQ=14)
                     agccttttGCTGTTCGCGTTTGCGCGCGTTGGCTCGGGCAATAGCCGCTTCTACGGCCGCTTTGCGCGGATCTACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGAcg                                                                                                                         <  5:399803/145‑1 (MQ=14)
                     aGCCTGGTGCTGGTCCCGTTTGGTCGCTTTGGCGCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCCACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGAcg                                                                                                                         <  7:1017105/151‑1 (MQ=14)
                     aGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCTGCTTCGACGGCAGATTTACGAGGATCAGCAACTTGTTCTAGATCGG‑‑AATAACCAAGTTGATGTTAAAACTTTCGTGCTGTGGCTAGCTCTATATCGGCTTCGaat                                                                                                                         >  4:302233/1‑149 (MQ=255)
                      cccttctgctgggcgcGTTTGCGGGCGTTCGCGCGGGCAATAGCCGCTTCCACGGCCGCTTTGCGCGGATCTACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGAcgg                                                                                                                        <  3:1004418/136‑1 (MQ=14)
                      cccGGTCGCTGTTCCATTTTTCGCGCTTTGGCTCGGGCAATAGCCGCTTCTACGGCGTGACTGCCGTGAACCCGCTGTTTCTCCCGATCTG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGAcgg                                                                                                                        <  3:447294/150‑1 (MQ=9)
                        cggttgcggcggccGTGTGCGGGCTTCGGCGGGGGCATTAGCGGCTGCGACGGCCGCTTTGCGCGGATCGACCTGTTCTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGAcggcg                                                                                                                      <  5:320800/139‑1 (MQ=11)
                           ttGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGAGCTCGTAGGCGGCGTCTTCTGCAGGC‑‑CAGTGACGCGGTAGCTGTTCCGGCTTGCGGGCTTTGGCACGGTCAATAGAGGCTTGGCcggggggg                                                                                                                   >  1:804624/1‑149 (MQ=2)
                             gCTGTTCCGGTTTTCGGGCGTTTGCGCGGGCAATAGCCGCGTCTACGGCCGCTTTGCGCGGATCTACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCttt                                                                                                                 <  8:150467/151‑1 (MQ=14)
                             gCTGTGCGCCTTTGCGGGCGTGCGCGCGGGCAATAGCCGCTTCGACGGCCGCTTTGCGCGGATCCACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCttt                                                                                                                 <  5:74546/151‑1 (MQ=14)
                               tGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGATTTGCGAGGATCGACCTGTTGTTCTGGGTCGG‑‑CATTAGCCTCTTGATGTTCAAGATTGCGTGCTTAGGCACGGGCAATAGCTGCTTCAACGGCGGAtttc                                                                                                                >  6:488783/1‑149 (MQ=14)
                                 gtCCCGCTTGCGCGCTTGGGCGCGGGCAATAGCCGCGTCGACTGCCGCTTTGCGCGGATCAACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTgcgc                                                                                                             <  7:879355/150‑1 (MQ=14)
                                     cgattGCGCGCTTTGGCTCGTGCAATAGCTGCTTCGACTGCGGCTTTGCGCGGATCTACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGAt                                                                                                         <  4:345321/148‑1 (MQ=255)
                                     cGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGTTCGG‑‑CATTAGCGTGTTGATGTTCCCGCTTGGGAGCTGTGGGCCGTGCGGTGGCGGGGTGGGCGGGGGCTTTGCGGGGAt                                                                                                         >  7:115520/1‑151 (MQ=11)
                                      gTTTGCGCGCTTTGGCTCGGGCAATAGCTACATCGCACGCGCAGTGGGCGGGAGCGAAGTCTGTAGGTGGATCCG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGAGAGATTTGGCGCGGGCGGTGGCGGGGTcggcggcggggggggggggggg                                                                                                        >  5:222662/1‑138 (MQ=2)
                                      gTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGATTTGCGCGGATCGCCCTGGTGTTCTGGGTCGG‑‑CATTAGCCTGTTGGTGTTCGAGCTGGCGTGCTATGGCCCGGCCCATAGCTGCTTCGACGGCGGGTTTCCGGGGATc                                                                                                        >  7:837755/1‑151 (MQ=14)
                                       gttGCGCGCGTGGGCGCGGGCAATAGCGGCTGCTACGGCGGCTTGGCGCGGAGCGACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCg                                                                                                       <  3:204066/150‑1 (MQ=14)
                                       attgcgaggtgttccTCGGTCATTACCTGCGTCGACGTCCGCTTTGCGCTTATCGACGTGCTTGTCTGGTACCG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATc                                                                                                        <  6:420804/136‑1 (MQ=2)
                                           cgcgcTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGGTTTGCGCGTATCGACCTGTGCGTGGGGATCTG‑‑CATTAGCGGCTTGCTGTTCCAGCTTGCGGGCTTTGGCTCGGGCAATAGCAGCTTCGACGGCGGGTTTGCggggcgggcgcg                                                                                                   >  1:69436/1‑143 (MQ=2)
                                           cgagggtgcgCTCGGGCTCTAGCGGCTTCGTCGGCTGCTGCGCTCTGATCCACCTGTTGTTCTGGCTCCG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCt                                                                                                   <  3:462891/142‑1 (MQ=11)
                                            tggcttgtgcgcGGGCAATAGCGGCTTCTACTGCCGCTTTGCGCGGATCCACCTGTTGTTCTGGTTCGG‑‑CATTAGCCTGTTGCTGTTCGAGCTTGCGTGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCt                                                                                                   <  8:1014758/139‑1 (MQ=11)
                                            gggcctgcgcgcGGGCAATAGCCGCGTCCACGGCCGCTTTGCGCGGATCGACCTGTTGTTCTGGCACGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCGTTGGCCCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTg                                                                                                  <  5:634399/140‑1 (MQ=11)
                                             tgcTTTTGCGCGGGCAATAGCCGCTGCTACTGCTGCTTTGCGCGTATCTCCCTGTTCTTCTGGGTCGG‑‑CATTAGCCTGTTGCTGTTCGAGGTTGCGTGCTTTTGCTCGGGCAATAGCGGCTTCGACTGCGGCTTTGCGCGGATCGACCTg                                                                                                  <  4:390159/149‑1 (MQ=11)
                                               ccgtGGCTCGGGCAATAGCGGCTTCTAGTGCGGCTTTGCGCGGAGCGACCTGGTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGttc                                                                                               <  5:958964/148‑1 (MQ=14)
                                               cTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGGTTTGCGCGGATCGCCATGTTTTTCTGGGTCGG‑‑CATTATCATGTAGCTGTTCCAGCTTTCTGGATTTGGCTCGCGCCATAGCTGCTTCGACGGCGGGTTTGCGAGGATcaaaaaattt                                                                                               >  3:326497/1‑142 (MQ=2)
                                                 ttcgcgagggCAATAGCCGCGTCCACTGCCGGTTTGCGCTGATCTATCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGGTGCTGTGCCTGCTTGCGTGCTTTGGCTCGGGCAATAGCGGCTTCTACGGCGGCTTTGCGCGGATCGACCTGttc                                                                                               <  2:625367/142‑1 (MQ=14)
                                                  tcgcgcGGGCAATAGGGGCTGCGACGTCGGCTTTGCTCTGATCTACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCTGCTTGCTTGCTTTGGCCCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGttctt                                                                                            <  3:923160/145‑1 (MQ=14)
                                                  tccctcgcccaagcgCGGCCTCGCCTGCCGCTTTGCGCGGATCTACCTGTTCTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCTCGTGCAATAGCTGCTTCGACGGCGGCTTTGCGCGGATCGACCTGttctt                                                                                            <  4:168751/136‑1 (MQ=11)
                                                  tGGCGCGGGCACGCGCGGCTTGGACGGCAGCCGTGGGCTGATCCCAGCGTTCTTCTGGCTCCG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGttctt                                                                                            <  5:608636/150‑1 (MQ=11)
                                                    gcgcGGGCAATAGCCGCGGCAACTGCGGCTTTGCGCGGATCTACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGGTGCTGTTCCAGCTTGCGTGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTtcttct                                                                                          <  2:1005225/147‑1 (MQ=11)
                                                    gCTCGGGCAATAGCGGCTTCGAGGGCGGCTGAGCGCGGATCGCCCGGTTGCTCTGGCTCGG‑‑GATTAGCCGGTTGCTGTTCCAGCTTGCGGGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTCTTCTg                                                                                        <  3:156618/151‑1 (MQ=14)
                                                     cgcGGGCAATAGCCGCGTCCACGGCCGCTTTGCGCGGATCCACCTGTTGTTCTGGCGCGG‑‑CATTCGCCGGTGGCCGGTCCAGCTTTCGGGCGTTGGCTCGGGCAATAGCGGCTTCCACGGCGGCTTTGCGCGGATCGACCTGTTCTTCTgg                                                                                       <  5:15173/149‑1 (MQ=11)
                                                      gcGGGCAATAGCCGCGTCAACGGCCGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTGAAGCTTGCTTGCTTTTGCTCGGGCAATAGCGGCTTCTACTGCGGCTTTGCGCGGATCGACCTGTTCTTCTGGt                                                                                      <  8:378585/150‑1 (MQ=11)
                                                       cGGGCAATAGCCGCTTCCACGGCCGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTCTTCTgg                                                                                       <  7:516703/149‑1 (MQ=255)
                                                         ggCAATAGCTGCTTCAACTGCCGCTTTGCGCGGAGCTACCTCTTGTTCTGGATCCG‑‑CATTAGCCTGTTGCGGTTCCAGCTTGCGTGCTTTCGCTCGGGCAATAGCGGCTTCTACTGCCGCTTTGCGCGGATCGACCTGTTCTTCTGGTTcc                                                                                   <  8:733571/151‑1 (MQ=11)
                                                          gcattttgcgCTTCTCCGGCTTCTTTTCGCTGATCGACGTGTTGTTGTGTTTGTG‑‑CATTTGCCTGTTGCTGGTCGTGTTTTTGTGCTTTTGCTCGGGCAATAGCGGCTTGTAGTGCGGCTTTGCGCGGATCGACCTGTTCTTCTGGTTCCg                                                                                  <  6:679653/142‑1 (MQ=2)
                                                           ccatccccgCTTCCACCGCCGCTTTGCGCTTATCTAAATTGTGTTGTGGACCCG‑‑CATTAGCCTGTTGCGGTTCCAGCTTGCTTGCTTTCTCTAGGGCAATAGCCGCTGCTACTGCGGCTTTGCGCGGATCGACCTGTTCTTCTGGTTCCGc                                                                                 <  4:142292/143‑1 (MQ=2)
                                                                gcggcgTCCACGGCGTGTTTGCGCGGATACACCTGTTTTTCTGTCCACG‑‑AATCTGGTGGTTGCTGGGCGTGGTTGTTGGCTTGTGCGCGGGCAATAGCGGCGTCTACTGCGGCTTTGCGCGGATCTACCTGTTCTTCTGGTTCCGCATTa                                                                             <  8:672367/150‑1 (MQ=255)
                                                                            gcggcTTTGCGCGGATCGACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTCTTCTGGTTCCGCATTAGCCTGTTGCTGtt                                                                >  4:535888/1‑151 (MQ=40)
                                                                                cTTTGCGCGGATCGACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCTATAGCGGCTTCGACGGCAGATTTGCGTGGATCACCCACTTCTTCTGCTTCCTTAGTAGCCTGAAGGTGTTcc                                                              >  7:394965/1‑149 (MQ=14)
                                                                                 tttGCGCGGATCGACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCCCGGGCGATAGCTGCTTCGACGGCGGGTTTGCGGGGATCCCCAACGTTTTCTGGTTCCTTAGTATCCTGATGATGTTCCAgg                                                           >  3:807043/1‑150 (MQ=14)
                                                                                  ttGCGCGGATCGACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGAGCTTTGGCACGGGCAATAGCGGCTTCCACGGCGGGTTTCCGCGGATCGCCAACTGTTTCTGGGTCCGGAGTAGCCTGATGTTGTGCCacg                                                           >  1:220210/1‑148 (MQ=11)
                                                                                         ggatctccgtgttgaGCTGCCTCCG‑‑GACTAGCCGGCTGCTGATCGCGCTTTGTTGCGTTGGCTCGGGCAATAGCCGCGGCTACGGCCGCTTTGCGCGGATCGACCTGTTCTTCTGGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTgcgcg                                                    <  7:952816/142‑1 (MQ=2)
                                                                                          gATCGACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGCCCTGTTCTTCTGGGTCCGGATTATGCTGTTGCTGTTCCAGCTTgcgagg                                                   >  7:829839/1‑148 (MQ=14)
                                                                                           aTCGACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCTGCTTCGACGGCGGGTTTGCGAGGATCGCCATCGTGTGCTGCGTCCGGAGTAGCCTGTTGCTGTTCCAGCTTgcgccgt                                                  >  3:300733/1‑148 (MQ=11)
                                                                                           aTCGACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGGGCTTTGGCGCGGGCAATAGCTGCTTCAACGGCAGCTGTGAACCGAACCAGGAGGACATCGGGTAACTCAGTAGCATGTTGCTGTTCCCGCTTGCGCGCt                                                  >  3:607454/1‑151 (MQ=2)
                                                                                           aTCGACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGAGCTTTGGCACGGGCAATAGCGGCTTCGACGGCAGGTTTACGAGGATCGACCTCTTGTTCTGCTTCCTAATTATCCTGATGATCTTCCAGGTTccgacat                                                  >  6:224074/1‑147 (MQ=11)
                                                                                              gACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACATGTTGTTCTGGTTCGGCATTATCCTGTTGCTGTTCCAGCTTGCGCCCtttt                                               >  7:5678/1‑150 (MQ=255)
                                                                                                 ccggtcgtCTGGCTCGT‑‑GATTAGCCGGGTGCTGGTCGCGCTTGCGGGCTCTGCGCAGGGCAATCCCGGCATAGTCGGCCGCTTCGGTCGTATCGACCTGTTGTTCTGGCTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGcg                                            <  5:641422/144‑2 (MQ=2)
                                                                                                 cTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCCCGGGCAATAGCTGCTTCGACGGCGGGTTTGCGGGGATCGCCCACTTGTTCTGGTTCCTGATTATCCTGTTGATGTTCCAGCttcccacatttgtgc                                            >  5:1022531/1‑138 (MQ=11)
                                                                                                 cTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGAGCTTTGGCGCGGGAAATAGCAACATCGACGTCTTCGGTAAGCGTGTCGACCGCTTGTTATGTTTCCGGGTTAGCCGGTTGATTTTAAAACTTGCGCGCGTGGGcg                                            >  4:904487/1‑150 (MQ=2)
                                                                                                   ggtgttCTGGCTCGG‑‑CATTAGCCTTGTGCTGGTCGTGGTTGGGGGCCTTCGCGCGGGCAATAGCGGCTGCGACGGCGGCTTTGCGCGGATCGACCTGTTCTTCTGGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACg                                          <  7:726215/149‑1 (MQ=14)
                                                                                                    ccccccccccTCCG‑‑CCTTCGCCGGGTGCTGTTCCAGCTTGCGTGCTTTGGCCCGGGCAATAGCCGCTTCTACGGCGGGATTGCAGGTATCACCGGTGCGCTCCGGATCGGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACgg                                         <  1:457718/142‑1 (MQ=2)
                                                                                                      gttCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCCCGGGCGATAGCGGCTTCCACGGCGGGTTTGCGGGGATCGGCCTCGTGTTCTGGTTCGGGAGTTGCCTGTTGGTGTTCCTGTTTgcgtccggggtgcggtgt                                       >  1:772249/1‑136 (MQ=2)
                                                                                                      gttCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGGTTTGCGAGGATCGGCCTGTTGTTCTGGTTCCTGATTTTCCTGTTGGTGTTCCTGATTGCGGGGTTTgggcggggg                                       >  7:696489/1‑144 (MQ=11)
                                                                                                         cTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCTGCTTCGACGGCGGGTTTGCGGGGATCGGCCAGGTGTTCTGGTTCGGGATTAGCCTGTTGGTGTTCCAGCTTGCGCGCGTTGGCCCGTGCAAt                                    >  7:821677/1‑151 (MQ=11)
                                                                                                         cTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGTCCTGTTGTTCTGCTTCGGCATTAGCCTGTTGCTGTTCCAGCTTCCGCGGTTTGGCCCGGGCact                                    >  7:367385/1‑149 (MQ=11)
                                                                                                          tGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGTTCCGGATTAGCCTGTTGGTGTTCCAGCTTGCGCGGTTTGGCCGGGGCAAtt                                   >  3:416423/1‑150 (MQ=11)
                                                                                                             cTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTCTTCTGGTTCCGGATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCCCGGGAAATAGct                                >  7:336341/1‑150 (MQ=255)
                                                                                                                    cATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGGTTTGCGCGGATCGCCCTGGTGTTCTGGTTCGGGATTAGCGTGTTGGGGTTCGAGGTAGCTCGCTTTGGCCCGGGGAATAGCGGCTTc                           >  3:313046/1‑151 (MQ=11)
                                                                                                                     aTTAGCCTGTTGCTGTTGCATCTTGGGTGCTTTGTCTCTTTCAATAGACGATGCGATTGCGTGGTTGCGCGGATCTACCTGTTCTTCTGGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCGGCTTCg                          <  8:42053/151‑1 (MQ=14)
                                                                                                                     aTTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCGCGTGCAATAGCCGCTTCGACGGCAGCTTTGCGCGGATCAACCTGTTCTTCTGGCTCCGCATTAGCCTGTTGCTGTTCGGCAGTGCGCGCGTGGGCACTGGcgtgaggggggtgc                          >  7:996830/1‑138 (MQ=14)
                                                                                                                        tAGCCTGTTTCTGTTCCAGCGTGCGAGCGTTCGCACGGCCTCTACCGCCATCGACGGCGTCCTTGCGCGGATCGAGATCTTCTTCTGGCTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCGGCTTCGa                         <  8:775088/150‑1 (MQ=14)
                                                                                                                        tAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTCTTCTGGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCGGCTTCgcc                        >  5:1032324/1‑149 (MQ=255)
                                                                                                                         aGCCTGTTGCTTTTCGAGCGTGCGGGCTGTGGCTCGGGCAATAGCGGCGGCGTGGGCGGCTTTGCGCGGTGCGAGCTGGCGATCTGGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCGGCTTCGAcg                       <  3:932558/151‑1 (MQ=11)
                                                                                                                         aGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCTGCTTCCACGGCCGGTTTGCGAGGATCGACATGTTGTTCTGGTTCCTCATTATCCTGATGATGTTCCAGATTCCGAGTTTTTTCACGTTGAATAGCTGCTTctgct                       >  4:795799/1‑147 (MQ=2)
                                                                                                                         aGCCGGTTGCGGTGCCCGCTTGCCGGCTTTGCCGCGGGCAATCGCGGCGTCGCGGGGGGCTTTCCGCGGTTCTACCTGTTCTTCTGGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTGGCACGGGCAATAGCGGCTTCGAcg                       <  7:348149/150‑1 (MQ=11)
                                                                                                                                ggCTGGTCCAGCTTGCGTGCTTTGGCGCGGGCAATAGCGGCGTCGACGGCTGCTTCGCGCGGAGCGCCCTGGTGTGCTGGGGGCGCCGTAGCGTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCtt                <  1:326636/150‑1 (MQ=14)
                                                                                                                                agctgtgcgctgctGCGCGCCTTGGCGCGGGCAATAGCGGCTTCTACGGCGGCTTTGCGCGGATCTACCTGTTCTTCTGGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCt                 <  3:295496/137‑1 (MQ=11)
                                                                                                                                 gCTGTTCCAGCTTGCGTGCTTTGGCACGGGTAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACATGTTCTTCTGGTTCCGCATTAGCCTGTTGCTGTTCCAGATTGCGCGCTATGGCCCGGGGAATCCCGGCATCAACGGCGGGttt               >  5:712990/1‑151 (MQ=14)
                                                                                                                                 cccgTTCCAGCTTTTTTTTTTTGCCACGGGCAATAGCGTCTTCTACCGAGGCTTTGAGAGTATTTACCTCCTCTACTCTCTTCCGATCAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCttt               <  2:694290/148‑1 (MQ=2)
                                                                                                                                    gagggccgacgCGTGCTCTGGCGCGCTCAATAGCCGCGGCCACTGCCGCTTTGCGCGGATCCACCTGTTCTTCTGGTACCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCCCGGGCAATAGCGGCTTCGACGGCGGCTTTgcg            <  1:829590/141‑1 (MQ=11)
                                                                                                                                        cAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGTCCCGGGCAATAGCTGCTTCGACGGCTGGTTTTCGAggt        >  7:18220/1‑150 (MQ=11)
                                                                                                                                         aGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTCTTGTTCTGCTTCCGGATTAGACTGTTGCTGTGCCAGCTTGCGCCCTTTCCCTCGGGCGATAGCGGCTTCGACGGCGGTTTTtcgttgat       >  3:913983/1‑146 (MQ=11)
                                                                                                                                            ttGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTCTTCTGGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCTGCTTCGACGGCTGGTTTgcggcgctcg     >  3:517905/1‑143 (MQ=14)
                                                                                                                                               cGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTCTTCTGGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGATTTTGCACGTGAAATAGCTGCTTCTCCTGCAGTTTTTCGAGGATctccc  >  8:502999/1‑146 (MQ=11)
                                                                                                                                                           |                                                                                                                                            
CTGTTCTTCTGGCTCCGCATTAGCCGGTTGCTGTTCCCGTTTGCGCGCTTTGGCTCGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTGTTCTGGCTCGG‑‑CATTAGCCTGTTGCTGTTCCAGCTTGCGTGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACCTGTTCTTCTGGTTCCGCATTAGCCTGTTGCTGTTCCAGCTTGCGCGCTTTGGCACGGGCAATAGCGGCTTCGACGGCGGCTTTGCGCGGATCGACC  >  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: