Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,945,5540GT21.5% 450.2 / 17.1 220P781Q (CCG→CAG) ftsKcell division protein FtsK
Reads supporting (aligned to +/- strand):  ref base G (81/87);  new base T (6/40);  total (88/132)
Fisher's exact test for biased strand distribution p-value = 8.27e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 229 total aligned reads displayed.

CGGCTGTCGCCATTACGCATCAACAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAGGTTCAACAATTGGCGTAA  >  DH1/2945404‑2945701
                                                                                                                                                      |                                                                                                                                                    
cggCTGTCGCCATTACGCATCAACAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGGGCCGAACCCGGGTGTTCCg                                                                                                                                                      >  3:1144677/1‑151 (MQ=25)
 ggCTGTCGCCATTACGCATCAACAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCgg                                                                                                                                                     <  7:929896/151‑1 (MQ=255)
  gCTGTCGCCATTACGCATCAACAGTGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGgc                                                                                                                                                    <  1:53515/151‑1 (MQ=255)
    tGTCGCCATTACGCATCAACAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATAGTGCTGCTGAGGCGCACCCGGTTGTTGCGgctc                                                                                                                                                  >  1:1087373/1‑150 (MQ=255)
     gTCGCCATTACGCATCAACAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGgct                                                                                                                                                   <  6:131060/149‑1 (MQ=255)
     gTCGCCATTACGCATCAACAGCGGATGAAGCAGGGTATCCTGCGGTTGGGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGgct                                                                                                                                                   <  2:989039/149‑1 (MQ=255)
      tCGCCATTACGCATCAACAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGAGCAACCGGTTGTTGCGgctgct                                                                                                                                                >  7:1252798/1‑151 (MQ=255)
       cGCCATTACGCATCAACAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGgctgct                                                                                                                                                <  4:1287167/150‑1 (MQ=255)
       cGCCATTACGCATCAACAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGctgctg                                                                                                                                               <  5:666558/151‑1 (MQ=255)
        gCCATTACGCATCAACAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGa                                                                                                                                              <  8:1063136/151‑1 (MQ=255)
          cATTACGCATCAACAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGATGTGGAGCAACAGGTTGTTGCGGCTGCTGaa                                                                                                                                             >  1:127854/1‑149 (MQ=255)
           aTTACGCATCAACAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGAGCAACCGGTTCTTGCGGCTGCTGATAATGCTGCAGAGGCGCCCCAGGTTCTTGCGCCTGCTGaca                                                                                                                                            >  2:1079250/1‑148 (MQ=25)
             tACGCATCAACAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGGGGATATTGCTGCTGTGGGGCAACCCGTAGGTGCggagcgtgtcatat                                                                                                                                         >  1:1098801/1‑139 (MQ=21)
               cGCATCAACAGGGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATtgtg                                                                                                                                       <  3:543227/151‑1 (MQ=255)
               cGCATCAACAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATtgtg                                                                                                                                       <  4:547222/151‑1 (MQ=255)
                  aTCAACAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGAGCAACAGGTTGTTGAGgctcatctaaatgttcat                                                                                                                                    >  3:314780/1‑136 (MQ=21)
                    cAACAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTg                                                                                                                                  <  8:255270/150‑1 (MQ=255)
                     aaCAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTgcg                                                                                                                                <  3:999724/151‑1 (MQ=255)
                          cGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCa                                                                                                                           <  7:872786/151‑1 (MQ=255)
                          cGGAGGAAGCAGGGCACCCTGCGGTTGTGCCGCAACTGGCTGTTGCGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCa                                                                                                                           <  6:175588/151‑1 (MQ=255)
                           ggATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCaa                                                                                                                          <  6:328755/151‑1 (MQ=255)
                              tGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGATGCTGTGGGGCAACCGGTTGTTCCGGCAGTATAGATAGTGGCTGGCGgcgcccgg                                                                                                                        >  3:271165/1‑147 (MQ=18)
                               gAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCg                                                                                                                       <  5:1088690/150‑1 (MQ=255)
                               gAAGCAGGGTATCCTGCGGGTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCg                                                                                                                       <  5:600435/150‑1 (MQ=255)
                                  gcggaggcacttccGGATGAGGCGCAAACGGTTGATGCGGCTGCGGATCGTGGGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTtg                                                                                                                   <  1:748943/138‑1 (MQ=18)
                                   cAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCtgt                                                                                                                  <  4:621974/151‑1 (MQ=255)
                                    aGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACAGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCACCCGGCtgat                                                                                                                 >  2:1150834/1‑149 (MQ=255)
                                       gTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTtgtg                                                                                                              >  8:80603/1‑151 (MQ=255)
                                         aTCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATAATGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGAAAATGTGCCTGCGGCGCATGCGTGTGTtgtgcg                                                                                                            >  5:753102/1‑149 (MQ=18)
                                           ccTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGgctg                                                                                                          <  1:1012696/151‑1 (MQ=255)
                                           ccTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGATGAGGAGCAACAGGATCTTGAGGCTGCTGATCATCTTCCTGCCGCGCAATCTGCtgttgt                                                                                                               >  5:269101/1‑146 (MQ=17)
                                           ccTGCGGTTGTGGCGCAACTGGCTGTTGGGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGgct                                                                                                           <  4:788113/150‑1 (MQ=255)
                                              gCGGGTGTGGCGCATCTGGCTGTGGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCGGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGctgctg                                                                                                       <  5:886833/151‑1 (MQ=255)
                                                ggTTGTGGCGCATGGGGCTGTTGCTGCTTATGATACTGCTGCTGCGGCGGAAGCGGTTGTTGCGGCTGCTGATTTTGCTGATTTGGCTCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGAt                                                                                                     <  6:1183034/151‑1 (MQ=21)
                                                ggTTGTGGCGCAACTGGCTGTTGTCGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCATCCGGCTGTGGTGGCTGCTGAt                                                                                                     >  6:1088675/1‑151 (MQ=255)
                                                    gttgcggacggggCTGTGGTGGCTTCTGATCCTGCGGCTGCGGCGGACGCGGTTGTTGCGGCGGCTGATATTGGGGCTGTGGCGCAACCGGCTGTTGCTGCTGCTGATAGTGGGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATAct                                                                                                  <  3:473998/139‑1 (MQ=255)
                                                    gtgGCGCTTCTGTTTTTTGTTGCTTATGATTATTCTGCTGCGGCGTATTCGGTTTTTGCTTCTGATGATTTTGGTGATGTGGCTCAACCTGTTGTTGAGGCTGCTTATATTGTGGGTGCGGCGCAACCGGCTGTTGTGGCTGCTGATAct                                                                                                  <  4:974768/150‑1 (MQ=255)
                                                    ggggCGCACGTGGCTGGTGGGGCTGCTGATCCTGCTGCTGCGGCGCACGCGGTTGATGCTGCTGCGGAGATTGCGGCTGTGGCGGAACCGGGTGTTGCGGCTGCGGATAGTGGGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATActg                                                                                                 <  7:1180117/149‑1 (MQ=255)
                                                         gcAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGAGGAGCAACAGGTTGTTGCGGCTGGTGAAAATGTTCCTGCCTCGCAATCGGCTGAGGTGCCTGCTGATaatcagtc                                                                                             >  3:636510/1‑143 (MQ=17)
                                                           atgtGGCTGTTGTTGCTTATGCTTCGTCTGCTGCGGCGTCTTCGGTTTTTGCTTCTGATGTTTTTGGTGATTTGGCTTAACCTGTTGTTGATGCTGGTTATATTGTGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATActgcggctg                                                                                           <  2:46715/147‑1 (MQ=9)
                                                           aagtGGCTGTTGGGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATActgcggctgc                                                                                          <  5:1160187/148‑1 (MQ=255)
                                                             tttttgtttgtGGTGGTGATTCTGATGTTGCGGCGTATTCGGTTTTTGCTGCTGAGGTTTTTGGTGATTTGGCTTAAACTGCTGTTGATGCTGGTTATATTGTGGTTGCGGCGCAACTGGCTGTTGTGGCTGCTGATActgcggctgcgg                                                                                        <  8:563867/143‑1 (MQ=9)
                                                             ggggCGGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGAGATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATActgcggctgcgg                                                                                        <  7:323495/149‑1 (MQ=255)
                                                              tGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATAATGCTGCTGAGGAGCAACCGGTTGTTGCTGCTGCTGATAATGTTCCTGAGGCGCATACATCTGGTGTGCGTGATGATAatgctgctgcaac                                                                                       >  3:1282408/1‑148 (MQ=14)
                                                               ggCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTTGCTGTCGAGCAAACGGTTCTTGCGCCTGCTGAAAATGTGCCTGCGCCGCAAACAGCTGAGGTGCCTCATGATActccggctccaccg                                                                                      >  5:1332887/1‑145 (MQ=12)
                                                                gCTGTGGTGGCTGCTGCTCCGGCTGCTGCGGCGCCCCCGGTTGTTGCGGCGGCTGCTAGTGCGGCTGTGGCGCAACCGGGTGTTGCGGCTGCTGATACTGGGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGcggcgg                                                                                     <  5:579516/151‑1 (MQ=37)
                                                                 cTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTCGGGCAACCGGTTGTTGCGGGTGCTGAAAATGTTCCTGCGGCGCAAGCAGCTGCTGTTCGTGCTGATCCTCCGGCTCccgcgg                                                                                     >  3:859634/1‑150 (MQ=17)
                                                                      gtgGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGAGGCGCAACCGGATGTTGCGCCTGCTGAGAATGTTCCTGACGCGCAATCTGATGCTGAGGCTGCTGATCCTGCGCCAGCGGCGGaccacg                                                                               >  7:489625/1‑146 (MQ=17)
                                                                           tgctgATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCTGCGCACCCGGCTGTTGTGGCTGCTGATACTGAGGCTGCGGCGGAACTGcgtgtt                                                                          >  3:886268/1‑151 (MQ=255)
                                                                            gctgctcctgaggctgcGGCGTAAGCGGTTGTTGCTGCTGCTGATTTTCCGGCTGTGGCTGAACCGGTTGTTGCTGCTGCTGATATTGCGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGgttgt                                                                           <  1:871507/137‑1 (MQ=25)
                                                                             atgATACTGATGCTGCGGCGGATGCGGTTGTTGCGGCTGATGTTATTGGTGCTGTGGCGCAACCTGTTGTTGATGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGttgttg                                                                         <  8:728512/149‑1 (MQ=37)
                                                                               gATACTGCTGCTGCGGCGCAAGCGGTTTTTGCTGCGGCGGAAATTCCGGCTGAGGCGGAACCGGTTGTTGAGGCTGCTGATACTGCGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCg                                                                       <  1:618324/150‑1 (MQ=17)
                                                                               gATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGAGCACCCGGTTGTTGCGGCTGCTGAAAATGTGCCTGCGGCGCACGCGGCTGGTGTTGCTGCGGAGACTGCTGCTGCGCCGGAACGCGTTGTTGCgg                                                                      >  7:392520/1‑151 (MQ=18)
                                                                               gATACTGCGGCTGCGGCGCAACCGGGTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACGGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCgg                                                                      <  5:152537/151‑1 (MQ=255)
                                                                               gATACTCCTGCTGCGGCGCAACCGGTTGATGCGGCTGCGGATATTGGTGATGTGGCGCAACCGGTTGTTGCGGCTGCGGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGttgttg                                                                         <  4:875901/148‑1 (MQ=21)
                                                                                attctaaggctgcGGCGTATTCGGTTTTTGCTTCTGCTGTTTTTGGGGATTTGGCTGAACCTGTTGTTGATGCTGGGTATAGTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCgg                                                                      <  6:1141668/142‑1 (MQ=2)
                                                                                aTACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATCCTGCGGCTGCGGCGGACCTGGTTGTTGCGgc                                                                     >  7:985550/1‑151 (MQ=255)
                                                                                  cctgctgctgGGGCTCAACCGGCTTTTGCTTCTGGTGTTTTTTCTGATGTGGCTTAACCTGGTGTTGATGCTGGTGATACTGCGGCTGTGGCGCAACAGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGgct                                                                    <  4:1327351/149‑1 (MQ=12)
                                                                                  aCTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGAGGCGCACCCGGTTGGTGCGGCTGCTGAGCTTGTTGCTGCGGCGCAAGCGTCTGCGGTGGCTGCTGATCATCCGGCTCCCCCGGATGCGGTGGTGGCGgctg                                                                   >  5:1018519/1‑151 (MQ=12)
                                                                                    tgctgctgcGGCGCAACCGGTTGTTGCGGCTGCTGATATTGGTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATAATGTTGCTGCGGCGCAAGCGTCTGGGGTGGCTGCTGATACTGCGGCTGCGGCGGATGTGGTGGATGCGggtgct                                                                 >  3:1233864/1‑151 (MQ=18)
                                                                                     gctgctgcGGCGCAAGCGGTTGATGCGGCTGCTGATATTGCCGATGTGGCGGAACCGGGTGTTGCTGCTGCTGATATTGCGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGgctgc                                                                  <  1:344891/149‑1 (MQ=37)
                                                                                     gcggctgcGGCGGATGCGGTTGTTGCTGCGTCTGATTTTGGGGCTGTGGCTCAACCGGTTGTTGCTGCTGGTGATATTGCGGTTGAGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGctgctg                                                                <  6:23841/148‑1 (MQ=18)
                                                                                       tgctgcGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTCCTGTGGAGCAACAGGATGTTGCGGCTGCTGAACTTCTTGCTGAAGAGCAAGCAGCTGTTGTGTGTGGTGATACTCCGCCTGCGACGGAAGTGcttgttcacgatcaggat                                                              >  7:82795/1‑139 (MQ=12)
                                                                                        gctgcGGCGCAACCGGTTGTTGCTGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGata                                                             <  8:536835/151‑1 (MQ=255)
                                                                                           ggggcgCACGCGGTTGTTGCGGCTGCTGATTTTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGGGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGatat                                                            <  1:468061/147‑1 (MQ=255)
                                                                                           gcggcgGAGGCGGTTGTTGCTGCTGCTGATATTGCTGATGTGGCTTAACCTGTTGTTGAGGCTGCGGATATTGCGGTTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATAtt                                                           <  8:342940/150‑1 (MQ=21)
                                                                                           gcggcgCATCCGGTTGTTGCGGCTGCTGATAGTCCGGATGAGGCGCAACCTGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATtg                                                          <  8:436924/151‑1 (MQ=21)
                                                                                                gcAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGGTGTGGCTGCTGATCCTGCGGCTGCGGCGGATGTGGTTGTTGCGGCTGCTGATATtgctgcg                                                     >  1:501412/1‑150 (MQ=255)
                                                                                                 gatgcGGTTGTTGCTGCTGATGAAGTTGCTGCTGTGGCGCCTTCTGTTGTTGCTGCTGCTGATATTGTGGCTGAGGCGCAATCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGCAACTGGTTGTTGCGGCTGCTGATATtgctgctg                                                    <  4:75630/147‑1 (MQ=18)
                                                                                                  cgcgggTTGATGCGTCTGCTGCCTTTGGGGCTGAGGCGGCACCGGTTGATGCGGCTGCGGATAGTGGGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATtgctgctg                                                    <  3:752953/146‑1 (MQ=18)
                                                                                                    gcGGTTGTTGCAGCTTATGATATTCCGGATGTGGCGGAACCGGTTGTTGCTGCTGGTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCgg                                                 <  1:147328/150‑1 (MQ=21)
                                                                                                       gaggtggccgctgatGATTTTCCTGCTGTGGCGGAAGCGGTTGATGCTGCTGCGGATATTGGGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGc                                              <  8:606989/142‑1 (MQ=21)
                                                                                                                gctgcgGATTTTGCGGCTGTGGCGGATCCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGgttg                                     <  8:947120/151‑1 (MQ=255)
                                                                                                                   gctgATATTGCGGCTGCGGCGGAACCGGTTGATGCTGCTGCTGATAGTGGGGATGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGgttgtt                                   <  1:370155/150‑1 (MQ=21)
                                                                                                                        tattcgtgctgTGGCGGAATCGGTTGTTGCTTCTGCTGATTGTGTGGCTGAGGCTCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGgc                              <  8:975154/144‑1 (MQ=18)
                                                                                                                        tatTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGCCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCTGCAGCGGCGCAACAGTTTGTTGCGGCTGCTGATAATCCTGCTGCGGAGCAAACGGTTGTAGCgtat                             >  2:702034/1‑148 (MQ=18)
                                                                                                                         tttGCCGCTGTGGCTGCACCTGTTGTTGCGGCTGCGGATAGTGGGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGgct                             <  7:1204054/149‑1 (MQ=18)
                                                                                                                          ttGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGAGGAGCAAGCGGCTGTTGTGGCTGCTGAGACTGCGGCTGCGGCGGAAGGGGTTGTTGCGGCTGATGAGATTGCTGCTGCGCTGCAGGCCGTTGTTGCGgctgc                           >  3:1305025/1‑151 (MQ=18)
                                                                                                                           tgctgcTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGTAGCAGGCGGCTGGTGTGGCTGGTGATACTGCGGCTGCGGCGGAAGTGGTTGTTGAGGCTGCTGATATTGCTGCTGCGGTGCAAGCGGTTGTTGCGggtgct                          >  3:1049813/1‑151 (MQ=21)
                                                                                                                             ctgctgTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTCCTGCTGCGGCGCAAATGGTTGTTGCTGCTGCTGATATTCCTGCTGCCGTGCAAACGGTTCTTGCGGATGCTGt                        >  4:144544/1‑151 (MQ=18)
                                                                                                                               gctgtcgctgcatctgttgttGCTGCTGAGGATTTTGGGGTTTTGGCGCAACGGGCGGTTGATGCTGGGGATACTGCGGTTGAGGCGGAACTGGGTGTTGTGGCTGCTGATAGTGCTGCTGCGGGGCAACCGGTTGTTGCGGCTGCTGTa                       <  8:216279/135‑1 (MQ=2)
                                                                                                                               gctgTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCATCCGGTTGTTGCGGCTGCTGTAc                      >  6:437889/1‑151 (MQ=255)
                                                                                                                               gctgTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGATGCTGCAACGGAACAAGTTGTTGAGGCTGATGATATTGCTGCTGCGGTGCAAACTGTTGTTGAGGCTGCTGTAc                      >  2:281713/1‑151 (MQ=18)
                                                                                                                               gctgTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGATGAGGAGCAACCGGCTCATGTGCCTGATGATAATGCGGCTGCCGCGCAACGAGTTGTTGAGGATCATGATATTCATGCTCAGGAGCAACCGCTTGATCAGGATGCTGAAc                      >  4:391353/1‑151 (MQ=12)
                                                                                                                                 tgtgGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAAACGGTTGTTGCGGCTGATGTAc                      >  8:1313605/1‑149 (MQ=255)
                                                                                                                                 tgtgGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTAGTTGCGGCTGCTGTAc                      >  8:191183/1‑149 (MQ=255)
                                                                                                                                  ttggctgaaCCGGTTGATGAGGCTGGTGATATTGTGGATGTGGCGCAACGTGCTGTTGTGGCTGGTGATACTGCGGGTGCGGCGCAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAg                    <  6:977896/143‑1 (MQ=17)
                                                                                                                                  ttcgcgcAACCGGTTGTTGCGGCTGCTGATAGTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAg                    <  3:203276/147‑1 (MQ=255)
                                                                                                                                   tgGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGGTGGGGATACTGCTGCTGCGGCGGAAGCGGTTGTTGCGCGTGCTGATATTGCTGCTGCAGTGCATCCGGTTGTTGCGGGTGCTGTACATGt                  >  7:1294106/1‑151 (MQ=21)
                                                                                                                                   tgGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAGGt                  >  7:307360/1‑151 (MQ=255)
                                                                                                                                        gcccctgctgatgCCGCTGCGGATCGTGGGGATGTGGCGCAACGGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAGGTTcaac             <  4:1142211/140‑1 (MQ=17)
                                                                                                                                         cagctgttgttGCCGCTGCTGATATTGGGGCTGAGGCGGAACGGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAGGTTcaac             <  2:405632/145‑1 (MQ=21)
                                                                                                                                          aCCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGAGCAACCGGCTGTTGTGGCTGCTGATACTGAGGCTGCGGCGGAACTGGTTGTTGAGGCTGCTGATATTGCTGCTGCGGTGCAAGCGGTTGTTCCGGCTGCTGTACAGGTTcaaca            >  7:582403/1‑150 (MQ=255)
                                                                                                                                            cGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATCATGAGCCTGCGGCGGAGCTGGTTGTTGAGCCTGTTGATAATCCTGCTGCGGGGCAACCGGTTGTTGCCGCTGATGATCAGGATCAACTAtt         >  5:962965/1‑151 (MQ=18)
                                                                                                                                              gttgttGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGATGCTGCTGATATTGCTGCTGCGGTGCAAGCGGTTGTTGCGGGTGCTGTCCAGGTTcacctagtgg       >  5:307295/1‑143 (MQ=21)
                                                                                                                                              gctgatgctgcTGTTGATAGTGGGGCTTTGGCGTAACCTGCTGTTGTTGCTGCTGATACTGCGGGTGCGGCGGAACTGGTTGTTGAGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAGGTTCAACAATTg        <  6:47430/145‑1 (MQ=18)
                                                                                                                                               gtgttgtGGCTGCTGATAGTGGTGCTGAGGCTTAACCTGCTGTTGTGGCTGCTGATATTGCTGCTGTGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGTACAGGTTCAACAATTGGc      <  8:1252136/150‑1 (MQ=17)
                                                                                                                                                 gttgCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAAGTGGTTGTTGCGGCTGCTGATATTGCTCCTGCAGTGCAAACGGTTGTTGAGGCTGATGTACAGGATCAACAATTGGAGt    >  6:1127482/1‑151 (MQ=25)
                                                                                                                                                   tgctgctgcTGATAGTGCGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCAGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAGGTTCAACAATTGGCGTa   <  7:244558/150‑1 (MQ=255)
                                                                                                                                                    gCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAGGTTCAACAATTGGCGTaa  >  8:1129727/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CGGCTGTCGCCATTACGCATCAACAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAGGTTCAACAATTGGCGTAA  >  DH1/2945404‑2945701

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: