Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,948,824 C→T A567T (GCA→ACA)  aspS ← aspartate‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,948,8240CT97.6% 261.3 / NA 83A567T (GCA→ACA) aspSaspartate‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (43/38);  minor base A (1/1);  total (44/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

AGATGACCACTAAGATCGAAACGGGACGCTTATACACTTTATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAA  >  NC_000913/1948584‑1949065
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
aGATGAGCACTAAGAGCGAAACGGGACGATTAGACACTTTATACTGCACGCTTCCTCGTCGGTGACTGCCTTCGAGCAGCCCATTCACATAGTTATCTATGCTCATGGGAGTGCACTCAGTTGCCGACTCGATGCAACGCGAAGGATGTCGTGTATTTGAGTCAGATCAGTTATTCTCAGCCTTCTGCACAACCTGAATGCGCAGCGCAGCCAGTGGAGTCGGGTTAGCAAAGCTCGGTGTTTCAGttat                                                                                                                                                                                                                                           <  2:188851/250‑4 (MQ=255)
aGATGACCACTAAGATCGAAACGGGACGCTTATACACTTTATCCTTCACGCTGCCGCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGtcat                                                                                                                                                                                                                                           <  2:188852/250‑1 (MQ=255)
  aTGACCACTAAGATCGAAACGGGACGCTTATACACTTTATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCGCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCGCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTCCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAg                                                                                                                                                                                                                                        <  1:578570/251‑1 (MQ=255)
    gACCACTCAGATCGAAGAGAGACGCTTATACACTTTATCCTTCACGCTGCCTACTCGATGACTGCATTCGCTCATACCATTCACATAGTTATCAATGCTAATGGGCGTTCAATCAGATGACGCCACGATGCAAAGCGAATGTTTTCGTGTATTTGAGACATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAACGCTCGGTGTTTCAGTCATCAGac                                                                                                                                                                                                                                      >  2:586954/1‑251 (MQ=255)
      ccACTAAGATCGAAACGGGACGCTTATACACTTTATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGacac                                                                                                                                                                                                                                    <  1:220469/251‑1 (MQ=255)
         cTAAGATCGAAACGGGACGCTTATACACTTTATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCt                                                                                                                                                                                                                                 >  1:79592/1‑251 (MQ=255)
         cTAAGATCGAAACGGGACGCTTATACACTTTATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCt                                                                                                                                                                                                                                 >  2:185771/1‑251 (MQ=255)
          tAAGATCGAAACGGGACGCTTATACACTTTATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTg                                                                                                                                                                                                                                >  2:510698/1‑251 (MQ=255)
                cGAAACGGGACGCTTATACACTTTATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGccgccg                                                                                                                                                                                                                          >  1:422134/1‑251 (MQ=255)
                cGAAACGGGACGCTTATACACTTTATCCTTCACGCTGCCTCTTCGGTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATGTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGGCATCAGACACGCTGccgccg                                                                                                                                                                                                                          <  1:587008/251‑1 (MQ=255)
                    aCGGGACGCTTATAATCTTTATCCTTCACGCTGCCTCTCCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCACAGCTCGGTGTTTCACTCATCAGCCACGCTGCATCCGTGGt                                                                                                                                                                                                                      >  2:319205/1‑251 (MQ=255)
                     cgggacgctcagcctctatacccgtCACGCTGCCTCCACGTTGACTGTCGTCGCTCATCCCATTCACACAGTTATCTATGCTCATGGGAGTTCTCTCAGTTGCCGCCTCGATGCCCCGCGAATGATTTCGTGTATTTGATTCATATCATTTATTCTCAGCCTGCGTCACAACCGGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGATTCATCCATCAGACACGATGCCGCCGTGGtt                                                                                                                                                                                                                     <  2:246933/227‑1 (MQ=255)
                      gggACGCTTATACACTTTATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGttt                                                                                                                                                                                                                    >  1:61990/1‑251 (MQ=255)
                       ggACGCTTATACACTTTATCCTTCTCGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATGCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCTTGTATTTGAGTCATATCAGTTATGCGCAGCCGTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTATGAAAGCTCTGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGtttt                                                                                                                                                                                                                   <  2:529682/251‑1 (MQ=255)
                         aCGCTTATACACTTTATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCCCATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGGGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCg                                                                                                                                                                                                                 >  2:307635/1‑251 (MQ=255)
                           gCTTATACACTTTATCCTTCACGCTGCCTCTTCGTTGAATGCCTTCGCTCATACCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGACTGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCCAAGCTCGGGGTTTCAGTCATCAGACACGCTGCCGCCTTGGTTTTCggg                                                                                                                                                                                                               >  2:477702/1‑251 (MQ=255)
                              tataCACTTTATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGTTGCCGCCGTGGTTTTCGGGaa                                                                                                                                                                                                             <  2:312401/250‑1 (MQ=255)
                                   acTTTATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGAt                                                                                                                                                                                                       >  1:78727/1‑251 (MQ=255)
                                     tttATCATTCAAGCTGCCTCTTCGCTGCATGCCTTCGCTCATCCTATTAACCTAGTTATCTATGCTCCTGGGTGTTCACTCAGTTGCCGCCTCGATGCACAGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCCGTGCCGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGCAAGGCGATa                                                                                                                                                                                                      >  2:70090/1‑250 (MQ=255)
                                       tATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACg                                                                                                                                                                                                   >  2:65737/1‑251 (MQ=255)
                                             tCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCAc                                                                                                                                                                                               <  1:350769/249‑1 (MQ=255)
                                             tCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACg                                                                                                                                                                                              <  2:438358/250‑1 (MQ=255)
                                                cGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGat                                                                                                                                                                                           >  1:139141/1‑250 (MQ=255)
                                                   tGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTg                                                                                                                                                                                       >  1:527882/1‑251 (MQ=255)
                                                    gCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATATTTATCTATGCTCATGGGAGTTCCCTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTg                                                                                                                                                                                       <  2:422134/250‑1 (MQ=255)
                                                       tctTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCgg                                                                                                                                                                                  >  2:4849/1‑251 (MQ=255)
                                                         ttCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGAATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTTCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTg                                                                                                                                                                                >  1:565571/1‑251 (MQ=255)
                                                                    ccTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTcagcag                                                                                                                                                                     <  1:513026/251‑1 (MQ=255)
                                                                     cTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCagcagc                                                                                                                                                                    >  1:357257/1‑251 (MQ=255)
                                                                        cGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATACCGTCACGGATATTGTCGGTGCCGGTCAGCAGCATg                                                                                                                                                                 >  1:354465/1‑251 (MQ=255)
                                                                               cccATTCACGTAGTTCTCTATGCTCATCGGAGTTCACTCAGATGCCGCCACGATGAACCGCGAATGTTTTCGTAAATTAGAGTAATATCAGTTATTCTCAGCCTTCTTCCCAACCTGAATGATCAGCTCAGCCAGTGCCGTCGGGCTAGCAAAGCTCGATGTTTCCGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAAAGTCACGGATATTGTCGGGGTCGGTCCGCAGGATGGTCAGa                                                                                                                                                           >  2:239645/1‑250 (MQ=255)
                                                                               cccATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGAc                                                                                                                                                          <  1:555618/251‑1 (MQ=255)
                                                                                  aTTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGt                                                                                                                                                       >  1:217888/1‑251 (MQ=255)
                                                                                  aTTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTAAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGGTTTCGGGAAGGCGATAACGTCACGGATATAGTACGGGCCGGTCAGCAGCCAGGTCAGACGGt                                                                                                                                                       >  2:282277/1‑251 (MQ=255)
                                                                                   ttCACATAGTTATCTATGCTCATGGGAGCTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTATCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGGTTTCGGGAAGGCGATAAAGGCACGGATATTGTCGGTGCCGGTCAGCAGCATGGCCAGCCGGTc                                                                                                                                                      >  2:579880/1‑251 (MQ=255)
                                                                                          aGTTATCTATGCTCATGGTAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTGCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCg                                                                                                                                               <  1:105213/251‑1 (MQ=255)
                                                                                            ttATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATCTCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGAGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTAGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGa                                                                                                                                              >  1:35300/1‑250 (MQ=255)
                                                                                             tATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAAt                                                                                                                                            >  1:168298/1‑251 (MQ=255)
                                                                                               tctctGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGACGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGCATGc                                                                                                                                          >  2:382777/5‑251 (MQ=255)
                                                                                                        aTGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCtgc                                                                                                                                 >  2:69891/1‑251 (MQ=255)
                                                                                                           ggAGTTCCCACAGTTGCCGCCGCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCtgcgtg                                                                                                                              <  2:401376/251‑1 (MQ=255)
                                                                                                           ggAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCa                                                                                                                                                             <  1:424110/220‑1 (MQ=255)
                                                                                                           ggAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCa                                                                                                                                                             >  2:424110/1‑220 (MQ=255)
                                                                                                           ggAGTTCACTCAGTTGCCGCCGCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCGTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCtgcgtg                                                                                                                              <  1:425822/251‑1 (MQ=255)
                                                                                                             aGTTCACTCAGTTGACGCATCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCGGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCCGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCtgcgtgcg                                                                                                                            >  1:504226/1‑251 (MQ=255)
                                                                                                              gTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTgcgg                                                                                                                           <  1:510698/251‑1 (MQ=255)
                                                                                                                 cACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTgcggcg                                                                                                                         <  1:65737/250‑1 (MQ=255)
                                                                                                                  aCTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCAGCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATATCGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGa                                                                                                                       <  1:273216/251‑1 (MQ=255)
                                                                                                                      aGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTAc                                                                                                                   >  1:405408/1‑251 (MQ=255)
                                                                                                                      aGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAAGGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTAc                                                                                                                   <  1:429410/251‑1 (MQ=255)
                                                                                                                         ggccgtcTCGATGCAACGCGAATGATTTCGTGTATTTTGGTCATATCAGTTATTCTCAGCCTTGTGCACAACCTGTAGGCTGAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTGGGTGTTTCAGTCATCAGACACGCTGCCGCCGGTGTTTTCGGGAGGGCGATAAGGTCTCGGATAGTGTCGGTGCCGGGCAGTAGCATGGGCAGACGGACAAGACCGAATGCCTTACGGGCGTGCGGCGGGGGACCGt                                                                                                                <  2:30909/250‑1 (MQ=255)
                                                                                                                                    tGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGAc                                                                                                                                                          <  2:157439/198‑1 (MQ=255)
                                                                                                                                    tGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGAc                                                                                                                                                          >  1:157439/1‑198 (MQ=255)
                                                                                                                                      cAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTc                                                                                                   >  1:218756/1‑251 (MQ=255)
                                                                                                                                          gcgAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCGCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGc                                                                                               <  2:78727/251‑1 (MQ=255)
                                                                                                                                           cgAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCa                                                                                              <  1:62052/251‑1 (MQ=255)
                                                                                                                                            gAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAg                                                                                             <  2:139141/251‑1 (MQ=255)
                                                                                                                                                gATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGa                                                                                           >  1:208300/1‑249 (MQ=255)
                                                                                                                                                      gtgtATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGACAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAAt                                                                                   >  1:69611/1‑251 (MQ=255)
                                                                                                                                                         tATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGGCGGGTTAGCAAAGCTCGGTTTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATTGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAAttt                                                                                 <  1:345694/250‑1 (MQ=255)
                                                                                                                                                            ttGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTtctc                                                                              <  1:322946/250‑1 (MQ=255)
                                                                                                                                                                     atCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTcc                                                                    >  1:92982/1‑251 (MQ=255)
                                                                                                                                                                     atCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTcc                                                                    >  1:82120/1‑251 (MQ=255)
                                                                                                                                                                       cAGTTATTCTCAGCCTTCTTCCCAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCTTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATATCTCGCGCTGTTCctc                                                                  <  2:218756/251‑1 (MQ=255)
                                                                                                                                                                        aGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCctct                                                                 >  1:485198/1‑251 (MQ=255)
                                                                                                                                                                        aGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGAGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGCAGCCGAATTTCTCGCGCTGTTCctct                                                                 >  1:54035/1‑251 (MQ=255)
                                                                                                                                                                        aGTTAGTCTCAGCCGTCTTCCATACCTGAGTGTCCAGCGCAGTCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGGCATCAGACAGGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGGCAGACGGTCAAGAGCGAATGCCAGACCTGCGTGCGGCGGAGGACCGTACTTCAGAGCGTCGAGCAGGATGCCGAATTGCTCGCGCGGTTccgct                                                                 <  2:273227/251‑4 (MQ=255)
                                                                                                                                                                         gTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCCGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCtt                                                                >  2:408189/1‑251 (MQ=255)
                                                                                                                                                                           tATTCTCAGCCGTCTTCACAACCTGAATGCTCAGCGCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTTCTCTTc                                                               <  1:258599/250‑1 (MQ=255)
                                                                                                                                                                            aTTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTGCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTGCTCTTCg                                                              <  2:34384/250‑1 (MQ=255)
                                                                                                                                                                              tctcAGCCTTCTTCACAACCTGAATGCTCAGCTCAGACAGTGCAGTCGGGTTAGCAAATCTCGGTGTTTCAGTCATCAGACACGATGCCGACGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGCCGGTCAAAACCGAATGGCAGACCTGCGTGCGGCGGAGAACCGTCTTTCAGAGAGTCGAGCAGGAAGCCGAATTTCTCGCGCTGCTCCTCTTCGTTg                                                           >  2:280623/1‑251 (MQ=255)
                                                                                                                                                                                                    gAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACAcc                                      <  1:69891/251‑1 (MQ=255)
                                                                                                                                                                                                      aTGCGCAGCTCATCCAGTGCAGGCGGGTTAGCAAAGCTCGGTGTGTCATTCAGCAGACAGGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCCCGGATATTGTCTGCGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGCCCTTTGTGCGGGGGATGACCGTATTTCAGAGCGTCGAGCATGAAGCTGAATTTCGCTCGCCGTGCCTCTTCTTTGAttcgtagattgtcatataacgt                                    <  2:75523/251‑22 (MQ=255)
                                                                                                                                                                                                           cagctcagcCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTctgct                               <  1:426178/251‑1 (MQ=255)
                                                                                                                                                                                                             gctcagcCAGTGAAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGCTAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGGCCCGTATTTCAGAGCGTCGAGCAGGAAGCCGACTTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCtgctgc                             >  2:353584/1‑251 (MQ=255)
                                                                                                                                                                                                                cagcCAGTGCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGCTACCCAGAATACCAAACACCGTCTGCTGCat                           >  2:535889/1‑250 (MQ=255)
                                                                                                                                                                                                                cagcCAGTGCAGTCGGGTTAGCAAAGCTCGGCGTTTCAGTAATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGAAGGTCAAGACCGAATTCCAGACCTGCGTGCGGCGGAGTACCGAATTTCAGAGCGTCGAGCAGGAACCCGAATTTTTCGAGCTGTTCCTCTTCGTTGATACCCCGACAACCAAACACCTTCTGCTGCat                           >  2:444610/1‑250 (MQ=255)
                                                                                                                                                                                                                        gCAGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCAt                   <  1:4849/250‑1 (MQ=255)
                                                                                                                                                                                                                          aGTCGGGTTAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTAt                <  1:22701/251‑1 (MQ=255)
                                                                                                                                                                                                                              gggTTAGCAAAGCTCGGTGATTCAGTCCTCAGACACGCTGCCGCCGTTGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCCTGGCCAGACGGTCAAGACCGCATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCCAAGCGTCGAGCAGGAAGCCGACTTTCTCGCGCTGATCCTCTTCGTTGATAACCAGAATACCCAACACCGTCTGCCGCATATCACCATTATGGAt            >  2:27151/1‑251 (MQ=255)
                                                                                                                                                                                                                                 ttAGCAAAGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACg         >  1:205570/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tAGCAAAGCTCTGTTTTTCAGTCCTCAGACACTCTGCCGCCGGGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGTTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGGATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTTATACCCAGAAGACCAAACACCGTCTGCTGCATATCACCATTATGGATACg         <  2:25012/250‑1 (MQ=255)
                                                                                                                                                                                                                                        aGCTCGGTGTTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGaa  <  1:408189/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
AGATGACCACTAAGATCGAAACGGGACGCTTATACACTTTATCCTTCACGCTGCCTCTTCGTTGACTGCCTTCGCTCATCCCATTCACATAGTTATCTATGCTCATGGGAGTTCACTCAGTTGCCGCCTCGATGCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAA  >  NC_000913/1948584‑1949065

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