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.

CGGCTGTCGCCATTACGCATCAACAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAGGTTCAACAATTGGCGTAA  >  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)
                    cAACAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑g                                                                                                                                  <  8:255270/150‑1 (MQ=255)
                     aaCAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑gcg                                                                                                                                <  3:999724/151‑1 (MQ=255)
                          cGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCa                                                                                                                           <  7:872786/151‑1 (MQ=255)
                          cGGAGGAAGCAGGGCACCCTGCGGTTGTGCCGCAACTGGCTGTTGCGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCa                                                                                                                           <  6:175588/151‑1 (MQ=255)
                           ggATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCaa                                                                                                                          <  6:328755/151‑1 (MQ=255)
                              tGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGATGCTGTGGGGCAACCGGTTGTTCCGGCAGTATAGATAGTGGCTGGCGgcgcccgg                                                                                                                        >  3:271165/1‑147 (MQ=18)
                               gAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCg                                                                                                                       <  5:1088690/150‑1 (MQ=255)
                               gAAGCAGGGTATCCTGCGGGTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCg                                                                                                                       <  5:600435/150‑1 (MQ=255)
                                  gcggaggcacttccGGATGAGGCGCAAACGGTTGATGCGGCTGCGGATCGTGGGGCTGCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCT‑GCGGTGCAACCGGTtg                                                                                                                   <  1:748943/138‑1 (MQ=18)
                                   cAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCtgt                                                                                                                  <  4:621974/151‑1 (MQ=255)
                                    aGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACAGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCACCCGGCtgat                                                                                                                 >  2:1150834/1‑149 (MQ=255)
                                       gTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTtgtg                                                                                                              >  8:80603/1‑151 (MQ=255)
                                         aTCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATAATGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGAAAATGTGCCT‑GCGGCGCATGCGTGTGTtgtgcg                                                                                                            >  5:753102/1‑149 (MQ=18)
                                           ccTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGgctg                                                                                                          <  1:1012696/151‑1 (MQ=255)
                                           ccTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGATGAGGAGCAACAGGATCTTGAGGCTGCTGATCATCTTCCT‑GCCGCGCAATCTGCtgttgt                                                                                                               >  5:269101/1‑146 (MQ=17)
                                           ccTGCGGTTGTGGCGCAACTGGCTGTTGGGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGgct                                                                                                           <  4:788113/150‑1 (MQ=255)
                                              gCGGGTGTGGCGCATCTGGCTGTGGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCGGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGctgctg                                                                                                       <  5:886833/151‑1 (MQ=255)
                                                ggTTGTGGCGCATGGGGCTGTTGCTGCTTATGATACTGCTGCTGCGGCGGAAGCGGTTGTTGCGGCTGCTGATTTTGCTGATTTGGCTCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGAt                                                                                                     <  6:1183034/151‑1 (MQ=21)
                                                ggTTGTGGCGCAACTGGCTGTTGTCGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCATCCGGCTGTGGTGGCTGCTGAt                                                                                                     >  6:1088675/1‑151 (MQ=255)
                                                    gttgcggacggggCTGTGGTGGCTTCTGATCCTGCGGCTGCGGCGGACGCGGTTGTTGCGGCGGCTGATATTGGGGCTGTGGCGCAACCGGCTGTTGCTGCTGCTGATAGTGGGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATAct                                                                                                  <  3:473998/139‑1 (MQ=255)
                                                    gtgGCGCTTCTGTTTTTTGTTGCTTATGATTATTCTGCTGCGGCGTATTCGGTTTTTGCTTCTGATGATTTTGGTGATGTGGCTCAACCTGTTGTTGAGGCTGCTTATATTGTGGGT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATAct                                                                                                  <  4:974768/150‑1 (MQ=255)
                                                    ggggCGCACGTGGCTGGTGGGGCTGCTGATCCTGCTGCTGCGGCGCACGCGGTTGATGCTGCTGCGGAGATTGCGGCTGTGGCGGAACCGGGTGTTGCGGCTGCGGATAGTGGGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATActg                                                                                                 <  7:1180117/149‑1 (MQ=255)
                                                         gcAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGAGGAGCAACAGGTTGTTGCGGCTGGTGAAAATGTTCCT‑GCCTCGCAATCGGCTGAGGTGCCTGCTGATaatcagtc                                                                                             >  3:636510/1‑143 (MQ=17)
                                                           atgtGGCTGTTGTTGCTTATGCTTCGTCTGCTGCGGCGTCTTCGGTTTTTGCTTCTGATGTTTTTGGTGATTTGGCTTAACCTGTTGTTGATGCTGGTTATATTGTGGTT‑GTGGCGCAACTGGCTGTTGTGGCTGCTGATActgcggctg                                                                                           <  2:46715/147‑1 (MQ=9)
                                                           aagtGGCTGTTGGGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATActgcggctgc                                                                                          <  5:1160187/148‑1 (MQ=255)
                                                             tttttgtttgtGG‑TGGTGATTCTGATGTTGCGGCGTATTCGGTTTTTGCTGCTGAGGTTTTTGGTGATTTGGCTTAAACTGCTGTTGATGCTGGTTATATTGTGGTT‑GCGGCGCAACTGGCTGTTGTGGCTGCTGATActgcggctgcgg                                                                                        <  8:563867/143‑1 (MQ=9)
                                                             ggggCGGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGAGATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATActgcggctgcgg                                                                                        <  7:323495/149‑1 (MQ=255)
                                                              tGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATAATGCTGCTGAGGAGCAACCGGTTGTTGCTGCTGCTGATAATGTTCCT‑GAGGCGCATACATCTGGTGTGCGTGATGATAatgctgctgcaac                                                                                       >  3:1282408/1‑148 (MQ=14)
                                                               ggCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTTGCTGTCGAGCAAACGGTTCTTGCGCCTGCTGAAAATGTGCCT‑GCGCCGCAAACAGCTGAGGTGCCTCATGATActccggctccaccg                                                                                      >  5:1332887/1‑145 (MQ=12)
                                                                gCTGTGGTGGCTGCTGCTCCGGCTGCTGCGGCGCCCCCGGTTGTTGCGGCGGCTGCTAGTGCGGCTGTGGCGCAACCGGGTGTTGCGGCTGCTGATACTGGGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGcggcgg                                                                                     <  5:579516/151‑1 (MQ=37)
                                                                 cTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTCGGGCAACCGGTTGTTGCGGGTGCTGAAAATGTTCCT‑GCGGCGCAAGCAGCTGCTGTTCGTGCTGATCCTCCGGCTCccgcgg                                                                                     >  3:859634/1‑150 (MQ=17)
                                                                      gtgGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGAGGCGCAACCGGATGTTGCGCCTGCTGAGAATGTTCCT‑GACGCGCAATCTGATGCTGAGGCTGCTGATCCTGCGCCAGCGGCGGaccacg                                                                               >  7:489625/1‑146 (MQ=17)
                                                                           tgctgATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCTGCGCACCCGGCTGTTGTGGCTGCTGATACTGAGGCTGCGGCGGAACTGcgtgtt                                                                          >  3:886268/1‑151 (MQ=255)
                                                                            gctgctcctgaggctgcGGCGTAAGCGGTTGTTGCTGCTGCTGATTTTCCGGCTGTGGCTGAACCGGTTGTTGCTGCTGCTGATATTGCGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGgttgt                                                                           <  1:871507/137‑1 (MQ=25)
                                                                             atgATACTGATGCTGCGGCGGATGCGGTTGTTGCGGCTGATGTTATTGGTGCTGTGGCGCAACCTGTTGTTGATGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGttgttg                                                                         <  8:728512/149‑1 (MQ=37)
                                                                               gATACTGCTGCTGCGGCGCAAGCGGTTTTTGCTGCGGCGGAAATTCCGGCTGAGGCGGAACCGGTTGTTGAGGCTGCTGATACTGCGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCg                                                                       <  1:618324/150‑1 (MQ=17)
                                                                               gATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGAGCACCCGGTTGTTGCGGCTGCTGAAAATGTGCCT‑GCGGCGCACGCGGCTGGTGTTGCTGCGGAGACTGCTGCTGCGCCGGAACGCGTTGTTGCgg                                                                      >  7:392520/1‑151 (MQ=18)
                                                                               gATACTGCGGCTGCGGCGCAACCGGGTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACGGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCgg                                                                      <  5:152537/151‑1 (MQ=255)
                                                                               gATACTCCTGCTGCGGCGCAACCGGTTGATGCGGCTGCGGATATTGGTGATGTGGCGCAACCGGTTGTTGCGGCTGCGGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGttgttg                                                                         <  4:875901/148‑1 (MQ=21)
                                                                                attctaaggctgcGGCGTATTCGGTTTTTGCTTCTGCTGTTTTTGGGGATTTGGCTGAACCTGTTGTTGATGCTGGGTATAGTGCGGTT‑GTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCgg                                                                      <  6:1141668/142‑1 (MQ=2)
                                                                                aTACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATCCTGCGGCTGCGGCGGACCTGGTTGTTGCGgc                                                                     >  7:985550/1‑151 (MQ=255)
                                                                                  cctgctgctgGGGCTCAACCGGCTTTTGCTTCTGGTGTTTTTTCTGATGTGGCTTAACCTGGTGTTGATGCTGGTGATACTGCGGCT‑GTGGCGCAACAGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGgct                                                                    <  4:1327351/149‑1 (MQ=12)
                                                                                  aCTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGAGGCGCACCCGGTTGGTGCGGCTGCTGAGCTTGTTGCT‑GCGGCGCAAGCGTCTGCGGTGGCTGCTGATCATCCGGCTCCCCCGGATGCGGTGGTGGCGgctg                                                                   >  5:1018519/1‑151 (MQ=12)
                                                                                    tgctgctgcGGCGCAACCGGTTGTTGCGGCTGCTGATATTGGTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATAATGTTGCT‑GCGGCGCAAGCGTCTGGGGTGGCTGCTGATACTGCGGCTGCGGCGGATGTGGTGGATGCGggtgct                                                                 >  3:1233864/1‑151 (MQ=18)
                                                                                     gctgctgcGGCGCAAGCGGTTGATGCGGCTGCTGATATTGCCGATGTGGCGGAACCGGGTGTTGCTGCTGCTGATATTGCGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGgctgc                                                                  <  1:344891/149‑1 (MQ=37)
                                                                                     gcggctgcGGCGGATGCGGTTGTTGCTGCGTCTGATTTTGGGGCTGTGGCTCAACCGGTTGTTGCTGCTGGTGATATTGCGGTT‑GAGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGctgctg                                                                <  6:23841/148‑1 (MQ=18)
                                                                                       tgctgcGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTCCTGTGGAGCAACAGGATGTTGCGGCTGCTGAACTTCTTGCT‑GAAGAGCAAGCAGCTGTTGTGTGTGGTGATACTCCGCCTGCGACGGAAGTGcttgttcacgatcaggat                                                              >  7:82795/1‑139 (MQ=12)
                                                                                        gctgcGGCGCAACCGGTTGTTGCTGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGata                                                             <  8:536835/151‑1 (MQ=255)
                                                                                           ggggcgCACGCGGTTGTTGCGGCTGCTGATTTTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGGGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGatat                                                            <  1:468061/147‑1 (MQ=255)
                                                                                           gcggcgGAGGCGGTTGTTGCTGCTGCTGATATTGCTGATGTGGCTTAACCTGTTGTTGAGGCTGCGGATATTGCGGTT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATAtt                                                           <  8:342940/150‑1 (MQ=21)
                                                                                           gcggcgCATCCGGTTGTTGCGGCTGCTGATAGTCCGGATGAGGCGCAACCTGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATtg                                                          <  8:436924/151‑1 (MQ=21)
                                                                                                gcAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGGTGTGGCTGCTGATCCTGCGGCTGCGGCGGATGTGGTTGTTGCGGCTGCTGATATtgctgcg                                                     >  1:501412/1‑150 (MQ=255)
                                                                                                 gatgcGGTTGTTGCTGCTGATGAAGTTGCTGCTGTGGCGCCTTCTGTTGTTGCTGCTGCTGATATTGTGGCT‑GAGGCGCAATCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGCAACTGGTTGTTGCGGCTGCTGATATtgctgctg                                                    <  4:75630/147‑1 (MQ=18)
                                                                                                  cgcgggTTGATGCGTCTGCTGCCTTTGGGGCTGAGGCGGCACCGGTTGATGCGGCTGCGGATAGTGGGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATtgctgctg                                                    <  3:752953/146‑1 (MQ=18)
                                                                                                    gcGGTTGTTGCAGCTTATGATATTCCGGATGTGGCGGAACCGGTTGTTGCTGCTGGTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCgg                                                 <  1:147328/150‑1 (MQ=21)
                                                                                                       gaggtggccgctgatGATTTTCCTGCTGTGGCGGAAGCGGTTGATGCTGCTGCGGATATTGGGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGc                                              <  8:606989/142‑1 (MQ=21)
                                                                                                                gctgcgGATTTTGCGGCTGTGGCGGATCCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGgttg                                     <  8:947120/151‑1 (MQ=255)
                                                                                                                   gctgATATTGCGGCTGCGGCGGAACCGGTTGATGCTGCTGCTGATAGTGGGGAT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGgttgtt                                   <  1:370155/150‑1 (MQ=21)
                                                                                                                        tattcgtgctgTGGCGGAATCGGTTGTTGCTTCTGCTGATTGTGTGGCT‑GAGGCTCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGgc                              <  8:975154/144‑1 (MQ=18)
                                                                                                                        tatTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGCCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCTGCAGCGGCGCAACAGTTTGTTGCGGCTGCTGATAATCCTGCTGCGGAGCAAACGGTTGTAGCgtat                             >  2:702034/1‑148 (MQ=18)
                                                                                                                         tttGCCGCTGTGGCTGCACCTGTTGTTGCGGCTGCGGATAGTGGGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGgct                             <  7:1204054/149‑1 (MQ=18)
                                                                                                                          ttGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GAGGAGCAAGCGGCTGTTGTGGCTGCTGAGACTGCGGCTGCGGCGGAAGGGGTTGTTGCGGCTGATGAGATTGCTGCTGCGCTGCAGGCCGTTGTTGCGgctgc                           >  3:1305025/1‑151 (MQ=18)
                                                                                                                           tgctgcTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGTAGCAGGCGGCTGGTGTGGCTGGTGATACTGCGGCTGCGGCGGAAGTGGTTGTTGAGGCTGCTGATATTGCTGCTGCGGTGCAAGCGGTTGTTGCGggtgct                          >  3:1049813/1‑151 (MQ=21)
                                                                                                                             ctgctgTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTCCTGCTGCGGCGCAAATGGTTGTTGCTGCTGCTGATATTCCTGCTGCCGTGCAAACGGTTCTTGCGGATGCTGt                        >  4:144544/1‑151 (MQ=18)
                                                                                                                               gctgtcgctgcatctgttgttGCTGCTGAGGATTTTGGGGTT‑TTGGCGCAACGGGCGGTTGATGCTGGGGATACTGCGGTTGAGGCGGAACTGGGTGTTGTGGCTGCTGATAGTGCTGCTGCGGGGCAACCGGTTGTTGCGGCTGCTGTa                       <  8:216279/135‑1 (MQ=2)
                                                                                                                               gctgTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCATCCGGTTGTTGCGGCTGCTGTAc                      >  6:437889/1‑151 (MQ=255)
                                                                                                                               gctgTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGATGCTGCAACGGAACAAGTTGTTGAGGCTGATGATATTGCTGCTGCGGTGCAAACTGTTGTTGAGGCTGCTGTAc                      >  2:281713/1‑151 (MQ=18)
                                                                                                                               gctgTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGAT‑GAGGAGCAACCGGCTCATGTGCCTGATGATAATGCGGCTGCCGCGCAACGAGTTGTTGAGGATCATGATATTCATGCTCAGGAGCAACCGCTTGATCAGGATGCTGAAc                      >  4:391353/1‑151 (MQ=12)
                                                                                                                                 tgtgGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAAACGGTTGTTGCGGCTGATGTAc                      >  8:1313605/1‑149 (MQ=255)
                                                                                                                                 tgtgGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTAGTTGCGGCTGCTGTAc                      >  8:191183/1‑149 (MQ=255)
                                                                                                                                  ttggctgaaCCGGTTGATGAGGCTGGTGATATTGTGGAT‑GTGGCGCAACGTGCTGTTGTGGCTGGTGATACTGCGGGTGCGGCGCAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAg                    <  6:977896/143‑1 (MQ=17)
                                                                                                                                  ttcgcgcAACCGGTTGTTGCGGCTGCTGATAGTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAg                    <  3:203276/147‑1 (MQ=255)
                                                                                                                                   tgGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGGTGGGGATACTGCTGCTGCGGCGGAAGCGGTTGTTGCGCGTGCTGATATTGCTGCTGCAGTGCATCCGGTTGTTGCGGGTGCTGTACATGt                  >  7:1294106/1‑151 (MQ=21)
                                                                                                                                   tgGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAGGt                  >  7:307360/1‑151 (MQ=255)
                                                                                                                                        gcccctgctgatgCCGCTGCGGATCGTGGGGAT‑GTGGCGCAACGGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAGGTTcaac             <  4:1142211/140‑1 (MQ=17)
                                                                                                                                         cagctgttgttGCCGCTGCTGATATTGGGGCT‑GAGGCGGAACGGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAGGTTcaac             <  2:405632/145‑1 (MQ=21)
                                                                                                                                          aCCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGAGCAACCGGCTGTTGTGGCTGCTGATACTGAGGCTGCGGCGGAACTGGTTGTTGAGGCTGCTGATATTGCTGCTGCGGTGCAAGCGGTTGTTCCGGCTGCTGTACAGGTTcaaca            >  7:582403/1‑150 (MQ=255)
                                                                                                                                            cGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATCATGAGCCTGCGGCGGAGCTGGTTGTTGAGCCTGTTGATAATCCTGCTGCGGGGCAACCGGTTGTTGCCGCTGATGATCAGGATCAACTAtt         >  5:962965/1‑151 (MQ=18)
                                                                                                                                              gttgttGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGATGCTGCTGATATTGCTGCTGCGGTGCAAGCGGTTGTTGCGGGTGCTGTCCAGGTTcacctagtgg       >  5:307295/1‑143 (MQ=21)
                                                                                                                                              gctgatgctgcTGTTGATAGTGGGGCT‑TTGGCGTAACCTGCTGTTGTTGCTGCTGATACTGCGGGTGCGGCGGAACTGGTTGTTGAGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAGGTTCAACAATTg        <  6:47430/145‑1 (MQ=18)
                                                                                                                                               gtgttgtGGCTGCTGATAGTGGTGCT‑GAGGCTTAACCTGCTGTTGTGGCTGCTGATATTGCTGCTGTGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGTACAGGTTCAACAATTGGc      <  8:1252136/150‑1 (MQ=17)
                                                                                                                                                 gttgCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAAGTGGTTGTTGCGGCTGCTGATATTGCTCCTGCAGTGCAAACGGTTGTTGAGGCTGATGTACAGGATCAACAATTGGAGt    >  6:1127482/1‑151 (MQ=25)
                                                                                                                                                   tgctgctgcTGATAGTGCGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCAGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAGGTTCAACAATTGGCGTa   <  7:244558/150‑1 (MQ=255)
                                                                                                                                                    gCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAGGTTCAACAATTGGCGTaa  >  8:1129727/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CGGCTGTCGCCATTACGCATCAACAGCGGATGAAGCAGGGTATCCTGCGGTTGTGGCGCAACTGGCTGTTGTGGCTGCTGATACTGCTGCTGCGGCGCAACCGGTTGTTGCGGCTGCTGATATTGCTGCTGTGGCGCAACCGGTTGTTGCGGCTGCTGATATTGTGGCT‑GCGGCGCAACCGGCTGTTGTGGCTGCTGATACTGCGGCTGCGGCGGAACTGGTTGTTGCGGCTGCTGATATTGCTGCTGCGGTGCAACCGGTTGTTGCGGCTGCTGTACAGGTTCAACAATTGGCGTAA  >  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: ‑