Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,997,2800TG39.7% 100.3 / 11.0 59Y145D (TAT→GAT) xdhCxanthine dehydrogenase iron sulfur‑binding subunit XdhC
Reads supporting (aligned to +/- strand):  ref base T (7/28);  new base G (23/0);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 2.43e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTTGCTTATACCTTGCCGCCTGGGCTGAAGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTATGACGTGAAAGCGCCTCGCGGATCATTTGCTTTTCCATCTCCTCCAGCGCCGTGCCGCCCGC  >  NZ_CP009273/2997035‑2997519
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
gTTGCTTATACCTTGCCGCCTGGGCTGAAGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCa                                                                                                                                                                                                                                             <  1:101196/250‑1 (MQ=255)
  tGCTTATACCTTGCCGCCTGGGCTGAAGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGGAt                                                                                                                                                                                                                                               >  1:4288/1‑246 (MQ=255)
        taacTTGCCGCCTGGGCTGAAGGAAAAGAGATCAGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAA‑ATTTTTCGCTGCGAGGGGGATAAGGTGAATg                                                                                                                                                                                                                                     >  1:401984/4‑250 (MQ=255)
        taCCTTGCCGCCTGGGCTGAAGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGGAAT‑CTTTGTGGCTGCAAGGGGTATCAGATGATTg                                                                                                                                                                                                                                     >  1:401981/1‑250 (MQ=255)
          ccTTGCCGCCTGGGCTGAAGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGGAAT‑CTTTTTCGCTGCACGGGGGATCAGATGATTGt                                                                                                                                                                                                                                    >  1:334806/1‑249 (MQ=255)
          ccTTGCCGCCTGGGCTGAAGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGt                                                                                                                                                                                                                                    <  2:270172/249‑1 (MQ=255)
                      ggCTGAAGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTg                                                                                                                                                                                                                       <  1:166496/250‑1 (MQ=255)
                       gCTGAAGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTgg                                                                                                                                                                                                                      <  2:241831/250‑1 (MQ=255)
                         tGAAGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGGAATCCTTTGT‑GCCGCCCGGGGGATTAGATGGATGGAAATaaaagtttgggt                                                                                                                                                                                                                    >  2:543228/1‑239 (MQ=255)
                           aaGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGAtt                                                                                                                                                                                                                   <  1:253830/249‑1 (MQ=255)
                           aaGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTg                                                                                                                                                                                                                  <  1:148048/250‑1 (MQ=255)
                            aGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGGAAT‑CTTTGTTGCTGCACGGGGGATCAGATGATTGTAAATACCGTTTTGGATTGc                                                                                                                                                                                                                 >  1:122476/1‑250 (MQ=255)
                                   gagaTCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGGATCAGATGATTGTAAATACAGTTCTGGATTGCGAGaaa                                                                                                                                                                                                           >  1:10475/1‑249 (MQ=255)
                                        ccGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTTTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAg                                                                                                                                                                                                     >  1:343195/1‑250 (MQ=255)
                                                tGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCCGCAAGGGGTATCAGATGATTGTTAATACAGTTCTGGAGTGCGAGAAAACGAAGGaaa                                                                                                                                                                                                 >  1:438737/1‑246 (MQ=255)
                                                  gAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGGAAA‑CTTTTTTCGTTCACGGGGGATTAAATTATTGTAAAAACCAGTCTGGGTTTGGGGAAAACAAAAGAAAAAGat                                                                                                                                                                                            >  2:560511/1‑249 (MQ=255)
                                                     ggTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTTTTGCTGCACGGGGGATTCGATGATTGTAAATACAGTTCTTGGTTGCGAGAAAACGAAGGAAAAGGATATcc                                                                                                                                                                                        >  1:142001/1‑250 (MQ=255)
                                                     ggTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATcc                                                                                                                                                                                        <  2:223708/250‑1 (MQ=255)
                                                       tGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGGAAA‑CTTTTTCGCTTCACGGGGGATCAGATGTTTGTTAATACAATTTTGGGTTTCGGGGAAAAGAAGGAAAAGGATAACCgg                                                                                                                                                                                      >  2:309544/1‑250 (MQ=255)
                                                       tGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCgg                                                                                                                                                                                      <  2:13077/250‑1 (MQ=255)
                                                       tGAAGCGAAAGGCGGAAAACTGTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCgg                                                                                                                                                                                      <  2:338516/250‑1 (MQ=255)
                                                               aaGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGGATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAAAGAAGTAAAAGGATATCCGGCCTGAAtt                                                                                                                                                                              >  1:162949/1‑250 (MQ=255)
                                                                  gCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGGAAT‑CTTTGTCGCTGCACGGGGGATCAGATGATTGTAAATAAAGTTCTTGATTGCGAGAAAAAGAAGTAAAAGGATATCCCGCCCGAATTAAg                                                                                                                                                                           >  2:205566/1‑250 (MQ=255)
                                                                        aaCTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGa                                                                                                                                                                     <  2:577120/250‑1 (MQ=255)
                                                                            ttCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTTTTGCTGCACGGGGGATCAGATGATTGTAAATACAGTTCTGGGTTGGGAGAAAAAGAAGGAAAAGGATATCCGGCCTGAAATCCGACCggggtct                                                                                                                                                                 >  2:319493/1‑245 (MQ=255)
                                                                                 aTGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGa                                                                                                                                                             <  2:124636/249‑1 (MQ=255)
                                                                                  tGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAgg                                                                                                                                                           <  2:563297/250‑1 (MQ=255)
                                                                                     tCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGGATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTa                                                                                                                                                        >  1:273100/1‑250 (MQ=255)
                                                                                          agGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTTTCGCTGCACGGGGGATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTTAGGCCGGATTCACTGAGGGTTTgtgt                                                                                                                                                   >  1:138277/1‑250 (MQ=255)
                                                                                                   cGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACt                                                                                                                                          <  1:503349/250‑1 (MQ=255)
                                                                                                      aaTCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGGAAA‑ATTTTTTGCTGCACGGGGGGTCAGATGATTGTAAATAAAAGTCTGGGTTTCGAGAAAAAGAAGAAAAAAGATATTCCGCCTGAATTTAGGCCGGGTTCACTGAGGGTTTGTGttttaaaaatcc                                                                                                                                        >  2:347861/1‑240 (MQ=255)
                                                                                                                      gtgcGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAA‑CTTTTTCGCTGCACGGGGGATCAGATGATTGTAAATAAAGTTCTGGATTGCGAGAAAACGAAGGAAAAGGATATCCGGCCTGAATTTAGGCCGGGTTCACTGAGGGTATGTGTTTTACAAATCATATTTCTTAATCTTTCg                                                                                                                       >  1:554676/1‑250 (MQ=255)
                                                                                                                       tgcGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCAGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGa                                                                                                                      <  1:406082/250‑1 (MQ=255)
                                                                                                                                    cGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCaa                                                                                                         <  1:12522/250‑1 (MQ=255)
                                                                                                                                     gcctggcctgATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTTTCTCTGCACGGGGTATCAGATGGTTGTAAAAACCGATCTGGGTTGCGAGAAAACCAAGTAAAAAGAAATTCCGCCTGAATTTTGGCCCGATTCACTGAGGTTATTTGTTTTACAACTCATATTTTTTTATCTTTCGCTCGAGCGTTGCAAt                                                                                                        >  2:518600/1‑250 (MQ=255)
                                                                                                                                      cctggcctgATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGGATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTTACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATg                                                                                                       >  1:200479/1‑250 (MQ=255)
                                                                                                                                                  aTGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCCCGGGGGATCAGATGATTGTAAATACAGTTTTGGATTGCGAGAAAACGAAGTAAAAGGATATTCGGCCTGAATTCAGGCCGGATTCACTGAGGGTATGTGTTTAACAAATCATATTTCTTAATTTTTCGATAGAGCGTAGCAATGCCGATGCCCAGt                                                                                           >  1:324136/1‑250 (MQ=255)
                                                                                                                                                     gCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTTTTGCTGCACGGGGGATCAGATGATTTTAAATAAAGTTCTGGGTTGGGGGGAAACGAAGGAAAAAGATATTCCGCCCGAATTTAGGCCGGATTCACTGAGGGTATGTGGTTAACAACTCATATTTCTTAATCTTTCGATAGAGCGTAGCAAATGCGATGCCCAAttta                                                                                        >  2:370825/1‑248 (MQ=255)
                                                                                                                                                          cACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGACGCCCAGTTCATCAGc                                                                                   <  1:96202/250‑1 (MQ=255)
                                                                                                                                                            cGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCaa                                                                                 <  2:115099/250‑1 (MQ=255)
                                                                                                                                                                aaTGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTg                                                                             <  1:245323/250‑1 (MQ=255)
                                                                                                                                                                         gAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCt                                                                    <  2:268751/250‑1 (MQ=255)
                                                                                                                                                                           aaCCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGGAAT‑CTTTTTCCCCGCACGGGGGATTAGATGGTTGTTAATACAGTTCTGGGTTGCGGGAAAAAGGAGTAAAAGGATATTCCGCCCGAATTTAGGCCGGGTTCACTGAGGGTTTGTGTTTTAAAACACATATTTTTTTATTTTTCGCTAGAGAGTAGCAATGACGCTGCCCCGGTCATCCACCACTTGtttctttttgt                                                                  >  2:384259/1‑245 (MQ=255)
                                                                                                                                                                                  cgAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTTTTGCCGCACGGGGTATCAGATGGTTGTAAATAAAGTTCTTGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGGTTTGGGTTTAACAACTCCTATTTTTTAATCTTTCGATAGAACGTAGCAATGCCGGTGCCCCGTTTATCAGCAACTTGCTTCTTGCTTTTTTggag                                                           >  2:157729/1‑247 (MQ=255)
                                                                                                                                                                                  cgAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTATGACg                                                           <  1:233763/250‑1 (MQ=255)
                                                                                                                                                                                  cgAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTATGACg                                                           >  1:84545/1‑250 (MQ=255)
                                                                                                                                                                                         ccATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTTTCGCTGCACGGGGGATCAGATGATTGTAAATAAAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTATGACGTGAAAgc                                                    >  1:162895/1‑250 (MQ=255)
                                                                                                                                                                                           aTTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTTTCTCCGCACGGGGTTTGCGATTATTTTAAATACAGGTCTGGGTTTCGAGAAAAAGAAGTAAAAAGATATTCGGCCTTAATTCAGGCCGGATTCACTTAGGTTATGTGTTTTACAAATCATATTTCTTTATCTTGGGGTAGAGCGGAGCAATGCCGCAGCACCGTTCATCAGCAACATGCTTTTTGTGGTTATTTAGTGAAAgcgc                                                  >  2:309325/1‑250 (MQ=255)
                                                                                                                                                                                           aTTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGGATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGGTATGTGTTTAACAAATCATATTTTTTAATCTTTCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTTCTGTTATGACGTGAAAgcgc                                                  >  1:603631/1‑250 (MQ=255)
                                                                                                                                                                                               aCCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTATGACGTGAAAGCGCCTcg                                              <  1:142246/250‑1 (MQ=255)
                                                                                                                                                                                                     aCGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTATGACGTGAAAGCGCCTCGCGGAt                                         <  1:319493/249‑1 (MQ=255)
                                                                                                                                                                                                                   cgGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATGCAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAGCTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTATGACGTGAAAGCGCCTCGCGGATCATTTGCTTTTCCAt                          <  1:197249/250‑1 (MQ=255)
                                                                                                                                                                                                                   cgGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTATGACGTGAAAGCGCCTCGCGGATCATTTGCTTTTCCAt                          <  2:266699/250‑1 (MQ=255)
                                                                                                                                                                                                                     gACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATAGCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTATGACGTGAAAGCGCCTCGCGGATCATTTGCTTTTCCAtct                        <  2:478676/250‑1 (MQ=255)
                                                                                                                                                                                                                        tGGCGGGAAAT‑CTTTGTCGCTGCACGGGGGATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTATGACGTGAAAGCGCCTCGCGGATCATTTGCTTTTCCATctcct                     >  2:384558/1‑250 (MQ=255)
                                                                                                                                                                                                                        tGGCGGGAAAT‑CTTTGTCGCTGCACGGGGCATCCGATGATTGTAACTACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGCAGCAATGCCGATGCCCAGTTCATCAGCAACCTGCTTCTTGCTGTTATGACGTGAAAGCCCCCCGCGGCTCATTTGCTTTTCCCTctcct                     >  2:581962/1‑250 (MQ=255)
                                                                                                                                                                                                                          gCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTATGACGTGAAAGCGCCTCGCGGATCATTTGCTTTTCCATctcct                     <  1:422538/248‑1 (MQ=255)
                                                                                                                                                                                                                              gAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTATGACGTGAAAGCGCCTCGCGGATCATTTGCTTTTCCATCTCCTCCAgcg               <  1:491593/250‑1 (MQ=255)
                                                                                                                                                                                                                                            cTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTATGACGTGAAAGCGCCTCGCGGATCATTTGCTTTTCCATCTCCTCCAGCGCCGTGccgcccgc  <  1:422735/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
GTTGCTTATACCTTGCCGCCTGGGCTGAAGGAAAAGAGATCCGCACGCTGGAAGGTGAAGCGAAAGGCGGAAAACTTTCTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAAT‑CTTTGTCGCTGCACGGGGTATCAGATGATTGTAAATACAGTTCTGGATTGCGAGAAAACGAAGTAAAAGGATATCCGGCCTGAATTCAGGCCGGATTCACTGAGGTTATGTGTTTAACAACTCATATTTCTTAATCTTGCGATAGAGCGTAGCAATGCCGATGCCCAGTTCATCAGCAACTTGCTTCTTGCTGTTATGACGTGAAAGCGCCTCGCGGATCATTTGCTTTTCCATCTCCTCCAGCGCCGTGCCGCCCGC  >  NZ_CP009273/2997035‑2997519

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