Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,941,2100CT36.2% 111.6 / 30.2 69E40K (GAA→AAA) hdfRHTH‑type transcriptional regulator HdfR
Reads supporting (aligned to +/- strand):  ref base C (12/32);  new base T (25/0);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 5.41e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCGGA  >  NZ_CP009273/3940967‑3941453
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
aGCTTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGtt                                                                                                                                                                                                                                                    <  2:126352/249‑1 (MQ=255)
   cTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCt                                                                                                                                                                                                                                                <  1:79474/250‑1 (MQ=255)
      tATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGAt                                                                                                                                                                                                                                             <  2:234667/250‑1 (MQ=255)
         gggCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGG‑G‑GAAAGGGTTACCACCAG‑ATGATTTTTTAGtttttttgtttt                                                                                                                                                                                                                                          >  2:104169/1‑240 (MQ=255)
           gCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑‑AAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGaaa                                                                                                                                                                                                                                       <  1:164635/250‑1 (MQ=255)
              tCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGG‑G‑AAAAGGGTATCACCCAG‑TTGAGTTTTTCGTTGTCTGATTTGAAAg                                                                                                                                                                                                                                      >  2:447021/1‑249 (MQ=255)
                                    ccAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑GAAAAGTTCACACCCCC‑GTGTTTTTTCACTTTTTTTATTTGAAAACTCCCCTCTGACTGGGTCAGa                                                                                                                                                                                                                 >  1:281845/1‑248 (MQ=255)
                                      aGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATaga                                                                                                                                                                                                             <  1:644760/250‑1 (MQ=255)
                                         cACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCg                                                                                                                                                                                                          <  1:292437/250‑1 (MQ=255)
                                                tAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGG‑GAAAAAGGTTCA‑ACCCCG‑CTGTATTTTCTCTTGTTTGATTCGAAAACACACGGTTGCATGGGTGTGCTAGAGCGATTCTG‑c                                                                                                                                                                                                  >  2:426364/1‑250 (MQ=255)
                                                    cATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCt                                                                                                                                                                                              <  2:422353/250‑1 (MQ=255)
                                                        cATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGG‑G‑AAAAAGTTTACACCCAG‑GTGATTTTTTCGTTTTCTGATTCGAAAGATCACCGCTGGCTGGGTTAGGAAGAGCGAGTCAG‑CCGCTCGAc                                                                                                                                                                                          >  1:32128/1‑250 (MQ=255)
                                                               acaACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGG‑G‑GAAAAGTTTACACCCCGATTGTTTTTCC‑GTTTTCTGTTTTTGAAACTCACTGCTGGCTGGGGCAGATAGAGCGAGTCAG‑ccccccggcccaa                                                                                                                                                                                      >  2:100865/1‑236 (MQ=255)
                                                                    gAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGG‑G‑AAAAAGGTCACACCCAG‑TTGATTTTTCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGGCAGATAGAGGGATTCAG‑CCGCTCGACCAAAAAGAc‑gc‑g                                                                                                                                                                            >  1:362354/1‑250 (MQ=255)
                                                                         gcTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTGG‑AAAAGGTT‑ACACCCAC‑TTGATTTTTCAGTTGTCTGTTTCGAAAGCTCACCGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCCCGCCCAAAGAGAA‑GC‑GTTCgg                                                                                                                                                                       >  1:735310/1‑250 (MQ=255)
                                                                         gcTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGG‑G‑AAAAGGGTTACACCCCG‑GTGATTTTTCAGTTGTTTTATTCGAAAACTCACTGCCGACTTGGTCAGATAGAGGGATTCCA‑CCGCCCGACCAAAAGGAC‑GG‑GGTCgg                                                                                                                                                                       >  1:21503/1‑250 (MQ=255)
                                                                              cACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGTGG‑AAAAGGTTC‑CACCCAC‑TTGGTTTTTCTGTTTTCTTTTTTGGAAAGTCACTGCCGTCTGGGGCAGATAAAGAGATTCAG‑GCGCCCCGCCAAAAAGAC‑GC‑GATTGGCGAAc                                                                                                                                                                  >  2:441233/1‑250 (MQ=255)
                                                                               aCCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACt                                                                                                                                                                 <  2:295147/250‑1 (MQ=255)
                                                                                      gAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAgg                                                                                                                                                           <  2:225408/249‑1 (MQ=255)
                                                                                      gAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAgg                                                                                                                                                           <  1:247307/249‑1 (MQ=255)
                                                                                        aaaCTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGaaa                                                                                                                                                        <  1:249505/250‑1 (MQ=255)
                                                                                                 aTGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑GAAAGGGTCACACCCCG‑TTTATTTTTCCGTTGTTTGATTTGAAAACTCACTGCTGACTGGGGCAGATAGAGCGATTTAG‑CCGCTCGACCAAAGTGAC‑GC‑GGTCGGCTAACTTCCAGGGAAGtttttta                                                                                                                                                >  1:293396/1‑247 (MQ=255)
                                                                                                 aTGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGG‑G‑AAAAAGGTTACACCCAG‑TTGATTTTTCAGTTGTCTGGTTCGAAAACTCACTGCTGGCTGGGGCAGATAGAGCGATTTAG‑CCGCTCGGCCCAAAGGAA‑GC‑GGTCGGCTAACTTCCAGGGAAGgttttta                                                                                                                                                >  1:725924/1‑247 (MQ=255)
                                                                                                                cgcCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCcac                                                                                                                                 <  1:279516/249‑1 (MQ=255)
                                                                                                                 gcCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAAGGTCACACCCCG‑TTGATTTTTCCGTTGTCTGATTCGAAAACTCACTGCTGACTGGGGCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGAGAC‑GC‑GTTCGGCTAAATTCCAGGAAAGTTTTTAACAATTCCGTATCCACAg                                                                                                                               >  1:24292/1‑250 (MQ=255)
                                                                                                                    aCCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAAGGTCACACCCCG‑TTGGTTTTTCCGTTTTCTGGTTTGAAAACTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCCCGCCCCAAGGGAC‑GC‑GTTCGGGTCACTTCCAGGAAAATTTTTTACAATTTCGTATCCACAGAGc                                                                                                                            >  1:510197/1‑250 (MQ=255)
                                                                                                                     cctcctTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGct                                                                                                                           <  1:725258/250‑1 (MQ=255)
                                                                                                                     cctcctTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGct                                                                                                                           <  1:725220/250‑1 (MQ=255)
                                                                                                                         ctTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑GAAAAGGTCACACCCCG‑TTGTTTTTTCCGTTTTCTGATTTGAAAGCTCACTGCCGACTGGGGCAGATAGAGCGGTTTAG‑CCGCTCGCCCAAACAGACACC‑GTTTGCTCAATTCCCAGGACAG‑TTTTAAAAATTTCGTATCTCCAGTGCTcc‑cc                                                                                                                      >  2:544974/1‑250 (MQ=255)
                                                                                                                             cGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTa                                                                                                                                                                    <  1:216032/201‑1 (MQ=255)
                                                                                                                             cGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCaa                                                                                                                  <  1:505369/250‑1 (MQ=255)
                                                                                                                             cGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAAGGTTACACCCAG‑TTGATTTTTCTGTTGTCTTATTTTGAAACTCACCGCTGACTGGGGCAGATAGAGCGATTCAG‑CCGCCCGACCCAAAAGAC‑GC‑GTTCGGCGTACTTCTACGAAAGTTTTTTACAATTCCCTATCCCCAGTGCTCT‑CCGcc                                                                                                                   >  2:21401/1‑248 (MQ=255)
                                                                                                                              gggCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCaaa                                                                                                                 <  2:663603/250‑1 (MQ=255)
                                                                                                                                gcggcCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTTACACCCAG‑TTGATTTTTCAGTTGTCTGATTCGAAAACTCACTGCTGACTGGGGCAGATAGAGCGATTTAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTTCGTATCCACAGTGCTCT‑CCGCAAAat                                                                                                               >  1:321019/1‑250 (MQ=255)
                                                                                                                                  ggcCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑GAAAAGGTCACACCCCG‑TTGATTTTTCCGGTGTCTGATTCGAAAACTCACTTCTTACTGGGTCAGATAGAGCGATTTAG‑CCGCTCGACCCAAATGAC‑GC‑GGTTGGGTAACTTCCAGGGAAGTTTTTTAAAAATTCGTATCCACCGTGCTCT‑CCGCAAAattt                                                                                                             >  1:153145/1‑249 (MQ=255)
                                                                                                                                   gcCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGG‑G‑AAAAGGGTCACACCCAG‑GTGATTTTTCAGTTGTCTGATTTGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGAGAC‑GC‑GTTCGGCGAACTTCCCGGAAAGTTTTTTACAATTCCGTATCCACAGTGCTCT‑CCGCAAAatttt                                                                                                            >  1:389352/1‑248 (MQ=255)
                                                                                                                                        ccAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGG‑G‑AAAAAGGTCACACCCAG‑ATGATTTTTCAGTTGTCTGATTCGAAAAATCACTGCTGACTGGGGCAGATAGAGCGATTCAG‑CCGCCCGACCCAAAAGAC‑GC‑GGTCGGCTAACTTCCAGGGAAGTTTTTTACAATTACGTATCCACAGTGCTCT‑CCGCAAAATTATTTGtt                                                                                                       >  1:595005/1‑249 (MQ=255)
                                                                                                                                           ggTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGtt                                                                                                    <  1:375627/250‑1 (MQ=255)
                                                                                                                                            gTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGG‑G‑AAAAAGATTACACCCCGTTTGTTTTTCCATTTTTTTTGTT‑GCAAAAGCACTGCTGGCTGGGGGAGATAAAGGGATTCTGACCCCCCGCCCAAA‑TCAA‑GA‑GGTTGGCTCACTGTACGGGAAGGATTTAACAATTCCCGATCCCCAATGCTCT‑CCGCCAAATAATTTTTCGttt                                                                                                   >  2:517233/1‑249 (MQ=255)
                                                                                                                                                 cATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAAGGTCACACCCAG‑TTGGTTTTTCAGTTGTTTGGTTTGAAAGCTCACTGCTGACTGGGGCAGATAGAGCGATTTCG‑CCGCCCGCCCAAAAAGAC‑GCGGTTGGCCAAATTTCCGGGAAGGTTTTTAA‑AATTCCGTATCCACAGTGCTCT‑CCCCAAAATTATTTTTCGTTTTGGtt                                                                                              >  1:323926/1‑250 (MQ=255)
                                                                                                                                                     aGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGG‑G‑GAAAGGGTTAAACCCAG‑TTGATTTTTCCTTTTTTTGTTTTGAAAGCACACTGCTGACTGGGTTAGATAAAGGGATTTAG‑CCGCCCGCCCAAAAGGAC‑GC‑GTTCGGGGAACTTCCAGGAAAATTTTTTAAAATTCCGTTTCCACAGGGTTCT‑CCGCAAAATTTTTTTTCtttttttttttaa                                                                                          >  2:192232/1‑240 (MQ=255)
                                                                                                                                                      gCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGG‑G‑AAAAAGGTCACACCCCG‑TTTGTTTTTCCGTTTTTTGGTTTGAAAACTCACTGCTGGCTGGGGCAGGTAGAGCGGTTCAA‑CCGCCCGACCCAACAGAA‑GC‑CGTCGGGTAACCTCCAGGAAAGATTTTTACAATTTCGTATTCCCAATGCTCC‑CCCCCAAAATATTTTTTGTTATTGTTTTAAt                                                                                         >  2:211318/1‑250 (MQ=255)
                                                                                                                                                         tttCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAAGTTCACCCCCAG‑TTGATTTTTCCGTTTTCTGATTCGGAAACTCACTGCTGACTGGGTCCGAAAGAGCGATTCAG‑ACGCTCGACCCAAATGAC‑GC‑GGTCGGCGAAACTCCAGGGAAGATTTTTACAATTCCGTATCCACAGTGCTCT‑CCCCCCAAATTTTTTTTGTTTTGATTTTACtgtt                                                                                      >  2:749793/1‑250 (MQ=255)
                                                                                                                                                           tCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAAtgttt                                                                                     <  1:685891/249‑1 (MQ=255)
                                                                                                                                                           tCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGG‑G‑AAAAGGGTTAAACCCAG‑TTGATTTTTCAGTTGTCTGATTCGAAAACTCACTGCTGACTGGGGCAGATAGAGCGATTCAG‑CCGCCCGACCCAAGGGAT‑GC‑GTTCGGCTAACTTTCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCCAAAATATTTTTCGGTTTGATTTTAAtttttt                                                                                    >  1:144433/1‑250 (MQ=255)
                                                                                                                                                                   aGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGTAAAGTTTTTAACAATTCCGTATCCTCAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTac                                                                             <  2:605603/249‑1 (MQ=255)
                                                                                                                                                                       agtagtTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACActct                                                                        <  1:661867/250‑1 (MQ=255)
                                                                                                                                                                       agtagtTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACActct                                                                        <  2:422243/250‑1 (MQ=255)
                                                                                                                                                                             ttttCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGG‑G‑AAAAGGATCTCCCCCAC‑CTTATTTTCTCTTTTTTTGTTTTGAAAGCTCCCTTCTGGCTGGGGCCGATTTAGGGGTTTAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTTGGCCACCTTCCACGGAAGTTTTTTACACTTCCGTTTTCACACCGCTCTCCCGCCACATT‑TTTTTTGTTATTTTTTTTATGTTTTGTTTTTCAATCTTTctct                                                                  >  2:117553/1‑247 (MQ=255)
                                                                                                                                                                                    cGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACt                                                           <  2:168641/250‑1 (MQ=255)
                                                                                                                                                                                    cGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACt                                                           <  2:167192/250‑1 (MQ=255)
                                                                                                                                                                                    cGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGG‑G‑AAAAAGGTCACACCCCC‑GTGATTTTTCAGTTTTCTGATTTGAAAACTCACCGCTGACTGGGGCAGATAGAGAGATTCCG‑CCCCCCGAACAAAAAGAC‑GC‑GGTCGGCTAACTTCCCGGAAAGTGTTTTACAAATCCCTATCCACACAGCTCT‑CCGCAAAATTATTTTTCGTTATTATTTATATGTTTTTTTTTATCTTCTGTCAAGCGttaac                                                           >  2:375978/1‑246 (MQ=255)
                                                                                                                                                                                        aGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATa                                                       <  1:654752/250‑1 (MQ=255)
                                                                                                                                                                                           ggTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGGTTACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAACTCACTGCTGACTGGGTCAGATAGAGCGATTTAG‑CCGCCCGACCAAAAGGAC‑GC‑GTTCGGCTCACTTCCAGGAAAGTTTTTTACAAATTCGTATCCACAGAGCTCT‑CCGCCAAATTATTTTTTGTTTTGGTTTTAATGTTTTGTTTTTCACTCTGTCAAGCGTAAATAATACTc                                                    >  1:292258/1‑250 (MQ=255)
                                                                                                                                                                                            ggtAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTcc                                                   <  2:735310/248‑1 (MQ=255)
                                                                                                                                                                                            gTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTcc                                                   <  1:415913/250‑1 (MQ=255)
                                                                                                                                                                                                 aCGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGcc                                              <  1:100865/250‑1 (MQ=255)
                                                                                                                                                                                                    gATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGTCTGGGGCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGAGAC‑GC‑GTTCGGCTAACTTTCAGGAAAGATTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCCAAATTATTTGTCGTTATGTTTTTAATGTTTTGTTTTTCACTCTGTCTAGCGTAACTAATACTTCGCGCCAt                                            >  1:65054/1‑249 (MQ=255)
                                                                                                                                                                                                     atatTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATATTTTGTTTTACACTCTGTCAAGCTTAACTAATACTCCGCGCCATaa                                          <  2:673233/250‑1 (MQ=255)
                                                                                                                                                                                                      tatTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAAc                                         <  2:398866/250‑1 (MQ=255)
                                                                                                                                                                                                                 tgGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTc                               <  2:692242/249‑1 (MQ=255)
                                                                                                                                                                                                                  ggCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCa                              >  1:6117/1‑249 (MQ=255)
                                                                                                                                                                                                                          g‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAAt                        >  1:555214/1‑248 (MQ=255)
                                                                                                                                                                                                                                   tCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTg              >  1:655074/1‑250 (MQ=255)
                                                                                                                                                                                                                                    cacaCCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTg              >  2:621886/1‑249 (MQ=255)
                                                                                                                                                                                                                                      cacCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAg           <  1:114874/250‑1 (MQ=255)
                                                                                                                                                                                                                                       acCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAg           >  1:141254/1‑249 (MQ=255)
                                                                                                                                                                                                                                              ttGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCg    >  1:502566/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  ttttCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCGGa  >  2:472711/1‑248 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
AGCCTGTATGGGCATCCTGATTTTGATACAAGCGTCCCAGCCACTGATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGT‑G‑AAAAGGTTCACACCCAG‑TTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAG‑CCGCTCGACCAAAGTGAC‑GC‑GTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCT‑CCGCAAAATTATTTGTCGTTATGATTTAAATGTTTTGTTTTACACTCTGTCAAGCGTAACTAATACTCCGCGCCATAACTAGCTCGGTCAAAGAATTAGGAGCGTGCAGGATGGCGGA  >  NZ_CP009273/3940967‑3941453

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