Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,590,063 A→G 16.8% E824G (GAG→GGG)  acrD → aminoglycoside/multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,590,0630AG16.8% 155.0 / 3.5 71E824G (GAG→GGG) acrDaminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base A (23/33);  new base G (11/4);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 4.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACACTGCAAACCGCCTGGGGTTCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTA  >  NC_000913/2589831‑2590289
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                  
acacTGCAAACCGCCTGGGGTTCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGGGAttgt                                                                                                                                                                                                                              >  2:643167/1‑239 (MQ=255)
  acTGCAAACCGCCTGGGGTTCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGGGAttgttg                                                                                                                                                                                                                            >  2:54364/1‑239 (MQ=255)
  acTGCAAACCGCCTGGGGTTCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGAttgttg                                                                                                                                                                                                                            >  2:697399/1‑239 (MQ=255)
    tGCAAACCGCCTGGGGTTCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTCTTCTGCGG‑TTGGGATTGTTggg                                                                                                                                                                                                                          >  2:764676/1‑239 (MQ=255)
    tGCAAACCGCCTGGGGTTCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGGCTGGAACGCTATAACGGTTATTCTGCGG‑TTGGGTTTGTTggg                                                                                                                                                                                                                          >  2:435318/1‑239 (MQ=255)
            gCCTGGGGTTCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGGGCCCTTCTCTGCTTTCGCGACCTCCCGCTGGGGAACAGG‑CTGCCCGCGTCTGGAAGCCCATAACGGTTTTTTTGGGG‑TTGGGGTTTTTGGGGAAGccc                                                                                                                                                                                                                   >  1:736728/1‑237 (MQ=255)
             ccTGGGGTTCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGc                                                                                                                                                                                                                  >  2:790878/1‑238 (MQ=255)
                ggggTTTGGGTTTTGTGAATGACTTTTTGGATCGCGGTCGCGTGAAAAAAGTCTATGTGCGGGCAGCTGCGCCGTACCGCATGCTGCCGGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCg                                                                                                                                                                                                              <  1:435318/239‑1 (MQ=255)
                ggggTTCGAGCTATGTGAATGGATTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTAGCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCg                                                                                                                                                                                                              <  2:443053/239‑1 (MQ=255)
                 gggTTCGGGTTTTTTGAATGTCTTTATGGATGGCGGTCGCGTGAAAAAAGTCTTTGTGCAGGCACCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCGTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCg                                                                                                                                                                                                              <  1:591438/238‑1 (MQ=255)
                 gggTTCGAGTTTTGTGAATGACTTTTTGGATCGGGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCg                                                                                                                                                                                                              <  1:581401/238‑1 (MQ=255)
                     ttcgtgtgtgtgAAGGACTTTAGGGATCGCGGGGGCGTGAAGAAAGTCTGTGTGCAGGCAGCTGGGCCGTATCGCATGCTGCCAGATGACATCAATCCCTGGTATGTCCGAAATAAAGAGGGCGGCAGGGTGCCCTTCTCTGGTTTCGCGACCTCCCGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGt                                                                                                                                                                                                         <  1:774073/231‑1 (MQ=255)
                     tCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGt                                                                                                                                                                                                         <  2:579939/239‑1 (MQ=255)
                     tCGAGCTATGTGAATGACTTTAGGGATCGGGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACCTCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGt                                                                                                                                                                                                         <  2:361013/239‑1 (MQ=255)
                         gCTATGTGAATGGCTTTAGGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGAATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAg                                                                                                                                                                                                      <  2:533260/238‑1 (MQ=255)
                            aTGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTa                                                                                                                                                                                                    >  2:581401/1‑237 (MQ=255)
                             tgtgAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCg                                                                                                                                                                                                 >  2:4191/1‑239 (MQ=255)
                              gtgAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAACTAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGGGATTGTTGGGGAAGCCGCCCCGGGGGTTCGTACCgg                                                                                                                                                                                                >  2:341343/1‑239 (MQ=255)
                                gAATGACTTTTTGGATCCGGGTCGCGGTGAGAAAATCTTGTTGACGGCAGCCGCCCCGTTTCGCAGGCTGCCAGATGACAACAATCTCTGGTGTGTCCGAAATAAAGAGGGCGGCATGGTGCCCTTTTCTGCTTTCTCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTa                                                                                                                                                                                              <  1:544297/239‑1 (MQ=255)
                                gAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTTCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTa                                                                                                                                                                                              <  2:876759/239‑1 (MQ=255)
                                      tttgtggtgtCGCGGTTGGGTGAAGGAAGTTTTTGGTAAGGCGGATGCGCCGTATCGCATGCTGCCAGATGACACCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTGCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCg                                                                                                                                                                                         <  1:535562/229‑1 (MQ=255)
                                       ttttgggttCGGGGTCGCGTGAATAAAGTTTATGTGAAGGCGGCTGCGCGGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGAt                                                                                                                                                                                       <  1:63872/231‑1 (MQ=255)
                                       tttATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAACCAGG‑CTCCCCGCGTCTGGAACGCTATACCGGTTATTCTGCGG‑TTGGGATTGTTGGGGAAGCGGCACCGGGGGCCAGGACAGGGACggg                                                                                                                                                                                          >  1:211370/1‑235 (MQ=255)
                                          gttgtgtcggtggtcgTTAAAAAAGTTTCTGTGCGGGCAGTGGCACCGTATCCCATGCTGCCAGATGACATCAACCTCTGGTATGTCCGAAATAAAGAGTGCGGCATGGGGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACATG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGTTCAGTACCGGTACGGCGATGGa                                                                                                                                                                                    <  1:420787/225‑1 (MQ=255)
                                               tCGCGGTCCCGTGAAGAAAGGCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGTAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGAtatt                                                                                                                                                                                <  2:27174/238‑1 (MQ=255)
                                                cgcgGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGGGATTGTGGGGGAAGCCGACCCGGGGGTCCGTACCGGCACGGCGATGGAttttat                                                                                                                                                                              >  1:125655/1‑239 (MQ=255)
                                                cgcgGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCCGTACCGGTACGGCGATGGAtattat                                                                                                                                                                              >  2:776575/1‑239 (MQ=255)
                                                cgcgGTCCGGTGAAAAAAGCTTTTGTGCGGGAAGCTGCGCCGTATCGCAGGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCCCGGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGAtattat                                                                                                                                                                              <  1:357758/239‑1 (MQ=255)
                                                   ggTCGCGTGAAGAAAGTCTATGTGGAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGa                                                                                                                                                                           <  2:124073/239‑1 (MQ=255)
                                                       gcgTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCCGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATc                                                                                                                                                                        <  2:460676/238‑1 (MQ=255)
                                                                       tgtgCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCGGGAACCCTATAACGGTTCTTCTGCGG‑TTGGGATTGTTGGGGAAGCGGCACCGGGGGTCGGAACCGGTACGGGGAGGGAGATTTTGGAATAGGTAGTGAAGAAGGTGc                                                                                                                                                       >  1:777889/1‑239 (MQ=255)
                                                                           catggcgCTTCGCCGTTTCGCAGTGTGACAAATGACAACAATCTCTGGTATGTCCTAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGTAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGTTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCaaa                                                                                                                                                   <  1:597565/233‑1 (MQ=255)
                                                                           caggcatgggcgcGGAATCGAACGCTGCCAGATGACCTCAATCTGTGGTAGGGCCGAAATAAAGATGGCGGCAGGGGGCCCATCTTTGCTTTCGGGACCGCCCGCGGGGAAACAGG‑CTCGCCGCGTCTGGAACGGTATAACGGTTATTCTGCGG‑TTGAGAGTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCaaa                                                                                                                                                   <  1:822913/230‑1 (MQ=255)
                                                                              gcagCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTTCGG‑TTGGGATTGTTGGGGAAGCCGCACCGGGGGAGCGGACCGGTCAGGCGATGGGTATTTGTGAATCGTTTGAGAAGCAGCCGCCAAACg                                                                                                                                                 >  1:442949/1‑238 (MQ=255)
                                                                                ctctggcgcgtagtgcttgctgcCAGAGGACATCAATCTCTGGTAGGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGCCCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCt                                                                                                                                              <  1:170261/222‑1 (MQ=255)
                                                                                     cccctgttCGCCTGCTGCCAGATGACCTCAATCTCTGTTACGGCCGAAATAATGACGGCGGCAGGGCGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGt                                                                                                                                         <  1:780238/232‑1 (MQ=255)
                                                                                                aTGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGa                                                                                                                               <  2:736728/238‑1 (MQ=255)
                                                                                                  gctgcCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACgg                                                                                                                            <  2:777889/239‑1 (MQ=255)
                                                                                                    tgcCAGATGACATAAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCg                                                                                                                          <  2:207000/239‑1 (MQ=255)
                                                                                                                        tGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGGGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGAt                                                                                                                                                                                       <  1:186601/158‑1 (MQ=255)
                                                                                                                        tGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGGGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGAt                                                                                                                                                                                       >  2:186601/1‑158 (MQ=255)
                                                                                                                        tGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGAAATCGTTAGTGACGCAGGTGCCAACGGGCTTTGGTGTGGGTGGGCCGGCAGTGTGGTATCAGGAGGGGCt                                                                                                      >  1:763537/1‑239 (MQ=255)
                                                                                                                         ggTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGCTTGGAGGGGACGGCGGTTTCGTATCAGGACCGGCtt                                                                                                     >  2:385780/1‑239 (MQ=255)
                                                                                                                             tGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGCTGGGTATTATGGAATCGTTAGTGAAGCAGCTGCCAAAAGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCCGGAGCGGCTTTCCg                                                                                                 >  2:293966/1‑239 (MQ=255)
                                                                                                                              gTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGCTATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCg                                                                                                 >  2:35595/1‑238 (MQ=255)
                                                                                                                                   aaaTAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGAtatt                                                                                                                                                                                >  1:546997/1‑154 (MQ=255)
                                                                                                                                   aaaTAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGAtatt                                                                                                                                                                                <  2:546997/154‑1 (MQ=255)
                                                                                                                                    aaTAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGAGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCACACGGCTTTGGACTGGAGTGGACGGCGATTTGGTATAGGGAGGGGCTTTACGgggcgcc                                                                                          >  1:481076/1‑238 (MQ=255)
                                                                                                                                              gcggcATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTTTGGAGTGGACGGCGATGGCGGAACAGGAGAGGCTTTCCGGCGCGCAGGCCCGGgcg                                                                                >  1:336931/1‑239 (MQ=255)
                                                                                                                                                       tGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGTTCCGGCGCTGtttg                                                                         >  2:758724/1‑234 (MQ=255)
                                                                                                                                                                    tgtCGCGACCCTACGTTGGGAAACCGGCCGCGCGGCG‑GTGGAACGCTTGCACGGTTATTCTGTGG‑TGGAGATGTTGGGGGAAGCCGCACCGGGGGTCAGTCCCGGTACGGCGATGGATTTTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGTTCCGTCGCTGTACGCCATTTCCTTGCtgg                                                          <  1:529589/237‑1 (MQ=255)
                                                                                                                                                                      tCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCCTTTCCTTGCtgggg                                                        >  2:119343/1‑237 (MQ=255)
                                                                                                                                                                        gcgACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTAACGGTCCGGCGATGGATATTTTGGAATCGTTAGTGAAGCAGCTTCCAAACGGCTTTGGGCTTGAGGGGAGGGCGAGGTCGAATCAGGAGCGGCTTTCCGGCGGGGGGGCACCGGGGCGGGCCGCCAATTCCtttgggggggg                                                      >  1:618243/1‑230 (MQ=255)
                                                                                                                                                                          gACCTCACGCTGGGAAGCGGG‑CTCGCCGCGTCGGGAACGCTATATCGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCTTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTa                                                     <  2:723449/238‑1 (MQ=255)
                                                                                                                                                                             cTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTTCGCCATTTCCTTGCTGGGGGGAtt                                                   >  2:134035/1‑237 (MQ=255)
                                                                                                                                                                             cTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATAAGGAGCGGCTTTCCGGCAACAAGGCTCCGGCGCTGTAGGCCATTTCCTTGCTGGTGGTCtt                                                   >  2:589535/1‑237 (MQ=255)
                                                                                                                                                                                     ggggaacgcggcgcgcggctttgggAACGCTATAACGGTTATTAGGCGG‑TGGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTc                                           <  1:756062/215‑1 (MQ=255)
                                                                                                                                                                                      ggAAACAGG‑CTCGCTCGGTGTGGAACCCTTTAACGGTTTTTTGGCGG‑TTGAGATTGTTGGGGAAGCCGCCCCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTgg                                        <  1:854153/239‑1 (MQ=255)
                                                                                                                                                                                       gAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTgg                                        <  2:200040/238‑1 (MQ=255)
                                                                                                                                                                                             gggtctgccgcGTCTGGAACGCTATAACGGTTATTCTGGGG‑TTGGGATTGTTGGGGAAGCCGGCCCGGGGGTCAGTTCCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGGCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGt                                  <  1:119343/233‑1 (MQ=255)
                                                                                                                                                                                             ggggtccgtgtttgtGGAACGCTATGAGGGTTATTTGGGGG‑TGGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTgtgt                                  <  1:776575/225‑1 (MQ=255)
                                                                                                                                                                                              g‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGGATTCCTGGGTCTGGCTGCggtt                                >  2:756062/1‑236 (MQ=255)
                                                                                                                                                                                                cgggccgcGTCGGGAACTGTATAAGTGTT‑TTCTGGGGTTTGGGATTGTGGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAAGCGGTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGGGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGt                                  <  1:568912/233‑1 (MQ=255)
                                                                                                                                                                                                 tCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTa                              >  2:577139/1‑239 (MQ=255)
                                                                                                                                                                                                            ggAACGCTATAACGGTTAGTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGt                   <  2:211370/239‑1 (MQ=255)
                                                                                                                                                                                                                  cTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGGTGTATGAAAGCTGGTCGGTg              >  2:885136/1‑238 (MQ=255)
                                                                                                                                                                                                                      caCGGTTATTGTGGGG‑GTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGtt         <  1:54364/238‑1 (MQ=255)
                                                                                                                                                                                                                      aaCGGTTATTCTGCGG‑TTGGGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGtt         <  2:618243/239‑1 (MQ=255)
                                                                                                                                                                                                                        cGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGGGTGGACGGCGAGGTCGTATCAGGAGCGGCTTTCCGGCGCCCGGGGCCCGACGCGTTACGCACTTTCCTTGCTGGTGGTCTTCCTGTGCTTGGCGGCGTGTTATCAAACCTGGTCGGTGCCGt          >  1:619238/1‑236 (MQ=255)
                                                                                                                                                                                                                             aTTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGGTGGCGGATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGGATTCCGGTGTCTGGCTGTGTTGTATGAAAGCTGGTCGGGGGCGTTCTCggga  >  1:372315/1‑237 (MQ=255)
                                                                                                                                                                                                                             aTTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTa  >  2:825896/1‑239 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                  
ACACTGCAAACCGCCTGGGGTTCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGG‑CTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGG‑TTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTA  >  NC_000913/2589831‑2590289

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: ‑