Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,090,968 G→A E208K (GAG→AAG)  endA → DNA‑specific endonuclease I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,090,9680GA98.8% 292.3 / ‑4.3 83E208K (GAG→AAG) endADNA‑specific endonuclease I
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (35/47);  total (35/48)
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.81e-01

GGTGAGGTGAATGGCGATCGCGGCAACTTTATGTACAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCC  >  NC_000913/3090725‑3091205
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
ggtgaggtgaATGGCGATCGCGGCAACTTTATGTACAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGa                                                                                                                                                                                                                                        <  2:412057/251‑1 (MQ=255)
    gggtgaATGGCGGTCGCGGCATCTTTAAGTACAGACAGTGGAATGGCGGTGAAGGCCAGTACGTTCAATGCGCCATGAAGGTCGATTCAAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGTCAGTGCGCGCACCTACTTCTATATGCGCGACCAAGACAACCTGATACTCTCTCGCCAGCAAACGCAGCTGATCAACGCATGGAACAAGATGTATCCGGTTACCGTCTGGAAGTGCGAgcgc                                                                                                                                                                                                                                    <  2:628247/250‑1 (MQ=255)
       tgagtGGCGATCGCGGCAACTTTATGTACAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGGCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCCTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAAAAAGATGTATCCCGTTACCGACTGGAAGGGCGGGCGCGa                                                                                                                                                                                                                                  >  1:645789/5‑250 (MQ=255)
         aaTGGCGATCGCGGCAACTTTATGTACAGCCAGTGGAATGGCGGTGAAGGCCAGCACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACTCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGa                                                                                                                                                                                                                               <  2:531133/251‑1 (MQ=255)
          aTGGCGATCGCGGCAACTGTATGTACAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGaa                                                                                                                                                                                                                              <  1:551229/251‑1 (MQ=255)
                      gCAACTTTATGTACAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAg                                                                                                                                                                                                                  <  2:645789/251‑1 (MQ=255)
                        aaCTTTATGTACAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGt                                                                                                                                                                                                                <  1:408287/251‑1 (MQ=255)
                        aaCTGTATGTACAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCTCCAGCGCGTGCACGCGGTGCCATTGCGCGCTCCTGCTTCTATATGCGCGACCAATACAACCTGACACACTCTCGCCAGCAAACGCAGCGGTTCAACGCATGGAACAAGATGGATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGt                                                                                                                                                                                                                <  1:701618/251‑1 (MQ=255)
                            ttaggtccAGCCAGTGGAATGGCGATGAAGGCCAGTACGGTCAACGCGCCAGTGCGGTCGATTCCAAAGAAAAAGCTTCCTAACCACCAGCGCGTGCACGCGGTGCCATTCCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCCAACGCAGCTGTGCAACGCATGGAACAAGATGGAGCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCa                                                                                                                                                                                                             <  2:138468/243‑1 (MQ=255)
                            tgaTGTACAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGc                                                                                                                                                                                                              <  2:454590/247‑1 (MQ=255)
                              aTGTACAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAggg                                                                                                                                                                                                          <  2:400900/251‑1 (MQ=255)
                                gTACAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGACAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCCATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGGTACCGAATGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGGGCAGGGc                                                                                                                                                                                                         >  1:496351/1‑250 (MQ=255)
                                  aCAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCaa                                                                                                                                                                                                       <  1:252710/250‑1 (MQ=255)
                                  aCAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAAt                                                                                                                                                                                                      <  2:568614/251‑1 (MQ=255)
                                        aGTGGAATGGCGGTGAAGGCCCAGACGGTCAATGCGCCATGAAGGGCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAAACTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAAATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATaa                                                                                                                                                                                                 >  2:31378/1‑250 (MQ=255)
                                             aaTGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTa                                                                                                                                                                                           <  1:161595/251‑1 (MQ=255)
                                             aaTGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATAGGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTa                                                                                                                                                                                           <  1:506387/251‑1 (MQ=255)
                                                ggCGGTGCAGGCCAGTACGGTCAATGCGCCATGACGGTCGCTTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGACCGCATCGCGAAGGTGCAGGGCAATCATAACCCGTAtgt                                                                                                                                                                                        >  2:89949/1‑251 (MQ=255)
                                                ggCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTAtgt                                                                                                                                                                                        >  1:457136/1‑251 (MQ=255)
                                                ggCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTAtgt                                                                                                                                                                                        >  2:258400/1‑251 (MQ=255)
                                                  cGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATGGCGCGCACCTTCTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGc                                                                                                                                                                                      <  1:252945/251‑1 (MQ=255)
                                                    gTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCa                                                                                                                                                                                     >  1:595686/1‑250 (MQ=255)
                                                         ggCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACgcgc                                                                                                                                                                               >  2:246686/1‑251 (MQ=255)
                                                         ggCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACgcgc                                                                                                                                                                               >  2:87596/1‑251 (MQ=255)
                                                         ggCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTGCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCATGGCAATCATAACCCGTATGTGCAACgcgc                                                                                                                                                                               <  2:687045/251‑1 (MQ=255)
                                                                 aCGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGCAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGcc                                                                                                                                                                          >  2:501121/1‑249 (MQ=255)
                                                                 aCGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGcc                                                                                                                                                                          >  1:544746/1‑249 (MQ=255)
                                                                  cGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAgg                                                                                                                                                                       <  1:31378/251‑1 (MQ=255)
                                                                                  gAAGGTCGATGTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTGCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTaa                                                                                                                                                         <  2:672203/250‑1 (MQ=255)
                                                                                   aaGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCTCCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAGCGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGGTACCGACTGGAAGGGCGAGCGCGATGAACGCATTGCGAAGGTGCAGGGCAATCATAACCGGTATGGGCATCGCGGATGGCAGGTGTGAAAGAGCTAAcc                                                                                                                                                       <  2:35329/251‑1 (MQ=255)
                                                                                    aGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCTCGAAAGAGCTAACCt                                                                                                                                                      <  1:258400/251‑1 (MQ=255)
                                                                                       tCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTACCCTac                                                                                                                                                    >  1:387796/1‑250 (MQ=255)
                                                                                        cGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGAAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTacac                                                                                                                                                  >  2:218833/1‑251 (MQ=255)
                                                                                              cAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCgg                                                                                                                                            <  1:606925/251‑1 (MQ=255)
                                                                                                 aGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGAt                                                                                                                                         <  2:338734/251‑1 (MQ=255)
                                                                                                     aaaGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGACTCtt                                                                                                                                     >  2:438624/1‑251 (MQ=255)
                                                                                                          tGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGCCTGGAAGTGCGAGCGCGATGAACGCATAGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGt                                                                                                                                >  2:302999/1‑251 (MQ=255)
                                                                                                          tGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGt                                                                                                                                >  1:636739/1‑251 (MQ=255)
                                                                                                               aaCCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAAcc                                                                                                                           >  1:223516/1‑251 (MQ=255)
                                                                                                                   accaGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCt                                                                                                                        >  2:70146/1‑250 (MQ=255)
                                                                                                                     caGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAATTGCGAGCGCGATGAACGCATCGCGAAGTTGCAGGGCAATCATAACCCGTAGGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTAcc                                                                                                                     <  1:441682/251‑1 (MQ=255)
                                                                                                                         gcgTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCAc                                                                                                                 <  1:120257/251‑1 (MQ=255)
                                                                                                                         gcgTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCAc                                                                                                                 <  1:618030/251‑1 (MQ=255)
                                                                                                                          cgTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACTCTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCAcg                                                                                                                <  1:620189/251‑1 (MQ=255)
                                                                                                                           gTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCAcgc                                                                                                               >  1:500668/1‑251 (MQ=255)
                                                                                                                                cgcgGTGCCCATGCACGCCCCTAATTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGGATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCACAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTGGTTAACCCCTACCCCACGCGTACa                                                                                                          <  2:497669/251‑1 (MQ=255)
                                                                                                                                cgcgGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTAc                                                                                                           >  1:589781/1‑250 (MQ=255)
                                                                                                                                        cATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACGCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTg                                                                                                  <  2:118891/251‑1 (MQ=255)
                                                                                                                                          ttGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTggg                                                                                                >  1:331638/1‑251 (MQ=255)
                                                                                                                                           tGCGCGCACCTACTTCTATAGGCGCGACCAATACAACCTGACACCCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTGTTGTTAACCCCTACCCCACGCGTACAACCGCGTggg                                                                                                <  2:55933/250‑1 (MQ=255)
                                                                                                                                           ggcgctccccTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCGTGGAACAAGATGTATCCGGTTACCGACTGTAAGTGCGAGCGCGATTAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTTCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGGTTTTTTTTGTTAACCCCTGCCCCGCTCGTACAACCTCTTgggg                                                                                               <  2:730475/243‑1 (MQ=255)
                                                                                                                                            gcgcgcACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTgggg                                                                                               <  2:387796/250‑1 (MQ=255)
                                                                                                                                               cgccccTACTGCTATATGCGCGACCAATACAACCTGACTCTCTCGCGCCATCAAACGCAGCTGGTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGATCGCGATGAACGCATCGCGGAGGTGCAGGGCAAGCATAACCCGTATGTTCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACCCTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAgac                                                                                           <  1:382748/247‑1 (MQ=255)
                                                                                                                                                gcACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCa                                                                                                                  <  2:556944/227‑1 (MQ=255)
                                                                                                                                                gcACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCa                                                                                                                  >  1:556944/1‑227 (MQ=255)
                                                                                                                                                 cACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAGCCGCGTGGGGAgacg                                                                                          >  2:133257/1‑250 (MQ=255)
                                                                                                                                                  aCCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAgacgac                                                                                        <  1:488217/251‑1 (MQ=255)
                                                                                                                                                      aCTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCgg                                                                                    <  1:416901/251‑1 (MQ=255)
                                                                                                                                                       cTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGCTTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCgg                                                                                    >  1:34189/1‑250 (MQ=255)
                                                                                                                                                        ttCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGAt                                                                                  >  1:391054/1‑251 (MQ=255)
                                                                                                                                                            atatGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGAttttt                                                                              >  2:202978/1‑251 (MQ=255)
                                                                                                                                                              atGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTa                                                                             >  1:617518/1‑250 (MQ=255)
                                                                                                                                                                 cgcgACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTa                                                                         <  1:202978/251‑1 (MQ=255)
                                                                                                                                                                  gcgACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTAt                                                                        <  1:406061/251‑1 (MQ=255)
                                                                                                                                                                    gACCAATACACCCTGACACTCTCTCGCCAGCACACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCAAGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGc                                                                      >  2:28942/1‑251 (MQ=255)
                                                                                                                                                                     aCCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCg                                                                     <  2:616271/251‑1 (MQ=255)
                                                                                                                                                                         aTACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTCACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATc                                                                 <  1:705840/251‑1 (MQ=255)
                                                                                                                                                                           aCAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATccc                                                               <  2:180972/251‑1 (MQ=255)
                                                                                                                                                                           aCAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATccc                                                               <  1:501121/251‑1 (MQ=255)
                                                                                                                                                                                       tctctcGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTatc                                                    >  1:598037/1‑250 (MQ=255)
                                                                                                                                                                                        ctctcGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCtt                                                                                                                                                                             >  1:633862/1‑128 (MQ=255)
                                                                                                                                                                                        ctctcGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCtt                                                                                                                                                                             <  2:633862/128‑1 (MQ=255)
                                                                                                                                                                                           tcGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCt                                               >  1:642786/1‑251 (MQ=255)
                                                                                                                                                                                            cGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCt                                               <  1:13292/250‑1 (MQ=255)
                                                                                                                                                                                              ccAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCTAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGaa                                            <  2:41477/251‑1 (MQ=255)
                                                                                                                                                                                                 gCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCa                                         <  1:87596/251‑1 (MQ=255)
                                                                                                                                                                                                      cGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGAc                                    >  1:94428/1‑251 (MQ=255)
                                                                                                                                                                                                           gtGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATACCCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTACCCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGATCCACTGACCAGcc                               <  1:197903/250‑1 (MQ=255)
                                                                                                                                                                                                             gTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCAt                             <  1:528795/251‑1 (MQ=255)
                                                                                                                                                                                                                      aTGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATc                    >  2:544325/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ccGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATg    >  1:242630/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATgc   <  2:223516/251‑1 (MQ=255)
                                                                                                                                                                                                                                        ggTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATgcc  >  2:105195/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
GGTGAGGTGAATGGCGATCGCGGCAACTTTATGTACAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCC  >  NC_000913/3090725‑3091205

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