Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,716,3680GA39.3% 167.2 / 75.6 112L38F (CTC→TTC) thiKthiamine kinase
Reads supporting (aligned to +/- strand):  ref base G (65/3);  new base A (0/44);  total (65/47)
Fisher's exact test for biased strand distribution p-value = 5.41e-28
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

CATCCAGTCACGGAGATATAAATGCGGCTTCGGTGCAATGCATGCGGGTAGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGCACCTTTACCTGACCAGATAATTTCGCCCGTCTGCACCAGC  >  DH1/2716219‑2716517
                                                                                                                                                      |                                                                                                                                                      
cATCCAGTCACGGAGATATAAATGCGGCTTCGGTGCAATGCATGCGGGTAGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCAGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATgaga                                                                                                                                                       <  2:648386/151‑1 (MQ=255)
 aTCCAGTCACGGAGATATAAATGCGGCTTCGGTGCAATGCATGCGGGTAGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATgaga                                                                                                                                                       >  1:1050120/1‑150 (MQ=255)
  tCCAGTCACGGAGATATAAATGCGGCTTCGGTGCAATGCATGCGGGTAGATGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGaa                                                                                                                                                      >  2:2170532/1‑150 (MQ=255)
   ccAGTCACGGAGATATAAATGCGGCTTCGGTGCAATGCATGCGGGTAGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGaaa                                                                                                                                                     >  1:2187966/1‑150 (MQ=255)
            gagaTATAAATGCGGCTTCGGTGCAATGCATGCGGGTAGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCa                                                                                                                                          <  1:86365/151‑1 (MQ=255)
                    aaTGCGGCTTCGGTGCAATGCATGCGGGTAGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGccc                                                                                                                                  >  1:443964/1‑151 (MQ=255)
                         ggtttCGGTGCAATGAGGGGGGGTTTTTGTGAAAACCCCCGATAGTGGCGCAAAACCGCGGCCCGCGGGAAATCGGGACCGGGCGGCTACCGCAAAAAAAACCGCTGCCCTTTA‑TGTAAAATAAAAAAC‑CTCCCGCCACTAAGCCCACTAt                                                                                                                             <  2:2138285/148‑1 (MQ=255)
                            ttCGGTGCAATGCATGCGGGTAGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAg                                                                                                                                     >  2:1328147/1‑140 (MQ=255)
                             ggggggaaaTGCAGGGGGGTAGTTGTAAAAAGCCCCGATCCGGGCGAAAAAACGGGGCCGGCGGCGAATAAGGATGGGGCGGCTGACGCAAAAAAAACCGCGGGCCTGGA‑TGTAAAAAGAAAAAC‑CTCCCGCCACTAAGCCCACTATTaa                                                                                                                          <  2:1167456/143‑1 (MQ=255)
                              cGGTGCAATGCATGCGGGTAGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACt                                                                                                                               >  2:371421/1‑144 (MQ=255)
                              cGGTGCAATGCATGCGGGTAGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAAc                                                                                                                        >  2:366223/1‑151 (MQ=255)
                                gggCAATGCATGGGGGTGGTTGTAATAAAGCCCGATCCGGGCGAAAGACCGGGGCCTGCGGCGCATCAGGATCGTGCGGCTGCCCCACAACAAACCGCTGGCCTTTT‑TGTTAAATAAAAAAC‑CTCCCGCCACTAAGCCCACTATTAAAcgt                                                                                                                      <  2:2187975/149‑1 (MQ=255)
                                gTGCAATGCATGCGGGTAGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAAcgt                                                                                                                      >  1:1993112/1‑151 (MQ=255)
                                gTGAAATGCGGGCGGGTGGTTGTAAAAAGCCCCGATACGGGAGAAAGACCGCGGACGGCGGCGCATCAGGATCGTGCGGCGGACCAAAAAAAAACCGCTGGCCTTGA‑TTTAAAAAAAAAAAC‑CCCCCGCCACTAAGCCCACTATTAAAcgt                                                                                                                      <  2:1275571/151‑1 (MQ=255)
                                  gaaatgaaggggggTAGTTGTAAAAAAGCCCGATACGGGCGAAAGAACGCGGCCTGCGGCGAATCAGGATCGGGCGGCTGCCGCAAAAAAAACCGCGGGCCTTGA‑TTTTAAAGAAAAAAA‑CTCCCGCCACTAAGCCCACTATTAAAcgtc                                                                                                                     <  2:2113711/139‑1 (MQ=255)
                                    aaTGCATGCGGGTAGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTg                                                                                                                  >  1:2278448/1‑151 (MQ=255)
                                     atgcatgcGGGTAGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTg                                                                                                                  >  2:1162472/1‑150 (MQ=255)
                                      ggcaggggggTGGTTGTAAAAAACCCCAATGCGGGCGAAAGAGCGCGGCCGGCGCCGAAGCAGGAGCGGGCGGCTCACACAAAAAAAAACGCTGGCTCTGA‑TGTTAAAGAAAAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGAc                                                                                                                <  1:2037190/144‑1 (MQ=255)
                                       gcatgcGGGTAGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACg                                                                                                               >  2:1031545/1‑151 (MQ=255)
                                       gcatgcGGGTAGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACg                                                                                                               >  2:337788/1‑151 (MQ=255)
                                       gcatgcGGGTAGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACg                                                                                                               >  1:1226143/1‑151 (MQ=255)
                                        ggggcGGTTGGTTGTAAAAAACCCCGATACGGGCAAAAAAACGCGCCCTGCGGCGCATCGGAATCGGGCGGCTGCCCCAAAAAAAACCGCTGGCCTGGT‑TTTTAAATGAAAAAC‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACgg                                                                                                              <  1:572873/148‑1 (MQ=255)
                                        catgcGGGTAGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACgg                                                                                                              >  2:358121/1‑151 (MQ=255)
                                        atgggggttAGTTGTGAAAACGCCCGATGCGGGCAAAAAACCGGGGTCGGCGGCACATGAGGATGGGGGGGGTGACGCACAAAAAACCGCGGGCCGTGT‑TTTTAAAGGAAAAAC‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACgg                                                                                                              <  2:2058155/143‑1 (MQ=255)
                                          tgcGGGTAGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACggg                                                                                                             >  1:1340635/1‑150 (MQ=255)
                                          tgcGGGTAGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACggg                                                                                                             >  2:83549/1‑150 (MQ=255)
                                                 aGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACt                                                                                                     >  1:49759/1‑151 (MQ=255)
                                                  gtttaaaaaaaCCCAGATACTGGAGCAAAAAGGGGGACGGCGGCGCATCAGGTCCGGGCGGCTCACCCAAAAAAAACCGGTGGCCCTGA‑TTTAAAATAAAAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTg                                                                                                    <  1:198881/143‑1 (MQ=255)
                                                  cccaatcaaacccccGATACGGGCGAAAGACCGCGGACGGGGGCGCATCAGGTTCGGGCGGCTGACCCAAAAAAAACCGCGGCCCCTTA‑TGTTAAATGAAAAAC‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATAc                                                                                                      <  2:394493/137‑1 (MQ=255)
                                                     gtaaaaaaCCCCGATACGGGCGCAAGAGCGCGACCCGCGGCACAACAGGATCGTGCGGCTCCCCCACAAAAAACCGTTGGCCTTGT‑TGTAAAATGAAAAAC‑CCCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCg                                                                                                  <  1:750283/145‑1 (MQ=255)
                                                       gATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCggg                                                                                                >  2:1978157/1‑150 (MQ=255)
                                                        aaacaccccGATACGGGGGAAATACCGGGGCCGGCGCCGAAGAAGGATCGGGGGGCGGACGCACAAAAAACCGCCGCCCCTGA‑TGTAAAATGAAAAAC‑CCCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGaa                                                                                              <  1:358122/145‑1 (MQ=255)
                                                         aaaaGCCCAATCCTGGAGAAAGACGGCGGACGGGGGCACATAAGGATCGGGCGGCTGACCAAAAAAAAACCGTTGGCCTTG‑‑TTTTAAAAGAAAAAACCCCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGaaa                                                                                             <  1:1031544/150‑1 (MQ=255)
                                                            cgCCCAATCCGGGCGAAAGACCGGGGCCGGGGGCGCACAAGGTTCGTGCGGCTAACCCACAAAAAACCGTGGGCATTGA‑TTTTAAAAGAAAAAC‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAg                                                                                          <  1:337789/150‑1 (MQ=255)
                                                            aGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGaaaa                                                                                            >  1:2427826/1‑149 (MQ=255)
                                                             gCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACggg                                                                                                             >  1:2554988/1‑131 (MQ=255)
                                                               ccGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAgcgc                                                                                       >  2:864931/1‑151 (MQ=255)
                                                               ccGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAgcgc                                                                                       >  2:2272937/1‑151 (MQ=255)
                                                                 gctACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGcgcg                                                                                      >  2:1719136/3‑150 (MQ=255)
                                                                 gATACGGGCAAAAAAACGCGGACGGCGGCGCATCAGGTTCGGGCGGCTGACCCACAAAAAACCGCGGCCCTTGA‑TGTAAAATAAAAAAA‑CCCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGaaaa                                                                                            <  2:127434/144‑1 (MQ=255)
                                                                  aTACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGAc                                                                                    >  2:1181742/1‑151 (MQ=255)
                                                                  aTACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGAc                                                                                    >  1:1258074/1‑151 (MQ=255)
                                                                    aCTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGAc                                                                                    >  2:2494790/1‑149 (MQ=255)
                                                                    aCTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGAc                                                                                    >  2:337299/1‑149 (MQ=255)
                                                                     cctgcacAAGACCGCGGCCGGCGGCCCACCAGGATCGGGCGGCTGCCACAAAACAAACCGCTGCCCTTGA‑TGTTAAAAAAAAAAA‑CCCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGc                                                                                 <  2:1839847/145‑1 (MQ=255)
                                                                       ggcaaaaGAACGCGGACTGCGCCGCAACAGGCTCGTGCGGCTGACCCACAAAAAAACGCGGCCCCTGA‑TTTAAAAAGAAAAAC‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCaa                                                                               <  2:443963/146‑1 (MQ=255)
                                                                       ggCGAAAGAAGGCGGCGTGCGGCGCATCAGGTTGGGGCGGCTGCCGCACAAAAAAACGCGGGGCTTGA‑TTTTAAATGAAAAAC‑CCCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGc                                                                                 <  2:1191692/149‑1 (MQ=255)
                                                                        gcgcAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGc                                                                                 >  1:1839842/1‑148 (MQ=255)
                                                                          gcAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAgc                                                                                         >  1:139639/1‑138 (MQ=255)
                                                                          gcAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATtc                                                                            >  1:1238175/1‑151 (MQ=255)
                                                                            aaGAACGGGGACGGCGGCGAATCAGGTTCGGGCGGCTGACCCAAAAAAAACCGCGGCCCCTAA‑TGTTAAAGAAAAAAC‑CCCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATtcgt                                                                          <  2:611261/151‑1 (MQ=255)
                                                                            aaGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATtcgt                                                                          >  2:508335/1‑151 (MQ=255)
                                                                            aaGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATtcgt                                                                          >  1:82940/1‑151 (MQ=255)
                                                                            aaGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATtcgt                                                                          >  1:611260/1‑151 (MQ=255)
                                                                             aGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATtcgtc                                                                         >  2:74918/1‑151 (MQ=255)
                                                                              gAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATtcgtcg                                                                        >  1:1784115/1‑151 (MQ=255)
                                                                                gcggggACGGCGGCGCATCAGGTGCGGGCGGCTGACAAAAAAAAAACCGCGGCCCTTGA‑TGTAAAATAAAAAAC‑CCCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTcgcg                                                                      <  1:1181741/147‑1 (MQ=255)
                                                                                aCGCGGACGGCGGCGCATGAGGATCGTGCGGCTGACGCAAAAAAAACCGCTGCCCCTGA‑TGTTAAAAAAAAAAC‑CCCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTcgcg                                                                      <  1:864930/151‑1 (MQ=255)
                                                                                 cgcgGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTcgcgc                                                                     >  2:1621648/1‑151 (MQ=255)
                                                                                 cctatcttGCGGCGCATCAGATTGGGGCGGC‑GAACCAAAAAAAAACCGCGGCCCTTATTGTTCAAAAAAAAAC‑CCCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTcgcg                                                                      <  2:686330/143‑1 (MQ=255)
                                                                                  ggggccgggggcgCATAAGGATGGTGCGGCTGACCCAAAAAAAACCGCTGCCCGTTA‑TGTAAAAGAAAAAAC‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTcgcgcg                                                                    <  1:1621643/142‑1 (MQ=255)
                                                                                    ggACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTg                                                                  >  2:2490202/1‑151 (MQ=255)
                                                                                      aCTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGAt                                                                >  2:27928/1‑151 (MQ=255)
                                                                                       cgggggcgAATCAGGATCGGGGGGCTCACGCAAAAAAAACCGCGGCCCTTGT‑TTTAAAATAAAAAAC‑CCCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATg                                                               <  2:2427838/147‑1 (MQ=255)
                                                                                        tGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTcgcg                                                                      >  1:686329/1‑143 (MQ=255)
                                                                                          cggcgCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATggg                                                             >  1:1157953/1‑150 (MQ=255)
                                                                                           ggcgCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGat                                                           >  1:2391078/1‑151 (MQ=255)
                                                                                           ggcaCATAAGGATCGGGGGGCTGACACAAAAAAAAACGCGGGCCTTGA‑TTTAAAATAAAAAAC‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGat                                                           <  2:1157954/151‑1 (MQ=255)
                                                                                             ccatagagggtggTGCGGCTGACCCAAAAAAAACCTCGGCCCCTGT‑TATAAAAAGAAAAAT‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGatta                                                         <  2:1238177/139‑1 (MQ=255)
                                                                                              gcATCTGGATGGGGGGGCTGCCCCAAAAAAAACCGCTGGCCCTAT‑TTTAAAATAAAAAAC‑CCCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGattat                                                        <  2:1993120/151‑1 (MQ=255)
                                                                                              gcATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGat                                                           >  1:526904/1‑148 (MQ=255)
                                                                                              gaaTCAGGATCGTGCGGCGGACGAAAAAAAAACCGCGGCCCTTGA‑T‑TTAAAAGAAAAAACCTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGat                                                           <  2:526904/146‑1 (MQ=255)
                                                                                                aTCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATtg                                                      >  1:1537569/1‑151 (MQ=255)
                                                                                                  aaGGATCGTGCGGCTCACACACAAAAAAACGCTGCCCTTTA‑TGTTAAATGAAAAAC‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATtgct                                                    <  1:625469/150‑1 (MQ=255)
                                                                                                     gATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGc                                                                                 >  1:1065234/1‑119 (MQ=255)
                                                                                                       tCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGa                                                >  1:1781863/1‑150 (MQ=255)
                                                                                                        cGTGCGGCGGACGCACAAAAAAGCGCTGGCTCTGT‑TTTAAAATAAAAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGaaa                                              <  1:1341199/151‑1 (MQ=255)
                                                                                                          ggcgcCTCACCCACAAAAAACCGCGGGCTCTGA‑TTTAAAATGAAAAAC‑CCCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACg                                            <  2:1668951/150‑1 (MQ=255)
                                                                                                             ggCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGc                                          >  1:1701779/1‑150 (MQ=255)
                                                                                                               cTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGc                                          >  1:570783/1‑148 (MQ=255)
                                                                                                               cTGACCCAAAAAAAACCGCTGCCCCTGA‑TTTAAAATGAAAACC‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGc                                          <  2:570784/148‑1 (MQ=255)
                                                                                                                tGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGattat                                                        >  1:1911395/1‑133 (MQ=255)
                                                                                                                 gACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATtgct                                                    >  1:911252/1‑136 (MQ=255)
                                                                                                                 gACGAAAAACAAAACGCGGCCCCTGA‑TTTTAAATAAAAAAC‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGCACCtt                                     <  2:2546379/151‑1 (MQ=255)
                                                                                                                    gCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGCACCTTTAc                                  >  2:2339988/1‑151 (MQ=255)
                                                                                                                            aaaCGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGCACCTTTACCTGAcc                            >  2:1370864/1‑149 (MQ=255)
                                                                                                                              ccgcgGGCCTTGA‑TTTTAAATAAAAAAC‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGCACCTTTACCTGACCAGAt                        <  2:1537574/150‑1 (MQ=255)
                                                                                                                              aggcgggCCATGA‑TGTTAAATAAAAAAC‑CCCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGCACCTTTACCTGACCAGAt                        <  1:366224/146‑1 (MQ=255)
                                                                                                                                 cTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGCACCTTTACCTGACCAGATAAt                     >  2:2364277/1‑151 (MQ=255)
                                                                                                                                  tGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGCACCTTTACCTGACCAGATAAtt                    >  1:2310288/1‑151 (MQ=255)
                                                                                                                                   ggCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGCACCTTTa                                   >  2:312122/1‑135 (MQ=255)
                                                                                                                                    gCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGCACCTTTACCTGACCAGATAATTTc                  >  2:513523/1‑151 (MQ=255)
                                                                                                                                       cTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGCACCTTTACCTGACCAGATAATTTCGcc               >  2:1235502/1‑151 (MQ=255)
                                                                                                                                          a‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGCACCTTTACCTGACCAGATAATTTCGCCCGt            >  1:369097/1‑151 (MQ=255)
                                                                                                                                                   tGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGCACCTTTACCTGACCAGATAATTTCGCCCGTCTGCAcc     >  1:1955269/1‑150 (MQ=255)
                                                                                                                                                    gagaAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGCACCTTTACCTGACCAGATAATTTCGCCCGTCTGCAcc     >  1:1834278/1‑149 (MQ=255)
                                                                                                                                                     agaAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGCACCTTTACCTGACCAGATAATTTCGCCCGTCTGCACCAGc  >  1:631038/1‑151 (MQ=255)
                                                                                                                                                     agaAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGCACCTTTACCTGACCAGATAATTTCGCCCGTCTGCACCAGc  >  2:1853136/1‑151 (MQ=255)
                                                                                                                                                      gaAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGCACCTTTACCTGACCAGATAATTTCGCCCGTCTGCACCAGc  >  1:1185352/1‑150 (MQ=255)
                                                                                                                                                      gaAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGCACCTTTACCTGACCAGATAATTTCGCCCGTCTGCACCAGc  >  1:1143183/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CATCCAGTCACGGAGATATAAATGCGGCTTCGGTGCAATGCATGCGGGTAGTTGTGATAAAGCCCGATACTGGCGCAAGAACGCGGACTGCGGCGCATCAGGATCGTGCGGCTGACGCACAACAAAACGCTGGCCCTGA‑TGTTCAATGAGAAAA‑CTCCCGCCACTAAGCCCACTATTAAACGTCGTGACGGGATGATACTGCGGGAAAAAGCGCGACAGCAATTCGTCGCGCGTGATGGGATTATTGCTGCGAAACGGCACCTTTACCTGACCAGATAATTTCGCCCGTCTGCACCAGC  >  DH1/2716219‑2716517

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