Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,156,590 C→A 5.8% intergenic (+72/‑33) yqhD → / → dkgA aldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,156,5900CA 5.8% 360.2 / 10.9 104intergenic (+72/‑33)yqhD/dkgAaldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A
Reads supporting (aligned to +/- strand):  ref base C (39/59);  new base A (3/3);  total (42/62)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-04

Only 100 of 105 total aligned reads displayed.

TATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAGCTGTGGAACGACG  >  NC_000913/3156320‑3156863
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                  
cttgtagtccccgATTTCCGCAACCCTCAATTTCTTGGAGCACTGAGGCGTTCCGACCCTCCTCTCCAACTACGTTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGTCCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCGGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCt                                                                                                                                                                                                                                                                                <  4:161353/264‑1 (MQ=255)
     aCGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCgg                                                                                                                                                                                                                                                                            <  1:127156/274‑1 (MQ=255)
      cgccgcGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCCGACTTAAGTTCCCAACACCTCAC‑CGGAGCCTGCTCCggcg                                                                                                                                                                                                                                                                          >  1:150585/1‑273 (MQ=255)
       gcagcgATGGCCGCAACCCTCAATTTCTGTGAGCAATTAGGCGTGCCGACCCACCTCACCGACTCAGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACGGGAAAAGAACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGAGGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCACAC‑CGGAGCCTGCTCCGGTg                                                                                                                                                                                                                                                                          <  2:479598/271‑1 (MQ=255)
        ccgcGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATTACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAAACTTAATTTCCCAAACCCCCCC‑CGGAACCCGCCCCGGGGaa                                                                                                                                                                                                                                                                        >  4:468160/1‑274 (MQ=255)
         cgggATGGCCGCAACCCGCAATTTCTTTGAGCAATGAGGCGTGCAGAACCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGGATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGt                                                                                                                                                                                                                                                                       <  2:171246/272‑1 (MQ=255)
         cgcgATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGt                                                                                                                                                                                                                                                                       <  2:450964/275‑1 (MQ=255)
         cgcgATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGt                                                                                                                                                                                                                                                                       <  2:275005/275‑1 (MQ=255)
          gcgATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGt                                                                                                                                                                                                                                                                       <  2:502025/274‑1 (MQ=255)
          gcgATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGt                                                                                                                                                                                                                                                                       <  2:349594/274‑1 (MQ=255)
          gcgATTGCCGCAACCCGCAATTGCTTTGGGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGt                                                                                                                                                                                                                                                                       <  2:527584/274‑1 (MQ=255)
          gcgATTGCCGCAACCCGCAAATTCTTTGAGCAAATAGGCGGTGCGACCCACCTCTCCGACTACGGTCTGGACCGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGt                                                                                                                                                                                                                                                                       <  2:190861/274‑1 (MQ=255)
                      aCCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAaggag                                                                                                                                                                                                                                                          <  4:311511/275‑1 (MQ=255)
                       cccGCATTTTCTTTGGGAAATTAGGCGTGCCGACCCCCCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAaggagg                                                                                                                                                                                                                                                         <  4:58193/275‑1 (MQ=255)
                                    tGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCCTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGCTCAAAACACCCCAA‑CGGAGCCTGCTCCGGTGAGGTCATATAAAGGAGGAACGTATGGCTaa                                                                                                                                                                                                                                            >  2:459314/1‑275 (MQ=255)
                                       gCAATTAGGCGTGCCGACCCACCTCCCCGACGACGGTCTGGACGGAAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATcc                                                                                                                                                                                                                                         <  1:71041/275‑1 (MQ=255)
                                          aTTAGGGGGTCCGACCCAACCACCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTGGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCAGTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAc                                                                                                                                                                                                                                      <  4:370157/275‑1 (MQ=255)
                                                 gTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttatt                                                                                                                                                                                                                               <  2:156338/275‑1 (MQ=255)
                                                   gCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttatt                                                                                                                                                                                                                               <  2:314680/273‑1 (MQ=255)
                                                   gCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAACGttatt                                                                                                                                                                                                                               >  4:314679/1‑273 (MQ=255)
                                                       aCCCACCTCTCCGACTACTGTCGGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTTGAAGAGCACGGCATGACCCAACTGGTCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTAGATACGAAGCCGCCCGCTAAGCTTTTTACTCATCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCAGATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGc                                                                                                                                                                                                                           <  2:516599/273‑1 (MQ=255)
                                                                 ccGACTACGTGCTGGACGGCAGTCCCATCCCGGCTGTGCTGAAAAAACTGGAAGAGAACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCa                                                                                                                                                                                                               <  4:24199/275‑1 (MQ=255)
                                                                  cGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTGCACAACACCTCAC‑CGGAGCCTGTTCCGGTGAGTTCATATAAAGGAGGAAAGTAGGGCTAATCTAACCGTTATTAAGCTACAGGATGGc                                                                                                                                                                                                                >  1:150819/1‑273 (MQ=255)
                                                                         ggTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGCCTTAAGTTCACAACACCTCAC‑CGGAGCCTGCCCCGGCGAGTTCATATAACGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGCCAATGTAAtc                                                                                                                                                                                                       >  4:77998/1‑274 (MQ=255)
                                                                         ggTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCATAACACCCCAC‑CGGCGCCTGCTCCGGTGAGTTCATATAAAGGTGGACAGTATCGCTAATCCAACCNTTATGNNACTAAAGGATGGaaaatttcat                                                                                                                                                                                                       >  3:4112/1‑268 (MQ=255)
                                                                           gttGGACGGAAGCTCCACCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGAAATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGcc                                                                                                                                                                                                     <  4:321673/273‑1 (MQ=255)
                                                                             tGGACGGCAGCTCCACCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGcc                                                                                                                                                                                                     <  1:4112/273‑1 (MQ=255)
                                                                               gcctgCAGCTCCATCCCGGCTTTTCTGCAACAACCGGAAGCGCCCCGCATGACCCCACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAATCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAg                                                                                                                                                                                                 <  4:486260/271‑1 (MQ=255)
                                                                                                 gCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGCTGGCCATGTCATGCCCCAGCTGGGAGTGGGCGTTTgg                                                                                                                                                                               >  4:301294/1‑275 (MQ=255)
                                                                                                      tgtGAAAAAACTGGAAAAGCACGGCATGACCCAACTGGGGGAAAATCATGACATTAAGTGGGATGTCAACCGCCCTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTGGGTCAGAATTCAGTTCACAACACCTCAC‑CGGAGGCTGCTCCGGTTAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGAAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaagc                                                                                                                                                                          <  4:417268/273‑1 (MQ=255)
                                                                                                         gAAAAAACTGGAAGGGCACGGCAGGCCCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGcaagcaag                                                                                                                                                                       <  4:181267/275‑1 (MQ=255)
                                                                                                             aaaCTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGACGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAAt                                                                                                                                                                   <  4:223621/275‑1 (MQ=255)
                                                                                                                   gAAGAGCACGGCATGACCCAACTGGGCGAAAATCAGGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGaa                                                                                                                                                             <  2:301295/275‑1 (MQ=255)
                                                                                                                         cACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATc                                                                                                                                                       >  4:8483/1‑275 (MQ=255)
                                                                                                                                  aCCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCAtt                                                                                                                                              >  2:317745/1‑275 (MQ=255)
                                                                                                                                       aCTGGGCGAAAATCATGACATTTCGTTGAATTTCAGCCGCCGTTTATACGAAGCCGCCCACAAAGCTTTTTACCCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaaa                                                                                                                                         <  2:296164/275‑1 (MQ=255)
                                                                                                                                        cTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGGGGACGTAATCACCGCCCTTTaaaaa                                                                                                                                        >  4:425262/1‑275 (MQ=255)
                                                                                                                                             cGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTaa                                                                                                                                                                                                                             >  2:112512/1‑185 (MQ=255)
                                                                                                                                             cGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGt                                                                                                                                    <  2:8483/274‑1 (MQ=255)
                                                                                                                                             cGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGt                                                                                                                                    <  2:425263/274‑1 (MQ=255)
                                                                                                                                             cGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCAGGGACTGGGCGTCTGGCAAGCAAGTAATGCGGACGTAATCACCGCCATTCACAAAGcatg                                                                                                                                   >  2:508708/1‑272 (MQ=255)
                                                                                                                                                             aCGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGc                                                                                                                  <  2:83566/275‑1 (MQ=255)
                                                                                                                                                              agTTGGATGGTCGCCTCCGTATATCCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTGTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCt                                                                                                                 <  1:108152/274‑1 (MQ=255)
                                                                                                                                                                   gATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGAt                                                                                                                                                                                                                   >  2:25405/1‑173 (MQ=255)
                                                                                                                                                                       tCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCa                                                                                                                                                                                                               <  4:25404/173‑1 (MQ=255)
                                                                                                                                                                        cAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAAGGAGGAAGTAATCACCGCCATTCAAAAAACGTTAGAAGGTGGGTTA‑CGCCCGATTGAAAc                                                                                                       >  4:521783/1‑275 (MQ=255)
                                                                                                                                                                          gccgccGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCgg                                                                                                                                                                                                                                                                            <  7:767083/109‑1 (MQ=255)
                                                                                                                                                                          gccgccGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCgg                                                                                                                                                                                                                                                                            >  8:767083/1‑109 (MQ=255)
                                                                                                                                                                              ccggttataCGAAGCCGCCCGCTAAGATTTTTTCGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATAccgccg                                                                                                  <  4:293100/269‑1 (MQ=255)
                                                                                                                                                                              ccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATACcgccgc                                                                                                 >  2:85957/1‑275 (MQ=255)
                                                                                                                                                                               cgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCCACGGTGAGTTCATCAAAAAGAGGAAAGTATGGCTAATCNAACCGTGATTAAGCTACAGGATGGAAATGTCATGCCCCACCAGGGACTGGGCGTCGGGCAAGCAAGAAATGAGGAAGTAANCTCCGCCATCAAAAAACAGTTAGAA‑NGGGGTGATCGCCACACTGATAACGCagcg                                                                                                >  3:95109/1‑275 (MQ=255)
                                                                                                                                                                               cgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTACCAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATACCGCcgcg                                                                                                <  2:372930/275‑1 (MQ=255)
                                                                                                                                                                                      aCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACa                                                                                         <  2:468161/275‑1 (MQ=255)
                                                                                                                                                                                         aaGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCAGGCTCCGGTg                                                                                                                                                                                                                                                                          >  8:273374/1‑96 (MQ=255)
                                                                                                                                                                                         aaGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCAGGCTCCGGTg                                                                                                                                                                                                                                                                          >  5:48368/1‑96 (MQ=255)
                                                                                                                                                                                         aaGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCAGGCTCCGGTg                                                                                                                                                                                                                                                                          <  7:273374/96‑1 (MQ=255)
                                                                                                                                                                                         aaGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCAGGCTCCGGTg                                                                                                                                                                                                                                                                          <  6:48368/96‑1 (MQ=255)
                                                                                                                                                                                           gccgccCGCTAAGCTTTTTACGCCTCACACTTTCGTTGTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGCGCCTGCTCCGGTGAGTTCCTATAAAGGAGGAACGTATGCCTAATCCAACCGTTCTTAAGCTACAGGATGTCAATGTCATGCTCCAGCTGGTACTGGGCGTCTGGCAAGCAAGTCATGCGGAAGTAACAACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATGGCTCGATGTAACCCACCGCGGCCGACAAGaaa                                                                                    >  4:388922/1‑274 (MQ=255)
                                                                                                                                                                                           gccgccCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGaa                                                                                     >  2:494210/1‑274 (MQ=255)
                                                                                                                                                                                                    caaGCCTTTTACGCCTCGACCTTTCGTTTAAGGCAATTTCGGCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCCTATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCTATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGgt                                                                           <  2:117117/274‑1 (MQ=255)
                                                                                                                                                                                                             tACGCCCCAAACTTTCGTTTTCGGGCACTTCGGCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCCCCGGTGAGTGCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCaaa                                                                  <  4:317744/275‑1 (MQ=255)
                                                                                                                                                                                                                   tCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCaa                                                                                                                                                                                                                                       >  6:1014393/1‑105 (MQ=255)
                                                                                                                                                                                                                   tCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCaa                                                                                                                                                                                                                                       <  5:1014393/105‑1 (MQ=255)
                                                                                                                                                                                                                   tCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg                                                                                                                                                                                            >  5:552772/1‑148 (MQ=255)
                                                                                                                                                                                                                    cAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCg                                                                                                                                                                                                                                                                             <  5:779561/66‑1 (MQ=255)
                                                                                                                                                                                                                    cAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCg                                                                                                                                                                                                                                                                             >  6:779561/1‑66 (MQ=255)
                                                                                                                                                                                                                      aaCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttat                                                                                                                                                                                                                                <  7:93007/109‑1 (MQ=255)
                                                                                                                                                                                                                           tCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGc                                                                                                                                                                                                >  6:1126000/1‑136 (MQ=255)
                                                                                                                                                                                                                           tCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGc                                                                                                                                                                                                <  5:1126000/136‑1 (MQ=255)
                                                                                                                                                                                                                            cGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGcc                                                   >  4:457189/1‑275 (MQ=255)
                                                                                                                                                                                                                               ttgtgggCATTTAGG‑CAGAATTAAGGTTANGACACCCCACGCGGAGCCGGCTCCGNGGAGTTCAAATAAAGGAGGAACGTATGGCTAATCCAACAGTTATTTAGCTACAGGATGGCAATGTCATGCCCGAGCGGGGNGTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCa                                                <  3:99109/271‑1 (MQ=255)
                                                                                                                                                                                                                                   gggCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAAt                                                                                                                                                                                                                                           >  7:688482/1‑85 (MQ=255)
                                                                                                                                                                                                                                   gggCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAAt                                                                                                                                                                                                                                           <  8:688482/85‑1 (MQ=255)
                                                                                                                                                                                                                                    ggCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGAt                                                                                                                                                                                                                   <  5:249912/108‑1 (MQ=255)
                                                                                                                                                                                                                                    ggCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGAt                                                                                                                                                                                                                   >  6:249912/1‑108 (MQ=255)
                                                                                                                                                                                                                                      cATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAc                                                                                                                                                                                                                                      >  6:196790/1‑87 (MQ=255)
                                                                                                                                                                                                                                      cATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAc                                                                                                                                                                                                                                      <  5:196790/87‑1 (MQ=255)
                                                                                                                                                                                                                                           cGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGagaa                                    <  4:494209/275‑1 (MQ=255)
                                                                                                                                                                                                                                              ccAGACTTAAGGTCACAACACCGCAC‑CGGCGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCAGGCCCCAGCTGGGACTGGGCGTCTGTCAAGAAAGTAATGAGGAAGTAATCACCGCCATTAAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAgaagaa                                 <  2:533354/275‑1 (MQ=255)
                                                                                                                                                                                                                                                aGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACt                               >  4:404910/1‑275 (MQ=255)
                                                                                                                                                                                                                                                 gACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGa                                                                                                                                                              <  6:552772/148‑1 (MQ=255)
                                                                                                                                                                                                                                                     tAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCa                                                                                                                                                                                                                                        >  8:170527/1‑70 (MQ=255)
                                                                                                                                                                                                                                                     tAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCa                                                                                                                                                                                                                                        <  7:170527/70‑1 (MQ=255)
                                                                                                                                                                                                                                                     tAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTg                                                                                                                                                                                              <  8:454589/112‑1 (MQ=255)
                                                                                                                                                                                                                                                     tAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTg                                                                                                                                                                                              >  7:454589/1‑112 (MQ=255)
                                                                                                                                                                                                                                                         gtCACAACACCTCCC‑CCGGGCCTGCTCCTGTTACTTCATATGACGGAGGAACGTATGGCCAATCCAACCGTCATTAAGCTACATGATGGCAAGGTCCAGCCCCAGCCGGGACTGGGCGTCTGGCAAGCAAGGAATGAAGACGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATcac                      <  4:459313/274‑1 (MQ=255)
                                                                                                                                                                                                                                                          tCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCAc                                                                                                                                                     >  8:442436/1‑148 (MQ=255)
                                                                                                                                                                                                                                                             caacaCCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAGCGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATCCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATTCCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACt                  <  2:75217/275‑1 (MQ=255)
                                                                                                                                                                                                                                                             caacaCCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaa                                                                                                                                                                            >  5:691273/1‑122 (MQ=255)
                                                                                                                                                                                                                                                             caacaCCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaa                                                                                                                                                                            <  6:691273/122‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAGAtgtg          >  1:68200/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                     cAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCGTAAAAATGCCCCAGTCAACAGAGAAAAACTGTTCACCACCACTAACAtgtg          >  4:364439/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                       c‑cGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCaa                                                                                                                                                                                                              <  5:437509/78‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cGGGGCCTGCTCCGGTGAGTTCAGAAAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAGCTGTGGaa       <  2:448517/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATg                                                                                                                                                                  >  7:415525/1‑121 (MQ=255)
                                                                                                                                                                                                                                                                         cGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATg                                                                                                                                                                  <  8:415525/121‑1 (MQ=255)
                                                                                                                                                                                                                                                                            aGCCTGCTCCGGTGAGTTCATATAAAGGAGGCACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATg                                                                                                                                                                  >  7:382963/1‑118 (MQ=255)
                                                                                                                                                                                                                                                                            aGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATg                                                                                                                                                                  <  8:382963/118‑1 (MQ=255)
                                                                                                                                                                                                                                                                             gCCTGCTCCTGTGAGTTCATATAAAGGGGGAACGTTTGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAGCTGTGGAacgac   <  4:85956/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                              ccTGCTCCGGTGAGTTCATAAAAGGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCAACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAGCTGTGGAacgacg  <  1:1171/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                  
TATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAGCACGGCATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCAC‑CGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAGAA‑GTGGGTTATCGCTCGATTGATACCGCCGCGGCCTACAAGAACGAAGAAGGTGTCGGCAAAGCCCTGAAAAATGCCTCAGTCAACAGAGAAGAACTGTTCATCACCACTAAGCTGTGGAACGACG  >  NC_000913/3156320‑3156863

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑