Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,221,7570GC35.0% 62.9 / 10.1 40pseudogene (481/1436 nt)mdtQmultidrug resistance outer membrane protein MdtQ
Reads supporting (aligned to +/- strand):  ref base G (16/10);  new base C (6/8);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCC  >  NZ_CP009273/2221509‑2222002
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        
gggTTTTGCTGGCATTAATAACGGTTTCCGCACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGc                                                                                                                                                                                                                                                         <  1:974726/250‑1 (MQ=255)
  gTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTGCCCGGGTTCGCCCCCCCGGCCTTTAACCGACCCAGGCGGGGCGCAACcccccccc                                                                                                                                                                                                                                                       >  1:1136509/1‑250 (MQ=255)
  gTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTTCCGGTCCGCCCCCCGGGCCTTTACCGTACCCCGGCGGGCCGGAACCTcccccc                                                                                                                                                                                                                                                       >  1:823239/1‑250 (MQ=255)
   ttttGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCTCGTTTCGCCCCCCCAGTTTTTACCGCATCCACGCGGGCAGCAACCTccccccc                                                                                                                                                                                                                                                      >  1:376129/1‑249 (MQ=255)
        cTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCt                                                                                                                                                                                                                                                 <  1:270076/250‑1 (MQ=255)
                      ggTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTTTCGCGGTTCGCCCCCCCCGCTTTTACCGCACCCCGGGGGGCGGCAAACCCCCCCCCCCTTTTTTccccggtattc                                                                                                                                                                                                                                    >  2:575647/1‑243 (MQ=255)
                       gTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTTGCGGGGTTTTCTCGCGTCCCCCGCCCCGCGCTTTACCCTCACCCGGCGGGGCGTAATCTCCCCCCCGCTCTTTtccccccccatca                                                                                                                                                                                                                                  >  2:797832/1‑241 (MQ=255)
                              cacCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCCGGCGGGCAGTAACCTCCCCCCGGTTTTTTCCCCCGATATCgaaaatccc                                                                                                                                                                                                                           >  1:1188483/1‑243 (MQ=255)
                                 cTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGcc                                                                                                                                                                                                                       <  1:549184/250‑1 (MQ=255)
                                 cTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGGTCCCCCCCCCCTGCTTTTACCGGACCCAGGCGGGCAGTAACCTCCGCCCCGTTCTTTTCCCCGATATCCAGATTCCAGcc                                                                                                                                                                                                                       >  1:376227/1‑250 (MQ=255)
                                    cAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCccg                                                                                                                                                                                                                     <  1:188129/249‑1 (MQ=255)
                                                             tattggcgggCGTTCGTTCGCGATAATGTTTTTCGGCTCTTTTTTTATTTGCTGCAAGCCCGTTTTTACCGCGGCCTGGGTTGCCCCCTCCCGGTACGGGCGGGCTCCCCTGCCAGCCAGAATTTGGCGGGTTTGCCCGCTTCCCGCCGCCCGTGCTTTACCCGTCCCCAGGGGGGCAGTACCCCCCGCCCGATTTTTCCCCCAGATATCGAGTTGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGgtg                                                                                                                                                                                           <  2:482572/242‑1 (MQ=255)
                                                                    tGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTCCGCCGCCCCCGCCTTTAACCGACCCCGGCGGGGAGTAACCCCCGCCCCGCTCTTTTCCCCCATATCCAGATGGCAGCCCCCCCCTTACACCAAAGTACCCTTTGTGTTCCCCg                                                                                                                                                                                    >  2:164539/1‑250 (MQ=255)
                                                                    tGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGATGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTCCCCCGCCCCCGCCTTTACCCGCACCACGGGGGGAGGGAACCCCCCCCCCTTTTTTCCCCCCATATCGCGAAGGCAGCCCCCCCGTTAGCCCAAAAAACCGGTGGTGTagcccg                                                                                                                                                                                    >  2:356713/1‑245 (MQ=255)
                                                                                 cgATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTGTCGCCGTCCGCCCCCCCTGCTTTTACCCTACCCCAGCGGGGCGGAAACCCCCCCCCCTTCTTTTCCCCACAATCCGGAAGGCCACCCCCCCCTTAACCCAAAAAACCCCTGGGGGTACCACAGGCgggcgtcac                                                                                                                                                                        >  2:1041759/1‑244 (MQ=255)
                                                                                                      tttCTATTTACTGCAAGACCGTGTTTAACGCCGCTTGGGTTTCCCACCCCCGGTACGGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGTTTTAACCGTCCCAAGGGGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCagag                                                                                                                                                  <  1:273065/249‑1 (MQ=255)
                                                                                                                     agcccGTGTTAACCGCCGCCGGGGTTTGCCACTCCCAGAACGGGGGGGCTCCGCGGCCAGCCAGCAGTGGGCGGTTTTGCCCGCGTCCCGCCGCCCTTGCTTTACCCGTACCCAGGCGGGCAGTACCCCCCGCCCTATTCTTCCCCCAGATATCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTaa                                                                                                                                  <  2:32964/247‑1 (MQ=255)
                                                                                                                      gACGGTGTTTACCGCCGCCTGGGTTGCCCACCCCCAGTCCAGGGGGGCGCCGCGGCCAGCCGGCAATGGGCGGGTTTGCTCGCTTCCCCCCCCCCGTTTTTTACCCGTCCCCGGGCGGGCAGTACCCCCCGCCCGTTCCTTCCCCCAGATATAGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAg                                                                                                                                 <  2:120078/250‑1 (MQ=255)
                                                                                                                                     cccccgggTTTGCCCCCCCCCAGGCCGGGCGGCTTCCGCTGCCGGCCGTCATTGGGGGGTTTGGCCCGCTTCCCCCCCCCCTTGTTTAACCGTCCCCCGGGGGGCAGGTACCCCCCCCCCTTTTTTTCCCCCAGTAAGGGAGATGCAGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAt                                                                                                                   <  2:605255/245‑1 (MQ=255)
                                                                                                                                          ggttttGCCGCCCCCAGTACAGGCGGGCTCCGCTGCCAGCCAGCAATTGGGGGGTTTGCCCGCGTCCCGCCCCCCGTGCTTTAACCGTCCCCAGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTg                                                                                                             <  1:410122/247‑1 (MQ=255)
                                                                                                                                                 gccccccccAGCACGGCCGGCCTCCGCTGCCAGCCAGCATTGGGGGGGTTGCCCCCCGTCCCGCCCCCCGTGTTTAACCCGTACCCGGGCGGGCAGCACCCCCCCCCCTATCCTTCCCCCAGATACCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTc                                                                                                       <  2:258276/244‑1 (MQ=255)
                                                                                                                                                   cACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCCCCCGGTTCGCCCCCCCCGCTTTTACCCTAACCCCGCGGGGCGGCACCCCCCCCCCCCTCCTTCCCCCCAAAATCAGACGGCCAGCCCCCCGTTAACCCAAAACACCCCTGGGGTCACCCAGGGCGGGGGGACCCCCGGCCCGCACCGTTCGGGCCAAACGGCCCCATTCAGCCTTCCGCCGCCCTTTTTTCCCCTCCa                                                                                                     >  2:64010/1‑250 (MQ=255)
                                                                                                                                                   cACCCCCAGTAAGGGCGGGCTCCGCTGCCAGCCAGCAATTGGCGGTTTGGCTCGCTTCCCGCCGCCCGTGTTTTAACCGTACCCAGGCGGGAAGTAACCCCCGCCCGTTTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCa                                                                                                     <  1:177637/250‑1 (MQ=255)
                                                                                                                                                                           cTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCg                                                                            >  1:810852/1‑250 (MQ=255)
                                                                                                                                                                           cTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCCGATATCGAGATGCCAGCCCCCCGTTAAGCCAAAAGTACCGTTGGGGTACCACGGGCCGGGCGTACCTGCGGCCGGATCGTTCCGAGCAAACGGCCCCACTACGCCCTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAATCGACCTGCg                                                                            >  2:35204/1‑250 (MQ=255)
                                                                                                                                                                                                 gTTTCCTCGCTTTCCCCCCCCCGTGTTTTACCCGTACCCAGGCGGGAAGTACCCTCCCCCCGTTTTTTCCCCCAAATATGGGGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAATGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTg                                                       <  1:606910/250‑1 (MQ=255)
                                                                                                                                                                                                  gtttcccTCGCTTCCCCCCCCCCGTCCTTAACCCTACCCCGGGGGGGCAGTCCCCCCCCCCCTTTTCTTCCCCAGGTTCCCGAGGTGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACt                                                        <  2:964753/244‑1 (MQ=255)
                                                                                                                                                                                                   gttgccccgcgTCCCGCCGCCGGTGCTTAAACCGTCCCCGGGCGGGCAGTAACCTCCCCCCGATTTTTCCCCCAGATAGCGAGAGCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTg                                                       <  1:378126/243‑1 (MQ=255)
                                                                                                                                                                                                     tGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCt                                                    >  1:480946/1‑250 (MQ=255)
                                                                                                                                                                                                        cTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTc                                                  >  1:769700/1‑250 (MQ=255)
                                                                                                                                                                                                         tCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCg                                                 >  1:313270/1‑250 (MQ=255)
                                                                                                                                                                                                                 cgccgcCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCt                                         >  2:192578/1‑250 (MQ=255)
                                                                                                                                                                                                                  cccgccCGTCCTTAACCCGTCCCCGGGCGGGCAGTACCCCCCCCCCGATTTTTCCCCCAGATACCGAGAGGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCtt                                        <  1:164539/249‑1 (MQ=255)
                                                                                                                                                                                                                   ccgccCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTc                                       >  2:636046/1‑250 (MQ=255)
                                                                                                                                                                                                                    cgccCGTGCTTTAACCGTCCCCAGGCGGGGAATACCCCCCCCCCGATTCTTCCCCCAGATACCGAGATCCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCa                                      <  1:606192/250‑1 (MQ=255)
                                                                                                                                                                                                                        cGTGTTTAACCCGTCCCCAGGGGGGAAGTACCCCCCCCCCGATCCTTCCCCCAGATATCGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta                                  <  1:390268/250‑1 (MQ=255)
                                                                                                                                                                                                                        cGTGCTTTCACCGTACCCAGGCGGGCAGCAACCTCCGCCCGATTCTTTCCCCAGATATCGTGACGCCACCCCGCCTTTACGCCAAAAGTACCGTTGGCGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCCAACGGCCCCCTTAAGCCTTCTGCCGACCTTTTTTGCCGCTCCATACCCTCCGGAAAGGCGATCCCGCGGCCCTCCTGCGTTGGCACTGCCGTCGCCtcgcactccgccc                                  >  2:470434/1‑238 (MQ=255)
                                                                                                                                                                                                                        cGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagc                                    >  2:645250/1‑248 (MQ=255)
                                                                                                                                                                                                                             cTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGCGTACCACTGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGTCCCATTAAGCCTTCTGCCCACCTTTTTTGCCCCTCCATAACTGCCGAAAAGTCGCTCTGCCGCACATCTTGAGTGGCCACTGCCTTCGACCTGGCTTTAGCTAGCTg                              >  2:658497/1‑250 (MQ=255)
                                                                                                                                                                                                                                tAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCCGGGCTTCAGCTAGCTGAAt                           >  1:844599/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    cAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCCTATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCCACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTg         >  2:997035/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     aGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGACCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTAGATCTGCGGACCATCCTGAGAGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTg         >  1:1031719/1‑249 (MQ=255)
                                                                                                                                                                                                                                                             cTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCCTATcc  >  1:987460/1‑250 (MQ=255)
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        
GGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTGTTCAGCCACCTGCATATCC  >  NZ_CP009273/2221509‑2222002

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: