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

ACACTGCAAACCGCCTGGGGTTCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTA  >  NC_000913/2589831‑2590289
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                  
acacTGCAAACCGCCTGGGGTTCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGGGAttgt                                                                                                                                                                                                                              >  2:643167/1‑239 (MQ=255)
  acTGCAAACCGCCTGGGGTTCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGGGAttgttg                                                                                                                                                                                                                            >  2:54364/1‑239 (MQ=255)
  acTGCAAACCGCCTGGGGTTCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGAttgttg                                                                                                                                                                                                                            >  2:697399/1‑239 (MQ=255)
    tGCAAACCGCCTGGGGTTCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTCTTCTGCGGTTGGGATTGTTggg                                                                                                                                                                                                                          >  2:764676/1‑239 (MQ=255)
    tGCAAACCGCCTGGGGTTCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGGCTGGAACGCTATAACGGTTATTCTGCGGTTGGGTTTGTTggg                                                                                                                                                                                                                          >  2:435318/1‑239 (MQ=255)
            gCCTGGGGTTCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGGGCCCTTCTCTGCTTTCGCGACCTCCCGCTGGGGAACAGGCTGCCCGCGTCTGGAAGCCCATAACGGTTTTTTTGGGGTTGGGGTTTTTGGGGAAGccc                                                                                                                                                                                                                   >  1:736728/1‑237 (MQ=255)
             ccTGGGGTTCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGc                                                                                                                                                                                                                  >  2:790878/1‑238 (MQ=255)
                ggggTTTGGGTTTTGTGAATGACTTTTTGGATCGCGGTCGCGTGAAAAAAGTCTATGTGCGGGCAGCTGCGCCGTACCGCATGCTGCCGGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCg                                                                                                                                                                                                              <  1:435318/239‑1 (MQ=255)
                ggggTTCGAGCTATGTGAATGGATTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTAGCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCg                                                                                                                                                                                                              <  2:443053/239‑1 (MQ=255)
                 gggTTCGGGTTTTTTGAATGTCTTTATGGATGGCGGTCGCGTGAAAAAAGTCTTTGTGCAGGCACCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCGTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCg                                                                                                                                                                                                              <  1:591438/238‑1 (MQ=255)
                 gggTTCGAGTTTTGTGAATGACTTTTTGGATCGGGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCg                                                                                                                                                                                                              <  1:581401/238‑1 (MQ=255)
                     ttcgtgtgtgtgAAGGACTTTAGGGATCGCGGGGGCGTGAAGAAAGTCTGTGTGCAGGCAGCTGGGCCGTATCGCATGCTGCCAGATGACATCAATCCCTGGTATGTCCGAAATAAAGAGGGCGGCAGGGTGCCCTTCTCTGGTTTCGCGACCTCCCGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGt                                                                                                                                                                                                         <  1:774073/231‑1 (MQ=255)
                     tCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGt                                                                                                                                                                                                         <  2:579939/239‑1 (MQ=255)
                     tCGAGCTATGTGAATGACTTTAGGGATCGGGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACCTCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGt                                                                                                                                                                                                         <  2:361013/239‑1 (MQ=255)
                         gCTATGTGAATGGCTTTAGGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGAATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAg                                                                                                                                                                                                      <  2:533260/238‑1 (MQ=255)
                            aTGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTa                                                                                                                                                                                                    >  2:581401/1‑237 (MQ=255)
                             tgtgAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCg                                                                                                                                                                                                 >  2:4191/1‑239 (MQ=255)
                              gtgAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAACTAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGGGATTGTTGGGGAAGCCGCCCCGGGGGTTCGTACCgg                                                                                                                                                                                                >  2:341343/1‑239 (MQ=255)
                                gAATGACTTTTTGGATCCGGGTCGCGGTGAGAAAATCTTGTTGACGGCAGCCGCCCCGTTTCGCAGGCTGCCAGATGACAACAATCTCTGGTGTGTCCGAAATAAAGAGGGCGGCATGGTGCCCTTTTCTGCTTTCTCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTa                                                                                                                                                                                              <  1:544297/239‑1 (MQ=255)
                                gAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTTCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTa                                                                                                                                                                                              <  2:876759/239‑1 (MQ=255)
                                      tttgtggtgtCGCGGTTGGGTGAAGGAAGTTTTTGGTAAGGCGGATGCGCCGTATCGCATGCTGCCAGATGACACCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTGCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCg                                                                                                                                                                                         <  1:535562/229‑1 (MQ=255)
                                       ttttgggttCGGGGTCGCGTGAATAAAGTTTATGTGAAGGCGGCTGCGCGGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGAt                                                                                                                                                                                       <  1:63872/231‑1 (MQ=255)
                                       tttATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAACCAGGCTCCCCGCGTCTGGAACGCTATACCGGTTATTCTGCGGTTGGGATTGTTGGGGAAGCGGCACCGGGGGCCAGGACAGGGACggg                                                                                                                                                                                          >  1:211370/1‑235 (MQ=255)
                                          gttgtgtcggtggtcgTTAAAAAAGTTTCTGTGCGGGCAGTGGCACCGTATCCCATGCTGCCAGATGACATCAACCTCTGGTATGTCCGAAATAAAGAGTGCGGCATGGGGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACATGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGTTCAGTACCGGTACGGCGATGGa                                                                                                                                                                                    <  1:420787/225‑1 (MQ=255)
                                               tCGCGGTCCCGTGAAGAAAGGCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGTAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGAtatt                                                                                                                                                                                <  2:27174/238‑1 (MQ=255)
                                                cgcgGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGGGATTGTGGGGGAAGCCGACCCGGGGGTCCGTACCGGCACGGCGATGGAttttat                                                                                                                                                                              >  1:125655/1‑239 (MQ=255)
                                                cgcgGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCCGTACCGGTACGGCGATGGAtattat                                                                                                                                                                              >  2:776575/1‑239 (MQ=255)
                                                cgcgGTCCGGTGAAAAAAGCTTTTGTGCGGGAAGCTGCGCCGTATCGCAGGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCCCGGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGAtattat                                                                                                                                                                              <  1:357758/239‑1 (MQ=255)
                                                   ggTCGCGTGAAGAAAGTCTATGTGGAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGa                                                                                                                                                                           <  2:124073/239‑1 (MQ=255)
                                                       gcgTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCCGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATc                                                                                                                                                                        <  2:460676/238‑1 (MQ=255)
                                                                       tgtgCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCGGGAACCCTATAACGGTTCTTCTGCGGTTGGGATTGTTGGGGAAGCGGCACCGGGGGTCGGAACCGGTACGGGGAGGGAGATTTTGGAATAGGTAGTGAAGAAGGTGc                                                                                                                                                       >  1:777889/1‑239 (MQ=255)
                                                                           catggcgCTTCGCCGTTTCGCAGTGTGACAAATGACAACAATCTCTGGTATGTCCTAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGTAAACAGGCTCGCCGCGTCTGGAACGCTATAACGTTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCaaa                                                                                                                                                   <  1:597565/233‑1 (MQ=255)
                                                                           caggcatgggcgcGGAATCGAACGCTGCCAGATGACCTCAATCTGTGGTAGGGCCGAAATAAAGATGGCGGCAGGGGGCCCATCTTTGCTTTCGGGACCGCCCGCGGGGAAACAGGCTCGCCGCGTCTGGAACGGTATAACGGTTATTCTGCGGTTGAGAGTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCaaa                                                                                                                                                   <  1:822913/230‑1 (MQ=255)
                                                                              gcagCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTTCGGTTGGGATTGTTGGGGAAGCCGCACCGGGGGAGCGGACCGGTCAGGCGATGGGTATTTGTGAATCGTTTGAGAAGCAGCCGCCAAACg                                                                                                                                                 >  1:442949/1‑238 (MQ=255)
                                                                                ctctggcgcgtagtgcttgctgcCAGAGGACATCAATCTCTGGTAGGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGCCCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCt                                                                                                                                              <  1:170261/222‑1 (MQ=255)
                                                                                     cccctgttCGCCTGCTGCCAGATGACCTCAATCTCTGTTACGGCCGAAATAATGACGGCGGCAGGGCGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGt                                                                                                                                         <  1:780238/232‑1 (MQ=255)
                                                                                                aTGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGa                                                                                                                               <  2:736728/238‑1 (MQ=255)
                                                                                                  gctgcCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACgg                                                                                                                            <  2:777889/239‑1 (MQ=255)
                                                                                                    tgcCAGATGACATAAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCg                                                                                                                          <  2:207000/239‑1 (MQ=255)
                                                                                                                        tGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGGGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGAt                                                                                                                                                                                       <  1:186601/158‑1 (MQ=255)
                                                                                                                        tGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGGGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGAt                                                                                                                                                                                       >  2:186601/1‑158 (MQ=255)
                                                                                                                        tGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGAAATCGTTAGTGACGCAGGTGCCAACGGGCTTTGGTGTGGGTGGGCCGGCAGTGTGGTATCAGGAGGGGCt                                                                                                      >  1:763537/1‑239 (MQ=255)
                                                                                                                         ggTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGCTTGGAGGGGACGGCGGTTTCGTATCAGGACCGGCtt                                                                                                     >  2:385780/1‑239 (MQ=255)
                                                                                                                             tGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGCTGGGTATTATGGAATCGTTAGTGAAGCAGCTGCCAAAAGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCCGGAGCGGCTTTCCg                                                                                                 >  2:293966/1‑239 (MQ=255)
                                                                                                                              gTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGCTATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCg                                                                                                 >  2:35595/1‑238 (MQ=255)
                                                                                                                                   aaaTAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGAtatt                                                                                                                                                                                >  1:546997/1‑154 (MQ=255)
                                                                                                                                   aaaTAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGAtatt                                                                                                                                                                                <  2:546997/154‑1 (MQ=255)
                                                                                                                                    aaTAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGAGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCACACGGCTTTGGACTGGAGTGGACGGCGATTTGGTATAGGGAGGGGCTTTACGgggcgcc                                                                                          >  1:481076/1‑238 (MQ=255)
                                                                                                                                              gcggcATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTTTGGAGTGGACGGCGATGGCGGAACAGGAGAGGCTTTCCGGCGCGCAGGCCCGGgcg                                                                                >  1:336931/1‑239 (MQ=255)
                                                                                                                                                       tGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGTTCCGGCGCTGtttg                                                                         >  2:758724/1‑234 (MQ=255)
                                                                                                                                                                    tgtCGCGACCCTACGTTGGGAAACCGGCCGCGCGGCGGTGGAACGCTTGCACGGTTATTCTGTGGTGGAGATGTTGGGGGAAGCCGCACCGGGGGTCAGTCCCGGTACGGCGATGGATTTTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGTTCCGTCGCTGTACGCCATTTCCTTGCtgg                                                          <  1:529589/237‑1 (MQ=255)
                                                                                                                                                                      tCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCCTTTCCTTGCtgggg                                                        >  2:119343/1‑237 (MQ=255)
                                                                                                                                                                        gcgACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTAACGGTCCGGCGATGGATATTTTGGAATCGTTAGTGAAGCAGCTTCCAAACGGCTTTGGGCTTGAGGGGAGGGCGAGGTCGAATCAGGAGCGGCTTTCCGGCGGGGGGGCACCGGGGCGGGCCGCCAATTCCtttgggggggg                                                      >  1:618243/1‑230 (MQ=255)
                                                                                                                                                                          gACCTCACGCTGGGAAGCGGGCTCGCCGCGTCGGGAACGCTATATCGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCTTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTa                                                     <  2:723449/238‑1 (MQ=255)
                                                                                                                                                                             cTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTTCGCCATTTCCTTGCTGGGGGGAtt                                                   >  2:134035/1‑237 (MQ=255)
                                                                                                                                                                             cTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATAAGGAGCGGCTTTCCGGCAACAAGGCTCCGGCGCTGTAGGCCATTTCCTTGCTGGTGGTCtt                                                   >  2:589535/1‑237 (MQ=255)
                                                                                                                                                                                     ggggaacgcggcgcgcggctttgggAACGCTATAACGGTTATTAGGCGGTGGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTc                                           <  1:756062/215‑1 (MQ=255)
                                                                                                                                                                                      ggAAACAGGCTCGCTCGGTGTGGAACCCTTTAACGGTTTTTTGGCGGTTGAGATTGTTGGGGAAGCCGCCCCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTgg                                        <  1:854153/239‑1 (MQ=255)
                                                                                                                                                                                       gAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTgg                                        <  2:200040/238‑1 (MQ=255)
                                                                                                                                                                                             gggtctgccgcGTCTGGAACGCTATAACGGTTATTCTGGGGTTGGGATTGTTGGGGAAGCCGGCCCGGGGGTCAGTTCCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGGCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGt                                  <  1:119343/233‑1 (MQ=255)
                                                                                                                                                                                             ggggtccgtgtttgtGGAACGCTATGAGGGTTATTTGGGGGTGGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTgtgt                                  <  1:776575/225‑1 (MQ=255)
                                                                                                                                                                                              gCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGGATTCCTGGGTCTGGCTGCggtt                                >  2:756062/1‑236 (MQ=255)
                                                                                                                                                                                                cgggccgcGTCGGGAACTGTATAAGTGTTTTCTGGGGTTTGGGATTGTGGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAAGCGGTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGGGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGt                                  <  1:568912/233‑1 (MQ=255)
                                                                                                                                                                                                 tCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTa                              >  2:577139/1‑239 (MQ=255)
                                                                                                                                                                                                            ggAACGCTATAACGGTTAGTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGt                   <  2:211370/239‑1 (MQ=255)
                                                                                                                                                                                                                  cTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGGTGTATGAAAGCTGGTCGGTg              >  2:885136/1‑238 (MQ=255)
                                                                                                                                                                                                                      caCGGTTATTGTGGGGGTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGtt         <  1:54364/238‑1 (MQ=255)
                                                                                                                                                                                                                      aaCGGTTATTCTGCGGTTGGGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGtt         <  2:618243/239‑1 (MQ=255)
                                                                                                                                                                                                                        cGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGGGTGGACGGCGAGGTCGTATCAGGAGCGGCTTTCCGGCGCCCGGGGCCCGACGCGTTACGCACTTTCCTTGCTGGTGGTCTTCCTGTGCTTGGCGGCGTGTTATCAAACCTGGTCGGTGCCGt          >  1:619238/1‑236 (MQ=255)
                                                                                                                                                                                                                             aTTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGGTGGCGGATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGGATTCCGGTGTCTGGCTGTGTTGTATGAAAGCTGGTCGGGGGCGTTCTCggga  >  1:372315/1‑237 (MQ=255)
                                                                                                                                                                                                                             aTTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTa  >  2:825896/1‑239 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                  
ACACTGCAAACCGCCTGGGGTTCGAGCTATGTGAATGACTTTATGGATCGCGGTCGCGTGAAGAAAGTCTATGTGCAGGCAGCTGCGCCGTATCGCATGCTGCCAGATGACATCAATCTCTGGTATGTCCGAAATAAAGATGGCGGCATGGTGCCCTTCTCTGCTTTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTATTCTGCGGTTGAGATTGTTGGGGAAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTA  >  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: