Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,980,4760AC32.4% 148.4 / 13.0 71Y24S (TAC→TCC) ftnBnon‑heme ferritin‑like protein
Reads supporting (aligned to +/- strand):  ref base A (34/14);  new base C (0/23);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 2.41e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTGATGATTTAAGCATTTTCTTATACCCGTTCAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATG  >  NZ_CP009273/1980228‑1980724
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                         
ttGATGATTTAAGCATTTTCTTATACCCGTTCAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTAc                                                                                                                                                                                                                                                          >  1:763507/1‑250 (MQ=255)
    tgatTTAAGCATTTTCTTATACCCGTTCAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTc                                                                                                                                                                                                                                                     <  1:356591/250‑1 (MQ=255)
       tttAAGCATTTTCTTATACCCGTTCAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGACTGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCAcc                                                                                                                                                                                                                                                  >  2:49007/1‑250 (MQ=255)
          aaGCATTTTCTTATACCCGTTCAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTg                                                                                                                                                                                                                                                <  2:224304/249‑1 (MQ=255)
            gCATTTTCTTATACCCGTTCAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAg                                                                                                                                                                                                                                              >  2:25837/1‑249 (MQ=255)
             cATTTTCTTATACCCGTTCAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGt                                                                                                                                                                                                                                             >  2:367973/1‑249 (MQ=255)
              aTTTTCTTATACCCGTTCAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTaa                                                                                                                                                                                                                                           >  2:476013/1‑250 (MQ=255)
                   cTTATACCCGTTCAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGgt                                                                                                                                                                                                                                      >  2:570379/1‑250 (MQ=255)
                       taCCCGTTCAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGtt                                                                                                                                                                                                                                   <  1:61066/249‑1 (MQ=255)
                       taCCCGTTCAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGtt                                                                                                                                                                                                                                   <  2:434072/249‑1 (MQ=255)
                             ttCAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGaa                                                                                                                                                                                                                              >  1:261441/1‑248 (MQ=255)
                               cAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACag                                                                                                                                                                                                                           <  2:388483/249‑1 (MQ=255)
                                aGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACagag                                                                                                                                                                                                                         <  2:374229/250‑1 (MQ=255)
                                 gACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGt                                                                                                                                                                                                                        <  1:339578/250‑1 (MQ=255)
                                     ttATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTg                                                                                                                                                                                                                     >  2:474664/1‑249 (MQ=255)
                                      tATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGaa                                                                                                                                                                                                                   >  2:576867/1‑250 (MQ=255)
                                        ctCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCCGCAGAAAGGAATCATTCTGCAAACAGTAATTCTGGTGATTTGATTTATCTTGCACCTCTCCCCTTCTGGATATAAGGATATTAGGCACGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTATGAACAGAGTCTGAACg                                                                                                                                                                                                                 >  2:711941/2‑250 (MQ=255)
                                                 cAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCa                                                                                                                                                                                                        >  2:549348/1‑250 (MQ=255)
                                                      caccgtaaaaaTTGACTCCACGATCACTTTGCGGACCGGAAAAAAGGAATTATTCTGAAAAAAGTATTTATGGTTTTTTGATTTATCTTCCCCCTCCCCCTTTCGGGATATAAGGATATTAGGTAGGCCACCCGCTGAAATGCTTCTAAACCTCACCTCTCAAATGACCCGCGATTTTTACGCATCCAATCTCTCCCTTCCCCTGAGTAACTGGTGTTCTGAACAGAGTCGGACCGGCACCGCCACTTTc                                                                                                                                                                                                   <  2:221679/241‑1 (MQ=255)
                                                         cGTCAAAATTGACTCCACGATCCCATTTCGGACCGGCAAAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTCCCCCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCACCCGCTGAAATGCTTCTCAAACTCAACTCTCAAATGACCCGCGATTTTTACGCATCCAATCTCTCCCTTCCCCTGAGTAACTGTTTTTCTAAAAAGAGTCTGAACGGCACCGCCACTTTCCtt                                                                                                                                                                                                <  2:216093/250‑1 (MQ=255)
                                                              gAATGGACTCCACGATACCATTCCGGACGGGAAAAAAGGATTTATTCTGCAAACATTATTTAGGGTTTTTTGATTTATCTTGCACCTCCCCACTTCGGGATATAAGGATATTAGGTATGGCACCCGCTGGAATGCTTCTCAAACTCACCTCTAAAAGGACCCGCGATTTTTACGCATCCATTCTCTCCCTCCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCgcgc                                                                                                                                                                                           <  2:621258/250‑1 (MQ=255)
                                                               aaTTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGcc                                                                                                                                                                                          >  1:352852/1‑250 (MQ=255)
                                                                    aCTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTATTTATGGTGTTTTGATTTATCTTGCCCCTCTCCCCTTCGGGATAAAAGGATATTAGGTATGGCACCCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGACCCGCGAGTTTTACGCATCCAATCTCTCCCTTCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTGACCGGCACCGCCACTTTCCTTCGCGCCCAgg                                                                                                                                                                                      <  1:4995/249‑1 (MQ=255)
                                                                        cACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGca                                                                                                                                                                                    >  1:460889/1‑247 (MQ=255)
                                                                        cACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTAGGGTGTTTTGTTTTATCTTGCACCTCTCCTCTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGAAATGCTTCTCAAACACACCTCTCAAATGAACCGCGATTTTTACGCACCCAATCTCTCCCCTCACCTGAGTACCTGTTGTTCTGAACAGATTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGcac                                                                                                                                                                                   <  1:502647/248‑1 (MQ=255)
                                                                             tCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTaa                                                                                                                                                                            >  1:65703/1‑250 (MQ=255)
                                                                               acaTTTCGGACCGGAAGAAAGGAATTTTTCTGCAAACAGTAATTTTGGTTTTTTGATTTATCTTGCCCCTCTCCACTTCTGAATATAAGAATATTAGGTATGGCAACCGCTGGAATGCTTCTCAACCTCAACTCTCAAATGAACCGCGAGTTTAACGCATCCAATCTCTCCCTTCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTGACCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAAtg                                                                                                                                                                          <  1:97214/250‑1 (MQ=255)
                                                                                 aTTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTTTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTAAAATGACCCGCGATTTTTACGCATCCAATCTCTCCCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAAtgtg                                                                                                                                                                        <  1:651444/250‑1 (MQ=255)
                                                                                        aCCGGTAAAAAGGAATTATTCTGCAAACAGTAATTATGTTGTTTTGATTTTCTTTCCCCCCCTCCACTTTTGGATAAAGGGATATTAGGTAGGCCACCCGCGGTAATGCTTCTCAAACTCACCTCCAAAATGACCCGCGAGTTTTACGCATCCAATCTCTCCCTTCCCCTGAGTACCTGGTGTTCTAACCAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCaaa                                                                                                                                                                 <  2:485076/250‑1 (MQ=255)
                                                                                              aGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATgc                                                                                                                                                           >  1:347635/1‑250 (MQ=255)
                                                                                                   ggAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATg                                                                                                                                                      >  1:830548/1‑250 (MQ=255)
                                                                                                    gAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGt                                                                                                                                                     >  2:834999/1‑250 (MQ=255)
                                                                                                      tttattCTGCAAAAAGTAATTATGTTGTTTTGATTTATCTTGCCCCTCCCCTTTTCTGAAAATAAGGATATTGGGTAGGGCAACCGCTGAAATGCTTCTCAAACCCAACCCAAAAATGACCCGCGAGTTTTACCCATCCATTCTCTACCTCCCCCTGAGAAACTGGTTTTCTGAACAGAGTCTGACCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGttt                                                                                                                                                   <  2:713316/249‑1 (MQ=255)
                                                                                                           tCTGCAAACAGTAATTATGGTGTTTTGATTTTTCTTGCCCCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCACCCGCTGGAATGCTTCTCAAACCCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTCCCTTCCCCTGAGTAACTGGTTTTCTGAACAGAGTCGGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACt                                                                                                                                               <  1:793998/249‑1 (MQ=255)
                                                                                                            cTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACtt                                                                                                                                              >  1:369019/1‑249 (MQ=255)
                                                                                                                aaaCAATAATTATGGTTTTTTAAATAATAGTGACCCTCTCGAGGTATAAAAAGAAGGTAATAGGGATGGCGAACCGGTGGAATGCTTCTAAAACTCAACTCTAAAGTGAACCGAGGGTTTTACACATCTTATCTCTCCCTTCCCCGGAGAAACGGTTGTTCGGAACAGGGTCTGAACGGCGCCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAAGGTGACCCAAATGATGCGCATGTTTAACTTTATGa                                                                                                                                         <  2:327768/250‑1 (MQ=255)
                                                                                                                   cAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGaa                                                                                                                                        >  2:268197/1‑248 (MQ=255)
                                                                                                                    tgTAATTATGGTGTTTTGATTTTTCTGACACCTCTCCACTACTGAATATAAGGATATTAGTTATGCCAACGCCGGAAATGTTCCTCAAACCAACCTCAAAAAGGACCCGCAATTTTACCCAAACCATTCCCTCCCTTCCCCTGGGTAACTGGTGTTCGAAACAGATTCTGAAGGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAagag                                                                                                                                     <  2:743031/249‑1 (MQ=255)
                                                                                                                    aGTAATTATGTTTTTTTAATTTATCTTCCCCCTCTCCTCTTCTGGATATAAGGATATTAGGTATGGCACCCGCTGGAATGCTTCTCAAACTCAACTCTCAAATAAACCGCGAGTTTTACGCATCCAATCTCTCCCTTCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAagag                                                                                                                                     <  1:387192/250‑1 (MQ=255)
                                                                                                                    aGTAATTATGGTGTTTTGATTTATCCTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAagag                                                                                                                                     >  2:406785/1‑250 (MQ=255)
                                                                                                                     gTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGAAAAGAAGGATATTAGGAAGGGCAACCGCTGGAATGCTTCCAAAACTAAACTCTAAAATGAACCGAGAGTTTTCCGCAGCCAATCTCTACCTTCACCGGAGTAACTGGTGTTCGGAAAAGAGTCTGAACGTCTCCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAgt                                                                                                                                    <  1:410273/250‑1 (MQ=255)
                                                                                                                        ttttgggTTTTTTGATTTATCTTCCCCCTTCCCACTTTTGTATAAAAGTATATTAGGTAGGGCACCCGCTGAAATGCTTCTCAAACTCAACTCTCAAATGACCCGCGAGTTTTACGCATCCAATCTCTCCCTCCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAgtgt                                                                                                                                  <  2:279082/244‑1 (MQ=255)
                                                                                                                             ggTGTTTTTATTTATCTTCCACCTCTCCACTTCGTGATATAAGGATATTAGGTATGGCCACCGCTGGAATGCTTCTCAACCTCACCTCTCAAATGACCCGCGAGTTTTACGCATCCAATCTCTCCCTTCCCCTGAGTAACTGGTGTTCTAAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGgcg                                                                                                                             <  2:412046/249‑1 (MQ=255)
                                                                                                                                gtTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTAc                                                                                                                         >  1:36639/1‑250 (MQ=255)
                                                                                                                                    tGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTAcccc                                                                                                                      >  1:114090/1‑249 (MQ=255)
                                                                                                                                       tttATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATc                                                                                                                  >  1:452084/1‑250 (MQ=255)
                                                                                                                                       tttATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATc                                                                                                                  >  2:42530/1‑250 (MQ=255)
                                                                                                                                            cTTGCCCCCCTCCCCTTCTGAATAAAAGGATATTAGGTATGGCACCCGCTGGAATGCTTCTCAAACTCACCCCCCAAATGACCCGCGATTTTTACGCATCCAATCTCTCCCTCCCCCTGAGTAACTGGTGTTCTAAACAGATTCTAAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTa                                                                                                              <  1:339618/249‑1 (MQ=255)
                                                                                                                                               gCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCACCTCTCAAAGGAACCGCGAGTTTTACGCATCCAATCTCTCCCTTCACCTGATTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAg                                                                                                           <  1:698996/249‑1 (MQ=255)
                                                                                                                                                ccccTCTCCCCTTCTGGATATAAGGATTTTAGGTATGGCACCCGTGGAAATGTTTCCAAAACTCACCTCTAAAATGAACCGGGATTTTTCCGCATCCAATCTCCCCCTCCCCCGGATTACCGGGTGTCCTAAACAGAGCCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGc                                                                                                          <  2:221538/247‑1 (MQ=255)
                                                                                                                                                cACCTCTCCTCTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGACCCGCGAGTTTTACGCATCCAATCTCTCCCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGc                                                                                                          <  1:337748/249‑1 (MQ=255)
                                                                                                                                                  ccTCTCCCCTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAAGGAACCGCGATTTTTACGCATCCAATCTCTCCCTCCACCGGAGTAACTGGTGTTCTGAACAGAGTCGGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCAt                                                                                                       <  1:389013/250‑1 (MQ=255)
                                                                                                                                                         aCTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGtt                                                                                                >  2:281442/1‑250 (MQ=255)
                                                                                                                                                           gtcgggATATAAGGATTGAAGATAGGACACCCGGTGGACGGCTTCTCAACCGCACCTCACAAATGACCCGCAATTTATACCCACCAAATCTCCCCCTTCCCCAGAGAACCTGGTGTTCTTAACAGGTTCTGACCGGCTCCGCCACTTTCCGTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCAGCGTTAAAGCCATTGATGtt                                                                                                <  2:743798/244‑1 (MQ=255)
                                                                                                                                                                     aaGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGaaaaa                                                                                    >  1:405159/1‑250 (MQ=255)
                                                                                                                                                                     aaGAATATTAGTAAGGGCACCCGCTGAAACGCTTCCAAACCCAACCCCCAAAATAACCCGCGATTTTTACGCCTCCAACCCCCACCTCCCCCTGAGTACCGGGTGTTCGGAACAGAGTCTGACCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGaaaaa                                                                                    <  2:69482/250‑1 (MQ=255)
                                                                                                                                                                       gaatatAAGTTATGCCACCCGCTGAAATGCTTCCCAAACCCAACCCTAAAAGGAACCGCGATTTTTACCCATCCATTCTCCCCCTTCACCGGAGTACCTGGTGTTCTGAACAGAGTCTGAACGCCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACt                                                                                  <  2:487051/248‑1 (MQ=255)
                                                                                                                                                                                           cGCTGAAATGCTTCCAAAACTCACCCCTAAAATGACCCGCGATTTTTCCGCACCCAATCTCCCCCTCCCCCGGATTAACGGGTGTCCGAAACAGAGTCTAACCGCCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTg                                                               <  2:607633/250‑1 (MQ=255)
                                                                                                                                                                                                 aaTGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTACAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCa                                                          >  2:507878/1‑249 (MQ=255)
                                                                                                                                                                                                            aactcacctcGAAAATAAACCGCGAGTTTTACGCATCCATTCTCTCCCCCCCCCTGATTAATTGTGGTCCTGAACAGAGTCTGACCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGga                                              <  2:155397/250‑1 (MQ=255)
                                                                                                                                                                                                               tcaactcTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCGCAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGgaag                                            >  1:382207/1‑249 (MQ=255)
                                                                                                                                                                                                                caactcTCAAATGAACCGCGATTTTTACGCATCCAATCTCTACCTTCCCCTGAGTAACTGGTGTTCTAAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGgaagaa                                          <  1:720515/250‑1 (MQ=255)
                                                                                                                                                                                                                     ctcAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAAGATGGGAGAATaacag                                    >  1:516688/1‑246 (MQ=255)
                                                                                                                                                                                                                     cccaaaaTGAACCGCGGTTTTTACCCCTCAAACTCCAACTTCCCCCTGAGTAACTGTTTTTCTGAAAAGAGTCTGACCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACg                                      <  2:791537/246‑1 (MQ=255)
                                                                                                                                                                                                                         aaaTGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGc                                   >  1:512337/1‑249 (MQ=255)
                                                                                                                                                                                                                         aaaTGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGc                                   >  1:672964/1‑249 (MQ=255)
                                                                                                                                                                                                                           ataaaCCGCGAGTTTTACGCACCCAATCTCTACCTTCCCATGGTTAACTGGTGTTCTGAAAAGAGTCTGACCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAAc                                <  1:691902/247‑1 (MQ=255)
                                                                                                                                                                                                                                       ttttACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGt                    >  1:87481/1‑250 (MQ=255)
                                                                                                                                                                                                                                        tttACGCATCCAATCTCTACCTTCACCTGAGTATCTGGTGTTCTGAATAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGACATTGATGTTCCAGGTGAAGAACTGCACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACgat                   >  2:541023/1‑248 (MQ=255)
                                                                                                                                                                                                                                         ttACGCATCCAATCTCTACCTTCACCGGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGtt                   <  1:549348/249‑1 (MQ=255)
                                                                                                                                                                                                                                             gCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAACAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGTca               >  2:752782/1‑249 (MQ=255)
                                                                                                                                                                                                                                                         ccaCCTTCCCCGGAGTACCTGGTGTCCTGAACAGAGCCTGAACGGCCCCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGAt   <  2:104321/248‑1 (MQ=255)
                                                                                                                                                                                                                                                           aCCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATg  >  1:515954/1‑249 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                         
TTGATGATTTAAGCATTTTCTTATACCCGTTCAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATG  >  NZ_CP009273/1980228‑1980724

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