Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,221,7540TC43.9% 57.9 / 18.1 57pseudogene (484/1436 nt)mdtQmultidrug resistance outer membrane protein MdtQ
Reads supporting (aligned to +/- strand):  ref base T (22/10);  new base C (8/17);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 7.98e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATC  >  NZ_CP009273/2221517‑2222001
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                
cTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACCCTTCCAACCAGCAATTGGCGGGGTTGGCCCCGTTCCGCCGCCCGCGCTTTCACCGTACCCAACCCGGCAGCAAACTccccccccttcc                                                                                                                                                                                                                                                     >  1:159378/1‑244 (MQ=255)
 tGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCCGGCTTTTAACCTAACCCGGCGGGCGGTAACCTCCCCCCGGTTCtt                                                                                                                                                                                                                                                   >  1:597713/1‑250 (MQ=255)
  ggCATAAAAATCGGCTTCCACTCCGTCAACGAAAGAAGTAATCCCGTTGTGAGATTGCTGGCGATCGGTCGCGATAAGGGTGTTCGGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTGGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTGCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCtt                                                                                                                                                                                                                                                   <  2:887030/249‑1 (MQ=255)
   gCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCttt                                                                                                                                                                                                                                                  <  2:225947/249‑1 (MQ=255)
   gCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTc                                                                                                                                                                                                                                                 <  2:766030/250‑1 (MQ=255)
          tatCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTGTCCGGGTTCGCCCCCCCCGCCTTTACCCTGCCCAGGGGGGGCGTAACCCCCGCCCCGCTTTTTTcccccat                                                                                                                                                                                                                                           >  2:1407198/1‑247 (MQ=255)
            tCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTTCCCGCCAGCAATTGGCGGGGTTTCTCGCGCTCCGCCCCCCGCGCTCTTACCCTACCCCCGCGGGCAGCATACCCCCCCCCCCTCTTTCCCCCCatat                                                                                                                                                                                                                                       >  2:1269976/1‑250 (MQ=255)
                 ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTGTCGCGGTCCGCCGCCCCGGCTTTTACCCTAACCCGGCGGGCCGGAACCCCCCCCCGATTCTTTCCCCCGATATc                                                                                                                                                                                                                                      >  1:289262/1‑246 (MQ=255)
                      cacCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTATTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCa                                                                                                                                                                                                                            <  1:950089/250‑1 (MQ=255)
                            cAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGccc                                                                                                                                                                                                                       <  2:574965/248‑1 (MQ=255)
                            cAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTCCGCCCCCCCTGCTTTTACCCGACCCAGGCGGGCCGTAAACCCCCCCCGCTTCTTTCCCCCGATATCGGGATGCccgcccgc                                                                                                                                                                                                                     >  1:9675/1‑250 (MQ=255)
                            cAACTGAAGAAGTAATCCCGTTCTGATATAGCGGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGccc                                                                                                                                                                                                                       <  2:153170/248‑1 (MQ=255)
                                                  cTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGGTCGCGGTCCGCCCCCCCTGCTTTTACCCTACCCCGGCGGGCCGTAAACCCCGCCCGGTTTTTTTCCCCGATATCGGGATTGCCGCCCGCCGTTAAGCCCAAAATAACCTTg                                                                                                                                                                                             >  1:642812/1‑250 (MQ=255)
                                                              gCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCGCCGCCCGGGCTTTAACCGTACCCAGGCGGGCCGTAAACTCCGCCCGATTCTTTTCCCCGATATCGAGATGCCCACCCGCCCGTAAGCCAAAAGCACCGTTGGGGTACCACagg                                                                                                                                                                                 >  1:438950/1‑250 (MQ=255)
                                                                 gtcggtcggTATAATGGTGTTCGTCTTTTTTTCTATTGGCTGCAAGACCGTTTTAACCGCCGCGGGGTTTTGCCCCCCCCAGACCGGGCGGGCCACGCTCCCGGCCAGCATTGGGGGGTTTTGCCCGCGTCCCGCCCCCCGTCCTTAACCCGTCCCAGGGGGGGCGGAACCCCCCGCCCGATTTTTCCCCCAGAAATCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCg                                                                                                                                                                              <  2:601068/249‑1 (MQ=255)
                                                                    ggTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGGTCGCGGTCCGCCCCCCCCGGTTTTACCGTAACCCAGCGGGCCGCAAACTCCCCCCCGTTCTTTCCCCCGCTATACAGAAGCCCGACCCCCACTTAACCCAAAAAAACCTGTGGGGACCACGGGCCgggc                                                                                                                                                                           >  2:1209340/1‑248 (MQ=255)
                                                                      tgccgttAATGTTGTTCTGCTCTTTTCCTTTTGGTGGAAAACCGTGGTTAACCGCCGCTGGGTTGGACCCCCCCGAAGAGGGCGGGTCTAGGCGGCCGGCCGGCAGTGGGGGGGTTTGCCCGCGTCCCGCCGCCCGTGCTTAACCCTTACCCAGGCGGGGAATACCCCCCGCCCGATTCTTCCCCCAAAAATCGGGATGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGt                                                                                                                                                                        <  2:562387/244‑1 (MQ=255)
                                                                                     tCTGCTTTTTTTCTTTTGCCGACAACCCGGTTGTTACCGCCCCTGGGTTGTCCCTCCCCCAGTACGGGCGGTCTCCGCTGCCACCCCGAAATTGGCGGGTTTGCCCGCTTCCCGCCGCCCGTTTTTTACCCGTACCCGGGCGGGCAGTACCCTCCCCCCGTTTCTTCCCCCATATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGAt                                                                                                                                                          <  2:1497607/249‑1 (MQ=255)
                                                                                                            ccaacacggttttAACCCCCCCTGGGGTTGCCCCCCCCCAGTACGGGGGGTCTCCGCCCCCGGCCAGAATTGGGCGGTTTTGCCCGCCTCCCGCCCCCCGTTTTTTACCCGTCCCCGGGCGGGCAGTACCCTCCCCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCAt                                                                                                                                   <  2:22885/240‑1 (MQ=255)
                                                                                                            cAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCGCGCGTTCCCCCCCCCCCGCGTTTACCGCACCCCAGCGGGCCGTAACCCCCCCCCCAGTTTTTTCCCCAATAACGAGATGCCCACCCCCCCGTAAGCCAAAAGTATCGCTCGGTGACCCCGGGCCGGGCGGTCCTTGCGGCGGGTCCGTTCGTGCAGAACGCCCCCAt                                                                                                                                   >  2:685144/1‑250 (MQ=255)
                                                                                                                acctgttttcagcgccgcCTGGGTTGCCCACCCCCAGTACGGGCGGCTTCCGCTGCCCGCCAGAAATGGGGGGGTTTCCTCCCGTCCCCCCCCCCGTGCTTTAACCGCACCCAGGCGGGCAGAACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAg                                                                                                                               <  2:759048/238‑1 (MQ=255)
                                                                                                                 ccGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCCACGCCGCCAGCCAGCAAATGGCGGGGTTTCTTTCGGTCCCCCCCCCCCGCTTTAACCCCACCCCAGCGGGGCGGAACCCCCCCCCCCATTTTTTCCCCAGGAATGCGATGACAACCCCCCGCTTAACCCAAAAAATCCCTGGGGTACCCCGGGCCGGTCCCACCCTCTTCCCGCGCGTTCCGAGCAACCGGCCCCACTAAgac                                                                                                                             >  2:960047/1‑248 (MQ=255)
                                                                                                                                      ggttgcgcccccccccGTCCAGGGGGGCTCCGCGCCCCGCCGGCATTGGGCGGGTTTCCTCGCGTCCCGCCCCCGTTGCTTTACCCCCCCCGGGGCGGGAAGTCCCCCCCCCCCTGTTCTTTCCCCAGATATGGAGGGGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAtttttt                                                                                                            <  2:1507852/237‑1 (MQ=255)
                                                                                                                                             ccACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTTCGCCGCCCGGGGTTTTAACGTAACCCGGGGGGGAGGAACCTCCCCCCCGATCTTTCCCCCGATAACGGGATGGCCAGCCCCCGGTTAGCCAAAAAGAACGGTGGGGGTCCACCGGGCGGGGCGAACCCGCGCCGGGGCGGTTAGAGGACAACGGGCCCCTAAACCTTCCTGCGGCCAtttttgtcgcccc                                                                                                     >  2:499426/1‑242 (MQ=255)
                                                                                                                                                                ggCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCCCCGCCCGCTTCTTTCCCCCGATATCGAGATGCCAGACCGCCCTTAAGCCCAAAATACCGTTGGTGTACCACGGGCCGGGCGGACCCACGGCCGGATCGGTCAGAGCAAACAGCGCCATTTAGCCCTCTGCCCGCATTTTTTGTCGCCCCATCTCCGCCGCAAAgac                                                                                   >  2:320890/1‑247 (MQ=255)
                                                                                                                                                                      cGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGTTCCGCCGCCCGTGCTTTTACCGTACCCCGGCGGGCAGTAACCTCCGCCCGCTTCTTTCCCCAGATATCGAGATGCCAGCCCCCCCTTAAGCCCAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGGTCGTTCCGAGCAAACGGCCCCCTTAAGCCTTCTGCCGCCATTTTTTTCCCCCCCCTCTCCCCCCAAAAGTCGATCtt                                                                             >  2:976221/1‑249 (MQ=255)
                                                                                                                                                                          gccagccaACAATTGGCGGGTTTGCTCGCGTCCCGCCGCCCGTGTTTTACCCGCCCCCAGGCGGGAAGAACCCCCCGCCCGATTCTTCCCCCAGATATCGAGAGGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGa                                                                         <  1:856397/250‑1 (MQ=255)
                                                                                                                                                                                 agcAATGGGCGGTTTTGCGCGCGTCCCCCCGCCCGTGCTTTACCCGTCCCCAGGCGGGCAGTACCCCCCGCCCAATTCTTTCCCCAGATATCGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCt                                                                  <  1:75299/250‑1 (MQ=255)
                                                                                                                                                                                   cccttGGCGGGTTTCGCCCGCTTCCCCCCCCCCTTTTTTAACCCGACCCCGGGGGGGCAGTACCCCCCCCCCTTTTTTTCCCCCAATTTCCGGGTGCCCGCCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTg                                                                 <  2:289262/247‑1 (MQ=255)
                                                                                                                                                                                      ttGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCCGCCCGCCGTTAAGCCAAAAATACCGGTTGTGTACCACCGGCCGGGCGTACCTGCGGGCGGAACGGTCAGAGCAAAAGGCCCCATTTAGCCTTCTGCCGACATTTTTTGCCGCTCCATATTCCCCGCAAAGTCGATCTGCGGACCATCCTGCGTg                                                             >  2:526205/1‑250 (MQ=255)
                                                                                                                                                                                          ggctgtttgccgcgcgTCCCGCCCCCCGCTTTTTACCCGCCCCCGGGCGGGCAGTACCCCCCCCCCTTTTCTCTCCCCATATCTGGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTgg                                                            <  2:329817/237‑1 (MQ=255)
                                                                                                                                                                                              gTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTg                                                      >  1:735847/1‑250 (MQ=255)
                                                                                                                                                                                              gTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTg                                                      >  1:602011/1‑250 (MQ=255)
                                                                                                                                                                                                            ccgccgCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGGCGATCTGCGGACCATCCTGAGTTGCAACTGCCTTCGCCTTGGc                                         >  1:1536632/1‑250 (MQ=255)
                                                                                                                                                                                                            ccgccgCCCGTCTTTTACCCGTCCCCAGGCGGGAAGTACCCCCCCCCCGTTTCTTCCCCCAGATTCCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGc                                         <  1:362212/250‑1 (MQ=255)
                                                                                                                                                                                                               ccgccCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCtt                                       >  1:32922/1‑249 (MQ=255)
                                                                                                                                                                                                                  cccGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagc                                   >  1:930438/1‑250 (MQ=255)
                                                                                                                                                                                                                     gTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagctag                                >  2:1401924/1‑250 (MQ=255)
                                                                                                                                                                                                                          ttACCCGTCCCCAGGCGGGCAGTACCCCCCCCCCGATTCTTCCCCCAAATAGCGAGACGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGa                            <  1:1025046/249‑1 (MQ=255)
                                                                                                                                                                                                                            aaCCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATg                         >  1:1176599/1‑250 (MQ=255)
                                                                                                                                                                                                                                  cgtccccGGGGGGCAAGAACCCCCCCCCCGTTTTTTCCCCAAGATATGGAGAGGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATgcgc                      <  1:163534/247‑1 (MQ=255)
                                                                                                                                                                                                                                    cccagggggggCAGTACCCTCCCCCCGTTTCTCCCCCAAGAAGCCGACATGCCGCCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGt                   <  2:422598/242‑1 (MQ=255)
                                                                                                                                                                                                                                     cccAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGtt                  >  1:1418594/1‑250 (MQ=255)
                                                                                                                                                                                                                                       ccggCGGGTAGTAACCCCCCGCCATTTTCTTCCCCCAAATATGTAGTTGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTc                 <  1:1296918/248‑1 (MQ=255)
                                                                                                                                                                                                                                       cAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGtt                  >  1:1266206/1‑248 (MQ=255)
                                                                                                                                                                                                                                       cAGGCGGGAAGTACCCCCCCCCCGATTCTTTCCCCAGATATCGAGATCCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCa                <  1:924321/250‑1 (MQ=255)
                                                                                                                                                                                                                                         ggCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGc              >  1:172255/1‑250 (MQ=255)
                                                                                                                                                                                                                                         ggCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGc              >  1:1414574/1‑250 (MQ=255)
                                                                                                                                                                                                                                         ggCGGGATGTACCCCCCCCCCGATTTTTCCCCCAGATGCCGAGATGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGc              <  1:272619/250‑1 (MQ=255)
                                                                                                                                                                                                                                           aggcaaGTAACCCCCGCCCGATCTTTCCCCCAGATACCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCa            <  1:1429358/245‑1 (MQ=255)
                                                                                                                                                                                                                                            ggggAGGAACCCCCGCCAGTTTCTTCCCCCAGGTATGGGGATACCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCAc           <  1:1239787/250‑1 (MQ=255)
                                                                                                                                                                                                                                             ggCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCAcc          >  1:1217233/1‑250 (MQ=255)
                                                                                                                                                                                                                                             agCAGTACCCTCCGCCTGATTCTTGCCCCAGATAACGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCAcc          <  1:1539345/249‑1 (MQ=255)
                                                                                                                                                                                                                                               cAGTACCCCCCGCCCGTTTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTg        <  1:1407198/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  tAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCa      >  1:425888/1‑249 (MQ=255)
                                                                                                                                                                                                                                                   aaCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCata    >  1:759273/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      ccTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGCGTGGCCACTGCCTTCGCCTGGTCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCGGCatacc  >  2:161295/1‑248 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                
CTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATC  >  NZ_CP009273/2221517‑2222001

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