Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,590,065 A→G 18.2% 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,0650AG18.2% 143.8 / 3.9 67E824G (GAG→GGG) acrDaminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base A (23/29);  new base G (11/4);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 7.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

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

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