Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,515,348 A→C L444V (TTG→GTG)  fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,515,3480AC97.1% 349.9 / ‑4.9 105L444V (TTG→GTG) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base C (47/55);  total (50/55)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 105 total aligned reads displayed.

AGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCAT  >  NC_000913/4515104‑4515590
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
aGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATTCCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCg                                                                                                                                                                                                                                               <  1:351167/251‑1 (MQ=255)
  cGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGAt                                                                                                                                                                                                                                             <  1:158417/251‑1 (MQ=255)
  cGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGAt                                                                                                                                                                                                                                             <  1:543449/251‑1 (MQ=255)
  cGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGGGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGAt                                                                                                                                                                                                                                             >  1:199668/1‑251 (MQ=255)
  cGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGGGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGAt                                                                                                                                                                                                                                             >  1:102849/1‑251 (MQ=255)
     tGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑CCGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAg‑c                                                                                                                                                                                                                                         >  1:631706/1‑251 (MQ=255)
          tataGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAg‑cccacc                                                                                                                                                                                                                                    >  1:669701/1‑250 (MQ=255)
             aGCTCACTTCTACGTGCGCTCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAg‑cccacgccc                                                                                                                                                                                                                                 <  1:72580/251‑1 (MQ=255)
              gCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAg‑cccacgccc                                                                                                                                                                                                                                 >  1:520590/1‑250 (MQ=255)
               cTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATAACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAg‑cccacgcccac                                                                                                                                                                                                                               >  1:601130/1‑251 (MQ=255)
                 cccttcttCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACtt                                                                                                                                                                                                                             <  1:117374/249‑1 (MQ=255)
                    ttcttcGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCCTGCGTTGTTTCATTCACATAGCGATAG‑CCCACGCCCACTTCGt                                                                                                                                                                                                                          >  1:298572/1‑251 (MQ=255)
                     tcgtcGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCg                                                                                                                                                                                                                           <  2:596004/246‑1 (MQ=255)
                       ttcGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGGAGTAACGCATTTCCTGCGTTGATTCATTAACATAGCGATAG‑CCCACGCCCCCCTCGTgcg                                                                                                                                                                                                                       >  2:171759/1‑251 (MQ=255)
                        tcGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTgc                                                                                                                                                                                                                        <  1:557836/249‑1 (MQ=255)
                         cggtcgtgtCTGGGATGGCGGTGTTCTGGTATGACTCGATTTTTTCGAAGCGCATCCCCGGCGTGATGGGCCAGTTGCCGATGTCGAGTTTGTCATCCAGATACCAGGCTTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCATCTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGcgcg                                                                                                                                                                                                                     <  2:276853/247‑1 (MQ=255)
                          gtgggtgcCTGGGATGGCTTTGTTCTGGTAGGACTCGATATGTGCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCgg                                                                                                                                                                                                                    <  2:567166/247‑1 (MQ=255)
                          gtgcgtgcCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCAGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCACGCGTTGATTCATTAACATAGAGGTAg‑cccacgcccacatagttcgcgg                                                                                                                                                                                                                    >  2:501544/1‑240 (MQ=255)
                            gcgggcCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGTAACTCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGaa                                                                                                                                                                                                                  <  2:351212/247‑1 (MQ=255)
                             cgtgcCTGTGAGGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGGTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGGATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACGTCGTGCGCGGAAg                                                                                                                                                                                                                 <  2:30445/251‑1 (MQ=255)
                              gtgccggggATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCACACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑TCCACGCCCACTTCGTGCGCGGAAgg                                                                                                                                                                                                                <  2:473633/243‑1 (MQ=255)
                              gtgcCTGTTATGGCGTTTTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑CCGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCCTTTCCGGCGTAGATTCATTCACCTAGCGATAG‑CCCACCCCCACTTCGTGCGCGGAAgg                                                                                                                                                                                                                >  2:651045/1‑251 (MQ=255)
                              gtgcCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAgg                                                                                                                                                                                                                <  1:42231/251‑1 (MQ=255)
                               tgcctgggctgtctgtgttCTGGTATGCCTCGATATGTTCGAAACGCATTCTCTGCGTGTTGGTCCAGTTGCCGAGGGCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGTATCGCGTATCGCGGTCGTAAGGGCTTGTGCCGG‑ACGTCAACTGCCCGCTGCTGGTGGCTGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGAGAT‑CCCACGCCCACTTTGTGCGCGGAAg                                                                                                                                                                                                                 <  2:393614/234‑1 (MQ=255)
                                gcCTGGGATGGCGTTGTGTTGGTATGACTCGATATGTTCGAAACGCAGACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGGATCGCGGGCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATGCATTCACATAGCGATAG‑CCCACGCCCACTGCGTGCGCGGAAggg                                                                                                                                                                                                               <  2:603831/250‑1 (MQ=255)
                                 ccTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGcc                                                                                                                                                                                                             <  1:305243/251‑1 (MQ=255)
                                         ggCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATaa                                                                                                                                                                                                     <  1:328122/251‑1 (MQ=255)
                                         ggCGTTGTTCTGGTATGACGCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATaa                                                                                                                                                                                                     <  1:384024/251‑1 (MQ=255)
                                                ttcagggataaCTCGATATGTTAGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTAGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAAAGCGTATCGCGGTAGTAAGGGCTTGAGCCGGCACG‑CAACTGCCCGCTGCTGGTAGCGGTGGAGTAACGACGTTCCTGCGTTAATTAATTCAAGTAGCCATAGGCCCA‑TACCAATTCGTGCGCGGGCGGACCGATCCTAAAGAACTg                                                                                                                                                                                              >  2:217724/11‑251 (MQ=255)
                                                   ccgaactaCTCGATATGGTCGAAACGCATACCCGGCGTGATGGTCCAGCTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCTG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCt                                                                                                                                                                                           <  2:157315/244‑1 (MQ=255)
                                                    ggTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGCGCCGG‑CCGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACACATTTCATGGGTTGATTCAATCACATGGCGATAC‑CCCACGCCCCCTTCATGCGCGGGAGGCCAGATCATAAAGATCTGGCTg                                                                                                                                                                                          >  2:207638/1‑251 (MQ=255)
                                                      tATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGCTCATAAAGATCTGGCTGt                                                                                                                                                                                         >  1:152457/1‑250 (MQ=255)
                                                       aTGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGt                                                                                                                                                                                         >  1:199421/1‑249 (MQ=255)
                                                          aCTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTTGCt                                                                                                                                                                                           >  1:157315/1‑244 (MQ=255)
                                                             cGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGc                                                                                                                                                                                 >  1:651321/1‑251 (MQ=255)
                                                              gATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCt                                                                                                                                                                                <  1:193739/251‑1 (MQ=255)
                                                              gATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAAATGCCCGCTGCTGGTGGCGGTGTAGTAACGCCTTTCATGGGTTGATTCATTCACATAGCGATAGACCAACACCC‑CTTCGTGCGCGGAAAGGCCTATCATAAAGATTTGGCttgagcagggcc                                                                                                                                                                                >  2:594204/1‑240 (MQ=255)
                                                                 atGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATCCCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCa                                                                                                                                                                              <  2:619792/250‑1 (MQ=255)
                                                                   gTTCGAAACGCATACCCGGCGTGATGGTCCAGTGGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAAt                                                                                                                                                                            <  1:665963/250‑1 (MQ=255)
                                                                    ttCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGTTTCATTCCCATAGCGATAA‑CCCCCGCCCACTTCGTGCGCGGAAGGGCCGATCCAAAAGATATGGCTGTAGCGTGGCCCAATAc                                                                                                                                                                          >  2:475468/1‑251 (MQ=255)
                                                                     tCGAAACGCATACCCGGCGTGATGGTTCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGGCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATAc                                                                                                                                                                          <  1:150349/250‑1 (MQ=255)
                                                                     gcGAAACGCATCACCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCTGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGGTCATAAAGATCTGGCTGTAGCGTGGCTCAATAc                                                                                                                                                                          <  1:576217/249‑1 (MQ=255)
                                                                            gCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAAAGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGAGGCTCAATACCGCGCAcc                                                                                                                                                                  >  1:193324/1‑251 (MQ=255)
                                                                              aTACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGGTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCAccc                                                                                                                                                                 >  1:139968/1‑250 (MQ=255)
                                                                                  ccGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACt                                                                                                                                                                                                                              >  1:385601/1‑185 (MQ=255)
                                                                                  ccGGCGGGATGGGCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACt                                                                                                                                                                                                                              <  2:385601/185‑1 (MQ=255)
                                                                                    ggcttGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCCCCCAGTAGtt                                                                                                                                                          <  2:32807/247‑1 (MQ=255)
                                                                                    ggCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACcg                                                                                                                                                                        >  1:521965/1‑237 (MQ=255)
                                                                                    ggCGTGATGGGCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACcg                                                                                                                                                                        <  2:521965/237‑1 (MQ=255)
                                                                                     gCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTa                                                                                                                                                         <  1:590392/251‑1 (MQ=255)
                                                                                      cGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTAc                                                                                                                                                        <  1:171759/251‑1 (MQ=255)
                                                                                       gTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑CCGGCAACTGCCCGCTGCTGGTGGCGGTGGAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGCTAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAACGTGGCCCAATACCGCGCACCCCGTAGTTAcg                                                                                                                                                       >  2:325073/1‑251 (MQ=255)
                                                                                          aTGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACgcgg                                                                                                                                                    >  1:168636/1‑251 (MQ=255)
                                                                                          aTGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGAGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACgcgg                                                                                                                                                    >  1:593432/1‑251 (MQ=255)
                                                                                           gggTCCAGTGGCCGATGTCGATTTTGTCATCCAGATCCCAGGCGTGCGCCTCGGGGCCGGAACGCGTATCTCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCAGGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACgcggc                                                                                                                                                   <  2:651321/250‑1 (MQ=255)
                                                                                            ggTCCCGTCGCCGATGCCGATTTTGTCATCCAGCTACCAGGCGTGCGCCTCGGTGCCGGAACGCGTCTCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCCTGCGTTGATTCATTCACATAGCGATAG‑CCCACCCCCACTTCGTGCGCGGAAGGGCCGAACATAAAGATCTGGCTGTCGCGTGGCCCAATAccccgacccccggagttaaccggcg                                                                                                                                                  >  2:635124/1‑228 (MQ=255)
                                                                                               ccAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCgaga                                                                                                                                               <  2:252057/251‑1 (MQ=255)
                                                                                                cAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCCTGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTTGCTGTAGCCGGGCTCAATTCAGCGCACCCAGCAGTTAcgccggggagg                                                                                                                                              >  2:644217/1‑243 (MQ=255)
                                                                                                        gATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGc                                                                                                                                      >  1:157275/1‑251 (MQ=255)
                                                                                                         aTGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCg                                                                                                                                     >  1:619973/1‑251 (MQ=255)
                                                                                                           gTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGt                                                                                                                                    <  1:475468/250‑1 (MQ=255)
                                                                                                              gATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGc                                                                                                                                <  1:207868/251‑1 (MQ=255)
                                                                                                                ttttGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCt                                                                                                                              >  1:72941/1‑251 (MQ=255)
                                                                                                                    gTCATCCATATACCAGGCGTGCGCCTCGGTGCCGGATCGCGTATCGCGTTCTTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATTCCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGc                                                                                                                           <  2:157275/250‑1 (MQ=255)
                                                                                                                            gtttccAGGCGGGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAATGGCTTGTGCCGT‑ACGGCAACTGCCGGCTGCTTGTGGTGGTGTAGTACCGCATTTCACGCGTTGATTCATTCACTTAGCGATAG‑CCCACGCCCTCTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATCCCTCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAg                                                                                                                  <  2:162200/247‑1 (MQ=255)
                                                                                                                                ccAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGGGATGCGTTTGCCTTGCTCCAGGTAGCCGc                                                                                                              >  1:256916/1‑251 (MQ=255)
                                                                                                                                    gCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTgcg                                                                                                          >  1:230465/1‑251 (MQ=255)
                                                                                                                                     cGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTgcgc                                                                                                         <  2:168217/251‑1 (MQ=255)
                                                                                                                                         cgcCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAggg                                                                                                     <  1:170149/251‑1 (MQ=255)
                                                                                                                                          gccccggggccggagtgcgattcgtggTCGTAAGGGCTTGGGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGGTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTTGCCGCTGCGCAGGGt                                                                                                    <  2:229260/226‑1 (MQ=255)
                                                                                                                                          gcCTCGGTGCCGGAACGCGTATCGCGGGCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTGCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGt                                                                                                    <  1:619016/251‑1 (MQ=255)
                                                                                                                                           ccTCGGTGCCGGAACGCGTATCGCGGTCGAAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGCAGTAACGGATTTCATGCGTTGATTCATTCACATAACGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGGAGCGTGGCTCAATACCTCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCAGAGCAGGGtt                                                                                                   >  1:634399/1‑251 (MQ=255)
                                                                                                                                            cccggggcctgAACGCGTATCGCGTTCGTAATGGCTTGAGCCGC‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGCAGTAACCCATTTCATGCTTTGATTCATTCACATATCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAACACCGCGCACCCAGTAGTTACGCGGCGAGATGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAggggtt                                                                                                  <  2:218013/241‑4 (MQ=255)
                                                                                                                                                      gAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGCTTCATTCACATAGCGCTAG‑CCCCCGCCCACTTCGTGCGCGGGAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCCCCCAGTAGTTAAGCGGCGCGAGGGGGATGCGGTTGCCTTGCTCCAGGTAGACGCTGCGCAGGGGTTGGGTGTAgc                                                                                         >  2:665177/1‑249 (MQ=255)
                                                                                                                                                              attgcgGTCTTCAGGGCTTGAGCAGT‑GCGGCAACTTCCCGCTGCTGGTGGCGGTGTATTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGTTCATAAAGATCTTGCTGTAGTGTGGCTCAATACCGCGCACCCAGTATTTACGCGCCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTTCGCAGTGTTTGGTTGTAGATCCCCTGa                                                                                 <  2:449412/247‑1 (MQ=255)
                                                                                                                                                                       gTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGGGGGTGATGCGTTTGCCTTGCTCCCGGGAGCCGCTGCGCAGGGGTTGGGTGTGGAACCCCTTGATGTTGAAtt                                                                       >  2:202996/1‑251 (MQ=255)
                                                                                                                                                                                 gAGCCGG‑ACGGCCACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTGCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCAGATCATAAAGATCGGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTGTGCCTTGCTCCAGGTATCCGCTGCGCAGGGGTTGGGTGTAGAACCGCTGAATGTTGAATTTATgctggc                                                              <  2:385517/250‑1 (MQ=255)
                                                                                                                                                                                 gAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTTCGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATgctggc                                                              <  1:46477/250‑1 (MQ=255)
                                                                                                                                                                                   gCCGG‑ACGGCATCTGCCCGCTGCTGGTGGCGGGGTATTAACGCAGTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGt                                                           <  2:653607/251‑1 (MQ=255)
                                                                                                                                                                                   gCCGG‑ACGGCAACTGCCCGCTGCTGGTGTCGGTGTAGTAACGCATTTCATGCTTTGATTCAGTCACATAGCGATAT‑CCCACGCCCACTTCGTGCGCGGAATGGCCGATCATAAAGATCTGGCTTTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGTTGAGTCGTTTGCCTTGCTCCAGTTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGt                                                           <  1:435237/251‑1 (MQ=255)
                                                                                                                                                                                    ccGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTc                                                          <  1:389066/251‑1 (MQ=255)
                                                                                                                                                                                     cGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAg                                                                                                       <  2:369891/205‑1 (MQ=255)
                                                                                                                                                                                     cGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAg                                                                                                       >  1:369891/1‑205 (MQ=255)
                                                                                                                                                                                          cGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTctggct                                                     >  1:293203/1‑251 (MQ=255)
                                                                                                                                                                                           ggCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGCATGTTGAATTTATGCTGGCTGTctggctg                                                    >  1:295334/1‑251 (MQ=255)
                                                                                                                                                                                           ggCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGGAGTTAAGCGGCGAGAGGGGGATGCGTTTGCCTTGCTCCAGGGAGCCGCTGCGCAGGGTTTGGGGGTTGAACCCCGGGATGTTGAATTTATgctggcctgccggctg                                                    >  2:65192/1‑241 (MQ=255)
                                                                                                                                                                                                cTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACt                                               >  1:505088/1‑251 (MQ=255)
                                                                                                                                                                                                      gctgctGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGc                                         >  1:336640/1‑251 (MQ=255)
                                                                                                                                                                                                         gctgGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCa                                      <  1:450222/251‑1 (MQ=255)
                                                                                                                                                                                                              tggCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCCCTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATGCCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAAccc                                                                                     <  1:286632/199‑1 (MQ=255)
                                                                                                                                                                                                              tggCGGTGTAGGAACGCATGTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGGGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAg                                                                                                       <  1:76552/181‑1 (MQ=255)
                                                                                                                                                                                                                gCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCCCTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGc                               <  1:207638/251‑1 (MQ=255)
                                                                                                                                                                                                                 gggTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCCCGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGGGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGcc                              <  2:364273/250‑1 (MQ=255)
                                                                                                                                                                                                                      tagtaACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCtt                         >  1:217235/1‑251 (MQ=255)
                                                                                                                                                                                                                       agtaACGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTg                        >  1:538294/1‑251 (MQ=255)
                                                                                                                                                                                                                          aaCGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGAGCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGTCTCGCCAGCTTGcga                     <  1:152412/251‑1 (MQ=255)
                                                                                                                                                                                                                            cGCATTTCATGCGTTGATTCCTTCACATAGCGATAG‑CCCACGCCCAATTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGCAGTTACGCGGCGAGAGGGCGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGCGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCCCGcccgtcttcgcac                   >  2:301287/1‑240 (MQ=255)
                                                                                                                                                                                                                            cGCATTTCATGCGTTGATTCATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGcgacg                   >  1:253889/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ttGATACATTCACATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCCTCAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGGTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGAAACTGGTAGACCAGGCTCGACAGCTTGCGACGACCCCAGATGCgt      >  2:358824/1‑250 (MQ=255)
                                                                                                                                                                                                                                              tcattcaCATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGTGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCAt  >  1:584123/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
AGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGG‑ACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAG‑CCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCAT  >  NC_000913/4515104‑4515590

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