Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 625,930 T→C N369D (AAC→GAC)  unknown ← N‑acetylneuraminate epimerase (EC 5.1.3.24)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1625,9300TC100.0% 182.9 / NA 67N369D (AAC→GAC) unknownN‑acetylneuraminate epimerase (EC 5.1.3.24)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/31);  total (36/31)

ATTTTTTATTAAACGGAAGAGAGCATATGACTATGCGATTTCCTTTGGCAAATAAAAA‑CTTCCGTTTTTGTGAGAAGTATCACGCAACGGCGTGGTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTTGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGTAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTCTTTGTGCCAGGCCTTGCTTAAGCCTTTGTGGGCGTGAAGCATGCCTTCTTTGAATTGTTTTACGGAGCCCGGGAAGTTTGCGC  >  1/625690‑626172
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
aTTTTTTATTAAACGGAAGAGAGCATATGCCTATGCGATTTCCTTTGGCAAATAAAAA‑CTTCCGTTTTTGTGAGAAGTATCACGCAACGGCGTGGTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATg                                                                                                                                                                                                                                          >  1:313912/1‑251 (MQ=255)
aTTTTTTATTAAACGGAAGAGAGCATATGCCTATGCGATTTCCTTTGGCAAATAAAAA‑CTTCCGTTTTTGTGAGAAGTATCACGCAACGGCGTGGTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATg                                                                                                                                                                                                                                          >  1:393508/1‑251 (MQ=255)
 ttttttATTAAACGGAAGAGAGCATATGCTTATGCGTTTTCCTTGGGAAAAAAAAAA‑CTTCTTTTTTTGTGAAAAGTATCCCGCACCGGCGTGTTATTGGCGGTGGTTTTTTTAAAAACGG‑TCCC‑TTTTT‑TCGGAAATTTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGc                                                                                                                                                                                                                                         <  1:572629/251‑1 (MQ=255)
 tttaaaagagacaggAAGAGAGAATATGCCTATGC‑TTTTTCTTTGGAAAAAAAAAACCTTCCTTTTTTGGGAAAAGTATCACGCACCGGCGTGTTTTTGGGGGTGTTTTATTTAAAACCGGTTCCC‑TTGTT‑CCATAAT‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGc                                                                                                                                                                                                                                         <  1:842816/238‑1 (MQ=255)
   ttttATAAAAGGGAAGAGAGCATATGCTTATGCGTTTTCTTTTGGAAAAAAAAAA‑CTTCCTTTTTTGGGAGAATTACCCCGCACCGGCGTGTTATGGGGGGTGG‑TTTTTTTAAAACCGTTCCC‑TTTTTCCCGTAAT‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCtt                                                                                                                                                                                                                                       <  1:165539/251‑1 (MQ=255)
            aCGGAAGAGAGCATATGCCTATGCGATTTCCTTTGGCAAATAAAAA‑CTTCCGTTTTTGTGAGAAGTATCACGCAACGGCGTGGTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGa                                                                                                                                                                                                                               >  1:898767/1‑250 (MQ=255)
            aCGGAAGAGAGCATATGCCTATGCGATTTCCTTTGGCAAATAAAAA‑CTTCCGTTTTTGTGAGAAGTATCACGCAACGGCGTGGTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGa                                                                                                                                                                                                                               >  1:814266/1‑250 (MQ=255)
                 agagagCATATGCCTATGCGATTTCCTTTGGCAAATAAAAA‑CTTCCGTTTTTGTGAGAAGTATCACGCAACGGCGTGGTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTaa                                                                                                                                                                                                                         >  1:109466/1‑251 (MQ=255)
                      gCATATGCCTATGCGATTTCCTTTGGCAAATAAAAA‑CTTCCGTTTTTGTGAGAAGTATCACGCAACGGCGTGGTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGc                                                                                                                                                                                                                                         >  1:733766/1‑230 (MQ=255)
                             ccTATGCGATTTCCTTTGGCAAATAAAAA‑CTTCCGTTTTTGTGAGAAGTATCACGCAACGGCGTGGTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTgccgc                                                                                                                                                                                                              >  1:157128/2‑250 (MQ=255)
                              cTATGCGATTTCCTTTGGCAAATAAAAA‑CTTCCGTTTTTGTGAGAAGTATCACGCAACGGCGTGGTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGccgccg                                                                                                                                                                                                            >  1:673750/1‑251 (MQ=255)
                                    gATTTCCTTTGGCAAATAAAAA‑CTTCCGTTTTTGTGAGAAGTATCACGCAACGGCGTGGTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGtt                                                                                                                                                                                                      >  1:40757/1‑251 (MQ=255)
                                      tttCTTTTGGAAAAAAAAAT‑CTTCCTTTTTTGTGAGAAGTATCCCGCAGCGGCGTGTTATTGGCGGTGTTTTTTTTTAAAACGTTTCCT‑TTGTT‑CCGTAAT‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTc                                                                                                                                                                                                    <  1:748948/251‑1 (MQ=255)
                                         cttttgcaaaaaaaaaa‑CTTCTTTTTTTGTGAAATGTACCACGCAACGGCGTGGTATTAGGGGTGTTTTTTTTAAAAACGGTTCCC‑TTGTT‑CCATAAT‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTcacc                                                                                                                                                                                                 <  1:761246/239‑1 (MQ=255)
                                         ccTTTGGCAAATAAAAA‑CTTCCGTTTTTGTGAGAAGTATCACGCAACGGCGTGGTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCt                                                                                                                                                                                                                                        >  1:16593/1‑212 (MQ=255)
                                                  aaTAAAAA‑CTTCCGTTTTTGTGAGAAGTATCACGCAACGGCGTGGTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTaa                                                                                                                                                                                        >  1:66438/1‑251 (MQ=255)
                                                                           aattATCCCGCAGCGGCGTGTTATTAGGGTTGGTTTA‑TTTAAAACCGTTCCC‑TTTTTCCCTTAAT‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGac                                                                                                                                                                  <  1:423742/246‑1 (MQ=255)
                                                                           aaGTCTCACGCACCGGGGTGGTTTTGGGGGTGGTTTATTTAAAAACGTTCCCC‑TTGTC‑CCGTAAT‑TTTGGGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGac                                                                                                                                                                  <  1:417991/249‑1 (MQ=255)
                                                                                          gCGTGGTATTGGGGGTGTTTTATTTAAAAACGGTCCCT‑TTGTT‑CCGTAAT‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGc                                                                                                                                                  <  1:817787/250‑1 (MQ=255)
                                                                                           tgtgtttttgggggtgttttATTTAAAACCGGTTCCC‑TTGTT‑CCGTAAT‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGcc                                                                                                                                                 <  1:528068/232‑1 (MQ=255)
                                                                                           ggtggttttgggggtgttttATTTTAAAACGGTTCCT‑TTTTT‑CCATAAT‑TTTGTCCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGcc                                                                                                                                                 <  1:77662/232‑1 (MQ=255)
                                                                                              ggTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAacac                                                                                                                                                        >  1:390169/1‑240 (MQ=255)
                                                                                               gTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCt                                                                                                                                                                                                                                              >  1:57393/1‑153 (MQ=255)
                                                                                               gTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTc                                                                                                                                                                                                                                             >  1:388603/1‑154 (MQ=255)
                                                                                               gTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGc                                                                                                                                                                                                                      >  1:878187/1‑177 (MQ=255)
                                                                                               gTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGtt                                                                                                                                                                                                      >  1:168215/1‑193 (MQ=255)
                                                                                               gTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATa                                                                                                                                             >  1:218071/1‑250 (MQ=255)
                                                                                               gTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATa                                                                                                                                             >  1:488626/1‑250 (MQ=255)
                                                                                                ttttggggggtgttttttttAAAAACGTTCCCTTTTTT‑CCTTAAT‑TTTGTCCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATAc                                                                                                                                            <  1:406660/235‑1 (MQ=255)
                                                                                                tATTGGCGGTGGTTTTTTTAAAAACGGTCCCT‑TTTTT‑CATTAAT‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATAcc                                                                                                                                           <  1:406347/251‑1 (MQ=255)
                                                                                                          ggCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAAt                                                                                                                                 >  1:298658/1‑251 (MQ=255)
                                                                                                               aTTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTCTcaccac                                                                                                                                                                                               >  1:178219/1‑184 (MQ=255)
                                                                                                                ttttAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCg                                                                                                                           >  1:55181/1‑251 (MQ=255)
                                                                                                                  ttAAAACGGGTCCC‑TTTTT‑TCCTAAATTTTGTCCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGa                                                                                                                          <  1:526309/251‑1 (MQ=255)
                                                                                                                                   gTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCg                                                                                                         >  1:317904/1‑251 (MQ=255)
                                                                                                                                       ccGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCgttg                                                                                                      >  1:154378/1‑251 (MQ=255)
                                                                                                                                               tttGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAgg                                                                                               <  1:295454/251‑1 (MQ=255)
                                                                                                                                                         ccACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGacac                                                                                                                                                                >  1:845880/1‑176 (MQ=255)
                                                                                                                                                            caATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAg                                                                                                                                                                                     <  1:21489/152‑1 (MQ=255)
                                                                                                                                                            caATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGc                                                                                                                                                  >  1:142019/1‑187 (MQ=255)
                                                                                                                                                            caATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATAc                                                                                                                                            >  1:609764/1‑193 (MQ=255)
                                                                                                                                                            caATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATaa                                                                                                                       >  1:237675/1‑214 (MQ=255)
                                                                                                                                                                 aCAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTCTTTGTGcc                                                                             <  1:689489/251‑1 (MQ=255)
                                                                                                                                                                  cAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAacac                                                                                                                                                        >  1:57089/1‑175 (MQ=255)
                                                                                                                                                                  cAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTCTTTGTGCCa                                                                            >  1:859850/1‑251 (MQ=255)
                                                                                                                                                                    aTGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGac                                                                                                                                                                  <  1:241966/163‑1 (MQ=255)
                                                                                                                                                                         ggATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAg                                                                                                                                                                                                                             <  1:109070/99‑1 (MQ=255)
                                                                                                                                                                            ttAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTATTTTTTTGCCTTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTCTTTGTGCCAGGCCTTGCtt                                                                  <  1:340846/251‑1 (MQ=255)
                                                                                                                                                                            ttAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTCTTTGTGCCAGGCCTTGCtt                                                                  <  1:283336/251‑1 (MQ=255)
                                                                                                                                                                            ttAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTCTTTGTGCCAGGCCTTGCtt                                                                  <  1:882232/251‑1 (MQ=255)
                                                                                                                                                                            ttAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTCTTTGTGCCAGGCCTTGCtt                                                                  <  1:357351/251‑1 (MQ=255)
                                                                                                                                                                              acgccaatgacgtttAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGcc                                                                                                                                                 <  1:16116/156‑1 (MQ=255)
                                                                                                                                                                               aaCGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGac                                                                                                                                                                  >  1:305348/1‑152 (MQ=255)
                                                                                                                                                                                 cGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTCTTTGTGCCAGGCCTTGCTTAAGc                                                              >  1:891321/1‑250 (MQ=255)
                                                                                                                                                                                   tCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGc                                                                                                                                                                                    >  1:627744/1‑130 (MQ=255)
                                                                                                                                                                                   tCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTCTTTGTGCCAGGCCTTGCt                                                                   <  1:891820/243‑1 (MQ=255)
                                                                                                                                                                                   tCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTCTTTGTGCCAGGCCTTGCTTAAGCCtt                                                           >  1:601195/1‑251 (MQ=255)
                                                                                                                                                                                         gCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTCTTTGTGCCAGGCCTTGCTTAAGCCTTTGTggg                                                     >  1:108390/1‑251 (MQ=255)
                                                                                                                                                                                                ccTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTc                                                                                     >  1:229851/1‑212 (MQ=255)
                                                                                                                                                                                                   tCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTCTTTGTGCCAGGCCTTGCTTAAGCCTTTGTGGGCGTGAAGCAt                                           >  1:651962/1‑251 (MQ=255)
                                                                                                                                                                                                      ttttttACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTCTTTGTGCCAGGCCTTGCTTAAGCCTTTGTGGGCGTGAAGCATGcc                                        <  1:196792/251‑1 (MQ=255)
                                                                                                                                                                                                        ttttACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCgttg                                                                                                      <  1:460827/187‑1 (MQ=255)
                                                                                                                                                                                                                 tttGACAGGTTATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTCTTTGTGCCAGGCCTTGCTTAAGCCTTTGTGGGCGTGAAGCATGCCTTCTTTGAAtt                             <  1:661654/251‑1 (MQ=255)
                                                                                                                                                                                                                          ttATTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTCTTTGTGCCAGGCCTTGCTTAAGCCTTTGTGGGCGTGAAGCATGCCTTCTTTGAATTGTTTTACgg                    <  1:710336/251‑1 (MQ=255)
                                                                                                                                                                                                                            aTTCGACACTTAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTCTTTGTGCCAGGCCTTGCTTAAGCCTTTGTGGGCGTGAAGCATGCCTTCTTTGAATTGTTTTACgg                    <  1:176206/249‑1 (MQ=255)
                                                                                                                                                                                                                                     ttAATTTTTTGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTCTTTGTGCCAGGCCTTGCTTAAGCCTTTGTGGGCGTGAAGCATGCCTTCTTTGAATTGTTTTACGGAGCCCGGGAAg         <  1:811226/251‑1 (MQ=255)
                                                                                                                                                                                                                                            tttGCCGTCGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGCAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTCTTTGTGCCAGGCCTTGCTTAAGCCTTTGTGGGCGTGAAGCATGCCTTCTTTGAATTGTTTTACGGAGCCCGGGAAGTTTgcgc  <  1:634410/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
ATTTTTTATTAAACGGAAGAGAGCATATGACTATGCGATTTCCTTTGGCAAATAAAAA‑CTTCCGTTTTTGTGAGAAGTATCACGCAACGGCGTGGTATTAGCGGTGGCTTATTTTAAAACGGTTCCC‑TTGTT‑CCGTAAA‑TTTGTGCCCGCCACAATAACAATGAAGGATTAAACGTCATTGGCGTAATCCTTCATTTTTTACGTCTTTGACAGGTTATTCGACACTTAATTTTTTGCCGTTGTAGCTCATGCTTTGTACGGAAGTTAAGGCTTTGCCGCCGTTGGTTTCACCACCGATTAAAAGCACTTGATTGTTGTAAGACACCGTAACACCATAGCCGATACCGGCAGGCAATTCACCGATAATTTGCCATTTACCGTTGTTTAAGGTGTAAATGTCTTTGTGCCAGGCCTTGCTTAAGCCTTTGTGGGCGTGAAGCATGCCTTCTTTGAATTGTTTTACGGAGCCCGGGAAGTTTGCGC  >  1/625690‑626172

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