Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 721,604 G→A R44C (CGT→TGT)  kdpE ← DNA‑binding transcriptional activator KdpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913721,6040GA97.4% 229.9 / ‑4.2 78R44C (CGT→TGT) kdpEDNA‑binding transcriptional activator KdpE
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (35/39);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

TGGTGGCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCGTGTACAT  >  NC_000913/721367‑721845
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
gggtggcaggggtgcggcgTAATGCGCCGCGAAGACGGGCCTGCAATTCGCCAATGCCGAACGGCTTACGCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGTCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTcc                                                                                                                                                                                                                                       <  2:624436/238‑1 (MQ=255)
  gtggCAGAGTGGCGGCGTAATGCGACGCGCAGACTGGCCGGCACTTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAg                                                                                                                                                                                                                                     <  2:518571/251‑1 (MQ=255)
    ggCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCCATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTTTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATAACCGGCACCGCGCCCCACTGGCGCCGGTCGCGGATAAACACAATCCCATCCCCATCCGGCAGGCCGAGATCGGAAATAATCAAATCTGGCTTACCGGTTGCCGATTCCatc                                                                                                                                                                                                                                    >  2:365895/1‑248 (MQ=255)
                         gACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCgg                                                                                                                                                                                                              >  1:156682/1‑251 (MQ=255)
                                aGACGGTCCTGCTATACTCCAATGCCAAACGGCTTACGCAGATAATCATCCGCGCCTGCAACCAGCGCGGCGATTTTGTCGCTCTCGTCGCGGCGTGTGGAAAGCCCAATCACCGGCACCGCGCTGCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGTGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCaaa                                                                                                                                                                                                       <  2:320242/251‑1 (MQ=255)
                                  aCGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGa                                                                                                                                                                                                     <  1:365895/251‑1 (MQ=255)
                                    gggCCTGCACTTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGc                                                                                                                                                                                                                                >  2:91619/1‑222 (MQ=255)
                                    gggCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGc                                                                                                                                                                                                                                <  1:91619/222‑1 (MQ=255)
                                    ggaccTGCACTTCGCAAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGAATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGCAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATAGGGCAGGCCGAGATCGAGACTAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCATTGCAGCGTTTCGGACTCCAAGAcc                                                                                                                                                                                                   >  2:339489/4‑250 (MQ=255)
                                          tcaaatCGCCAATGCCAAAAGGCTTACTCAGATAAGCACCCGCTCCGGCATCCAGCGCGGCGATTTTGGCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAAGCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTTAAAGACGCGCAt                                                                                                                                                                                              <  2:116987/245‑1 (MQ=255)
                                           cAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGACATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATcc                                                                                                                                                                                            >  1:258815/1‑251 (MQ=255)
                                               tCGCCAATGCCAAACGGCTTACTCAGATAACCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGt                                                                                                                                                                                         <  2:466358/250‑1 (MQ=255)
                                                cGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTTGGCCTCAAAGACGCGCATCCCGTCg                                                                                                                                                                                       <  1:207860/251‑1 (MQ=255)
                                                cGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCCCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTc                                                                                                                                                                                        <  2:208178/250‑1 (MQ=255)
                                                     aTGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGc                                                                                                                                                                                                                                >  1:318538/1‑205 (MQ=255)
                                                     aTGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTc                                                                                                                                                                                  >  2:407983/1‑251 (MQ=255)
                                                     aTGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGCAGCCCTc                                                                                                                                                                                  >  2:14760/1‑251 (MQ=255)
                                                        ccAAACGGCTTACCCAGATAATCATCCGCGCCTGCATCCAGCGCGGCGATTTTGTCGCTCGCGTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATTAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCTGTGCAGCGTGTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAg                                                                                                                                                                               <  2:36219/251‑1 (MQ=255)
                                                         cacaCGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGATCGCGGATAAACTCAATCCCATCCCCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGATGACGCTTCCAGCACGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGcccccc                                                                                                                                                                                 >  2:638806/4‑245 (MQ=255)
                                                            aCGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGcc                                                                                                                                                                           <  1:577575/251‑1 (MQ=255)
                                                                         aTAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCACTCTTCGCTGCGTGCGGCAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCAGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTACAGCGTTTCGGCCCCACAGACGCGACTCCCGCCGCCCTCCAGCGccagtcgcaaaaaac                                                                                                                                                              >  2:657494/1‑238 (MQ=255)
                                                                         aTAATCAGCCGCTCCGCCATCCAGCGCGGCGAGTTTTGCGCTCTCTCCGAGGCATGCGGAAAGCACAAGCACCGGCACAGCGCTCCACTGGCGCAGGTTGCGGATAAACTTAATCCCACCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTTCCGCTGCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGACGACCTGCAGCGCCGTGCGCAGAAAg                                                                                                                                                               <  2:670170/250‑1 (MQ=255)
                                                                               atcCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCCGAAAGcgacga                                                                                                                                                         >  2:411984/1‑250 (MQ=255)
                                                                                  cGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGACAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGCGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCCAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAgcacgcgcagaca                                                                                                                                                       >  2:173172/1‑237 (MQ=255)
                                                                                   gCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGc                                                                                                                                                    >  1:209987/1‑251 (MQ=255)
                                                                                   gCTCCGACATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGCAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGCTAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCACAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGcgacggctagc                                                                                                                                                    >  2:429992/1‑245 (MQ=255)
                                                                                    cTCCGGCATCCAGCGCGGCGATATTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGCCGCGCATCCCGTCGCCCTCCAgcgcagcgagcagacatcgccgtatagcc                                                                                                                                                   >  2:101711/1‑226 (MQ=255)
                                                                                         gCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTc                                                                                                                                              >  1:157459/1‑251 (MQ=255)
                                                                                          cATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGTGCATAAAGCGACGAATAGCCTGTTCa                                                                                                                                             <  2:93434/251‑1 (MQ=255)
                                                                                               ggcggggcgATTTTGTCGCTCTCTGCGCGGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTc                                                                                                                                        <  2:477238/246‑1 (MQ=255)
                                                                                                gcgcGGCGATTTTGTCGCGCTCTTCGCTTCGTGCGGAAAGCACAATCACCGTCACCGCGCTCCACTGGCGCAGGGCGCTGATAAACTCAATCCCATCACCATCGGGCAGGCCGGGATCGAGAATAATCAAAGCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCGCCAGCGCCGTGCGCGGAAAGTGACGAATAGCCTGTTCATCTTca                                                                                                                                       <  2:30865/251‑1 (MQ=255)
                                                                                                  gcggcgATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGCGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAACGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTAATCTTcaac                                                                                                                                     >  2:143726/1‑251 (MQ=255)
                                                                                                  gcggcgATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTAGCGGATAAACTCAATACCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCATCAAAGACGCGCATCCCGTCGCCCTCCAGCGACGTGCGCAGAAAGCGCCGAATAGCCTGTTCATCTTcaac                                                                                                                                     >  2:443661/1‑251 (MQ=255)
                                                                                                     gcgATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAAt                                                                                                                                  <  1:51313/251‑1 (MQ=255)
                                                                                                     gcgATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAAt                                                                                                                                  <  1:630983/251‑1 (MQ=255)
                                                                                                     gcgATTTTGTCGCTCTCTTCGCTGCGGGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAAt                                                                                                                                  <  1:429397/251‑1 (MQ=255)
                                                                                                       gATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATc                                                                                                                                 >  1:19214/1‑250 (MQ=255)
                                                                                                        aTTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATc                                                                                                                                 >  1:201405/1‑249 (MQ=255)
                                                                                                           tttGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGa                                                                                                                              <  1:411984/250‑1 (MQ=255)
                                                                                                            ttgtagccctctTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCACGCCGCGTTGCAGCGTTTCGGCCTCAAAGTCGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAAc                                                                                                                            >  2:662734/9‑251 (MQ=255)
                                                                                                            ttGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGaa                                                                                                                             <  1:662848/250‑1 (MQ=255)
                                                                                                                    ccctTCGCGGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACGGGCGCAGGTCGCGGATAAACTCAATCCCATCACCAGCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTca                                                                                                                    <  2:213976/249‑1 (MQ=255)
                                                                                                                        tCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCCCTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAAGCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACat                                                                                                                 <  2:615451/250‑1 (MQ=255)
                                                                                                                        tCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTcaca                                                                                                                  <  1:173172/249‑1 (MQ=255)
                                                                                                                         cGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACata                                                                                                                <  1:14760/250‑1 (MQ=255)
                                                                                                                                 gcgGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGa                                                                                                       <  1:194822/251‑1 (MQ=255)
                                                                                                                                     aaaGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCACATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAAtt                                                                                                   >  1:278327/1‑251 (MQ=255)
                                                                                                                                     aaaGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAAtt                                                                                                   <  1:501410/251‑1 (MQ=255)
                                                                                                                                          acaATCACCGGCACAGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAAttcttc                                                                                               >  1:509259/1‑250 (MQ=255)
                                                                                                                                               cACCGGCACCGCTCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCTCCATCGGGCAGGTCGAGATCGAGAATATTCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGGGCGCAGAAAGCGACGAATAGCCTTTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCGTCAAGtt                                                                                          <  2:256646/250‑1 (MQ=255)
                                                                                                                                                aCCGGCACCGCGCTCCACTGGCGCAGGTCGCGTATAAACTCAACCCCAACAACATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCATCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCa                                                                                        <  2:135620/251‑1 (MQ=255)
                                                                                                                                                    gCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCCCCATCGGGCAGGCCGAGATCGAGAATAATCAAAGCGGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAgggg                                                                                    <  2:498249/251‑1 (MQ=255)
                                                                                                                                                        cgctctccacgtgcgcAGTTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACTTTTGTCACATAGCCTCATGAAATTCTTCAAGTTCTGGGGCAGt                                                                                <  2:415719/239‑1 (MQ=255)
                                                                                                                                                          cgcTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGttt                                                                              >  1:58797/1‑251 (MQ=255)
                                                                                                                                                           gcTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCt                                                                                                                                                  >  1:20020/1‑182 (MQ=255)
                                                                                                                                                           gcTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCt                                                                                                                                                  <  2:20020/182‑1 (MQ=255)
                                                                                                                                                           gcTCCACGGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCGCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTg                                                                             <  2:437733/251‑1 (MQ=255)
                                                                                                                                                                  tGGCGCAGGTCGCGGATCAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGAGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACCATCCGaaag                                                                                                                           >  2:405981/1‑196 (MQ=255)
                                                                                                                                                                  tGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACg                                                                                                                           <  1:405981/198‑1 (MQ=255)
                                                                                                                                                                                     aCTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCCTCTTCAACAATCAGAACGTTTGTCACATCTCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGACGTGTAACACGAAAACAgg                                                   >  2:3712/1‑251 (MQ=255)
                                                                                                                                                                                       ccAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGca                                                 <  1:339489/250‑1 (MQ=255)
                                                                                                                                                                                       ccAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGCCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGca                                                 <  1:398735/250‑1 (MQ=255)
                                                                                                                                                                                          aTCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGcacc                                               >  1:239247/1‑250 (MQ=255)
                                                                                                                                                                                          aTCCCAGCACCCTCGGGCCGGCCGATATCGAGAATAATCAACTCTGGCTTACAGCTGGCCGCTTCCAGCACGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTGGCTGGGGAAGTGTAACACGAAAACAGGctcca                                             <  2:514822/251‑1 (MQ=255)
                                                                                                                                                                                             ccATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCt                                          >  1:346400/1‑251 (MQ=255)
                                                                                                                                                                                                tCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAATAAGGCAACACCTTcc                                       >  2:204446/1‑251 (MQ=255)
                                                                                                                                                                                                             ggCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAAcg                          <  1:527111/251‑1 (MQ=255)
                                                                                                                                                                                                              gCCGAGATCGAGAATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCAGGTCGGTTGAAcgg                         >  2:513496/1‑250 (MQ=255)
                                                                                                                                                                                                              accGAGATCGAGAATCATCAAATCTGGCTTACAGGTTGCCGCTTCCAGAAAGCCGCGTTGCAGCGTTTCGGCATCACAGACGCGCATCCCGTCGCCCTCCCGCGCCGAGCGCAGAAAGCGACGAATAGACAGATCATCTTCAACAATCAGACCGTTTGTCACATATCCTCATGAAATTGTTCAAGGTCAGGTGCCGCTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTccagtagattggacgg                         >  2:513503/2‑237 (MQ=255)
                                                                                                                                                                                                                   gaTCGAGAATAACTAAATCTGGCTTACAGATTGCCGCTTCCAGAAAGACGCGTTAAAGCGTTTAGGCATCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAACGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTGACATATCCTCATGACATTCTTCACATTCACGTACAGATTTCTGGATAAGGGTAACACCAAAACCACACCCACCTTCAGGTCGGTTGAACGCgggat                    >  2:235736/1‑248 (MQ=255)
                                                                                                                                                                                                                       gagaATAATCTACTCTGGCTTAAAGGTTGCCGCGTCCAGCAAGCCGCGTTGCTGCATTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCTTGCGTAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCATATATCCTCATGAAATTCTTCAAGTGCAgg                                                                                      <  2:213003/182‑1 (MQ=255)
                                                                                                                                                                                                                       gagaATAATCAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCCGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAAAGTTTGTCGCATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGt                >  1:52941/1‑251 (MQ=255)
                                                                                                                                                                                                                       gagaATAATCAAATCTGGCTTACAGGATGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGGAAGCGCCGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAgg                                                                                      >  1:213003/1‑182 (MQ=255)
                                                                                                                                                                                                                              atCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCAGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGGGCGCAGAAAACGACGAATAGCCTGTTCATCTTCAACAATAAGAACGTGTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGTAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGACCCCg         >  2:369407/1‑251 (MQ=255)
                                                                                                                                                                                                                              atCAAATCTGGAATACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTAGCAAACATCCTCATGAAATTCTTACAGTTCAGGGACAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCCACAACTGACGGTCGGTTGAACGCGGTAATAGTGCCCCCg         >  2:411840/1‑251 (MQ=255)
                                                                                                                                                                                                                                cAAATCTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGTCGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAACTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCgt        <  1:638806/250‑1 (MQ=255)
                                                                                                                                                                                                                                 acaTCTGGCTTACCCGGTTCCGCTTCCAGCAAGCCGCGGTGCCGCGTATCTGCCTCCAAAACGCGACTCCCGTCAACCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCCACAATCAGAACGTTTGTCACATATCCTCATGCAATTCTTAAAGTTCAGGGGAAGGTTGCTGGGGAAGGGCAAAACGGAAAAAGGCACCACCTTCCGGCCGTTTGAACGCGTTAATAGTGCCCCCgggt      >  2:107133/3‑248 (MQ=255)
                                                                                                                                                                                                                                     cTGGCTTACAGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGACATTCTTCACGTTCTGGGGCAGTTTGCTGGGGCAGTGTAACACGAAACCAGGCACCACCTTCCGGTCGGTTGCACGCGGAACTAGTACACCCGTGTACAt  >  2:220263/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
TGGTGGCAGAGTGGCGGCGTAATGCGACGCGCAGACGGGCCTGCAATTCGCCAATGCCAAACGGCTTACTCAGATAATCATCCGCTCCGGCATCCAGCGCGGCGATTTTGTCGCTCTCTTCGCTGCGTGCGGAAAGCACAATCACCGGCACCGCGCTCCACTGGCGCAGGTCGCGGATAAACTCAATCCCATCACCATCGGGCAGGCCGAGATCGAGAATAATCAAATCTGGCTTACGGGTTGCCGCTTCCAGCAAGCCGCGTTGCAGCGTTTCGGCCTCAAAGACGCGCATCCCGTCGCCCTCCAGCGCCGTGCGCAGAAAGCGACGAATAGCCTGTTCATCTTCAACAATCAGAACGTTTGTCACATATCCTCATGAAATTCTTCAAGTTCAGGGGCAGTTTGCTGGGGAAGTGTAACACGAAAACAGGCACCACCTTCCGGTCGGTTGAACGCGGTAATAGTGCCCCCGTGTACAT  >  NC_000913/721367‑721845

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