Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,221,7680TC36.6% 59.0 / 10.3 41pseudogene (470/1436 nt)mdtQmultidrug resistance outer membrane protein MdtQ
Reads supporting (aligned to +/- strand):  ref base T (21/5);  new base C (0/15);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 2.44e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGT  >  NZ_CP009273/2221523‑2222017
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                           
ttAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGGTCCGCCGCCCCGGCTTTTAACCGACCCCGGCGGGCAGGAAACTCCGCCCGATTCTTTtccc                                                                                                                                                                                                                                                         >  1:282536/1‑250 (MQ=255)
    tatCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGat                                                                                                                                                                                                                                                     <  2:212939/250‑1 (MQ=255)
      tCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCCCCGCCCGTGGCTTAACCCGACCCCGGCGGGCAGGAACCTCCGCCCGGTTCTTTTCCCAGatat                                                                                                                                                                                                                                                   >  1:6133/1‑250 (MQ=255)
        ggTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTTCGCCGCCCCCGCTTTTACCCTACCCCGGCGGGCAGTCACCTCCCCCCCCGTCTTTTCCCAGATATc                                                                                                                                                                                                                                                  >  2:259652/1‑249 (MQ=255)
         gTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCga                                                                                                                                                                                                                                                <  1:205075/250‑1 (MQ=255)
              cacaCCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATg                                                                                                                                                                                                                                           <  2:358904/249‑1 (MQ=255)
                  ccTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGc                                                                                                                                                                                                                                      <  2:24801/250‑1 (MQ=255)
                   cTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCCCGGTCCCCCCCCCCTGCCTTTACCGTACCCCGGCGGGGCGGAAACCCCCCCCCACTCCTTCCCCCGCTATCgcgctgaccacc                                                                                                                                                                                                                                      >  2:523276/1‑239 (MQ=255)
                          tGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTTCCCCGCCCCTGCTTTTAACCTACCCAGGCGGGCAGTAAACCCCCCCCCGTTCTTTTCCCCGATATCCAGATGCCACCCCCCCCtt                                                                                                                                                                                                                              >  1:75490/1‑250 (MQ=255)
                                               atatgtcTGGCGACCGGCCGCGATATTTGTTTTCTGTTTTTTTTCTATTTGCTGCAAGCCCGTTTTTACCGCCGCCTGGTTTGCCCCCCCCCATTCCGGGCGGGCTACGCTGCCCGCCAGCAATTGGCGGTTTTGCCCGCGTCCCGCCGCCCTTGCTTTAACCGTCCCCAGGGGGGCGGTACCCCCCGCCCTATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTg                                                                                                                                                                                                           <  2:132577/241‑1 (MQ=255)
                                                                    gATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCCGGCGGGCAGTAACCTCCGCCCGATTCTTTTCCCCGATATCCAGATGCCAGCCCGCCCTTTAGCCAAAAGTACCGCTTGTGTACCACCGGCCGGTCGTACCt                                                                                                                                                                                   >  1:64899/1‑250 (MQ=255)
                                                                                  gCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTCCCCCCCCCCCGCCTTTACCCTACCCCCGCGGGGCGGAAACCCCCCCCCCTTTTTTTCCCCCATATCCAGGAGGCCACCCGCCCCTTAGCCCAAAATACCGCTGGGGTTCCCCGGGGCGGGGGGACCCGCGGCCGGGTCggt                                                                                                                                                                     >  2:713605/1‑248 (MQ=255)
                                                                                   ctctTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTCCGCCCCCCCCGCTTTTACCCTACCCCAGGGGGGGGGGAACCCCCCCCCCCTTTTTTCCCCGATATTGGGATGGCCGCCCCCCCCTAAGCCCAAAAAAACCCTTGGGGTCCCACGGCCGGGCGGACCCCCGGGCGGGCCGTTc                                                                                                                                                                    >  2:52700/1‑250 (MQ=255)
                                                                                        tttCTATTTCCTGCAACACCGTTTTTACCCCCCCTGGGTTTGGCCCCCCCCAGTACGGGCGGCCTCCCCTGCCACCCAGCATTTGGCGGTTTTCCTCGCTTCCCCCCCCCCTTCTTTTACCCGTCCCCGGGCGGGCAGTACCCTCCCCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCagag                                                                                                                                                                <  2:260020/249‑1 (MQ=255)
                                                                                                ggcagcataccggtgtTCCACGCGGCCGTGGTTGGCCCCCCCAAGTCAGGGGGGGCTCCGCGGCCACCCAGCAGTGGGCGGGTTTGCTGGCGTCCCGCCGCCCTTTCTTTAACCGTCCCAGGGCGGGCGGTACCCCCCGCCCGATTTTTCCCCCAGATATCGGGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGc                                                                                                                                                        <  2:129358/237‑1 (MQ=255)
                                                                                                                     gccgccGGGTTTTGCCACCCCCAGTCCGGGGGGGCTCCGCTCCCAGCCAGCAATGGGCGGGTTTGCTCGCCTCCCGCCGCCCGTGCTTTAACCGTCCCCGGGCGGGCAGTACCCCCCCCCCGATTCTTCCCCCAAATATCGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACa                                                                                                                                  <  2:61076/250‑1 (MQ=255)
                                                                                                                                    ccACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTGTCGCGTTCCGCCCCCCCGGCTTTTACCCTACCCCGGGGGGGCGTAACCCCCCCCCCGTTTTTTCCCCCGATATTCCGATGGCCACCCCCCGTTCAGCCCAAAGAACCCGCGGTGTGACACGGGGCGGGCGTACCCCGGGCCGGGACGGTCCCAGGCAAAGGCCCCCCTAAGCCGTCGGCCCGCAtttttttgcgcccc                                                                                                                    >  2:375/1‑242 (MQ=255)
                                                                                                                                        ccccGGTACGGGGGGGCTCCGTGGCCCCCCCGCATTGGGGGGGTGTGCCCCCGTCCCGCCGCCCGTGCTTCACCCGCCCCCAGGGGGGCACTACCCCCCCCCCCTTTTTTCCCCCAGATAGGGAGCGGCCCGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCatat                                                                                                                <  2:575287/249‑1 (MQ=255)
                                                                                                                                                cAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCTCCGCCCGTGCTTTAACCGCACCCAGGCGGGGCGTTATCTCCGCCCGCCTCTTTCCCCAGATATCGAGATTCCAGCCCGCCGCTAAGCCCAAAGTACCGTCCGGGTACCACCGGCCCGTCGCACCTGCTGCCCGCTCGCTCAGAGCCAACAGACTCAGTAAGCCTTCCGCTGACATTTTTTGTCTCCCCATATccctccc                                                                                                         >  2:228711/1‑244 (MQ=255)
                                                                                                                                                     gggCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGCGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTc                                                                                                   >  1:319672/1‑250 (MQ=255)
                                                                                                                                                      tgctgcgctgcCACCCAGCATGGGGCGGTTTTCCCCGTCTCCCGCCCCCGTTGTTTAACCCGTCCCGGGGGGGGCAGCACCCCCCCCCCGTTTTTTCCCCCATATTTCGAGATCCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCg                                                                                                  <  2:645571/249‑1 (MQ=255)
                                                                                                                                                           cGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGGTTCTTTCCCCAGATATCGAGATGCCAGCCCCCCGTTAAGCCCAAAGTACCGTTGGTGTACCACCGGCCGGTCGTACCTGCGGCCCGAACGGTCAGAGCAAACGGCCCCACTAAGCCTTTTGCCGACATTTTTTGCCGCCCCATATCCGCCGGAAAGTCCATCTg                                                                                             >  2:175941/1‑250 (MQ=255)
                                                                                                                                                                      agcAATGGGCGGTTTTCCTCGCGTTCCGCCCCCCGTGCTTTACCCGTCCCCAGGCGGGCAGTAACCCCCCCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCt                                                                                  <  1:713605/250‑1 (MQ=255)
                                                                                                                                                                               cGGTTTTGCTCCCGTCCCGCCCCCCGTTCTTTACCCGTCCCCAGGCGGGCATTCCCCCCCCCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCaa                                                                         <  1:350026/250‑1 (MQ=255)
                                                                                                                                                                                 ggTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACt                                                                       >  1:490442/1‑250 (MQ=255)
                                                                                                                                                                                 ggTTCGCTCGCGTTCCCCCGCCCGTTCTTTACCCGTCCCCGGGGGGGCAGTACCCTCCCCCCGATTCTCTCCCCAGATAGCGAGCGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACt                                                                       <  2:19021/250‑1 (MQ=255)
                                                                                                                                                                                  gtttcgcgccgTTCCCCCCCCCCTTCTTTCACCCGACCCGGGGGGGGCAGTCCCCCCCCCCCGATTCTTCCCCCAGATAGCGAGATGCCCGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCAGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACt                                                                       <  2:258393/241‑1 (MQ=255)
                                                                                                                                                                                     tGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCt                                                                   >  1:35552/1‑250 (MQ=255)
                                                                                                                                                                                         cgcgTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGc                                                               >  1:549032/1‑250 (MQ=255)
                                                                                                                                                                                         cgcgTCCCCCCCCCCGTGTTTAACCCGTCCCCGGGCGGGCAGAACCCCCCGCCCGTTTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGc                                                               <  1:146357/250‑1 (MQ=255)
                                                                                                                                                                                           cgTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCt                                                             >  2:37039/1‑250 (MQ=255)
                                                                                                                                                                                                   cgccCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCtt                                                       >  1:684435/1‑248 (MQ=255)
                                                                                                                                                                                                     cccGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagc                                                   >  1:694690/1‑250 (MQ=255)
                                                                                                                                                                                                                   ccGTACCCAGGCGGGCAGTAACCTCCGCCCGCTTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGACAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATgcg                                       >  2:478003/1‑250 (MQ=255)
                                                                                                                                                                                                                      tACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTg                                    >  2:460680/1‑250 (MQ=255)
                                                                                                                                                                                                                         cccAGGCGGGCATAACCCTCCCCCCTATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGtt                                  <  1:606871/250‑1 (MQ=255)
                                                                                                                                                                                                                            cgggCGGGCAGTACCCCCCCCCCGATTCTTCCCCCAAATACCGGGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCa                                <  1:583754/248‑1 (MQ=255)
                                                                                                                                                                                                                                cGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCa                            >  2:678049/1‑250 (MQ=255)
                                                                                                                                                                                                                                          cccccccccgtaTTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCatat                   <  1:367930/239‑1 (MQ=255)
                                                                                                                                                                                                                                               cccccGTTTCTTCCCCCAGATATCGAGATGCCCGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGt              <  1:125245/248‑1 (MQ=255)
                                                                                                                                                                                                                                                        cTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTg      <  1:333118/249‑1 (MQ=255)
                                                                                                                                                                                                                                                           tCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTCCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCCTATCCGGTGCCTTCTGTAGt  >  2:125560/1‑250 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                           
TTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCCGGTGCGTTCTGTAGT  >  NZ_CP009273/2221523‑2222017

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