Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,221,7540TC40.9% 64.9 / 15.6 67pseudogene (484/1436 nt)mdtQmultidrug resistance outer membrane protein MdtQ
Reads supporting (aligned to +/- strand):  ref base T (21/18);  new base C (13/14);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGC  >  NZ_CP009273/2221508‑2221989
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  
cGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCg                                                                                                                                                                                                                                                 <  2:537128/250‑1 (MQ=255)
   gTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGccc                                                                                                                                                                                                                                              <  1:169844/250‑1 (MQ=255)
    ttttGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCGCCCCCCGGGCTTTTACCCTACCCCGGCGGGCAGTAACCTCCGcccc                                                                                                                                                                                                                                             >  1:345434/1‑249 (MQ=255)
     tttGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCGTGCTTTTACCGTACCCAGGCGGGCAGTAACCTCCGCCCgc                                                                                                                                                                                                                                            >  1:422778/1‑249 (MQ=255)
      ttGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTCCGCCGCCCGTGCTTTAAACCTAACCAAGCGGGGAGTAACCcccgcccggt                                                                                                                                                                                                                                           >  1:846654/1‑248 (MQ=255)
        gCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGTGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCGCCGCCCGAtt                                                                                                                                                                                                                                          <  2:331683/249‑1 (MQ=255)
        gCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGAtt                                                                                                                                                                                                                                          <  2:615111/249‑1 (MQ=255)
         cTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTc                                                                                                                                                                                                                                         <  2:133623/249‑1 (MQ=255)
         cTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCt                                                                                                                                                                                                                                        <  1:21201/250‑1 (MQ=255)
         cTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGATCGCGCTCCGCCGCCCCAGCCTTTACCCCAACCAAGCGGGCATAAAACTCCCCCCGAtttt                                                                                                                                                                                                                                        >  1:704866/1‑248 (MQ=255)
          tGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGTCTTGCGTCCGCCCCCCCCGCTTTAACCGCTACCCAGGCGGGAAGAATCACCCCCCCGAtatt                                                                                                                                                                                                                                        >  2:839410/1‑247 (MQ=255)
            gCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCttt                                                                                                                                                                                                                                      <  1:774025/249‑1 (MQ=255)
             cATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTcc                                                                                                                                                                                                                                    <  1:127745/250‑1 (MQ=255)
               ttAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGGTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTcccc                                                                                                                                                                                                                                  >  1:612760/1‑250 (MQ=255)
                     tCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGatat                                                                                                                                                                                                                          <  2:467673/250‑1 (MQ=255)
                         tttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGTTCCGCCCCCCCTGCTTTTAACCTACCCCGGCGGGCCGTAACCTCCCCCCCGTTCTTTTCCCCGAAATc                                                                                                                                                                                                                         >  1:666648/1‑247 (MQ=255)
                           tCCCCACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGGCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCgaga                                                                                                                                                                                                                     <  2:70460/249‑1 (MQ=255)
                           tCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTCCCCCCCCCCTGCCTTTACCCTACCCCAGCGGGGCGGAAACCCCCCCCCCTTTTTTCCCCCAgtatcccggt                                                                                                                                                                                                                      >  2:158319/1‑241 (MQ=255)
                                    tCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGGCGGAAACCCCGCCCGATTCTTTCCCCAGATATCGAGATgcccgcccc                                                                                                                                                                                                          >  1:53380/1‑249 (MQ=255)
                                         tGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTTTCGCGGTTCGCCCCCCCGTCTTTTAACCGTACCCAGCGGGGCGGTACCCCCCCCCCGCTTTTTTCCCCAGATACCGAGATCCAGcccccccct                                                                                                                                                                                                     >  2:932285/1‑247 (MQ=255)
                                           aagaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGTTCCGCCGCCCGTGCTTTTACCGTACCCAGGCGGGCCGGAACCTCCGCCCGATTCTTTCCCCCGATATCGAGATGCCAGCcccccggt                                                                                                                                                                                                     >  1:65674/1‑246 (MQ=255)
                                            agaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCGGGCTTTTACCGGACCCAGGGGGGGAGGAAACCCCCCCCGATTTTTTCCCCCAATATCGGGATGCCCACCCGCCGGTaaa                                                                                                                                                                                                  >  1:312147/1‑249 (MQ=255)
                                                  aaTCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTCCCCAGGCGGGAAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCaaa                                                                                                                                                                                             <  1:621666/249‑1 (MQ=255)
                                                          tCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGGTCCGCCGCCCGTGCTTTAAACGTACCCAGGCGGGGAGTAACCTCCGCCCGATTCTTTTCCCCGATATCGAGATGCCAGCCCGCCGTTAAGCCCAAAATACCGtt                                                                                                                                                                                    >  1:709872/1‑250 (MQ=255)
                                                               aataGCTGGCGATCGGTGGAGATAGTGGTTTTCTTCCCTTTTTCTTTTTGCTGAAAGACGGTGTTCAGCGCCGCCTGGGTTGGCCACCCCCAGTACGGGGGGCCTCCGCTGCCGCCCAGCAATTGGCGGTTTTGCCCGCTTCCCGCCGCCCGTGTTTTACCCGTACCCGGGGGGGCAGCACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGgtg                                                                                                                                                                                <  2:254032/248‑1 (MQ=255)
                                                                 taGCGGGCGTTCGGTCGCGATAATGTTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTCCCCCCCCCAGTACAGGCGGGCTACGCTGCCAGCCAGAAATTGGCGGTTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTCCCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTa                                                                                                                                                                              <  1:541182/249‑1 (MQ=255)
                                                                              gTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGGTCCGCCCCCCCTGCTTTTACCGGACCCAAGCGGGCAGTAACCCCCGCCCGATTTTTTCCCCAGATATCGGGATGCCCACCCCCCGGTTAACCCAAAAGAACGTTTGGGTACCACCGGGCGGTCg                                                                                                                                                               >  1:164765/1‑250 (MQ=255)
                                                                                 gcgATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCCCCGCCCGATTTTTTCCCCCGATATCGAGATGCCCGCCCGCCGTTTAACCCAAAATACCCTTTGTGTACCACCGGCCGGTCGTAc                                                                                                                                                            >  1:314653/1‑250 (MQ=255)
                                                                                  cgATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCCCCCCCCCTGCTTTTAACCTACCCCGGCGGGCGGTAACCTCCCCCCCCTTCTTTTCCCCGATATCTAGATGCCCGCCCCCCCTTTAACCCAAAGAACCCTTTGTGTACCACCGGCCGGGCGGAcc                                                                                                                                                           >  2:777367/1‑250 (MQ=255)
                                                                                    aTAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTTCGCCCCCCCCGGCTTTACCCGACCCCAGCCGGCCGGAACCCCCCCCCCCTTTTTTTCCCCAAATATCGGGAGCCCACCCCGCGCTACACCCACAAATACCGTGGTGTGCCACGGggcgggcggcacct                                                                                                                                                         >  2:458057/1‑244 (MQ=255)
                                                                                      aaTGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTGTCGCGTTTCGCCCCCCCGGCTTTTACCCCAACCCAGCGGGGCGTAACCCCCCCCCCAGTTTTTCCCCCGATATTCAGATGCCCCCCCGCCCGTTAGCCCAAAGTACCCCTTGGGTACCCCCGGgccggccggcccgcg                                                                                                                                                        >  2:679680/1‑243 (MQ=255)
                                                                                       gtGGTGTTCTGTTTTTTTTCTTTTTGCTGCAACACCTGTTTTACCGCCCCCTGGTTTTGCCCCCCCCAGTACAGGGGGGCTCCGCTGCCAGCCAGCAATGGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTACCCGTCCCCGGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGAGGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCg                                                                                                                                                      <  2:973337/248‑1 (MQ=255)
                                                                                           tgtTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAAACTCCCCCCGATTCTTTTCCCAGATATCGAGATGCCAGCCCGCCGTTTAACCAAAAATACCGTTGGGGTACCACGGGCCGGGCGTACCCGCGGCCgg                                                                                                                                                 >  1:43121/1‑250 (MQ=255)
                                                                                               cTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCCAAAGTACCGTTGGTGTACCACGGGCCGGGCGTACCTGCGGCCGGATCg                                                                                                                                             >  1:674953/1‑250 (MQ=255)
                                                                                                  ctctTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCGTGCTTTAAACGTACCCAGGCGGGGAGTAACCTCCGCCCGATTCTTTTCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGGGTACCACGGGGCGGGCGTACCTGCGGGCGGGTCGtt                                                                                                                                           >  1:189094/1‑249 (MQ=255)
                                                                                                     tttttCTATTTGCTGCAAGACCTTGTTTACCCCCGCCGGGTTTTGCCCCCCCCAGTACGGGGGGGCTCCGCTGCCAGCCAGCAATGGGCGGGTTTGCTCGCTTCCCGCCCCCCGTGTTTTACCCGTACCCGGGCGGGCAGTACCCTCCCCCCGTTTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTc                                                                                                                                          <  2:314653/247‑1 (MQ=255)
                                                                                                                       gACCGTGTTAACCGCCGCTGGGTTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATGGGCGGGTTTGCTCGCTTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGAAGTACCCCCCGCCCGATTCTTTCCCCAGATACCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAg                                                                                                                   <  1:199471/250‑1 (MQ=255)
                                                                                                                        cccGTTTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACGGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGTTTTACCCGTACCCAGGCGGGAAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGc                                                                                                                  <  1:164910/249‑1 (MQ=255)
                                                                                                                          cGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTCCCCCCCCCCCGCTTTTAACCTACCCCCGCGGGGGGGAACCCCCCCCCGCTTCTTTTCCCCGATATCAAGGTGCCCCCCCCCCCTTTAAACCACAACACCGCTGGTGTTCCACCGGCCGGGCGTACCCGCCGCCCGGGCGTTTACAGGACACGGACCCCTTTAAAcc                                                                                                                 >  2:237854/1‑250 (MQ=255)
                                                                                                                            tgtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCCTGCGGTTCGCCCCCCCCGGTTTTAACCTACCCCGGCGGGCGGGAAACCCCGCCCCCTTTTTTTCCCCAATATTCGGAATGCCGGCCCCCCTTTAGGCCAAAAGACCACTTGGGGAACACCGGGCGGGCGGGCCTGCGGGCGGGACCGTTCGAAGAAAAAGGCCCAATTAGACCtt                                                                                                               >  2:505369/1‑250 (MQ=255)
                                                                                                                             gtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTAAGCTGCCAGCCAGCAATTGGCGGGGTTGTTCGCCTTCCCCCCCCCCGGCCTTTACCCTATCCCGGGGGGGATGAACCCCCCCCCCCTTTTTTTCCCCGCTATCGCGATGGCCGGCCGCCGGTAAACCCACAGTAACCTTGGTGGTCCCCCGGCCGGGCCGACCCCCGGCCGGGTCGTTCAGAGCCAACCGGCCCACTacaacctct                                                                                                             >  2:31851/1‑242 (MQ=255)
                                                                                                                                tAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCCGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGcc                                                                                                          >  1:429318/1‑250 (MQ=255)
                                                                                                                                       ccgtGGGTTTGCCACCCCCATTAGGGGCGGGTTCGGCGCCCGCCCGTAAATGGGGGGGTTGGCCCCCGTCCCGCCCCCGTTGCTTTAACGGAACCCGGGGGGGCAGAACCCCCCCCCCTATTCTTCCCCAAAATATCGAGAGGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttt                                                                                                   <  2:128065/247‑1 (MQ=255)
                                                                                                                                           gtgTTTGCCCCCTCCCAGTACAGGGGGGTTCCTCGGTCAGCCAGCAATGGGCTGTTTTGCTCGCGTTCCCCCGCCAGTTTTCTAACCGCACCCAGGGGGGCATTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAtttttt                                                                                                <  2:592220/248‑1 (MQ=255)
                                                                                                                                              tttGCCACTCCCAGTACAGGCGGACTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCCAACGGCCCCATTAAGCCTTCTGCCGACCTTTTTTGCCg                                                                                            >  1:188256/1‑250 (MQ=255)
                                                                                                                                                      cccctgttccgggcggcCTACGCTGCCAGCCAGCAGTTGGGGGGTTTCCCCGCGTTCCCCCCCCCGTGCTTTACCCGCCCCCGGGCGGGCAGTACCCCCCCCCCGTTTCTTCCCCCAAATACGGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCata                                                                                     <  2:383607/239‑1 (MQ=255)
                                                                                                                                                       tCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCCGTCCGCCGCCCCCGCTTTTACCCGACCCCCGCGGGGCGTAACCCCCCCCCCCTTCTTTTCCCCCGTAACGAGGTGCCACCCCCCCCCTAAGCCCCAAATACCCTTGGTGTACCCCCGGGCCGGCGGACCCGCGGCCCGGGCGGTCCGAGCAAACGGCCCCCCTTAGCCCTCCGCCCCCCATTTTTGCCCCCCCataa                                                                                    >  2:463640/1‑249 (MQ=255)
                                                                                                                                                         ccAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTTCGCCGCCCCTGCTTTTACCGTACCCAGGCGGGGCGGAAACTCCCCCCCCTTTTTTCCCCCGATATCCAGATGCCCGCCCCCCCTTTAACCCAAAGAAACCGTGGGGGACCAAGGGGCCGGCGGACCTGCGGGCCGGGCGGTCCGAGCAAAACGGCCCCATTAACCCTCTGCCGCCATTTTTTTCCGCCCCCAATcc                                                                                  >  2:954632/1‑250 (MQ=255)
                                                                                                                                                                 ggCGGGCTCCGCTCCCACCCAGTATGGGGCGGGTTTCCGCCCTTCCCCCCCCCCGTGCTTTAACCCTCCCCGGGGGGGCATTCCCCCCCGCCCTTTTCTTCCCCCAGATCTCGAGAGGCCGCCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCTCTCCATATCCGCCGaaaa                                                                          <  2:831879/250‑1 (MQ=255)
                                                                                                                                                                     ggCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTCCCGCCCCCCTTGTTTTAACCGTACCCAGGGGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCg                                                                      <  1:35098/250‑1 (MQ=255)
                                                                                                                                                                       cTACGCGGCCAGCCAGAAATGGGCGGGTTTGCTCGCGTCCCGCCGCCCGTGCTTTACCCGTCCCCAGGCGGGCAGTACCCTCCGCCCGATTCTTCCCCAAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGAt                                                                    <  1:338237/250‑1 (MQ=255)
                                                                                                                                                                          cGCTGCCAGCCAGCATTGGGGGGGTTTGCTCGCGTCCCGCCGCCCGTGCTTTACCCGTACCCAGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTg                                                                 <  1:68234/250‑1 (MQ=255)
                                                                                                                                                                            cTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCg                                                               >  1:346503/1‑250 (MQ=255)
                                                                                                                                                                                  gccagcAATTGGCGGGTTTGCTCGCGTCCCGCCGCCCGTGCTTTAACCGTACCCGGGCGGGCAGTACCCTCCCCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCAt                                                         <  1:856154/250‑1 (MQ=255)
                                                                                                                                                                                           ttGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGAACGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTg                                                 >  2:219211/1‑250 (MQ=255)
                                                                                                                                                                                                          cggtgctCCGCCGCCGGTGCTAGACCCGTCCGGGGGCGCCAAGTCCCCCCCGCCGGATTCTCCCCCGAGAGCTCGAGAAGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCACATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCg                                    <  2:435771/244‑1 (MQ=255)
                                                                                                                                                                                                                     cgccCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGTGCCAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCTGCCGAACAGTCGGTCTGCGGATCATCCTGAGTGGCTACTGCCTTCGCCTGGGCtt                           >  2:252412/1‑248 (MQ=255)
                                                                                                                                                                                                                     cgccCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCtt                           >  2:224196/1‑248 (MQ=255)
                                                                                                                                                                                                                        ccGTGCTTAACCCGACCCAGGGGGGGAAGTACCCCCCCCCCTTTCTTTTCCCCATATATCGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagct                      <  1:974624/250‑1 (MQ=255)
                                                                                                                                                                                                                         ggtcctatAACCGCACCCGGGCGGTCTGCACCCCCCCCCCGTTTCTCCCCCCAGAACGCGAGTGCCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta                     <  2:749069/243‑1 (MQ=255)
                                                                                                                                                                                                                           cgctgcaaCGGTACCCGCGCGGGCTATCCCCCCCCCCCCATTCTGTCCCCGAGAAGCGAGACCCCCCCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta                     <  2:398962/242‑1 (MQ=255)
                                                                                                                                                                                                                              cttAACCCTTCCCAGGCGGGCAGTAACCCCCCCCCCGTTTCTTTCCCCAAATACCGGGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTg                 <  2:961913/249‑1 (MQ=255)
                                                                                                                                                                                                                                  tacccgaccccAGGCGGGCAGCACCCCCCGCCCGATTCTTCCCCCAGATATCGAGAGGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGaa               <  2:946384/241‑1 (MQ=255)
                                                                                                                                                                                                                                       gtccccGGGCGGGCAGTACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCt         <  1:857296/247‑1 (MQ=255)
                                                                                                                                                                                                                                           ccAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCGTATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGtt      >  1:325690/1‑249 (MQ=255)
                                                                                                                                                                                                                                           ccAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAGAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGtt      >  2:632570/1‑249 (MQ=255)
                                                                                                                                                                                                                                             cgggCGGGCAGTACCCCCCGCCCGATTCTTCCCCAAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCa    <  1:582264/248‑1 (MQ=255)
                                                                                                                                                                                                                                                gCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGc  >  1:889041/1‑249 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  
CGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGC  >  NZ_CP009273/2221508‑2221989

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