Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 689,840 T→C H72R (CAT→CGT)  unknown ← N‑acetylmuramoyl‑L‑alanine amidase (EC 3.5.1.28)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1689,8400TC100.0% 171.7 / NA 67H72R (CAT→CGT) unknownN‑acetylmuramoyl‑L‑alanine amidase (EC 3.5.1.28)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/37);  total (30/37)

CGAGCACGCTACCCGCTCCGCCCAACAGCTCCGAACGTTTTTCATGATCTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAATGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTACATCATCGTTATTTAAT  >  1/689595‑690078
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              
cGAGCACGCTACCCGCTCCGCCCAACAGCTCCGAACGTTTTTCATGATCTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGa                                                                                                                                                                                                                                           <  1:40681/251‑1 (MQ=255)
    cACGCTACCCGCTCCGCCCAACAGCTCCGAACGTTTTTCATGATCTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATc                                                                                                                                                                                                                                         >  1:893415/1‑249 (MQ=255)
      cGCTACCCGCTCCGCCCAACAGCTCCGAACGTTTTTCATGATCTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTa                                                                                                                                                                                                                                     >  1:721220/1‑251 (MQ=255)
        cTACCCGCTCCGCCCAACAGCTCCGAACGTTTTTCATGATCTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTAt                                                                                                                                                                                                                                    >  1:163814/1‑250 (MQ=255)
                   gCCCAACAGCTCCGAACGTTTTTCATGATCTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCtt                                                                                                                                                                                                                        >  1:192374/1‑251 (MQ=255)
                      cAACAGCTCCGAACGTTTTTCATGATCTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTa                                                                                                                                                                                                                     <  1:835339/251‑1 (MQ=255)
                       aaCAGCTCCGAACGTTTTTCATGATCTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTaa                                                                                                                                                                                                                    <  1:500321/251‑1 (MQ=255)
                                          cATGATCTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAgg                                                                                                                                                                                                 <  1:245241/251‑1 (MQ=255)
                                               tCTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTa                                                                                                                                                                                             >  1:679812/1‑250 (MQ=255)
                                                cTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACa                                                                                                                                                                                           >  1:517433/1‑251 (MQ=255)
                                                cTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACa                                                                                                                                                                                           >  1:174828/1‑251 (MQ=255)
                                                 ttCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACAt                                                                                                                                                                                          <  1:178048/251‑1 (MQ=255)
                                                 ttCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACAt                                                                                                                                                                                          <  1:463528/251‑1 (MQ=255)
                                                                       tCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGt                                                                                                                                                                     <  1:559203/250‑1 (MQ=255)
                                                                                      tttGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAgcgc                                                                                                                                                      <  1:302204/250‑1 (MQ=255)
                                                                                         gACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCggg                                                                                                                                                  >  1:231651/1‑251 (MQ=255)
                                                                                                cccATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGtt                                                                                                                                                                                              <  1:550068/200‑1 (MQ=255)
                                                                                                cccATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCttt                                                                                                                                            >  1:607068/1‑250 (MQ=255)
                                                                                                   aTACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGaaa                                                                                                                                                                                                   <  1:73/192‑1 (MQ=255)
                                                                                                   aTACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGa                                                                                                                                                                                                     <  1:82982/190‑1 (MQ=255)
                                                                                                   aTACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCg                                                                                                                                        >  1:897432/1‑251 (MQ=255)
                                                                                                                 cgTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGAt                                                                                                                           <  1:160546/250‑1 (MQ=255)
                                                                                                                 cgTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATc                                                                                                                          <  1:848267/251‑1 (MQ=255)
                                                                                                                    aaaTTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATGGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAAc                                                                                                                       <  1:515537/251‑1 (MQ=255)
                                                                                                                       ttCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATCCGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCGGGGGTCTTTCCCGCCATGACCGGGATCAACGGc                                                                                                                    <  1:782953/251‑1 (MQ=255)
                                                                                                                         cGGGTTTTCTGAGAAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTCCGCGCAATTTCGGACCGTTGCGGCACGAAAATAAAGTAATCGCCGTTGCGGGTAAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGTTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCaa                                                                                                                  <  1:23404/251‑1 (MQ=255)
                                                                                                                           gggTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAggg                                                                                                                                                                                                <  1:847361/171‑1 (MQ=255)
                                                                                                                             gTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGt                                                                                                              >  1:632981/1‑251 (MQ=255)
                                                                                                                                tCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTcc                                                                                                            >  1:211650/1‑250 (MQ=255)
                                                                                                                                        aTCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTAc                                                                                                   <  1:147716/251‑1 (MQ=255)
                                                                                                                                                 aaTGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCgggg                                                                                                                                                 >  1:466052/1‑196 (MQ=255)
                                                                                                                                                         ccAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGTAATAAAGTAATCGCCGTTGCGGGTTAATACAGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCCCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCa                                                                                    >  1:561739/1‑249 (MQ=255)
                                                                                                                                                                         gtgCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAgg                                                                                                                                                                                                 >  1:521719/1‑124 (MQ=255)
                                                                                                                                                                                cgcgAAATTTCGGAACGTTGCGGCACGAAAATAAAGTAATCGCCGTTGCGGGTTAATCCGGCATGAAAACGCGGATCCTTATCAAAAAGCGTTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGAc                                                           <  1:573731/251‑1 (MQ=255)
                                                                                                                                                                                   cggggagatgtgaaaaagagacaggAAAAAAAGTAATCGCCGTGGGGGGTTAATACGGCATGAAACCGCGGATCTTTATATAAAAGCGTTTTTAATTCCGGGGGGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCGGATGGCGCCGGGGTTTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAg                                                        <  1:894159/228‑1 (MQ=255)
                                                                                                                                                                                      aTTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGaaa                                                     >  1:900472/1‑251 (MQ=255)
                                                                                                                                                                                       tttCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAaat                                                    >  1:91357/1‑251 (MQ=255)
                                                                                                                                                                                       tttCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAACTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAaat                                                    >  1:313331/1‑251 (MQ=255)
                                                                                                                                                                                            gAACGTTGCGGCACGAAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGTTCCTTATCTAAAAGCGCTTTTAATTCCCGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAAtt                                                <  1:707271/250‑1 (MQ=255)
                                                                                                                                                                                                       aaagaaaaTAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTGGCGCCCGATGGCGCCGGGGTTTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGc                                    <  1:534376/246‑1 (MQ=255)
                                                                                                                                                                                                         cGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGa                                          <  1:676312/243‑1 (MQ=255)
                                                                                                                                                                                                           gAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTc                                                                                                                                               >  1:236998/1‑140 (MQ=255)
                                                                                                                                                                                                           gAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCa                                                                                                                                     >  1:626060/1‑150 (MQ=255)
                                                                                                                                                                                                           gAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAg                                >  1:817228/1‑251 (MQ=255)
                                                                                                                                                                                                               aaaaGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAGGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCGGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTc                            <  1:829495/250‑1 (MQ=255)
                                                                                                                                                                                                                   gTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAAttt                                               <  1:162973/228‑1 (MQ=255)
                                                                                                                                                                                                                   gTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATtg                        >  1:832392/1‑251 (MQ=255)
                                                                                                                                                                                                                   gTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATACGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATtg                        >  1:441856/1‑251 (MQ=255)
                                                                                                                                                                                                                    tAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGTTTTTAATTCCGGGGCGATAGAAAGGTTTACTTTTTTTTCATAATTTTTTAGGTTGCGCCCGATAGCGCGGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATtg                        <  1:245403/250‑1 (MQ=255)
                                                                                                                                                                                                                    tAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGTTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCGGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATtg                        <  1:662164/250‑1 (MQ=255)
                                                                                                                                                                                                                    tAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGTTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCGGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATtg                        <  1:198748/250‑1 (MQ=255)
                                                                                                                                                                                                                      aTCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCTTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAGGGGTCACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTgg                     <  1:424880/251‑1 (MQ=255)
                                                                                                                                                                                                                          ccGTTGCGGGTTAATCCGGCATGAAACCGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAAAAAGGTTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTAc                  <  1:642493/250‑1 (MQ=255)
                                                                                                                                                                                                                          ccGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAg                                                        >  1:724244/1‑212 (MQ=255)
                                                                                                                                                                                                                           cGTTGCGGGTTAATACGGCATAAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTAcat                <  1:796432/251‑1 (MQ=255)
                                                                                                                                                                                                                              tGCGGTTTAATACGGCATGAAACCGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTAcatcat             <  1:743525/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gggTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGc                                                                                                                                                                                                           >  1:206137/1‑58 (MQ=255)
                                                                                                                                                                                                                                 gggTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTAc                                                                                                   >  1:401383/1‑162 (MQ=255)
                                                                                                                                                                                                                                  ggTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAg                                                                                                                                                                       <  1:567668/93‑1 (MQ=255)
                                                                                                                                                                                                                                  ggTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAgcgc                                                                                                                                                      <  1:613256/110‑1 (MQ=255)
                                                                                                                                                                                                                                  ggTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCg                                                                                                       >  1:362065/1‑157 (MQ=255)
                                                                                                                                                                                                                                  ggTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAg                                           >  1:417177/1‑217 (MQ=255)
                                                                                                                                                                                                                                  ggTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTACATCATCGtt         >  1:369307/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ggTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTACATCATCGtt         >  1:137743/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ttAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTACATCATCGTTAt       <  1:197249/251‑1 (MQ=255)
                                                                                                                                                                                                                                         aCGGCATGAAAACGCGGTTCCTTATCTAAAAGCGTTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAATTTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTACATCATCGTTATTTAAt  <  1:727814/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gCATGAAAACGCGGATCCTTATCTAAAAGCGTTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCAAAAATTTTTGGGTGGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATtg                        <  1:374669/226‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              
CGAGCACGCTACCCGCTCCGCCCAACAGCTCCGAACGTTTTTCATGATCTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAATGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTACATCATCGTTATTTAAT  >  1/689595‑690078

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